Z変換

提供: miniwiki
2016/3/21/ (月) 20:58時点におけるja>Japanbirdによる版 ({{デジタル信号処理}})
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
移動先:案内検索

関数解析学において、Z変換(ゼッドへんかん、Z-transform)とは、離散群上で定義される、ローラン展開をベースにした関数空間の間の線形作用素。関数変換。

Z変換は離散群上でのラプラス変換とも説明される。なお、Z変換という呼び方は、ラプラス変換のことを「S変換」と呼んでいるようなものであり、定義式中の遅延要素であるzに由来する名前である。

定義

xnのZ変換は以下の式で定義される:

[math] \mathcal{Z}[x_n] = X(z) = \sum_{n=-\infty}^{\infty} x_n z^{-n}[/math]

ここでn整数z複素数である。なお後述の片側Z変換に対してこれを両側Z変換(two-sided Z-transform、bilateral Z-transform)と呼ばれる。

n<0 でxn=0のような場合は、総和の範囲を 0 〜 ∞ で計算できる:

[math] \mathcal{Z}[x_n] = X(z) = \sum_{n=0}^{\infty} x_n z^{-n}[/math]

これを元の定義と区別して片側Z変換(single-sided Z-transform、unilateral Z-transform)と呼ぶこともある。工学の分野などでは因果律を想定するので、こちらの式で定義することがある。

二次元信号(例えば画像)に対する二次元Z変換の定義は類似的である:

[math]\mathcal{Z}[x(n_1,n_2)] = X(z_1,z_2)=\sum_{n_1=-\infty}^{\infty}\sum_{n_2=-\infty}^{\infty} x(n_1,n_2) z_1^{-n_1}z_2^{-n_2}[/math]

収束領域

なお、Z変換の級数は一般には発散することがある。収束するzの領域(収束領域,Region of Convergence)を以下のように書ける:

[math]\mbox{ROC}=\left\{z : \sum_{n=-\infty}^{\infty}x_nz^{-n} \lt \infty\right\}[/math]

厳密にはこの収束領域内においてのX(z)を、xnのZ変換と定義する。

二次元Z変換の収束領域の定義は類似する:

[math]\mbox{ROC}=\left\{(z_1,z_2) : \sum_{n_1=-\infty}^{\infty}\sum_{n_2=-\infty}^{\infty} x(n_1,n_2) z_1^{-n_1}z_2^{-n_2} \lt \infty\right\}[/math]

逆Z変換

Z変換の逆変換である逆Z変換(inverse Z-transform)は次のようになる:

[math]x_n = \mathcal{Z}^{-1}[X(z)]=\frac{1}{2\pi i}\oint_C X(z)z^{n-1}\,dz[/math]

ここでi虚数単位で積分路CX(z)の極を全て含むような閉路である。

なおこの式は留数定理を用いて留数の和として計算することができる。しかし、手計算で計算するときは以下の方法がよく使われる:

  • X(z)が既に級数展開されている場合、z-kの係数をxkの値とすることで簡単に逆変換ができる。例えば、z+2-3z-1の逆変換は { ..., 0, x-1=1,x0=2,x1=-3, 0, ...} のように係数をならべるだけで得られる。
  • X(z)を部分分数分解し、各々の部分分数を変換表を用いて逆変換したものの和として逆変換を得る。

いずれにせよ、定義に示した積分計算そのものを直接計算することは稀である。

性質

線型性
Z変換は線型性を持ち、したがって特に重ね合わせの原理を用いて計算できる。したがって任意のxn,ynに対して
[math]\mathcal{Z}[a x_n + b y_n] = a \mathcal{Z}[x_n] + b \mathcal{Z}[y_n][/math]
が成立する。但し、a,bは定数。逆Z変換も同様に線型性を持つ。したがって、与えられた関数を部分分数分解できるとき、各因子が変換表にあるものに合致すれば、その変換が求められる。
シフト性
[math]\mathcal{Z}[x_{n-k}] = z^{-k}\mathcal{Z}[x_n][/math]
Z領域微分
[math]\mathcal{Z}[n x_n] = -z \frac{d}{dz} \mathcal{Z}[x_n][/math]
畳み込み
フーリエ変換のように畳み込み定理が成り立ち、畳み込みはZ変換によって積となる。
[math] \mathcal{Z}[x_n * y_n] = \mathcal{Z}[x_n] \mathcal{Z}[y_n][/math]
初期値定理
[math]f(0)=\lim_{z\to\infty}F(z)[/math]
最終値定理
[math]f(\infty)=\lim_{k\to\infty}f(k)=\lim_{z\to 1}\left\{\frac{z-1}{z}F(z)\right\}[/math]
時間領域の乗積
[math]\mathcal{Z}[x_n h_n]=\frac{1}{2\pi i}\oint_{C_1} X(v)H\left(\frac{z}{v}\right)v^{-1}\,dv=\frac{1}{2\pi i}\oint_{C_2} H(v)X\left(\frac{z}{v}\right)v^{-1}\,dv[/math]

積分路 [math]C_1[/math][math]X(v)[/math][math]H\left(\frac{z}{v}\right)[/math]ROCの共同区域にある閉路であり、 [math]C_2[/math][math]H(v)[/math][math]X\left(\frac{z}{v}\right)[/math]ROCの共同区域にある閉路である。

Parseval定理
[math]\mathcal{Z}\left[\sum_{n=-\infty}^{+\infty}x_nh_n^*\right]=\frac{1}{2\pi i}\oint_{C_1} X(v)H^*\left(\frac{1}{v^*}\right)v^{-1}\,dv=\frac{1}{2\pi i}\oint_{C_2} H^*(v)X\left(\frac{1}{v}\right)v^{-1}\,dv[/math]

積分路 [math]C_1[/math][math]X(v)[/math][math]H^*\left(\frac{1}{v^*}\right)[/math]ROCの共同区域にある閉路であり、 [math]C_2[/math][math]H^*(v)[/math][math]X\left(\frac{1}{v}\right)[/math]ROCの共同区域にある閉路である。

離散時間のLTIシステム

離散時間のLTIシステムは以下の定数係数の線形差分方程式としてモデル化されることができる:

[math]\sum_{i=0}^{N}a_{i}y(n-i)=\sum_{j=0}^{M}b_{j}x(n-j)[/math]

一般には、[math]a_0=1[/math]と認める。

方程式の両辺をZ変換すると、

[math]Y(z)\sum_{i=0}^{N}a_{i}z^{-i} = X(z)\sum_{j=0}^{M}b_{j}z^{-j}[/math]

を得られて、

[math]H(z)=\frac{Y(z)}{X(z)}=\frac{\displaystyle \sum_{j=0}^{M}b_{j}z^{-j}}{\displaystyle \sum_{i=0}^{N}a_{i}z^{-i}}[/math]

は、伝達関数と呼ばれ、その分母多項式は特性多項式と呼ばれる。

伝達関数を分析すれば、システム特性の解明に役立つ。

他の変換との関係性

ラプラス変換との関係

Z変換は両側ラプラス変換を離散化したものである。つまり離散化された関数

[math]f(t)\sum_{n=-\infty}^\infty \delta(t-nT)=\sum_{n=-\infty}^\infty f_n \delta(t-nT)[/math]

のラプラス変換

[math]\sum_{n=0}^{\infty} f_n \delta(t-nT) e^{-snT}[/math]

に対応する。但し、Tはサンプリング周期であり、esTがZ変換におけるzに対応する。

離散時間フーリエ変換との関係

Z変換は離散時間フーリエ変換(DTFT)の拡張である。DTFTはZ変換でz=eを代入したものと一致する。言い換えると、複素平面単位円上のZ変換がDTFTであると解釈できる。

変換表

元の関数 x(n) Z変換 X(z) 収束領域
δ(n) 1 複素数全体
u(n) [math] \frac{1}{1-z^{-1}}[/math] [math]|z| \gt 1[/math]
anu(n) [math]\frac{1}{1-a z^{-1}}[/math] [math]|z| \gt |a|[/math]
n an u(n) [math]\frac{az^{-1} }{ (1-a z^{-1})^2 }[/math] [math]|z| \gt |a|[/math]
an u(-n-1) [math]\frac{-1}{1-a z^{-1}}[/math] [math] |z| \lt |a|[/math]
n an u(-n-1) [math]\frac{az^{-1} }{ (1-a z^{-1})^2 }[/math] [math]|z| \lt |a|[/math]
cos(ω0n) u(n) [math]\frac{ 1-z^{-1} \cos(\omega_0) }{ 1-2z^{-1}\cos(\omega_0)+ z^{-2} }[/math] [math] |z| \gt 1[/math]
sin(ω0n) u(n) [math] \frac{ z^{-1} \sin(\omega_0) }{ 1-2z^{-1}\cos(\omega_0)+ z^{-2} }[/math] [math] |z| \gt 1[/math]
an cos(ω0n) [math] \frac{ 1-a z^{-1} \cos( \omega_0) }{ 1-2az^{-1}\cos(\omega_0)+ a^2 z^{-2} }[/math] [math] |z| \gt |a|[/math]
an sin(ω0n) [math] \frac{ az^{-1} \sin(\omega_0) }{ 1-2az^{-1}\cos(\omega_0)+ a^2 z^{-2} }[/math] [math] |z| \gt |a|[/math]

関連項目

テンプレート:デジタル信号処理