時間領域から周波数領域に信号を変換する数学的変換
数学 と 信号処理 において 、 Z変換は 実数 または 複素数 の 列 である 離散時間信号 を複素数値の 周波数領域 ( Z領域 または Z平面 )表現に変換します 。 [1] [2]
これはラプラス変換 ( s領域 または s平面 )の離散時間版と考えることができる 。 [3]この類似性は 時間スケール計算 の理論で探求されている 。
連続時間フーリエ変換は s 領域の垂直軸 (虚軸) 上で評価されます が、 離散時間フーリエ変換は z 領域の 単位円 に沿って評価されます。s 領域の左 半平面 は z 領域の単位円の内側の領域にマッピングされ、s 領域の右半平面は z 領域の単位円の外側の領域にマッピングされます。
信号処理において、デジタル フィルタ を設計する手段の 1 つは 、アナログ設計をs 領域から z 領域にマッピングする 双線形変換 にかけ、検査、操作、または数値近似によってデジタル フィルタを作成することです。このような方法は、複素数単位の近傍、つまり低周波数以外では正確ではない傾向があります。
歴史
現在ではデジタル制御システムの解析と設計の要となる Z 変換として認識されている基本概念は、20 世紀半ばに登場した当時はまったく新しいものではなかった。その萌芽的な原理はフランスの数学者 ピエール=シモン・ラプラスの研究に遡ることができる。ラプラスは、密接に関連する数学的手法である ラプラス変換 でよく知られている 。しかし、現在 Z 変換として理解されているものの明示的な定式化と応用は、1947 年に Witold Hurewiczと同僚によって大幅に進歩した。彼らの研究は、当時 レーダー 技術の文脈でますます重要になってきていたサンプルデータ制御システムによってもたらされた課題に 動機づけられた。Z 変換は、離散時間信号とシステムの解析で広く使用されている定数係数の線形差分方程式を解くための体系的かつ効果的な方法を提供した。 [4] [5]
この方法は、コロンビア大学のサンプルデータ制御グループの一員であったジョン・R・ラガッツィーニ と ロトフィ・A・ザデー の努力により、1952年にさらに改良され、「Z変換」という公式名称が付けられました 。彼らの研究は、Z変換の数学的枠組みを固めただけでなく、特に電気工学と制御システムの分野でその応用範囲を拡大しました。 [6] [7]
改良Z変換または 高度Z変換 として知られる注目すべき拡張は、後に Eliahu I. Jury によって導入されました。Juryの研究は、特に初期条件の処理とデジタル制御システムの解析のためのより包括的なフレームワークの提供において、Z変換の適用性と堅牢性を拡張しました。この高度な定式化は、離散時間制御システムの設計と安定性の解析において極めて重要な役割を果たし、デジタル信号処理の分野に大きく貢献しました。 [8] [9]
興味深いことに、Z 変換の概念的基礎は、組合せ論と確率論の強力なツールである 生成関数法として知られるより広範な数学的概念と交差しています。このつながりは、確率論の発展における先駆者である アブラハム・ド・モアブル によって 1730 年にはすでに示唆されていました 。ド・モアブルは生成関数を利用して確率の問題を解き、最終的に Z 変換へと発展するものの基礎を築きました。数学的な観点からは、Z 変換は ローラン級数の特定のインスタンスと見なすことができます。ローラン級数 では、調査対象の数の シーケンスが 解析関数 の (ローラン) 展開の 係数 として解釈されます 。この観点は、Z 変換の深い数学的ルーツを強調するだけでなく、数学と工学のさまざまな分野にわたるその汎用性と幅広い適用性を示しています。 [10]
意味
Z変換は片側 変換または 両側 変換として定義できます。( 片側ラプラス変換 と 両側ラプラス変換が あるのと同じです 。) [11]
離散時間信号の両側 Z 変換は、 次 の よう に定義される 正式なべき級数 です 。
x
[
ん
]
{\displaystyle x[n]}
バツ
(
ず
)
{\displaystyle X(z)}
バツ
(
ず
)
=
ず
{
x
[
ん
]
}
=
∑
ん
=
−
∞
∞
x
[
ん
]
ず
−
ん
{\displaystyle X(z)={\mathcal {Z}}\{x[n]\}=\sum _{n=-\infty }^{\infty }x[n]z^{-n}}
ここで は整数で、 は一般に 複素数 です。 極形式 では、 次のように表すことができます。
n
{\displaystyle n}
z
{\displaystyle z}
z
{\displaystyle z}
z
=
A
e
j
ϕ
=
A
⋅
(
cos
ϕ
+
j
sin
ϕ
)
{\displaystyle z=Ae^{j\phi }=A\cdot (\cos {\phi }+j\sin {\phi })}
ここで、 は の大きさ 、 は 虚数単位 、 は ラジアン単位の 複素引数 (角度 または 位相 とも呼ばれる ) です 。
A
{\displaystyle A}
z
{\displaystyle z}
j
{\displaystyle j}
ϕ
{\displaystyle \phi }
あるいは、 がに対してのみ定義されて いる場合 、 片側 Z 変換は 次 のように定義されます。
x
[
n
]
{\displaystyle x[n]}
n
≥
0
{\displaystyle n\geq 0}
X
(
z
)
=
Z
{
x
[
n
]
}
=
∑
n
=
0
∞
x
[
n
]
z
−
n
.
{\displaystyle X(z)={\mathcal {Z}}\{x[n]\}=\sum _{n=0}^{\infty }x[n]z^{-n}.}
信号処理 では、この定義を使用して 、離散時間 因果システムの 単位インパルス応答 の Z 変換を評価することができます 。
一方向 Z 変換の重要な例としては、 確率生成関数 が挙げられます。ここで、 は 離散確率変数が値を取る確率です。 Z 変換の特性 (§ 特性にリストされています) は、確率論の文脈で有用な解釈ができます。
x
[
n
]
{\displaystyle x[n]}
逆 Z 変換は次のとおりです
。
x
[
n
]
=
Z
−
1
{
X
(
z
)
}
=
1
2
π
j
∮
C
X
(
z
)
z
n
−
1
d
z
{\displaystyle x[n]={\mathcal {Z}}^{-1}\{X(z)\}={\frac {1}{2\pi j}}\oint _{C}X(z)z^{n-1}dz}
ここで、 は 原点を囲み、 収束領域 (ROC) 内に完全に含まれる反時計回りの閉じた経路です。ROC が因果関係にある場合 (例 2 を参照)、これは経路 が のすべての極を囲む必要があることを意味します 。
C
{\displaystyle C}
C
{\displaystyle C}
X
(
z
)
{\displaystyle X(z)}
この輪郭積分 の特殊なケースは、 が単位円である ときに発生します 。この輪郭は、ROC が単位円を含むときに使用できます 。これは、 が安定しているとき、つまりすべての極が単位円の内側にあるときに常に保証されます。この輪郭を使用すると、逆 Z 変換は、単位円の周りの Z 変換の周期値の 逆離散時間フーリエ変換 、または フーリエ級数 に簡略化されます。
C
{\displaystyle C}
X
(
z
)
{\displaystyle X(z)}
x
[
n
]
=
1
2
π
∫
−
π
+
π
X
(
e
j
ω
)
e
j
ω
n
d
ω
.
{\displaystyle x[n]={\frac {1}{2\pi }}\int _{-\pi }^{+\pi }X(e^{j\omega })e^{j\omega n}d\omega .}
の有限範囲と有限個の均一に間隔を置いた値を 持つ Z 変換は、 Bluestein の FFT アルゴリズム によって効率的に計算できます 。 離散時間フーリエ変換(DTFT) (離散フーリエ変換 (DFT)と混同しないでください) は、 を単位円上に
制限することによって得られる Z 変換の特殊なケースです。
n
{\displaystyle n}
z
{\displaystyle z}
z
{\displaystyle z}
逆変換の評価には、次の3つの方法がよく使用されます。
輪郭積分による直接評価
この方法では、コーシー留数定理を適用して逆 Z 変換を評価します。複素平面の閉じた輪郭の周りを積分することにより、ROC 内の Z 変換関数の極の留数が合計されます。この手法は、複素変数で表される関数を扱う場合に特に便利です。
変数の項のシリーズへの展開 ず そして ず -1
この方法では、Z 変換はべき級数に展開されます。このアプローチは、Z 変換関数が有理数である場合に役立ち、級数に展開して項ごとに信号係数を決定することで逆関数の近似が可能になります。
部分分数展開と表参照
この手法では、Z 変換を、それぞれが既知の Z 変換ペアに対応する、より単純な分数の合計に分解します。次に、Z 変換ペアの標準テーブルで各項を検索して、逆 Z 変換を決定します。この方法は、その効率性と単純さから広く使用されており、特に元の関数を認識可能なコンポーネントに簡単に分解できる場合によく使用されます。
例: [12]
A) 級数展開法によって次の逆Z変換を決定する。
X
(
z
)
=
1
1
−
1.5
z
−
1
+
0.5
z
−
2
{\displaystyle X(z)={\frac {1}{1-1.5z^{-1}+0.5z^{-2}}}}
解決:
ケース1:
中華人民共和国:
|
Z
|
>
1
{\displaystyle \left\vert Z\right\vert >1}
ROC は円の外部であるため、 因果関係があります (n≥0 の場合に信号が存在します)。
x
(
n
)
{\displaystyle x(n)}
X
(
z
)
=
1
1
−
3
2
z
−
1
+
1
2
z
−
2
=
1
+
3
2
z
−
1
+
7
4
z
−
2
+
15
8
z
−
3
+
31
16
z
−
4
+
.
.
.
.
{\displaystyle X(z)={1 \over 1-{3 \over 2}z^{-1}+{1 \over 2}z^{-2}}=1+{{3 \over 2}z^{-1}}+{{7 \over 4}z^{-2}}+{{15 \over 8}z^{-3}}+{{31 \over 16}z^{-4}}+....}
したがって、
x
(
n
)
=
{
1
,
3
2
,
7
4
,
15
8
,
31
16
…
}
↑
{\displaystyle {\begin{aligned}x(n)&=\left\{1,{\frac {3}{2}},{\frac {7}{4}},{\frac {15}{8}},{\frac {31}{16}}\ldots \right\}\\&\qquad \!\uparrow \\\end{aligned}}}
(矢印はx(0)=1の項を示す)
長除算プロセスの各ステップで、 の最小のべき乗項を削除することに注意してください 。
z
−
1
{\displaystyle z^{-1}}
ケース2:
中華人民共和国:
|
Z
|
<
0.5
{\displaystyle \left\vert Z\right\vert <0.5}
ROC は円の内部であるため、 反因果的です (n<0 の場合に信号が存在します)。
x
(
n
)
{\displaystyle x(n)}
長除算を実行すると、
X
(
z
)
=
1
1
−
3
2
z
−
1
+
1
2
z
−
2
=
2
z
2
+
6
z
3
+
14
z
4
+
30
z
5
+
…
{\displaystyle X(z)={\frac {1}{1-{\frac {3}{2}}z^{-1}+{\frac {1}{2}}z^{-2}}}=2z^{2}+6z^{3}+14z^{4}+30z^{5}+\ldots }
x
(
n
)
=
{
30
,
14
,
6
,
2
,
0
,
0
}
↑
{\displaystyle {\begin{aligned}x(n)&=\{30,14,6,2,0,0\}\\&\qquad \qquad \qquad \quad \ \ \,\uparrow \\\end{aligned}}}
(矢印はx(0)=0の項を示す)
長除算プロセスの各ステップで、 の最小のべき乗項を削除することに注意してください 。
z
{\displaystyle z}
注記:
信号が因果関係にある場合、 の正の累乗が得られ 、信号が反因果関係にある場合、 の負の累乗が得られます 。
z
{\displaystyle z}
z
{\displaystyle z}
z
k
{\displaystyle z^{k}}
は の期間を示し 、 は の期間を示します 。
x
(
−
k
)
{\displaystyle x(-k)}
z
−
k
{\displaystyle z^{-k}}
x
(
k
)
{\displaystyle x(k)}
B) 級数展開法によって次の逆Z変換を決定する。
負の累乗を消去し 、を で割ると 、
z
{\displaystyle z}
z
{\displaystyle z}
X
(
z
)
z
=
z
2
z
(
z
2
−
1.5
z
+
0.5
)
=
z
z
2
−
1.5
z
+
0.5
{\displaystyle {\frac {X(z)}{z}}={\frac {z^{2}}{z(z^{2}-1.5z+0.5)}}={\frac {z}{z^{2}-1.5z+0.5}}}
部分分数展開により、
X
(
z
)
z
=
z
(
z
−
1
)
(
z
−
0.5
)
=
A
1
z
−
0.5
+
A
2
z
−
1
A
1
=
(
z
−
0.5
)
X
(
z
)
z
|
z
=
0.5
=
0.5
(
0.5
−
1
)
=
−
1
A
2
=
(
z
−
1
)
X
(
z
)
z
|
z
=
1
=
1
1
−
0.5
=
2
X
(
z
)
z
=
2
z
−
1
−
1
z
−
0.5
{\displaystyle {\begin{aligned}{\frac {X(z)}{z}}&={\frac {z}{(z-1)(z-0.5)}}={\frac {A_{1}}{z-0.5}}+{\frac {A_{2}}{z-1}}\\[4pt]&A_{1}=\left.{\frac {(z-0.5)X(z)}{z}}\right\vert _{z=0.5}={\frac {0.5}{(0.5-1)}}=-1\\[4pt]&A_{2}=\left.{\frac {(z-1)X(z)}{z}}\right\vert _{z=1}={\frac {1}{1-0.5}}={2}\\[4pt]{\frac {X(z)}{z}}&={\frac {2}{z-1}}-{\frac {1}{z-0.5}}\end{aligned}}}
ケース1:
中華人民共和国:
|
Z
|
>
1
{\displaystyle \left\vert Z\right\vert >1}
どちらの用語も因果関係があり、したがって 因果関係があります。
x
(
n
)
{\displaystyle x(n)}
x
(
n
)
=
2
(
1
)
n
u
(
n
)
−
1
(
0.5
)
n
u
(
n
)
=
(
2
−
0.5
n
)
u
(
n
)
{\displaystyle {\begin{aligned}x(n)&=2{(1)^{n}}u(n)-1{(0.5)^{n}}u(n)\\&=(2-0.5^{n})u(n)\\\end{aligned}}}
ケース2:
中華人民共和国:
|
Z
|
<
0.5
{\displaystyle \left\vert Z\right\vert <0.5}
どちらの用語も反因果的であるため、 は 反因果的です。
x
(
n
)
{\displaystyle x(n)}
x
(
n
)
=
−
2
(
1
)
n
u
(
−
n
−
1
)
−
(
−
1
(
0.5
)
n
u
(
−
n
−
1
)
)
=
(
0.5
n
−
2
)
u
(
−
n
−
1
)
{\displaystyle {\begin{aligned}x(n)&=-2{(1)^{n}}u(-n-1)-(-1{(0.5)^{n}}u(-n-1))\\&=(0.5^{n}-2)u(-n-1)\\\end{aligned}}}
ケース3:
中華人民共和国:
0.5
<
|
Z
|
<
1
{\displaystyle 0.5<\left\vert Z\right\vert <1}
項の 1 つは因果関係があり (p=0.5 は因果関係の部分を提供)、もう 1 つは反因果関係があり (p=1 は反因果関係の部分を提供)、したがって 両側性です。
x
(
n
)
{\displaystyle x(n)}
x
(
n
)
=
−
2
(
1
)
n
u
(
−
n
−
1
)
−
1
(
0.5
)
n
u
(
n
)
=
−
2
u
(
−
n
−
1
)
−
0.5
n
u
(
n
)
{\displaystyle {\begin{aligned}x(n)&=-2{(1)^{n}}u(-n-1)-1{(0.5)^{n}}u(n)\\&=-2u(-n-1)-0.5^{n}u(n)\\\end{aligned}}}
収束領域
収束領域 ( ROC) は、Z 変換の合計 が収束する (つまり、大きさが無限大に爆発しない) 複素平面上の点の集合です。
R
O
C
=
{
z
:
|
∑
n
=
−
∞
∞
x
[
n
]
z
−
n
|
<
∞
}
{\displaystyle \mathrm {ROC} =\left\{z:\left|\sum _{n=-\infty }^{\infty }x[n]z^{-n}\right|<\infty \right\}}
例1(ROCなし)
間隔 を
広げる と
x
[
n
]
=
(
.5
)
n
.
{\displaystyle x[n]=(.5)^{n}\ .}
x
[
n
]
{\displaystyle x[n]}
(
−
∞
,
∞
)
{\displaystyle (-\infty ,\infty )}
x
[
n
]
=
{
…
,
(
.5
)
−
3
,
(
.5
)
−
2
,
(
.5
)
−
1
,
1
,
(
.5
)
,
(
.5
)
2
,
(
.5
)
3
,
…
}
=
{
…
,
2
3
,
2
2
,
2
,
1
,
(
.5
)
,
(
.5
)
2
,
(
.5
)
3
,
…
}
.
{\displaystyle x[n]=\left\{\dots ,(.5)^{-3},(.5)^{-2},(.5)^{-1},1,(.5),(.5)^{2},(.5)^{3},\dots \right\}=\left\{\dots ,2^{3},2^{2},2,1,(.5),(.5)^{2},(.5)^{3},\dots \right\}.}
合計を見ると
∑
n
=
−
∞
∞
x
[
n
]
z
−
n
→
∞
.
{\displaystyle \sum _{n=-\infty }^{\infty }x[n]z^{-n}\to \infty .}
したがって、この条件を満たす
の値は存在しません。
z
{\displaystyle z}
例2(因果ROC)
ROC (青)、| z | = .5 (破線の黒い円)、および単位円 (点線の灰色の円)。
(ここで は ヘヴィ サイドのステップ関数 )とする 。 区間を展開すると 、
x
[
n
]
=
(
.5
)
n
u
[
n
]
{\displaystyle x[n]=(.5)^{n}\,u[n]}
u
{\displaystyle u}
x
[
n
]
{\displaystyle x[n]}
(
−
∞
,
∞
)
{\displaystyle (-\infty ,\infty )}
x
[
n
]
=
{
…
,
0
,
0
,
0
,
1
,
(
.5
)
,
(
.5
)
2
,
(
.5
)
3
,
…
}
.
{\displaystyle x[n]=\left\{\dots ,0,0,0,1,(.5),(.5)^{2},(.5)^{3},\dots \right\}.}
合計を見ると
∑
n
=
−
∞
∞
x
[
n
]
z
−
n
=
∑
n
=
0
∞
(
.5
)
n
z
−
n
=
∑
n
=
0
∞
(
.5
z
)
n
=
1
1
−
(
.5
)
z
−
1
.
{\displaystyle \sum _{n=-\infty }^{\infty }x[n]z^{-n}=\sum _{n=0}^{\infty }(.5)^{n}z^{-n}=\sum _{n=0}^{\infty }\left({\frac {.5}{z}}\right)^{n}={\frac {1}{1-(.5)z^{-1}}}.}
最後の等式は無限等 比級数 から生じ、等式は の場合にのみ成立し、 は として書き直すこと が できます。 したがって、ROC は この場合、ROC は原点に半径 0.5 の円盤が「打ち抜かれた」複素平面です。
|
(
.5
)
z
−
1
|
<
1
,
{\displaystyle |(.5)z^{-1}|<1,}
z
{\displaystyle z}
|
z
|
>
(
.5
)
.
{\displaystyle |z|>(.5).}
|
z
|
>
(
.5
)
.
{\displaystyle |z|>(.5).}
例3(反因果ROC)
ROC (青)、| z | = .5 (破線の黒い円)、および単位円 (点線の灰色の円)。
(ここで は ヘヴィ サイドのステップ関数 )とする 。 区間を展開すると 、
x
[
n
]
=
−
(
.5
)
n
u
[
−
n
−
1
]
{\displaystyle x[n]=-(.5)^{n}\,u[-n-1]}
u
{\displaystyle u}
x
[
n
]
{\displaystyle x[n]}
(
−
∞
,
∞
)
{\displaystyle (-\infty ,\infty )}
x
[
n
]
=
{
…
,
−
(
.5
)
−
3
,
−
(
.5
)
−
2
,
−
(
.5
)
−
1
,
0
,
0
,
0
,
0
,
…
}
.
{\displaystyle x[n]=\left\{\dots ,-(.5)^{-3},-(.5)^{-2},-(.5)^{-1},0,0,0,0,\dots \right\}.}
合計を見ると
∑
n
=
−
∞
∞
x
[
n
]
z
−
n
=
−
∑
n
=
−
∞
−
1
(
.5
)
n
z
−
n
=
−
∑
m
=
1
∞
(
z
.5
)
m
=
−
(
.5
)
−
1
z
1
−
(
.5
)
−
1
z
=
−
1
(
.5
)
z
−
1
−
1
=
1
1
−
(
.5
)
z
−
1
{\displaystyle {\begin{aligned}\sum _{n=-\infty }^{\infty }x[n]\,z^{-n}&=-\sum _{n=-\infty }^{-1}(.5)^{n}\,z^{-n}\\&=-\sum _{m=1}^{\infty }\left({\frac {z}{.5}}\right)^{m}\\&=-{\frac {(.5)^{-1}z}{1-(.5)^{-1}z}}\\&=-{\frac {1}{(.5)z^{-1}-1}}\\&={\frac {1}{1-(.5)z^{-1}}}\\\end{aligned}}}
そして、再び無限等 比級数 を使用すると、等式は の場合にのみ成立し、 は として 書き直すことができます 。したがって、ROC は この場合、ROC は原点を中心とし、半径 0.5 の円盤です。
|
(
.5
)
−
1
z
|
<
1
{\displaystyle |(.5)^{-1}z|<1}
z
{\displaystyle z}
|
z
|
<
(
.5
)
.
{\displaystyle |z|<(.5).}
|
z
|
<
(
.5
)
.
{\displaystyle |z|<(.5).}
この例と前の例の違いは、 ROC
だけです。これは、変換結果だけでは不十分であることを示すために意図的に行われています。
例の結論
例 2 と 3 は、ROC を指定した場合のみ、 Z 変換 が一意であることを明確に示しています。因果関係のあるケースと反因果関係のあるケースの 極ゼロ プロット を作成すると、どちらのケースの ROC にも 0.5 の極が含まれていないことがわかります。これは、複数の極がある場合にも当てはまります。ROC に極が含まれることは ありません 。
X
(
z
)
{\displaystyle X(z)}
x
[
n
]
{\displaystyle x[n]}
例2では、因果関係のシステムではROCが生成され、 例3では反因果関係のシステムではROCが生成され、
|
z
|
=
∞
{\displaystyle |z|=\infty }
|
z
|
=
0.
{\displaystyle |z|=0.}
ROCは青いリングで表示されます 0.5 < | z | < 0.75
複数の極を持つシステムでは、ROCに極も 極も含まれない可能性がある 。ROCは円形のバンドを作成する。例えば、
|
z
|
=
∞
{\displaystyle |z|=\infty }
|
z
|
=
0.
{\displaystyle |z|=0.}
x
[
n
]
=
(
.5
)
n
u
[
n
]
−
(
.75
)
n
u
[
−
n
−
1
]
{\displaystyle x[n]=(.5)^{n}\,u[n]-(.75)^{n}\,u[-n-1]}
0.5と0.75に極があります。ROCは0.5 < | z | < 0.75となり、原点も無限大も含まれません。このようなシステムは因果項 と反因果項 を含むため、混合因果システムと呼ばれます。
(
.5
)
n
u
[
n
]
{\displaystyle (.5)^{n}\,u[n]}
−
(
.75
)
n
u
[
−
n
−
1
]
.
{\displaystyle -(.75)^{n}\,u[-n-1].}
システムの安定 性は 、ROC のみを知ることによっても判断できます。ROC に単位円が含まれている場合 (つまり、| z | = 1)、システムは安定しています。上記のシステムでは、因果システム (例 2) は、| z | > 0.5 に単位円が含まれているため安定しています。
ROC のないシステムの Z 変換 (つまり、あいまいな ) が提供されていると仮定します。 次の条件を満たす場合に、
一意の を決定できます。
x
[
n
]
{\displaystyle x[n]}
x
[
n
]
{\displaystyle x[n]}
安定性のためには、ROC に単位円が含まれている必要があります。因果システムが必要な場合は、ROC に無限大が含まれている必要があり、システム関数は右側のシーケンスになります。反因果システムが必要な場合は、ROC に原点が含まれている必要があり、システム関数は左側のシーケンスになります。安定性と因果性の両方が必要な場合は、システム関数のすべての極が単位円の内側にある必要があります。
そうすれば、ユニークなもの が見つかります。
x
[
n
]
{\displaystyle x[n]}
プロパティ
パーセバルの定理
∑
n
=
−
∞
∞
x
1
[
n
]
x
2
∗
[
n
]
=
1
j
2
π
∮
C
X
1
(
v
)
X
2
∗
(
1
v
∗
)
v
−
1
d
v
{\displaystyle \sum _{n=-\infty }^{\infty }x_{1}[n]x_{2}^{*}[n]\quad =\quad {\frac {1}{j2\pi }}\oint _{C}X_{1}(v)X_{2}^{*}({\tfrac {1}{v^{*}}})v^{-1}\mathrm {d} v}
初期値定理 : が 因果関係にある場合、
x
[
n
]
{\displaystyle x[n]}
x
[
0
]
=
lim
z
→
∞
X
(
z
)
.
{\displaystyle x[0]=\lim _{z\to \infty }X(z).}
最終値定理 : の極が 単位円の内側にある場合、
(
z
−
1
)
X
(
z
)
{\displaystyle (z-1)X(z)}
x
[
∞
]
=
lim
z
→
1
(
z
−
1
)
X
(
z
)
.
{\displaystyle x[\infty ]=\lim _{z\to 1}(z-1)X(z).}
ここ:
u
:
n
↦
u
[
n
]
=
{
1
,
n
≥
0
0
,
n
<
0
{\displaystyle u:n\mapsto u[n]={\begin{cases}1,&n\geq 0\\0,&n<0\end{cases}}}
は単位(またはヘヴィサイド)ステップ関数 であり 、
δ
:
n
↦
δ
[
n
]
=
{
1
,
n
=
0
0
,
n
≠
0
{\displaystyle \delta :n\mapsto \delta [n]={\begin{cases}1,&n=0\\0,&n\neq 0\end{cases}}}
は離散時間単位インパルス関数 です ( 連続時間版の
ディラックデルタ関数を参照)。2 つの関数は、単位ステップ関数が単位インパルス関数の累積(実行合計)になるように一緒に選択されます。
単位円 として知られる 領域 内 の の値に対して 、変換を 1 つの実変数の関数として 定義して表すことができます 。また、この双方向変換は フーリエ級数 に簡約されます。
z
{\displaystyle z}
|
z
|
=
1
{\displaystyle |z|{=}1}
ω
{\displaystyle \omega }
z
=
e
j
ω
.
{\displaystyle z{=}e^{j\omega }.}
これは、シーケンスの離散時間フーリエ変換 (DTFT) としても知られています 。この -周期関数は、 フーリエ変換 の 周期的な合計 であり、広く使用されている分析ツールとなっています。これを理解するために、 を、ある間隔でのサンプルがシーケンス に等しい任意 の関数 のフーリエ変換とします 。すると、シーケンスの DTFT は 次のように記述できます。
x
[
n
]
{\displaystyle x[n]}
2
π
{\displaystyle 2\pi }
X
(
f
)
{\displaystyle X(f)}
x
(
t
)
{\displaystyle x(t)}
T
{\displaystyle T}
x
[
n
]
{\displaystyle x[n]}
x
[
n
]
{\displaystyle x[n]}
ここで、 は秒単位、は ヘルツ 単位です 。2 つの系列を比較すると、 は サンプルあたりラジアン 単位の 正規化された周波数 であることがわかります 。値 は に対応します 。そして、 を代入すると、 式 1 は (フーリエ変換)
で表すことができます。
T
{\displaystyle T}
f
{\displaystyle f}
ω
=
2
π
f
T
{\displaystyle \omega {=}2\pi fT}
ω
=
2
π
{\displaystyle \omega {=}2\pi }
f
=
1
T
{\textstyle f{=}{\frac {1}{T}}}
f
=
ω
2
π
T
,
{\textstyle f{=}{\frac {\omega }{2\pi T}},}
X
(
ω
−
2
π
k
2
π
T
)
{\displaystyle X({\tfrac {\omega -2\pi k}{2\pi T}})}
パラメータ T が変化すると、 式 2 の個々の項は f 軸に沿って離れたり近づいたりします 。ただし、 式 3 では、中心は 2 π 離れたままですが、幅は拡大または縮小します。シーケンスが LTI システム の インパルス応答 を表す場合、これらの関数は 周波数応答 とも呼ばれます 。 シーケンスが周期的である場合、その DTFT は 1 つ以上の高調波周波数で発散し、他のすべての周波数では 0 になります。これは、高調波周波数で振幅可変の ディラック デルタ関数を使用することでよく表されます。周期性のため、一意の振幅は有限の数しか存在せず、これははるかに単純な 離散フーリエ変換 (DFT)によって簡単に計算できます。( 離散時間フーリエ変換 § 周期データ を参照してください 。)
x
(
n
T
)
{\displaystyle x(nT)}
x
(
n
T
)
{\displaystyle x(nT)}
双 線形変換を 使用すると、連続時間フィルタ (ラプラス領域で表現) を離散時間フィルタ (Z 領域で表現) に変換したり、その逆を行ったりすることができます。次の置換が使用されます。
s
=
2
T
(
z
−
1
)
(
z
+
1
)
{\displaystyle s={\frac {2}{T}}{\frac {(z-1)}{(z+1)}}}
ラプラス領域の 関数を Z領域の関数に変換する( タスティン変換 )、または
H
(
s
)
{\displaystyle H(s)}
H
(
z
)
{\displaystyle H(z)}
z
=
e
s
T
≈
1
+
s
T
/
2
1
−
s
T
/
2
{\displaystyle z=e^{sT}\approx {\frac {1+sT/2}{1-sT/2}}}
Z 領域からラプラス領域へ。双線形変換により、複素 s 平面 (ラプラス変換) は複素 z 平面 (z 変換) にマッピングされます。このマッピングは (必然的に) 非線形ですが、 s 平面の軸全体 を z 平面の 単位円 にマッピングするという点で便利です。このように、フーリエ変換 ( 軸上で評価されるラプラス変換) は離散時間フーリエ変換になります。これは、フーリエ変換が存在すること、つまり軸が ラプラス変換の収束領域にあることを前提としています。
j
ω
{\displaystyle j\omega }
j
ω
{\displaystyle j\omega }
j
ω
{\displaystyle j\omega }
時間サンプリングされた関数の 片側 Z 変換が与えられると、対応する 星印付きの変換 によってラプラス変換が生成され、 (サンプリング パラメータ) への依存性が復元されます。
X
(
z
)
{\displaystyle X(z)}
T
{\displaystyle T}
X
∗
(
s
)
=
X
(
z
)
|
z
=
e
s
T
{\displaystyle {\bigg .}X^{*}(s)=X(z){\bigg |}_{\displaystyle z=e^{sT}}}
逆ラプラス変換は、 インパルスサンプリング 関数として知られる数学的な抽象化です。
線形定数係数差分方程式
線形定係数差 (LCCD) 方程式は、自己回帰移動平均 方程式に基づく線形システムの表現です 。
∑
p
=
0
N
y
[
n
−
p
]
α
p
=
∑
q
=
0
M
x
[
n
−
q
]
β
q
.
{\displaystyle \sum _{p=0}^{N}y[n-p]\alpha _{p}=\sum _{q=0}^{M}x[n-q]\beta _{q}.}
上の式の両辺は ゼロでなければ で割ることができる。LCCDで正規化すると、次の 式が書ける。
α
0
{\displaystyle \alpha _{0}}
α
0
=
1
,
{\displaystyle \alpha _{0}{=}1,}
y
[
n
]
=
∑
q
=
0
M
x
[
n
−
q
]
β
q
−
∑
p
=
1
N
y
[
n
−
p
]
α
p
.
{\displaystyle y[n]=\sum _{q=0}^{M}x[n-q]\beta _{q}-\sum _{p=1}^{N}y[n-p]\alpha _{p}.}
この形式のLCCD方程式は、「現在の」出力が過去の出力、 現在の入力 、および以前の入力 の関数であることをより明確にするために有利である。
y
[
n
]
{\displaystyle y[n]}
y
[
n
−
p
]
,
{\displaystyle y[n-p],}
x
[
n
]
,
{\displaystyle x[n],}
x
[
n
−
q
]
.
{\displaystyle x[n-q].}
伝達関数
上記の式を Z 変換すると (線形性と時間シフトの法則を使用)、次のようになります。
Y
(
z
)
∑
p
=
0
N
z
−
p
α
p
=
X
(
z
)
∑
q
=
0
M
z
−
q
β
q
{\displaystyle Y(z)\sum _{p=0}^{N}z^{-p}\alpha _{p}=X(z)\sum _{q=0}^{M}z^{-q}\beta _{q}}
ここで 、およびはそれぞれ 、および の z 変換です 。(表記規則では通常、大文字を使用して、対応する小文字で表された信号の z 変換を参照します。これは、ラプラス変換の表記規則に似ています。)
X
(
z
)
{\displaystyle X(z)}
Y
(
z
)
{\displaystyle Y(z)}
x
[
n
]
{\displaystyle x[n]}
y
[
n
]
,
{\displaystyle y[n],}
システムの 伝達関数 を整理すると次のようになります。
H
(
z
)
=
Y
(
z
)
X
(
z
)
=
∑
q
=
0
M
z
−
q
β
q
∑
p
=
0
N
z
−
p
α
p
=
β
0
+
z
−
1
β
1
+
z
−
2
β
2
+
⋯
+
z
−
M
β
M
α
0
+
z
−
1
α
1
+
z
−
2
α
2
+
⋯
+
z
−
N
α
N
.
{\displaystyle H(z)={\frac {Y(z)}{X(z)}}={\frac {\sum _{q=0}^{M}z^{-q}\beta _{q}}{\sum _{p=0}^{N}z^{-p}\alpha _{p}}}={\frac {\beta _{0}+z^{-1}\beta _{1}+z^{-2}\beta _{2}+\cdots +z^{-M}\beta _{M}}{\alpha _{0}+z^{-1}\alpha _{1}+z^{-2}\alpha _{2}+\cdots +z^{-N}\alpha _{N}}}.}
ゼロと極
代数学の基本定理 から、 分子 には 根 (の零点に対応 ) があり、 分母に は根(極に対応) がある。 伝達関数を 零点と極の 観点から書き直すと、
M
{\displaystyle M}
H
{\displaystyle H}
N
{\displaystyle N}
H
(
z
)
=
(
1
−
q
1
z
−
1
)
(
1
−
q
2
z
−
1
)
⋯
(
1
−
q
M
z
−
1
)
(
1
−
p
1
z
−
1
)
(
1
−
p
2
z
−
1
)
⋯
(
1
−
p
N
z
−
1
)
,
{\displaystyle H(z)={\frac {(1-q_{1}z^{-1})(1-q_{2}z^{-1})\cdots (1-q_{M}z^{-1})}{(1-p_{1}z^{-1})(1-p_{2}z^{-1})\cdots (1-p_{N}z^{-1})}},}
ここで 、 は 零点、 は極 です 。零点と極は一般に複素数であり、複素平面 (z 平面) 上にプロットすると 極零点プロット と呼ばれます。
q
k
{\displaystyle q_{k}}
k
th
{\displaystyle k^{\text{th}}}
p
k
{\displaystyle p_{k}}
k
th
{\displaystyle k^{\text{th}}}
さらに、 および にも零点と極が存在する可能性があります 。 これらの極と零点、および多次の零点と極を考慮すると、零点と極の数は常に等しくなります。
z
=
0
{\displaystyle z{=}0}
z
=
∞
.
{\displaystyle z{=}\infty .}
分母を因数分解することで 部分分数 分解が可能になり、それを時間領域に戻すことができます。これにより、 システムの
インパルス応答と線形定数係数差分方程式が得られます。
出力応答
このようなシステムが 信号によって駆動される場合 、出力は の 部分分数 分解 を実行して から逆 Z 変換を行うことで、出力を見つけることができます。実際には、 その量を で乗算する前に 分数分解して、簡単に計算できる逆 Z 変換を持つ項を持つ の形式を生成することが役立つことがよく あります。
H
(
z
)
{\displaystyle H(z)}
X
(
z
)
{\displaystyle X(z)}
Y
(
z
)
=
H
(
z
)
X
(
z
)
.
{\displaystyle Y(z)=H(z)X(z).}
Y
(
z
)
{\displaystyle Y(z)}
y
[
n
]
{\displaystyle y[n]}
Y
(
z
)
z
{\displaystyle \textstyle {\frac {Y(z)}{z}}}
z
{\displaystyle z}
Y
(
z
)
{\displaystyle Y(z)}
参照
参考文献
^ Mandal, Jyotsna Kumar (2020). 「 Z 変換ベースの可逆エンコーディング」。 変換エンコーディングによる可逆ステガノグラフィーと認証 。 計算知能研究。 Vol. 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 変換はラプラス変換から派生します。
^ 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 変換は離散信号やシステムを扱うのに特に適しています。離散時間フーリエ変換よりもコンパクトで便利な表記法を提供します。
^ Palani, S. (2021-08-26). 「 離散時間信号とシステムの z変換分析」. 信号とシステム . チャム: Springer International Publishing. pp. 921–1055. doi :10.1007/978-3-030-75742-7_9. ISBN 978-3-030-75741-0 . S2CID 238692483。z 変換は、ラプラス変換の離散版です。z変換は 、 離散時間システムの差分方程式を代数方程式に変換し、離散時間システム解析を簡素化します。ラプラス変換と z 変換 は共通ですが、ラプラス変換は連続時間信号とシステムを扱う点が異なります。
^ ER カナセウィッチ (1981)。地球物理学における時系列解析。アルバータ大学。 186、249ページ 。ISBN
978-0-88864-074-1 。
^ ER Kanasewich (1981). 地球物理学における時系列分析(第3版). アルバータ大学. pp. 185–186. ISBN 978-0-88864-074-1 。
^ Ragazzini, JR; Zadeh, LA (1952). 「サンプルデータシステムの分析」. 米国電気学会誌、第2部: 応用と産業 . 71 (5): 225–234. doi :10.1109/TAI.1952.6371274. S2CID 51674188.
^ Cornelius T. Leondes (1996). デジタル制御システムの実装と計算技術。Academic Press。p. 123。ISBN 978-0-12-012779-5 。
^ Eliahu Ibrahim Jury (1958)。 サンプルデータ制御システム 。John Wiley & Sons。
^ Eliahu Ibrahim Jury (1973)。Z 変換法の理論と応用 。Krieger Pub Co. ISBN
0-88275-122-0 。
^ Eliahu Ibrahim Jury (1964)。Z 変換法の理論と応用 。John Wiley & Sons。p. 1。
^ 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 変換と呼ぶことにします。
^ Proakis, John; Manolakis, Dimitris. デジタル信号処理の原理、アルゴリズム、アプリケーション (第 3 版)。PRENTICE-HALL INTERNATIONAL, INC.
^ ボルツェルン、パオロ;スカトリーニ、リッカルド。スキアヴォーニ、ニコラ (2015)。 Fondamenti di Controlli Automatici (イタリア語)。 MCグラウヒル教育。 ISBN 978-88-386-6882-1 。
^ 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、 「Z変換に関する講義ノート」 、Lulu Press、ノースカロライナ州モリスビル、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変換に関する記事
Comp.DSP の Z 変換スレッド
ラプラス変換のs平面とZ変換のZ平面の関係を示すグラフ
エンジニア向けのZ変換のビデオによる説明
z 変換とは何ですか?