
No.4ベストアンサー
- 回答日時:
と言うことはご質問の趣旨は「Laplace変換のある伝達関数X(s)が与えられたとき、それをz変換の伝達関数X(z)に置き換えたい」ということでしょうか。
単に「SをZで表すにはどうしたら良いのでしょうか」と質問すれば、nubouさんや私のような回答が返ってくるのは無理からぬことです。
X(s)からX(z)に変換するにはいくつかの手法があります。
1. 標準z変換
2. 双1次z変換
3. 整合z変換
まず1.から。この手法は「X(z)のインパルス応答がX(s)のインパルス応答と等しくなるように、X(z)の各項の係数を決める」という設計思想に基づきます。ただしすべての時間tにおいて等しくなるわけでなく、サンプル時点での出力(x(kτ)、k=0, 1, 2, 3....)だけは少なくとも同じになるように決めます。この方法にはさらにやり方が
1-a z=exp(sτ)を代入する
1-b 部分分数に分ける
の二つあります。
1-aでz=exp(sτ)をそのまま代入すると実用的でないので、通常はさらに
z=exp(sτ)≒1+sτ (4)
の近似を行います。これは直ちに
s=(z-1)/τ (5)
と書き改められますから、これをX(s)に代入すればよいことになります。(delta-sさんのおっしゃる(1-z^(-1))と本質的に同じですが、サンプリング時間τが入ってくることをお忘れなく)
この近似は「サンプリング時間τが系の時定数に比べ十分に小さい」場合にのみ成り立ちますが、通常の信号処理回路ではそのように設計されているはずです(連続時間近似)。
もう一つ、1-bのX(s)を部分分数展開する方法では、X(s)が
∞
Σ R[k]/(s+α[k]) (6)
k=0
と部分分数に分解できたとして、それぞれの項で
1/(s+α[k])→1/(1-exp(-α[k]τ)(z^(-1)) (7)
と置き換えます。ここにR[k]はk番目の項の係数、α[k]はk番目の極です。
2.の双1次z変換では
z=exp(sτ)≒(1+sτ/2)/(1-sτ/2) (8)
と近似します。sについて書き下せば
s=(2/τ)(z-1)/(z+1) (9)
となります。これを代入します。
双1次z変換ではインパルス応答はX(s)と異なってきます。一方、折返し誤差は生じません(*1)
3.はすみません、私は名前しか知りません。
Laplace変換での微分要素sですが、離散時間の信号に微分という演算が存在しないのはご承知の通りです。ただしLaplace変換の伝達関数X(s)のsに相当する要素を考えることはできます。それが(5)や(9)です。
私の最初の回答の(1)→(2)の変形で、x(t)をx(kτ)で置き換えていますが、これはで標本化操作(離散化)に他ならないことにご注意ください。連続関数を直接z変換することはできないので、いったんx(kτ)に標本化し、それをz変換していることになります。
とりあえずキーワードと手がかりは上記に示しましたので、恐縮ですがあとはご自身で参考文献を探して下さいますよう。(当初の質問範囲を越える追加質問については、補足欄への記入でなく新しい質問を立ててください)
--------
*1 離散時間での信号処理では、高い周波数の成分が低い周波数に化けて出てくることがあります(エイリアス効果)。例えば、本来除去したい高い周波数の信号が低い周波数の信号に化けるとそれは除去できません。これを「折り返し誤差」などと呼びます。詳しくは自分で調べてみて下さい。
*2 z変換でkについての和を-∞から∞まで取るか(両側z変換)、0から∞まで取るか(片側z変換)の議論はありますが、信号処理では因果律(入力がある以前に出力があってはおかしい)を考慮して主に片側z変換が用いられるようです。
すごく丁寧でわかりやすいご回答、ありがとうございます。
Umadaさんの回答で、理解することが出来ました。
度々、丁寧なご回答、ありがとうございました。
No.3
- 回答日時:
「連続時間での微分は、S、離散時間での微分は、(1-z^(-1))と表せますが」:
1-z^(-1)は
{c(n/fs)}⇒{c(n/fs)-c((n-1)/fs)}
の操作であり「微分」と言うよりは「差分」なのです
「連続時間から離散時間に変換する時にSは、(1-z^(-1))としても良いのですか?」:
z=exp(s/fs)ですからs=(1-exp(-s/fs))と見なせる条件下ではいいですがそうでないときは駄目です
周波数帯域とfsの大きさから判断すべきです
s=j・2・π・fですからj・2・π・f=(1-exp(-j・2・π・f/fs))
と見なせる「fの範囲」においてはokでしょう
No.2
- 回答日時:
連続時間信号は実数全体を定義域とする実数値関数cで表現することができます
以下tを実数変数としnを整数変数とし
ある実数t0が存在してt<t0とのきc(t)=0とする
連続時間信号を表現する実数値関数c(t)をt=0を含む標本化周波数fsの標本化をしてできた離散時間信号は整数全体を定義域とする実数値関数dで表現することができます
ただしd(n)≡c(n/fs)です
dはcと同じように扱うことができないので離散時間信号を表現するためdの代わりに以下のように定義される実数全体を定義域とする離散時間関数pを使います
p(t)≡Σ(-∞<n<∞)・c(n/fs)・δ(t-n/fs)
c(t)の両側ラプラス変換をCs(s)とすると
Cs(s)≡∫(-∞<t<∞)dt・exp(-s・t)・c(t)
p(t)の両側ラプラス変換をPs(s)とすると
Ps(s)≡∫(-∞<t<∞)dt・exp(-s・t)・p(t)
=Σ(-∞<n<∞)・c(n/fs)・exp(-s・n/fs)
p(t)の両側z変換をPz(z)とすると
Pz(z)≡Σ(-∞<n<∞)・c(n/fs)/z^n
exp(-s・n/fs)=1/(exp(s/fs))^nであるから
Ps(s)=Pz(exp(s/fs))です
すなわち離散時間信号を表現する関数において次のことが言える
両側z変換においてz=exp(s/fs)とすれば両側ラプラス変換が得られる
連続時間信号を表現する関数にはz変換は定義されていませんのでこのような関係は存在しません
No.1
- 回答日時:
離散時間系の信号処理のお話か、あるいは制御工学のお話ですね。
sはLaplace変換のsのことですよね。zとsの変換は特に定理というほど大袈裟なものはありません。結論から申し上げますと通常は、サンプリングの周期をτとしてz=exp(sτ)とおきます。
以下はj簡単な背景の説明です。
ある信号x(t)があったとき、そのLaplace変換を
X(s)=∫exp(-st) x(t) dt (0→∞) (1)
とします。ごくごく一般的なLaplace変換です。
さて、周期τでサンプリングする離散時間系で考えるなら積分は数列の和に置き換えることができ、上記のX(s)に相当する概念として
∞
Σexp(-skτ) x(kτ) (2)
k=0
なる式を導入することができます。
このとき、簡単のためz=exp(sτ)とおけば(2)は
∞
Σ(z^(-k)) x(kτ) (3)
k=0
と書き直すことができます。これを数列x(kτ)のz変換と名付けたわけです。
X(z)などと表示されるのはご存じの通りです。
ご回答ありがとうございます。
信号処理の話に関してです。
連続時間での微分は、S、離散時間での微分は、(1-z^(-1))と表せますが
連続時間から離散時間に変換する時に
Sは、(1-z^(-1))としても良いのですか?
アナログのスイッチドキャパシタLPFの伝達関数をZを使って表す方法が
分からないので、ここで質問してみました。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 多変数関数の微分とテイラー展開について 5 2022/04/24 16:55
- 国産車 回答宜しくお願いします。 今度中古のHONDAステップワゴン購入考えておりますが、DVA-RKの型式 1 2022/03/30 12:11
- 数学 次の問題についてご教授願いします。 3点z=0 , z=-1 , z=i をそれぞれω=0 , ω= 4 2023/05/23 20:13
- 統計学 t分布導出時のヤコビ行列式について教えて下さい。 1 2022/07/04 21:36
- 数学 「f(z)=1/(z^2-1)に関して ローラン展開を使う場合、マクローリン展開を使う場合、テイラー 3 2022/08/27 19:56
- 数学 paythonを使用した周回積分に関する質問です。 2 2023/02/17 19:09
- Excel(エクセル) 【至急 詳しい方教えて下さい】Excelの表を変換したい 6 2022/04/21 18:35
- Excel(エクセル) Excelシフト表 固定シフトの自動変換化 1 2022/04/14 16:10
- 統計学 ある変数XとYの相関係数が0.8でZの影響を除いたXとYの偏相関係数が-0.8の時、XとZ,YとZの 3 2023/04/30 19:26
- 数学 画像の問題について質問です。問題式を楕円の式に変形して、積分範囲を0<=x<=a √(z^2-1) 3 2022/08/29 13:44
今、見られている記事はコレ!
-
企業へPTA業務をアウトソースする仕組みとは?PTA専用支援サービスに聞いてみた
新年度がはじまり、学校のPTAも仕事内容や役員決めなどが始まっている頃だろう。最近では共働きの家庭が増え、仕事との両立が難しくPTA活動に負担を感じる保護者が増えているためか、PTA業務代行サービスが話題にな...
-
中学時代の経験は、その後の人生に影響するのか?教育研究家に聞いた!
皆さんは、中学時代の過ごし方の人生への影響力について考えたことはあるだろうか。「教えて!goo」にも「人生は中学時代で決まりますか?」という質問が投稿され、閲覧者の反響を呼んだ。中学時代の経験は、その後...
-
今年も英検が始まる!英検の面接官はどんなポイントを評価する?英検コーチに聞いてみた
英検の2022年度第1回検定がはじまる。前年度の検定後、「教えて!goo」で「英検2級の二次面接、手応えなしで合格したのはなぜでしょう?」という投稿が注目を集めた。自分の手応えとは違う評価になることもある面接...
-
葬儀をテーマにした中学受験の入試問題を教育者ではなく葬儀社に聞いてみた
2月から中学受験のシーズンに突入する。そこで今回は中学受験に出題された問題の中でも、特に葬儀や死をテーマにした問題を紹介する。新型コロナウィルスによって、今まで以上に病気や死について考える機会が多かっ...
-
ファッション選びやプレゼン資料作成にも便利!?「色彩検定」って何だ
先日、久々に会った友人が「ワイン検定」の試験に向けて勉強していると言っていた。友人がワイン好きなのは知っていたが、まさかそこまで本腰を入れているとは思わなかったので驚いた。 試験というと大学の受験勉強...
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
緯度・経度の値をエクセルで変...
-
勝手に無変換になる
-
エクセル関数/10進法から60進法...
-
青の吹き出しの何をどう考えれ...
-
意外に?以外に??どちらが正...
-
Z変換について教えてください。
-
VB.NETで小数点以下の桁数を取...
-
Google 日本語入力の調子がおか...
-
CADソフトの拡張子変換
-
スクリーンショット:文字がぼ...
-
アマレココの動画をmp4に変換し...
-
ExcelでTEXT()計算後、SUMできない
-
エクセルで、小数点は、全て切...
-
MP3→WAVE変換 一体どれが一番...
-
iTunes でMP3へ変換するには?
-
ヤコビアンとはなんですか?
-
「送る」の謙譲語
-
vmxファイルをisoイメージファ...
-
レントゲン画像のプリントアウ...
-
スプレッドシートから、スプレ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
緯度・経度の値をエクセルで変...
-
エクセル関数/10進法から60進法...
-
勝手に無変換になる
-
意外に?以外に??どちらが正...
-
無料でも使えるSSDを、MBRからG...
-
つらそう…を丁寧に言うなら、お...
-
●Illustrator CMYKをRGB変換し...
-
Wordに囲い文字というのがない...
-
VB.NETで小数点以下の桁数を取...
-
カーナビにスマホ画面を表示さ...
-
フォルダ、ファイル名の一括変...
-
スクリーンショット:文字がぼ...
-
MP3を楽譜に変換するサイト、ア...
-
pdfの色データをCMYK→K+Cの2色...
-
土日祝 の正しい読みは
-
CADソフトの拡張子変換
-
エクセル DEC2BIN 範囲外数値の...
-
単位の リューベ って
-
Windows Media PlayerからMP3に...
-
Google 日本語入力の調子がおか...
おすすめ情報