信号処理について（１）：線形時不変システムと伝達関数

はじめに

初めに，これから扱う対象を明確にします．この記事で扱う対象は「信号」と「システム」，特に，関数$x \colon\ \zahl \to \complex$で表される離散時間の信号と，それに変形を加えるシステムです．

例として，マイクで声を録音することを考えましょう．マイクにより，音は電圧の時間変化に変換されます．これは時間的に連続な信号ですが，計算機で扱うために一定周期で値を記録（標本化）することで，離散時間の信号が得られます．

録音を始めてから，最初に得られたサンプルを$x[0]$とします．同様に，１つ目のサンプルを$x[1]$，２つ目のサンプルを$x[2]$，というふうにして，$n > 0$について$x[n]$を定義します．

ここで，ボイスチェンジャーなどの音声を加工する機器について考えます．ボイスチェンジャーは，信号$x$を加工後の信号$y$に写す写像として表現できます．広くこういうものを指してシステムといいます．

この例では$n < 0$のとき$x[n] = 0$としますが，一般には$n$が負のとき$x[n]$の値を$0$以外で定義することもあります．たとえば，録音済みの音声を加工するときは時刻$t$において既に$x[t-1]$以前の信号が分かっているので，ここを時刻$0$とすれば$x[-1]$が定義されます．

用語の定義

初めに，いくつかの用語を定義しておきます．

離散デルタ関数

遅延作用素

遅延作用素の逆は時間を$1$だけ進めることに注意してください．式で書くと（すべての$n$について）次の通りです．
$$ (\delay^{-1} x)[n] = x[n+1] $$

また，信号$x$に遅延作用素を$k$回作用させた信号を$\delay^k x$と書くことにします． たとえば，$k = 2$のときは
$$ (\delay^2 x)[n] = x[n-2] $$
です．２つを組み合わせて，$k < 0$のときは$\delay^k = (\delay^{-k})^{-1}$とします．

まとめると，$k \in \zahl$に対して，$\delay^k$は信号の時間を時間を$-k$だけ進める作用素です．

線形時不変システム

システムというとあまりに広いので，普通はシステムの中でも比較的扱いやすいものについて考察します．この記事では特に，次で定義される線形時不変システムについて考察します．

線形時不変システムは，現実のシステムを近似するのによく用いられます．その理由の１つは，次に述べる定理により性質を解析しやすいためです．

ここで，証明の最後に出てきた
$$ \sys x = \sum_{n=-\infty}^\infty x[n] \delay^n h $$
という式について少し考察しましょう．時刻$n$における$\sys$の出力$(\sys x)[n]$は次のようになります．
$$ (\sys x)[n] = \pqty{\sum_{k=-\infty}^\infty x[k] \delay^k h}[n] = \sum_{k=-\infty}^\infty x[k] h[n - k] $$

この右辺を$x, h$の「畳み込み」といいます．より正確に，畳み込みの定義を述べると以下の通りです．

これにより先述の定理を言い換えると，$\sys$が線形時不変システムなら常に式
$$ \sys x = x * h = x * (\sys \delta) $$
が成り立ちます．

z変換と伝達関数

前節によれば，線形時不変システムは畳み込みによって表現できるのでした．ここでは，システムを畳み込みによる表記よりも解析しやすい形で書くために，z変換という操作を導入します．

z変換が持つ重要な性質として，畳み込みを積に変換するというものがあります．

この定理により，$\sys$が線形時不変システムなら
$$ \zt (\sys x) = \zt (x * \sys \delta) = (\zt x)[\zt (\sys \delta)] $$
が成り立ちます．よって，$(\zt x)(z) \neq 0$なら
$$ \frac{[\zt (\sys x)](z)}{(\zt x)(z)} = [\zt (\sys \delta)](z) $$
という関係が成り立ちます．左辺を$H(z)$とおき$H$を定義すると，$H$は
$$ [\zt (\sys x)](z) = H(z) (\zt x)(z) $$
という関係を満たす関数です．この関数$H$を伝達関数といいます．

補足：厳密でない証明のフォロー

線形時不変システムとインパルス応答の関係

この証明について少し補足します．まず
$$ \sys \delay^{-1} x = \delay^{-1} \sys x $$
であることを確認します．時不変性の式
$$ \sys \delay x = \delay \sys x $$
について，入力を$\delay^{-1} x$に変えて
$$ \sys \delay \delay^{-1} x = \delay \sys \delay^{-1} x \quad \therefore \sys x = \delay \sys \delay^{-1} x $$
さらに，両辺に$\delay^{-1}$を作用させると
$$ \delay^{-1} \sys x = \sys \delay^{-1} x $$
が得られます．２つの式
$$ \delay \sys x = \sys \delay x \quad \delay^{-1} \sys x = \sys \delay^{-1} x $$
を繰り返し用いれば，すべての整数$n$について
$$ \delay^n \sys x = \sys \delay^n x $$
であることが確認できます．

z変換と畳み込みの関係

和の交換は本来自明ではないので，この証明はあくまでもイメージと言うべきでしょう．厳密に示したい場合，$n$が負になる項がすべて$0$であると仮定すれば整級数の問題に帰着できます．具体的には，$z^{-1}$を$z$と置きなおして，$z$が収束円の内側にあるときを考えれば，級数が絶対収束することから確認できます．この場合，z変換の収束領域はある円の外側（半径が$\infty$になることもある）なので，たまたま境界上にあるときなんかを除けば，$\zt a$と$\zt b$がともに有限な値を持つ範囲であれば成り立ちます．

一見仮定が人工的にも見えますが，伝達関数について考える場合，この仮定はシステムに因果性（後述）を課すということです．そんなに無茶な仮定でもありません．実際，正冪の項が無いz変換（片側z変換）でz変換を定義することもあります．