Z変換

数学信号処理においてZ変換は実数または複素数である離散時間信号を複素数値の周波数領域Z領域またはZ平面)表現に変換する[1] [2] [3]

これはラプラス変換s領域またはs平面)の離散時間版と考えることができる[4]この類似性は時間スケール計算理論で探求されている

連続時間フーリエ変換はs領域の垂直軸(虚軸)に沿って評価されますが、離散時間フーリエ変換はz領域の単位円に沿って評価されます。s領域の左半平面はz領域の単位円の内側の領域にマッピングされ、s領域の右半平面はz領域の単位円の外側の領域にマッピングされます。

信号処理において、デジタルフィルタを設計する手段の一つは、アナログ設計を双一次変換(s領域からz領域への写像)にかけ、検査、操作、あるいは数値近似によってデジタルフィルタを作成することである。このような方法は、複素数1の近傍、すなわち低周波数を除いて、正確ではない傾向がある。

歴史

デジタル制御システムの解析と設計の礎石として現在認識されているZ変換の基礎概念は、20世紀半ばに登場した当時、全く新しいものではなかった。その萌芽的な原理は、密接に関連する数学的手法であるラプラス変換でよく知られるフランスの数学者ピエール=シモン・ラプラスの研究に遡ることができる。しかし、現在Z変換として理解されているものの明示的な定式化と応用は、1947年にヴィトルド・ヒューレヴィッツとその同僚によって大きく進歩した。彼らの研究は、当時レーダー技術の文脈で重要性を増していたサンプル値制御システムの課題に触発されたものであった。Z変換は、離散時間信号およびシステムの解析で広く用いられる定数係数の線形差分方程式を解くための体系的かつ効果的な手法を提供した。[5] [6]

この手法は、コロンビア大学のサンプルデータ制御グループに所属していたジョン・R・ラガッツィーニロトフィ・A・ザデーの尽力により、1952年にさらに改良され、「Z変換」という正式な名称が付けられました。彼らの研究は、Z変換の数学的枠組みを確固たるものにしただけでなく、特に電気工学と制御システムの分野において、その応用範囲を拡大しました。[7] [8]

改良Z変換または高度Z変換として知られる注目すべき拡張は、後にエリアフ・I・ジュリーによって導入されました。ジュリーの研究は、特に初期条件の取り扱いにおいてZ変換の適用性と堅牢性を拡張し、デジタル制御システムの解析のためのより包括的な枠組みを提供しました。この高度な定式化は、離散時間制御システムの設計と安定性解析において極めて重要な役割を果たし、デジタル信号処理分野に大きく貢献しました。[9] [3]

興味深いことに、Z 変換の概念的基盤は、組合せ論と確率論における強力なツールである生成関数法として知られるより広範な数学的概念と交差しています。 このつながりは、確率論発展の先駆者であるアブラハム・ド・モアブルによって 1730 年にはすでに示唆されていました。 ド・モアブルは生成関数を利用して確率の問題を解き、最終的に Z 変換へと発展するものの基礎を築きました。 数学的な観点からは、Z 変換はローラン級数の特定の例と見なすことができます。ローラン級数では、調査対象の数のシーケンスが解析関数の (ローラン) 展開における係数として解釈されます。 この観点は、Z 変換の深い数学的ルーツを強調するだけでなく、数学と工学のさまざまな分野にわたるその汎用性と幅広い適用性を示しています。[3]

意味

Z変換は片側変換または両側変換として定義できます。(片側ラプラス変換両側ラプラス変換があるのと同じです。)[10]

双方向Z変換

離散時間信号の両側 Z 変換は、ように定義される正式なべき級数です

ここでは整数であり、は一般に複素数です。極形式では、は次のように表されます。

ここで、 は の大きさ虚数単位、 はラジアン単位の複素引数(角度または位相とも呼ばれるです

片側Z変換

あるいは、がに対してのみ定義されている場合片側Z 変換は次のように定義されます

信号処理では、この定義は離散時間因果システムの単位インパルス応答の Z 変換を評価するために使用できます

片側Z変換の重要な例として、確率生成関数が挙げられます。ここで、 は離散確率変数が値をとる確率です。Z変換の性質(§ 性質に列挙されている)は、確率論の文脈において有用な解釈をもたらします。

逆Z変換

Z 変換は次のとおりです

ここで、 は原点を囲み、収束領域(ROC)内に完全に収束する反時計回りの閉経路です。ROCが因果関係にある場合(例2を参照)、これは経路が のすべての極を囲む必要があることを意味します

この路面積分の特殊なケースは、が単位円であるときに発生します。この路面積分は、ROCが単位円を含む場合に使用できます。これは、 が安定しているとき、つまりすべての極が単位円の内側にあるとき、常に保証されます。この路面積分を用いると、逆Z変換は、単位円の周りのZ変換の周期値の逆離散時間フーリエ変換、またはフーリエ級数に簡略化されます。

有限範囲のと有限個の均一間隔の値を持つZ変換は、ブルースタインのFFTアルゴリズムを用いて効率的に計算できます離散時間フーリエ変換(DTFT)(離散フーリエ変換(DFT)と混同しないでください)は、 を単位円上に制限することによって得られるZ変換の特殊なケースです。

逆変換の評価には、以下の3つの方法がよく使用されます。

輪郭積分による直接評価

この手法では、コーシー留数定理を用いて逆Z変換を評価します。複素平面上の閉曲線の周りを積分することで、ROC内のZ変換関数の極における留数が合計されます。この手法は、複素変数で表される関数を扱う場合に特に有用です。

変数の項の列への展開zそしてz−1

この方法では、Z変換をべき級数に展開します。このアプローチは、Z変換関数が有理数の場合に有用であり、級数に展開して信号係数を項ごとに決定することで、逆変換の近似が可能になります。

部分分数展開と表参照

この手法では、Z変換をより単純な分数の和に分解します。各分数は既知のZ変換ペアに対応します。逆Z変換は、Z変換ペアの標準テーブルで各項を参照することで決定されます。この手法は、その効率性と簡便性から広く利用されており、特に元の関数が容易に認識可能な要素に分解できる場合に有効です。

[11]

A) 次の逆Z変換を級数展開法で決定します。

解決:

ケース1:

中華人民共和国:

ROCは円の外部であるため、因果関係があります(n ≥ 0に信号が存在します)。したがって、 (矢印はx(0)= 1の項を示します)。

長除算プロセスの各ステップで、 の最小のべき乗項を削除することに注意してください

ケース2:

中華人民共和国:

ROC は円の内部であるため、反因果的です ( n < 0 の場合に信号が存在します)。

長割り算をすると

(矢印はx (0)=0の項を示す)。

長除算プロセスの各ステップで、 の最小のべき乗項を削除することに注意してください

注記:

  1. 信号が因果関係にある場合、 の正の累乗が得られ、信号が反因果関係にある場合、 の負の累乗が得られます
  2. は の期間を示しは の期間を示します

B) 次の逆Z変換を級数展開法で決定します。

を で割ると負の累乗が消去され

部分分数展開により、

ケース1:

中華人民共和国:

どちらの用語も因果関係があり、したがって因果関係があります。

ケース2:

中華人民共和国:

どちらの用語も反因果的であるため、 は反因果的です。

ケース3:

中華人民共和国:

項の 1 つは因果関係があり (p=0.5 で因果関係の部分が示されます)、もう 1 つは反因果関係があり (p=1 で反因果関係の部分が示されます)、したがって両側性があります。

収束領域

収束領域ROC)は、Z変換の合計が絶対収束する複素平面上の点の集合である。

例1(ROCなし)

間隔を拡大する

合計を見ると

したがって、この条件を満たすの値は存在しません。

例2(因果ROC)

ROC (青)、| z | = 0.5 (破線の黒い円)、および単位円 (点線の灰色の円)。

(ここでヘヴィサイドの階段関数)とする。区間を展開する

合計を見ると

最後の等式は無限等比級数から生じ、等式が成り立つのは、 の場合のみであり、 はとして書き直すことができる。したがって、 ROC は である。この場合、 ROC は、原点に半径 0.5 の円盤が「打ち抜かれた」複素平面である。

例3(反因果的ROC)

ROC (青)、| z | = 0.5 (破線の黒い円)、および単位円 (点線の灰色の円)。

(ここでヘヴィサイドの階段関数)とする。区間を展開する

合計を見ると

そして、再び無限等比級数を使用すると、等式は の場合にのみ成立し、はとして書き直すことができます。したがって、ROC は です。この場合、ROC は原点を中心とし、半径 0.5 の円盤です。

この例と前の例の違いはROCだけです。これは、変換結果だけでは不十分であることを示すために意図的に加えられたものです。

例の結論

例2と例3は、ROCを指定した場合にのみ、Z変換が一意であることを明確に示しています。因果関係のあるケースと反因果関係のあるケースの極零点プロットを作成すると、どちらのケースのROCにも0.5の極が含まれていないことがわかります。これは複数の極を持つケースにも当てはまります。つまり、ROCに極が含まれることはありません

例2では、​​因果関係のシステムではROCが生成され、例3の反因果関係のシステムではROCが生成され、

ROCは青いリングで表示 0.5 < | z | < 0.75

複数の極を持つシステムでは、ROCがどちらの極も含まない可能性がある。ROCは円形の帯を形成する。例えば、

0.5と0.75に極を持つ。ROCは0.5 < | z | < 0.75となり、原点も無限大も含まない。このような系は因果項と反因果項を含むため、混合因果系と呼ばれる。

システムの安定性は、ROCのみを知ることでも判断できます。ROCが単位円(つまり| z | = 1)を含む場合システムは安定です。上記のシステムでは、因果システム(例2)は| z | > 0.5が単位円を含むため安定です。

ROCを持たない(つまり曖昧な)系のZ変換が与えられていると仮定しましょう。以下の条件を満たす場合、一意の を決定できます。

  • 安定性
  • 因果関係

安定性のためには、ROCは単位円を含まなければなりません。因果系が必要な場合は、ROCは無限大を含み、システム関数は右辺列になります。反因果系が必要な場合は、ROCは原点を含み、システム関数は左辺列になります。安定性と因果性の両方が必要な場合は、システム関数のすべての極が単位円内になければなりません。

そうすれば、ユニークなものが見つかります。

プロパティ

Z変換の特性

財産

時間領域Zドメイン証拠中華民国
Z変換の定義(Z変換の定義)

(逆Z変換の定義)

直線性ROC 1 ∩ ROC 2を含む
時間拡張

デシメーションohio-state.edu または ee.ic.ac.uk
時間遅延

ROC ただし
時間を進める

双方向Z変換:

片側Z変換: [12]

最初の違いは逆順

〜のため

ROC
最初の違いは前進
時間の逆転
Z領域でのスケーリング
複雑な活用
実部
虚数部
Z領域における分化ROC(有理数の場合)

ROCは無理数の場合、境界を除外する可能性がある[13]

畳み込みROC 1 ∩ ROC 2を含む
相互相関ROC
蓄積
乗算-

パーセバルの定理

初期値定理: が因果関係にある場合、

最終値定理: の極が単位円の内側にある場合、

一般的なZ変換ペアの表

ここ:

単位(またはヘヴィサイド)ステップ関数で あり、

は離散時間単位インパルス関数である(連続時間版であるディラックのデルタ関数を参照)。2つの関数は、単位ステップ関数が単位インパルス関数の累積(累積合計)となるように選択される。

信号、Z変換、中華民国
11すべてのz
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17 [14]、正の整数の場合[13]
18、正の整数の場合[13]
19
20
21
22

フーリエ級数とフーリエ変換との関係

単位円として知られる領域 内の の値については、変換を 1 つの実変数の関数として定義して表すことができます。また、この双対変換はフーリエ級数に簡約されます。

これは、シーケンスの離散時間フーリエ変換(DTFT)としても知られています。この周期関数は、フーリエ変換周期的な和であるため、広く用いられている解析ツールとなっています。これを理解するために、を任意の関数 のフーリエ変換とします。この関数の、ある区間におけるサンプル数がシーケンスに等しいとします。すると、シーケンスの DTFT は次のように表すことができます。

ここでは秒単位、はヘルツ単位です。2つの級数を比較すると、 はサンプルあたりラジアン単位の正規化周波数であることがわかります。この値は に対応します。そして、 を代入することで、式1は(フーリエ変換)で表すことができます。

パラメータTが変化すると、式 2の各項はf軸に沿って離れたり近づいたりします。ただし、式 3では、中心は 2 π離れたままですが、幅は拡大または縮小します。シーケンスがLTI システムインパルス応答を表す場合、これらの関数はその周波数応答とも呼ばれますシーケンスが周期的である場合、その DTFT は 1 つ以上の高調波周波数で発散し、他のすべての周波数では 0 になります。これは、高調波周波数で振幅可変のディラックのデルタ関数を使用することでよく表されます。周期性のため、一意の振幅は有限個しか存在せず、それらははるかに単純な離散フーリエ変換(DFT)によって簡単に計算できます。( 「離散時間フーリエ変換」の「§ 周期データ」を参照してください。)

ラプラス変換との関係

双一次変換

一次変換は、連続時間フィルタ(ラプラス領域で表現)を離散時間フィルタ(Z領域で表現)に変換するために、またその逆の変換に使用できます。以下の置換が用いられます。

ラプラス領域の関数をZ領域の関数に変換する(タスティン変換)、または

Z領域からラプラス領域へ。双線形変換により、複素s平面(ラプラス変換)は複素z平面(z変換)に写像されます。この写像は(必然的に)非線形ですが、s平面の軸全体をz平面の単位円に写像するという点で有用です。したがって、フーリエ変換(軸上で評価されるラプラス変換)は離散時間フーリエ変換になります。これは、フーリエ変換が存在すること、すなわち軸がラプラス変換の収束領域にあることを前提としています。

スター付き変換

時間サンプリングされた関数の片側 Z 変換が与えられると、対応する星印の付いた変換によってラプラス変換が生成され、(サンプリング パラメータ) への依存性が復元されます。

逆ラプラス変換は、インパルスサンプリング関数として知られる数学的な抽象化です。

線形定数係数差分方程式

線形定係数差 (LCCD) 方程式は、自己回帰移動平均方程式に基づく線形システムの表現です

上式の両辺は、ゼロでなければ割ることができる。LCCDで正規化すると、次の式が書ける。

このLCCD方程式の形は、「現在の」出力が過去の出力、現在の入力、および以前の入力の関数であることをより明確にするのに有利である。

伝達関数

上記の式を Z 変換すると(線形性と時間シフトの法則を使用)、次のようになります。

ここで、およびはそれぞれ、およびの Z 変換です。(表記規則では通常、大文字を使用して、対応する小文字で表された信号の Z 変換を参照します。これは、ラプラス変換の表記規則に似ています。)

システムの伝達関数を整理すると次のようになります。

零点と極点

代数学の基本定理より、分子(の零点に対応を持ち、分母は根(極に対応)を持つ。伝達関数を零点と極を用いて書き直すと、

ここで、 は番目の零点、は 番目の極です。零点と極は一般に複素数であり、複素平面(z平面)上にプロットすると極零点プロットと呼ばれます。

さらに、 および にも零点と極が存在する可能性があります。これらの極と零点、および多次の零点と極を考慮すると、零点と極の数は常に等しくなります。

分母を因数分解することで部分分数分解が可能になり、これを時間領域に戻すことができます。これにより、システムのインパルス応答と線形定数係数差分方程式が得られます。

出力応答

このようなシステムが信号によって駆動される場合、出力はの部分分数分解を行い、逆Z変換を行うことで求めることができます。実際には、を乗算する前に分数分解を行うことで、容易に計算可能な逆Z変換を持つ項を持つの形式を生成することがしばしば有用です。

参照

参考文献

  1. ^ Mandal, Jyotsna Kumar (2020). 「Z変換ベースの可逆符号化」.可逆ステガノグラフィーと変換符号化による認証. 計算知能研究. 第901巻. シンガポール: Springer Singapore. pp.  157– 195. doi :10.1007/978-981-15-4397-5_7. ISBN 978-981-15-4396-8. ISSN  1860-949X. S2CID  226413693. Zは複素変数です。Z変換は、離散空間領域信号を複素周波数領域表現に変換します。Z変換はラプラス変換から派生します。
  2. ^ Lynn, Paul A. (1986). 「ラプラス変換とZ変換」. 『電子信号とシステム』 . ロンドン: Macmillan Education UK. pp.  225– 272. doi :10.1007/978-1-349-18461-3_6. ISBN 978-0-333-39164-8ラプラス変換とZ変換はフーリエ変換と密接に関連しています。Z変換特に離散信号や離散システムを扱うのに適しています。離散時間フーリエ変換よりも簡潔で便利な表記法を提供します。
  3. ^ abc Jury, Eliahu Ibrahim (1964). Z変換法の理論と応用. ニューヨーク: John Wiley & Sons. pp. XIII, 330 s.
  4. ^ Palani, S. (2021-08-26). 「離散時間信号とシステムのZ変換解析」. Signals and Systems . Cham: Springer International Publishing. pp.  921– 1055. doi :10.1007/978-3-030-75742-7_9. ISBN 978-3-030-75741-0S2CID  238692483。z変換はラプラス変換の離散版です。z変換離散時間システムの差分方程式を代数方程式に変換し、離散時間システム解析を簡素化します。ラプラス変換とz変換は共通していますが、ラプラス変換は連続時間信号とシステムを扱う点が異なります。
  5. ^ ER カナセウィッチ (1981).地球物理学における時系列解析。アルバータ大学。 186、249ページ。ISBN 978-0-88864-074-1
  6. ^ ER Kanasewich (1981). 地球物理学における時系列解析(第3版). アルバータ大学. pp.  185– 186. ISBN 978-0-88864-074-1
  7. ^ Ragazzini, JR; Zadeh, LA (1952). 「サンプル値システムの解析」.アメリカ電気学会誌, 第2部: 応用と産業. 71 (5): 225– 234. doi :10.1109/TAI.1952.6371274. S2CID  51674188.
  8. ^ コーネリアス・T・レオンズ (1996). デジタル制御システムの実装と計算技術. アカデミック・プレス. p. 123. ISBN 978-0-12-012779-5
  9. ^ Eliahu Ibrahim Jury (1958).サンプルデータ制御システム. John Wiley & Sons.
  10. ^ Jackson, Leland B. (1996). 「Z変換」.デジタルフィルタと信号処理. ボストン, マサチューセッツ州: Springer US. pp.  29– 54. doi :10.1007/978-1-4757-2458-5_3. ISBN 978-1-4419-5153-3離散時間システムにおけるz変換は、連続時間システムにおけるラプラス変換に相当します。z複素変数です。これは両側z変換と呼ばれることもあり、片側z変換はn = 0から無限大までの総和を除いて同じです。片側変換の主な用途は因果シーケンスであり、この場合、2つの変換はいずれにせよ同じです。したがって、ここではこの区別は行わず、...を単にx ( n )のz変換と呼ぶことにします
  11. ^ Proakis, John; Manolakis, Dimitris.デジタル信号処理の原理、アルゴリズム、およびアプリケーション(第3版). PRENTICE-HALL INTERNATIONAL, INC.
  12. ^ パオロ・ボルツェルン;スカトリーニ、リッカルド。スキアヴォーニ、ニコラ (2015)。Fondamenti di Controlli Automatici (イタリア語)。 MCグラウヒル教育。ISBN 978-88-386-6882-1
  13. ^ abc AR Forouzan (2016). 「Z変換の微分の収束領域」. Electronics Letters . 52 (8): 617– 619. Bibcode :2016ElL....52..617F. doi :10.1049/el.2016.0189. S2CID  124802942.

さらに読む

  • Refaat El Attar, Lecture notes on Z-Transform , Lulu Press, Morrisville NC, 2005. ISBN 1-4116-1979-X
  • 緒方勝彦著『離散時間制御システム 第2版』Prentice-Hall Inc, 1995年、1987年。ISBN 0-13-034281-5
  • Alan V. Oppenheim、Ronald W. Schafer (1999). 離散時間信号処理 第2版、Prentice Hall信号処理シリーズ. ISBN 0-13-754920-2
  • 「Z変換」.数学百科事典. EMS Press . 2001 [1994].
  • Merrikh-Bayat, Farshad (2014). 「Z変換の数値逆変換のための2つの方法」arXiv : 1409.1727 [math.NA].
  • いくつかの一般的なラプラス変換のZ変換表
  • MathworldのZ変換に関する記事(2013年1月30日アーカイブ、Wayback Machine)
  • Comp.DSPのZ変換スレッドは2012年6月15日にWayback Machineにアーカイブされています
  • ラプラス変換のS平面とZ変換のZ平面の関係を示すグラフ
  • エンジニア向けのZ変換のビデオベースの説明
  • z 変換とは何ですか?
Retrieved from "https://en.wikipedia.org/w/index.php?title=Z-transform&oldid=1321489021"