第7章

フーリエ変換の性質(1): 時間シフトと変調

この章のねらい

  • 時間シフトが振幅スペクトルを変えず、位相を周波数に比例して傾ける(線形位相)ことを理解する
  • 変調がスペクトルの平行移動であり、実信号の cos 変調では ± に複製されることを理解する
  • 時間シフト⇔位相回転、変調⇔周波数シフトという双対の構図を持つ

7.1 時間シフト

やる夫

第3章でフーリエ変換そのものは習ったお。ここからしばらくは「性質」の章だって聞いたお。性質って、何のありがたみがあるんだお?

やらない夫

ありがたみは絶大だ。変換の式を毎回まじめに積分しなくても、「波形をこういじったらスペクトルはこう変わる」という対応表があれば、計算を一段すっ飛ばせる。今日はその対応表の最初の2行、時間シフト変調を作る。

やる夫

時間シフトってのは…波形を横にずらすだけかお? それだけで何か面白いことが起きるとは思えないお。

やらない夫

それが起きるんだ。信号 x(t)x(t) を時間的に t0t_0 だけ遅らせた信号は x(tt0)x(t - t_0) と書ける。これのフーリエ変換を計算してみよう。定義どおりに書くぞ。

F[x(tt0)]=x(tt0)ejωtdt(7.1)\mathcal{F}[x(t-t_0)] = \int_{-\infty}^{\infty} x(t - t_0)\, e^{-j\omega t}\, dt \tag{7.1}
やる夫

積分の中に tt0t - t_0 がいて、外に tt がいて、なんかややこしいお。

やらない夫

こういうときは置換積分だ。τ=tt0\tau = t - t_0 と置く。すると t=τ+t0t = \tau + t_0dt=dτdt = d\tau。積分範囲は τ\tau-\infty から \infty のままだ。代入する。

x(τ)ejω(τ+t0)dτ=ejωt0x(τ)ejωτdτ(7.2)\int_{-\infty}^{\infty} x(\tau)\, e^{-j\omega (\tau + t_0)}\, d\tau = e^{-j\omega t_0} \int_{-\infty}^{\infty} x(\tau)\, e^{-j\omega \tau}\, d\tau \tag{7.2}

ejωt0e^{-j\omega t_0}τ\tau に関係ない定数だから、積分の外に出せる。そして残った積分は、X(ω)X(\omega) そのものだ。

やる夫

あ、後ろの積分は元の X(ω)X(\omega) だお! つまり…

やらない夫

結論はこうなる。

x(tt0)    ejωt0X(ω)(7.3)x(t - t_0) \;\longleftrightarrow\; e^{-j\omega t_0}\, X(\omega) \tag{7.3}

**波形を時間で t0t_0 ずらすと、スペクトルには ejωt0e^{-j\omega t_0} という因子が掛かる。**これが時間シフトの性質だ。

やる夫

掛かるのが ejωt0e^{-j\omega t_0}…。これって絶対値が1の複素数だお。大きさが1なら、X(ω)|X(\omega)| には影響しないってことかお?

やらない夫

そこに気づけるなら大したものだ。そのとおり。ejωt0=1|e^{-j\omega t_0}| = 1 だから、

ejωt0X(ω)=X(ω)\left| e^{-j\omega t_0} X(\omega) \right| = |X(\omega)|

振幅スペクトルは時間シフトでまったく変わらない。変わるのは位相だけだ。元の位相を X(ω)\angle X(\omega) とすると、シフト後の位相は

 ⁣(ejωt0X(ω))=X(ω)ωt0\angle\!\left( e^{-j\omega t_0} X(\omega) \right) = \angle X(\omega) - \omega t_0

つまり位相が ωt0-\omega t_0 だけ足される。これは ω\omega に比例した傾きだ。だから「線形位相」と呼ぶ。

やる夫

待つお、直観が追いつかないお。波形を遅らせただけなのに、なんで「周波数に比例して位相がずれる」なんて妙なことになるんだお?

やらない夫

いい疑問だ。1個の正弦波で考えろ。cosωt\cos\omega t を時間で t0t_0 遅らせると cosω(tt0)=cos(ωtωt0)\cos\omega(t - t_0) = \cos(\omega t - \omega t_0) になるだろ。位相が ωt0\omega t_0 ぶん遅れた。ここで大事なのは、同じ時間 t0t_0 の遅れでも、周波数 ω\omega が高い波ほど位相のずれ ωt0\omega t_0 が大きいってことだ。

やる夫

あー…。速く振動してる波は、同じ時間ずらしても、何周期ぶんもずれちゃうってことかお。ゆっくりの波は1周期の途中までしかずれない。

やらない夫

その理解で完璧だ。「波形を遅らせる」を周波数の言葉に翻訳すると、「各周波数成分の位相を、周波数に比例して遅らせる」になる。下のデモで確かめろ。長さ5の矩形パルスを n0n_0 だけ遅らせる装置だ。スライダーを動かすと、上のパルスは右に動くが、真ん中の振幅スペクトルはぴくりとも動かない。下の位相だけが傾きを変える。

TIME SHIFT → LINEAR PHASE INTERACTIVE
シフト量 n0 を上げてみる。上段のパルスは右へ動くのに、中段の振幅スペクトルはぴくりとも動かず、下段の位相の傾きだけが急になる。これが「時間シフトは振幅を変えず位相を周波数に比例して傾けるだけ」=線形位相。情報の置き場所が位相に移ったことを示している。
やる夫

ほんとだお! 真ん中の山の形は固定で、下のギザギザの傾きだけがどんどん急になっていくお。位相だけが情報の置き場所なんだお。

やらない夫

そういうことだ。ついでに離散信号版も書いておく。連続の積分が和になるだけで、導出も結論も同じ形だ。

x[nn0]    ejωn0X(ejω)(7.4)x[n - n_0] \;\longleftrightarrow\; e^{-j\omega n_0}\, X(e^{j\omega}) \tag{7.4}

デモはこの離散版を描いている。パルス中心が n0+2n_0 + 2 にいるから、位相の傾きは (n0+2)-(n_0 + 2) だ。

7.2 変調

やらない夫

対応表の2行目、変調に行く。今度は波形を時間でずらすのではなく、波形に回転する複素指数を掛ける。x(t)ejω0tx(t) e^{j\omega_0 t} のフーリエ変換を計算しよう。

やる夫

掛け算かお。さっきは引数をずらしたけど、今度は外から掛ける。…なんかさっきの裏返しっぽい匂いがするお。

やらない夫

その匂いは正しい。定義に放り込むぞ。

F ⁣[x(t)ejω0t]=x(t)ejω0tejωtdt=x(t)ej(ωω0)tdt(7.5)\mathcal{F}\!\left[ x(t) e^{j\omega_0 t} \right] = \int_{-\infty}^{\infty} x(t)\, e^{j\omega_0 t}\, e^{-j\omega t}\, dt = \int_{-\infty}^{\infty} x(t)\, e^{-j(\omega - \omega_0) t}\, dt \tag{7.5}

掛けた ejω0te^{j\omega_0 t} と変換核の ejωte^{-j\omega t} が指数の肩で合体して、ej(ωω0)te^{-j(\omega - \omega_0)t} になった。これは XX の定義式で ω\omega のところを ωω0\omega - \omega_0 に置き換えたものだ。だから

x(t)ejω0t    X(ωω0)(7.6)x(t)\, e^{j\omega_0 t} \;\longleftrightarrow\; X(\omega - \omega_0) \tag{7.6}

波形に ejω0te^{j\omega_0 t} を掛けると、スペクトルが ω0\omega_0 だけ右に平行移動する

やる夫

おー、きれいだお! 時間シフトはスペクトルに ee^{} を掛けた、変調はスペクトルを引っ越しさせた。本当に裏返しだお。

やらない夫

その対称性は 7.3 で整理する。先に実用上の問題を片付けよう。ejω0te^{j\omega_0 t} は複素数だから、現実の信号(実数値)に掛けると複素信号になってしまう。実際の電波や音でやりたいのは、実数の信号 x(t)x(t)実数の cosω0t\cos\omega_0 t を掛けることだ。これを計算するとどうなるか。

やる夫

cos\cos はオイラーの公式で…12(ejω0t+ejω0t)\frac{1}{2}(e^{j\omega_0 t} + e^{-j\omega_0 t}) だったお。さっきのを2回使えばいいのかお?

やらない夫

自分で道筋が見えてるじゃないか。やってみろ。

やる夫

えーと、x(t)cosω0t=12x(t)ejω0t+12x(t)ejω0tx(t)\cos\omega_0 t = \frac{1}{2} x(t) e^{j\omega_0 t} + \frac{1}{2} x(t) e^{-j\omega_0 t} だお。前半は式 (7.6) で X(ωω0)X(\omega - \omega_0)、後半は ω0\omega_0 がマイナスだから X(ω+ω0)X(\omega + \omega_0) になるお。半分ずつ足して…

やらない夫

完璧だ。まとめるとこうだ。

x(t)cosω0t    12(X(ωω0)+X(ω+ω0))(7.7)x(t)\cos\omega_0 t \;\longleftrightarrow\; \frac{1}{2}\Big( X(\omega - \omega_0) + X(\omega + \omega_0) \Big) \tag{7.7}

実信号に cosω0t\cos\omega_0 t を掛けると、元のスペクトルが +ω0+\omega_0ω0-\omega_0 の2か所に、半分ずつコピーされる

やる夫

なんで2か所なんだお。1か所に引っ越すんじゃダメだったのかお。

やらない夫

ダメというか、できないんだ。1か所だけにずらすには複素指数1本が要るが、実数の信号を作るには第2章でやったとおり正負の回転をペアで使うしかない。cos\cos は正負の回転がペアになったものだから、スペクトルも正負の2か所に割れる。実数であることの代償だ。

やる夫

でも係数が 12\frac{1}{2} ってことは、高さが半分になってるお。情報も半分に減っちゃったってことかお? なんだか損した気分だお。

やらない夫

そこは安心しろ。減ったんじゃなくて、1か所だったものが2か所に分かれただけだ。第9章で「総量(パワー)」をきちんと測る道具を出すが、片方の山の高さが半分でも、それが正負の2つあるから足し合わせれば帳尻が合う。エネルギーで言えば、各山が (12)2\left(\frac{1}{2}\right)^2 になっても2つあるから、合計は元の半分——というのは「掛けた cos\cos の振幅が1で、その実効パワーが 12\frac{1}{2}」というだけの当たり前の話で、xx の中身(情報)は両方の山にちゃんと丸ごとコピーされている。だから受信側で片方を拾えば元に戻せる。分散しただけで消えてはいない。

MODULATION — SPECTRUM SHIFT INTERACTIVE
変調周波数 ω0 を上げると緑のスペクトルが右へ引っ越す(破線が元の X(ω))。チェックを入れて cos 変調にすると、山が +ω0 と −ω0 の2か所に半分の高さで割れる。これが「実信号を作るには正負の回転がペアで要る」代償であり、AMラジオが音声を高周波帯へ乗せ替える仕組みそのもの。
やる夫

チェックを外してると山がまるごと右へ滑っていって、入れると左右に分裂して半分の高さになるお! 物理の電波の話、これだったのかお。

やらない夫

そう、これがAMラジオの入り口だ。マイクで拾った音声 x(t)x(t) は数 kHz までの低い周波数しか持たない。そのままアンテナに流しても効率よく飛ばない。そこで搬送波 cosω0t\cos\omega_0 t(例えば 1000 kHz)を掛けてやると、音声のスペクトルが搬送波の周波数 ω0\omega_0 のまわりに引っ越す。電波として飛ばしやすい高い周波数帯に音声を乗せ替える――それが変調という言葉の由来だ。

やる夫

なるほどだお。スペクトルの引っ越し業者が変調だお。受け取った側は逆に引っ越し戻せばいいわけかお。

やらない夫

そういうことだ。受信側でもう一度 cosω0t\cos\omega_0 t を掛ければ、スペクトルの一部が原点に戻ってくる(復調)。詳細は省くが、変調が「スペクトルを動かす道具」だと掴めれば今は十分だ。

7.3 双対性の整理

やらない夫

今日の2つを並べて眺めよう。

x(tt0)ejωt0X(ω)(時間シフト)x(t - t_0) \longleftrightarrow e^{-j\omega t_0} X(\omega) \qquad\text{(時間シフト)}ejω0tx(t)X(ωω0)(変調)e^{j\omega_0 t}\, x(t) \longleftrightarrow X(\omega - \omega_0) \qquad\text{(変調)}
やる夫

うわ、見事に対称だお。

  • 時間でずらす ⇔ 周波数で位相を回す(ee^{} を掛ける)
  • 時間で位相を回す(ee^{} を掛ける)⇔ 周波数でずらす

時間の世界と周波数の世界で、「ずらす」と「ee^{} を掛ける」がそっくり入れ替わってるお。

やらない夫

それを双対性という。時間領域と周波数領域は鏡のような関係にあって、片方で「平行移動」にあたる操作は、もう片方では「位相の回転(複素指数の掛け算)」にあたる。この鏡の構図は、次の第8章の「たたみこみと積」でも、もっと劇的な形でまた出てくる。

やる夫

性質の章、思ったより面白いお。計算をサボるための表どころか、時間と周波数の世界観そのものだお。

この章のまとめ
  • 時間シフト x(tt0)ejωt0X(ω)x(t-t_0) \leftrightarrow e^{-j\omega t_0}X(\omega)。振幅スペクトルは不変、位相が ωt0-\omega t_0 と周波数に比例して傾く(線形位相)。離散版も x[nn0]ejωn0X(ejω)x[n-n_0]\leftrightarrow e^{-j\omega n_0}X(e^{j\omega}) で同形
  • 「波形を遅らせる」=「各周波数成分の位相を周波数に比例して遅らせる」。高い周波数ほど位相のずれが大きい
  • 変調 x(t)ejω0tX(ωω0)x(t)e^{j\omega_0 t} \leftrightarrow X(\omega-\omega_0)。スペクトルの平行移動
  • 実信号の cos 変調 x(t)cosω0t12(X(ωω0)+X(ω+ω0))x(t)\cos\omega_0 t \leftrightarrow \frac{1}{2}(X(\omega-\omega_0)+X(\omega+\omega_0))。±に半分ずつ複製される(AMラジオの原理)
  • 時間シフト⇔位相回転、変調⇔周波数シフトは双対の関係