The impulse response(脈衝響應)$h[n]$ completely describes an LSI (linear shift-invariant) system: the output is always the convolution $y[n] = x[n] \circledast h[n] = \sum_{m} h[m]\,x[n-m]$.脈衝響應 $h[n]$ 就完整描述了一個 LSI(線性平移不變)系統:輸出永遠是摺積 $y[n] = x[n] \circledast h[n] = \sum_{m} h[m]\,x[n-m]$。インパルス応答 $h[n]$ は LSI(線形シフト不変)システムを完全に記述する。出力はつねに畳み込み $y[n] = x[n] \circledast h[n] = \sum_{m} h[m]\,x[n-m]$ である。
Feed in a complex sinusoid(複數弦波)$x[n] = e^{j\omega_0 n}$ and watch the convolution factor:輸入一個複數弦波 $x[n] = e^{j\omega_0 n}$,看看摺積會變成什麼:複素正弦波 $x[n] = e^{j\omega_0 n}$ を入力して、畳み込みがどうなるかを見てみよう。
The complex scale factor is the DTFT of the impulse response evaluated at the input frequency. Writing it in polar form $H(e^{j\omega}) = |H(e^{j\omega})|\,e^{j\theta(\omega)}$, the output is這個複數比例因子,就是脈衝響應的 DTFT 在輸入頻率處的值。寫成極座標形式 $H(e^{j\omega}) = |H(e^{j\omega})|\,e^{j\theta(\omega)}$,輸出便是この複素倍率は、インパルス応答の DTFT を入力周波数で評価したものである。極形式 $H(e^{j\omega}) = |H(e^{j\omega})|\,e^{j\theta(\omega)}$ で書けば、出力は次のようになる。
What these three slides say is this: after a sinusoid passes through a linear shift-invariant (LSI) system, what comes out is still a sinusoid at the same frequency; only its size has been scaled and its phase shifted. And how much each frequency is scaled and shifted by is exactly the frequency response $H(e^{j\omega})$.這三張投影片說的是:正弦波通過線性非時變(LSI)系統後,出來的還是同頻率的正弦波,只是大小被縮放、相位被平移。而每個頻率被縮放和平移多少,就是頻率響應 $H(e^{j\omega})$。この三枚のスライドが言っているのはこうである。正弦波が線形時不変(LSI)システムを通ったあとも、出てくるのは同じ周波数の正弦波であり、大きさが拡大縮小され、位相がずらされるだけである。そして、各周波数がどれだけ拡大縮小され、ずらされるかこそが周波数応答 $H(e^{j\omega})$ である。
Any signal can be taken apart into sinusoids. The IDTFT writes $x[n]$ as任何訊號都能拆成正弦波。IDTFT 把 $x[n]$ 寫成どんな信号も正弦波に分解できる。IDTFT は $x[n]$ を次のように書く。
That is, a weighted sum of infinitely many complex sinusoids $e^{j\omega n}$, the weight being $X(e^{j\omega})$.也就是無數個複數正弦波 $e^{j\omega n}$ 的加權總和,權重是 $X(e^{j\omega})$。すなわち、無数の複素正弦波 $e^{j\omega n}$ の重み付き和であり、重みは $X(e^{j\omega})$ である。
A sinusoid is an eigenfunction of an LSI system. A sinusoid goes in, and what comes out is only "scaled"; the shape (the frequency) does not "change". This is the same idea as the eigenvector of linear algebra:正弦波是 LSI 系統的特徵函數(eigenfunction)。正弦波進去,出來只是「縮放」,形狀(頻率)不會「改變」。這和線性代數中特徵向量的概念相同:正弦波は LSI システムの固有関数(eigenfunction)である。正弦波が入り、出てくるのは「拡大縮小」されたものだけで、形(周波数)は「変わらない」。これは線形代数における固有ベクトルの考えと同じである。
| Linear algebra線性代數線形代数 | Signals and systems訊號系統信号とシステム |
|---|---|
| Matrix $A$矩陣 $A$行列 $A$ | LSI system $h[n]$LSI 系統 $h[n]$LSI システム $h[n]$ |
| Eigenvector $\mathbf{v}$: $A\mathbf{v} = \lambda\mathbf{v}$特徵向量 $\mathbf{v}$:$A\mathbf{v} = \lambda\mathbf{v}$固有ベクトル $\mathbf{v}$:$A\mathbf{v} = \lambda\mathbf{v}$ | Eigenfunction $e^{j\omega_0 n}$: output $= H(e^{j\omega_0})\,e^{j\omega_0 n}$特徵函數 $e^{j\omega_0 n}$:輸出 $= H(e^{j\omega_0})\,e^{j\omega_0 n}$固有関数 $e^{j\omega_0 n}$:出力 $= H(e^{j\omega_0})\,e^{j\omega_0 n}$ |
| Eigenvalue $\lambda$特徵值 $\lambda$固有値 $\lambda$ | $H(e^{j\omega_0})$ (a complex number)$H(e^{j\omega_0})$(一個複數)$H(e^{j\omega_0})$(一つの複素数) |
Know the scaling at every frequency and you know the whole system. Because the input can be taken apart into sinusoids, because each sinusoid is only scaled, and because the system is linear, writing down the amount of scaling at every frequency is enough to work out the output for any input at all. That set of scaling amounts is the frequency response.知道每個頻率的縮放,就知道整個系統。因為輸入可以拆成正弦波、每個正弦波只被縮放、而系統是線性的,所以把每個頻率的縮放量記下來,就能算出任何輸入的輸出。這組縮放量就是頻率響應。各周波数での拡大縮小を知れば、システム全体を知ったことになる。入力は正弦波に分解でき、各正弦波は拡大縮小されるだけであり、そしてシステムは線形である。だから各周波数での拡大縮小の量を書き留めておけば、どんな入力に対しても出力を計算できる。この拡大縮小の量の集まりが周波数応答である。
An LSI system is completely described by its impulse response $h[n]$, and the output is the convolution:LSI 系統由脈衝響應 $h[n]$ 完全描述,輸出是摺積:LSI システムはインパルス応答 $h[n]$ によって完全に記述され、出力は畳み込みである。
Feed in the complex sinusoid $x[n] = e^{j\omega_0 n}$:輸入複數正弦波 $x[n] = e^{j\omega_0 n}$:複素正弦波 $x[n] = e^{j\omega_0 n}$ を入力する。
The key is that $e^{j\omega_0(n-m)} = e^{-j\omega_0 m}\cdot e^{j\omega_0 n}$ can be split apart: the part that has to do with $m$ (the blue line) stays inside the sum, and the part that has to do with $n$ (the green line) is taken outside. The blue-line part is exactly the DTFT of $h[n]$ evaluated at $\omega_0$:關鍵在於 $e^{j\omega_0(n-m)} = e^{-j\omega_0 m}\cdot e^{j\omega_0 n}$ 可以拆開:和 $m$ 有關的部分(藍線)留在求和裡,和 $n$ 有關的部分(綠線)提到外面。藍線部分正好是 $h[n]$ 的 DTFT 在 $\omega_0$ 的值:要点は $e^{j\omega_0(n-m)} = e^{-j\omega_0 m}\cdot e^{j\omega_0 n}$ と分けられることである。$m$ に関わる部分(青い線)は和の中に残り、$n$ に関わる部分(緑の線)は外へ出す。青い線の部分はちょうど $h[n]$ の DTFT を $\omega_0$ で評価した値である。
where $H(e^{j\omega}) = |H(e^{j\omega})|\,e^{j\theta(\omega)}$. The output is a sinusoid at the very same frequency, its amplitude multiplied by $|H|$ and its phase with $\theta$ added to it.其中 $H(e^{j\omega}) = |H(e^{j\omega})|\,e^{j\theta(\omega)}$。輸出是同一個頻率的正弦波,振幅乘上 $|H|$,相位加上 $\theta$。ここで $H(e^{j\omega}) = |H(e^{j\omega})|\,e^{j\theta(\omega)}$ である。出力はまったく同じ周波数の正弦波であり、振幅には $|H|$ が掛かり、位相には $\theta$ が加わる。
What lets this derivation go through is the shift-invariance of the system: every term $h[m]$ applies the same delay $m$ to the input, and delaying a sinusoid changes only its phase, never its frequency.這個推導能成立,靠的是系統的非時變性:每一項 $h[m]$ 都對輸入做同樣的延遲 $m$,而延遲一個正弦波只會改變它的相位,不會改變頻率。この導出が成り立つのは、システムの時不変性のおかげである。どの項 $h[m]$ も入力に同じ遅延 $m$ を与え、正弦波を遅らせても位相が変わるだけで、周波数は変わらない。
| Name名稱名称 | Symbol符號記号 | What it does to a sinusoid對正弦波的作用正弦波への作用 |
|---|---|---|
| magnitude response大小響應(magnitude response)振幅応答(magnitude response) | $\lvert H(e^{j\omega})\rvert$ | gain: the factor the amplitude becomes增益:振幅變成幾倍利得:振幅が何倍になるか |
| phase response相位響應(phase response)位相応答(phase response) | $\theta(\omega) = \arg H(e^{j\omega})$ | a phase shift (equivalent to a delay in time)相位平移(等效於時間延遲)位相のずれ(時間遅延と等価) |
Real sinusoids: if $h[n]$ is real then $H(e^{-j\omega}) = H^*(e^{j\omega})$ (the conjugate symmetry met earlier). Substituting the two terms of $\cos\omega_0 n = \frac{1}{2}(e^{j\omega_0 n} + e^{-j\omega_0 n})$ separately gives實數正弦波:若 $h[n]$ 是實數,$H(e^{-j\omega}) = H^*(e^{j\omega})$(前面講過的共軛對稱)。把 $\cos\omega_0 n = \frac{1}{2}(e^{j\omega_0 n} + e^{-j\omega_0 n})$ 兩項分別代入,得到実正弦波:$h[n]$ が実数であれば $H(e^{-j\omega}) = H^*(e^{j\omega})$ である(さきに述べた共役対称性)。$\cos\omega_0 n = \frac{1}{2}(e^{j\omega_0 n} + e^{-j\omega_0 n})$ の二項をそれぞれ代入すると、次が得られる。
Again the same frequency, scaled and shifted.一樣是同頻率、被縮放與平移。やはり同じ周波数であり、拡大縮小とずれを受けるだけである。
$h[n] = [\,0.5,\ 0.5\,]$, that is, $y[n] = \frac{1}{2}(x[n] + x[n-1])$:$h[n] = [\,0.5,\ 0.5\,]$,即 $y[n] = \frac{1}{2}(x[n] + x[n-1])$:$h[n] = [\,0.5,\ 0.5\,]$、すなわち $y[n] = \frac{1}{2}(x[n] + x[n-1])$ である。
So $|H| = \cos(\omega/2)$ and $\theta = -\omega/2$ (a delay of half a sample).所以 $|H| = \cos(\omega/2)$、$\theta = -\omega/2$(延遲半個樣本)。したがって $|H| = \cos(\omega/2)$、$\theta = -\omega/2$(半標本ぶんの遅れ)である。
| Input frequency輸入頻率入力周波数 | Gain增益利得 | Phase相位位相 | The intuition直覺直観 |
|---|---|---|---|
| $\omega = 0$ (DC)$\omega = 0$(直流)$\omega = 0$(直流) | 1 | 0 | averaging a constant leaves it unchanged常數取平均,不變定数を平均しても変わらない |
| $\omega = \pi/2$ | 0.707 | $-\pi/4$ | partial cancellation部分抵消部分的に打ち消し合う |
| $\omega = \pi$ | 0 | — | $+1, -1, +1, \dots$ — the average of two neighbouring points is 0, complete cancellation$+1, -1, +1, \dots$ 相鄰兩點平均為 0,完全抵消$+1, -1, +1, \dots$ 隣り合う二点の平均は 0 であり、完全に打ち消し合う |
This is a lowpass filter.這是一個低通濾波器。これは低域通過フィルタである。
Only the result is quoted here; it is obtained step by step in Worked Example 10-A, continued — the $M=2$ case, in §10.1.4.此處只引用結論;完整的逐步推導在 §10.1.4 的〈範例 10-A(續)— $M=2$ 的情形〉。ここでは結論だけを引いている。一歩ずつの完全な導出は §10.1.4 の「例題 10-A(続き)— $M=2$ の場合」にある。
The $H_k(e^{j\omega})$ of the STFT filterbank earlier is the frequency response of the $k$-th channel filter: the window spectrum shifted to $\omega_k$. When some sinusoid passes through it, the gain is the value of $|H_k|$ at that frequency, so sinusoids whose frequency is near $\omega_k$ get through and those far from it are attenuated (leaving only the leakage of the sidelobes).前面 STFT 濾波器組中的 $H_k(e^{j\omega})$,就是第 $k$ 列濾波器的頻率響應:窗頻譜平移到 $\omega_k$。某個正弦波通過它時,增益就是 $|H_k|$ 在該頻率的值,所以頻率在 $\omega_k$ 附近的正弦波通過、遠離的被衰減(只剩旁瓣的洩漏量)。さきの STFT フィルタバンクにおける $H_k(e^{j\omega})$ とは、第 $k$ 列のフィルタの周波数応答である。窓のスペクトルを $\omega_k$ へ平行移動したものである。ある正弦波がそこを通るとき、利得はその周波数における $|H_k|$ の値であり、だから周波数が $\omega_k$ の近くにある正弦波は通り抜け、遠いものは減衰させられる(残るのはサイドローブの漏れだけである)。
That filterbank is §9.5 of Chapter 9.該濾波器組見第 9 章 §9.5。そのフィルタバンクは第 9 章 §9.5 にある。
Pick a filter, drag the input frequency $\omega_0$, and watch how the output is scaled and shifted.選一個濾波器、拖動輸入頻率 $\omega_0$,看輸出如何被縮放與平移。フィルタを一つ選び、入力周波数 $\omega_0$ をドラッグして、出力がどのように拡大縮小され、ずらされるかを見てほしい。
The blue dots in the figure are the output worked out by convolution directly, and the orange dashed line is what the frequency response predicts, $|H|\cos(\omega_0 n + \theta)$. Whichever system and whichever frequency are chosen, the blue dots land on the orange dashed line. That is what "a sinusoid is an eigenfunction" means: the frequency of the output is the same as that of the input, and only the amplitude and the phase change.圖中的藍點是直接用摺積算出的輸出,橘虛線是用頻率響應預測的結果 $|H|\cos(\omega_0 n + \theta)$。不論選哪個系統、哪個頻率,藍點都落在橘虛線上,這就是「正弦波是特徵函數」的意思:輸出的頻率與輸入相同,只有振幅和相位改變。図の青い点は畳み込みで直接求めた出力であり、オレンジの破線は周波数応答による予測 $|H|\cos(\omega_0 n + \theta)$ である。どのシステムを選んでも、どの周波数を選んでも、青い点はオレンジの破線の上に落ちる。これが「正弦波は固有関数である」ということの意味である。出力の周波数は入力と同じであり、変わるのは振幅と位相だけである。
Things to try:可以試試:試してみるとよい。
In practice signals are real. Decompose a real cosine into two complex sinusoids and apply the eigenfunction property to each:實務上的訊號是實數。把一個實餘弦拆成兩個複數弦波,再對各自套用特徵函數性質:実際の信号は実数である。実余弦を二つの複素正弦波に分解し、それぞれに固有関数の性質を適用する。
For a real impulse response $h[n]$ the frequency response is conjugate-symmetric, $H(e^{-j\omega}) = H^{*}(e^{j\omega}) = |H(e^{j\omega})|e^{-j\theta(\omega)}$, so the two terms recombine into a real cosine:若脈衝響應 $h[n]$ 是實數,頻率響應就具有共軛對稱性 $H(e^{-j\omega}) = H^{*}(e^{j\omega}) = |H(e^{j\omega})|e^{-j\theta(\omega)}$,於是兩項會重新合成一個實餘弦:インパルス応答 $h[n]$ が実数であれば、周波数応答は共役対称 $H(e^{-j\omega}) = H^{*}(e^{j\omega}) = |H(e^{j\omega})|e^{-j\theta(\omega)}$ となり、二つの項は再び実余弦にまとまる。
Most signals start at a finite time(訊號從有限時間開始), e.g. a gated complex sinusoid $x[n] = e^{j\omega_0 n}\,u[n]$. What does gating do? Split the convolution sum:大多數訊號都從某個有限時刻才開始(訊號從有限時間開始),例如被閘控的複數弦波 $x[n] = e^{j\omega_0 n}\,u[n]$。閘控會造成什麼?把摺積和拆開來看:ほとんどの信号は有限の時刻から始まる(信号は有限時間から始まる)。たとえばゲートされた複素正弦波 $x[n] = e^{j\omega_0 n}\,u[n]$ である。ゲートは何をもたらすのか。畳み込み和を分けてみよう。
Fig. 10-1 fixes one system at one frequency. The figure below puts the tail itself on screen: choose a system, slide $n$, and watch which part of $h[m]$ the gated input has not reached yet. Whatever is still missing is the transient — and the readout measures, rather than claims, that steady state plus transient rebuilds $y[n]$.圖 10-1 只固定了一個系統、一個頻率。下面這張圖把尾巴本身搬上螢幕:選一個系統、拖動 $n$,看看被閘控的輸入還沒碰到 $h[m]$ 的哪一段。還沒碰到的那一段就是暫態——而下方的讀數是把「穩態加暫態等於 $y[n]$」量出來,不是宣稱出來。図 10-1 は一つのシステムを一つの周波数で固定したものである。下の図は尾部そのものを画面に出す。システムを選んで $n$ を動かし、ゲートされた入力がまだ届いていない $h[m]$ の部分を見ればよい。まだ届いていない部分こそが過渡である。そして下の読み出しは、定常状態と過渡の和が $y[n]$ を復元することを、主張するのではなく測って示す。
Fig. 10-1 and the interactive figure above both watch the decomposition unfold in time. The identity behind it, though, is a statement about vectors: $y[n] = \sum_{m=0}^{n} h[m]\,e^{j\omega_0(n-m)}$ is a sum of rotating phasors, the steady state is the complete sum $H(e^{j\omega_0})e^{j\omega_0 n}$, and the transient is precisely the part of that vector sum the gated input has not yet supplied. The figure below draws the sum itself. Panel ① lays the terms $q_m = h[m]e^{-j\omega_0 m}$ tip to tail: the travelling end is the partial sum $P_n$, the far end is $H(e^{j\omega_0})$, and the arrow still separating them is the whole transient. Panel ② carries that same arrow over to the output. Where the figure above measured how big the tail is, this one shows which way it points — and why an FIR chain shuts exactly at $n = N-1$ while an IIR spiral only ever shrinks towards $H$.圖 10-1 與上面那張互動圖,都是看著這個分解隨時間展開。但它背後的恆等式講的其實是向量:$y[n] = \sum_{m=0}^{n} h[m]\,e^{j\omega_0(n-m)}$ 是一堆旋轉相量的和,穩態就是那個完整的和 $H(e^{j\omega_0})e^{j\omega_0 n}$,而暫態恰好是這個向量和裡、被閘控的輸入還沒補上的那一段。下面這張圖畫的就是這個和本身。面板 ① 把各項 $q_m = h[m]e^{-j\omega_0 m}$ 頭尾相接地排出來:會走動的那一端是部分和 $P_n$,另一端是 $H(e^{j\omega_0})$,而還隔在兩者之間的那支箭頭就是暫態的全部。面板 ② 再把同一支箭頭搬到輸出上。上面那張圖量的是尾巴有多大,這張圖看的是它指向哪裡——也因此說明了,為什麼 FIR 的折線會在 $n = N-1$ 恰好閉合,而 IIR 的螺旋卻只能一路縮向 $H$。図 10-1 と上のインタラクティブな図は、どちらもこの分解が時間とともに展開する様子を見るものである。しかし、その背後にある恒等式が語っているのはベクトルである。$y[n] = \sum_{m=0}^{n} h[m]\,e^{j\omega_0(n-m)}$ は回転するフェーザの和であり、定常状態はその完全な和 $H(e^{j\omega_0})e^{j\omega_0 n}$ であって、過渡とは、このベクトル和のうちゲートされた入力がまだ供給していない部分にほかならない。下の図はその和そのものを描く。パネル ① は各項 $q_m = h[m]e^{-j\omega_0 m}$ を頭と尾でつないで並べる。動く側の端が部分和 $P_n$ であり、反対の端が $H(e^{j\omega_0})$ であり、両者をなお隔てている矢印が過渡のすべてである。パネル ② は同じ矢印を出力の側へ移す。上の図が尾部の大きさを測ったのに対し、この図は尾部がどちらを向いているかを示す。FIR の折れ線が $n = N-1$ でちょうど閉じ、IIR のらせんは $H$ へ縮み続けるだけである理由も、そこから見えてくる。
The MA filter $y[n] = \frac{1}{M}\sum_{\ell=0}^{M-1} x[n-\ell]$ is convolution with $h[n] = \frac{1}{M}$ for $0 \le n \le M-1$ (zero otherwise).移動平均濾波器 $y[n] = \frac{1}{M}\sum_{\ell=0}^{M-1} x[n-\ell]$ 就是與 $h[n] = \frac{1}{M}$($0 \le n \le M-1$,其餘為零)做摺積。移動平均フィルタ $y[n] = \frac{1}{M}\sum_{\ell=0}^{M-1} x[n-\ell]$ は、$h[n] = \frac{1}{M}$($0 \le n \le M-1$、それ以外は 0)との畳み込みである。
The warning just above asks which kind of jump you are looking at; the figure below answers that question for this filter, and answers it by measurement. The interactive figure of §10.1.1 already carries a 5-point average among its filters, and the one in §10.1.3 carries an 8-tap FIR, but neither takes $H(e^{j\omega})$ apart: this is the only figure in the chapter that separates it into the signed real amplitude $A(\omega)$, the jump index $r$ and the phase $\theta$, and sets the sign change of $A$ beside the $\pi$ step of $\theta$ at the same $\omega$. Panel ③ draws the folded and the unfolded phase over one another and prints the largest gap it measures between them, so the second kind of jump the warning names — the cosmetic $2\pi$ wrap — is ruled out on screen rather than taken on trust. Park either frequency exactly on a zero and the readout refuses to name a sign or an angle there, and says why. Sliding $M$ down to $2$ arrives at the case worked out by hand below.上面那段警告要你分辨自己看到的是哪一種跳變;下面這張圖就針對這個濾波器把它分辨清楚,而且是用量測的。§10.1.1 的互動圖裡已經有一個 5 點平均,§10.1.3 的也已經有一個 8 抽頭的 FIR,但兩者都沒有把 $H(e^{j\omega})$ 拆開:全章只有這張圖把它拆成帶正負號的實數振幅 $A(\omega)$、跳變指標 $r$ 與相位 $\theta$,並把 $A$ 的變號和 $\theta$ 的 $\pi$ 階梯擺在同一個 $\omega$ 上並排。面板 ③ 把折回前與折回後的相位疊著畫,並把它量到的最大差距印出來,所以警告裡提到的第二種跳變——純屬繪圖造成的 $2\pi$ 捲繞——可以當場被排除,而不是被相信。把任一個頻率剛好停在零點上,讀數會拒絕在那裡指認正負號或角度,並說明理由。把 $M$ 一路拉到 $2$,落到的就是下面手算的那個情形。すぐ上の注意は、いま見ている跳びがどちらの種類かを問うている。下の図は、このフィルタについてその問いに測定で答えるものである。§10.1.1 の図にはすでにフィルタの一つとして 5 点平均があり、§10.1.3 の図には 8 タップの FIR があるが、どちらも $H(e^{j\omega})$ を分解してはいない。本章でこれを符号つきの実振幅 $A(\omega)$、跳びの指標 $r$、位相 $\theta$ に分け、$A$ の符号の反転と $\theta$ の $\pi$ の段差を同じ $\omega$ の上に並べて示すのは、この図だけである。パネル ③ は折り返す前と後の位相を重ねて描き、測定した最大の差を表示するので、注意が述べるもう一方の跳び、すなわち見かけだけの $2\pi$ の巻き戻しは、信じるのではなく画面上で排除できる。どちらかの周波数をちょうど零点の上に置くと、読み出しはそこで符号も角度も名指すことを拒み、その理由を述べる。$M$ を $2$ まで下げると、下で手計算した場合に行き着く。
Interactive — 互動圖 — インタラクティブ — The defaults are the example from the slides ($M = 5$, $\omega_0 = 0.1\pi$, $\omega_1 = 0.5\pi$). The readout above sets out the calculation of $A$, $r$, $\theta$ and $|H|$ term by term, in agreement with the written derivation above. The four panels correspond to the four steps of that derivation:預設值是投影片的例子($M = 5$、$\omega_0 = 0.1\pi$、$\omega_1 = 0.5\pi$)。上方讀數逐項列出 $A$、$r$、$\theta$、$|H|$ 的計算,與前面的文字推導一致。四個面板對應推導的四個步驟:既定値はスライドの例である($M = 5$、$\omega_0 = 0.1\pi$、$\omega_1 = 0.5\pi$)。上の読み出しは $A$、$r$、$\theta$、$|H|$ の計算を一項ずつ並べており、前の文章による導出と一致する。四つのパネルは導出の四つの段階に対応する。
① The impulse response and where $H$ comes from① 脈衝響應與 $H$ 的由來① インパルス応答と $H$ の由来
② The real part $A(\omega)$ and $r$② 實數部分 $A(\omega)$ 與 $r$② 実数部分 $A(\omega)$ と $r$
Green segments are $A \gt 0$, red segments are $A \lt 0$, the grey dashed lines are the zeros $\omega = 2\pi k/M$, and the $r$ of each segment is marked above that segment:綠色區段 $A \gt 0$,紅色區段 $A \lt 0$,灰色虛線是零點 $\omega = 2\pi k/M$,每一段上方標示該段的 $r$:緑の区間が $A \gt 0$、赤の区間が $A \lt 0$ であり、灰色の破線は零点 $\omega = 2\pi k/M$ である。各区間の $r$ はその区間の上に示してある。
③ Magnitude and phase③ 大小與相位③ 振幅と位相
④ The time-domain output④ 時域輸出④ 時間領域の出力
The blue input has a pronounced fast oscillation (the $0.5\pi$ component), the red output keeps almost nothing but the smooth low-frequency component, matching the picture of slide 3. The first 4 samples, on the orange background, are the transient; after that the red line and the grey dashed line (the steady-state formula $0.904\cos(0.1\pi(n-2)) + 0.2\cos(0.5\pi n)$) coincide completely.藍色輸入有明顯的快速振盪($0.5\pi$ 成分),紅色輸出幾乎只剩平滑的低頻成分,對應投影片 3 的圖。橘色底色的前 4 個樣本是暫態;之後紅線與灰色虛線(穩態公式 $0.904\cos(0.1\pi(n-2)) + 0.2\cos(0.5\pi n)$)完全重合。青い入力には速い振動($0.5\pi$ の成分)がはっきり見えるが、赤い出力にはなめらかな低周波の成分がほぼ残るだけであり、スライド 3 の図に対応する。橙色の地の上にある最初の 4 標本が過渡である。それ以降は赤い線と灰色の破線(定常状態の式 $0.904\cos(0.1\pi(n-2)) + 0.2\cos(0.5\pi n)$)が完全に重なる。
Worth trying可以試試試してみるとよいこと
The simplest moving average of all is the two-point one, $h[n] = [\,0.5,\ 0.5\,]$. Putting $M=2$ into Step 2 above gives a result clean enough to be worth deriving a second time from nothing:最簡單的移動平均就是兩點的那一個,$h[n] = [\,0.5,\ 0.5\,]$。把 $M=2$ 代進上面的步驟 2,會得到一個乾淨到值得從頭再推一次的結果:もっとも単純な移動平均は 2 点のもの、すなわち $h[n] = [\,0.5,\ 0.5\,]$ である。上のステップ 2 に $M=2$ を入れると、何もないところからもう一度導く価値があるほど簡潔な結果が得られる。
This formula is reached in three steps: write down $h[n]$, substitute it into the definition of the DTFT, then tidy up with the manoeuvre of factoring out the half-angle.這條式子分三步得到:寫出 $h[n]$、代入 DTFT 定義、再用「提出半角」整理。この式は三つの段階で得られる。$h[n]$ を書き下し、DTFT の定義に代入し、最後に「半角をくくり出す」という操作で整理するのである。
| Result結果結果 | What it means意義意味 | |
|---|---|---|
| Magnitude $\lvert H\rvert$大小 $\lvert H\rvert$振幅 $\lvert H\rvert$ | $\cos(\omega/2)$ | gain close to 1 at low frequency, 0 at $\omega = \pi$: lowpass低頻增益接近 1,$\omega = \pi$ 時為 0:低通低域では利得が 1 に近く、$\omega = \pi$ では 0 になる。すなわち低域通過である |
| Phase $\theta$相位 $\theta$位相 $\theta$ | $-\omega/2$ | linear phase; every frequency is delayed by 0.5 of a sample線性相位,所有頻率都延遲 0.5 個樣本線形位相であり、どの周波数も 0.5 標本 だけ遅れる |
Alternatively, take the z-transform of $y[n] = h[n] \circledast x[n]$: convolution becomes multiplication, $Y(z) = H(z)X(z)$. The transfer function(轉移函數) $H(z)$ can be obtained two equivalent ways:另一條路,是對 $y[n] = h[n] \circledast x[n]$ 取 z 轉換:摺積變成相乘,$Y(z) = H(z)X(z)$。轉移函數 $H(z)$ 可以用兩種等價的方式得到:別の道として、$y[n] = h[n] \circledast x[n]$ の z 変換をとる。畳み込みは積になり、$Y(z) = H(z)X(z)$ である。伝達関数 $H(z)$ は等価な二つの方法で得られる。
For an FIR filter with $h[n] = \{h_0, h_1, \dots, h_{M-1}\}$ the two views coincide trivially: $p_k = h_k$, $d_0 = 1$, and the difference equation is $y[n] = \sum_{k=0}^{M-1} h_k\,x[n-k]$.對於 $h[n] = \{h_0, h_1, \dots, h_{M-1}\}$ 的 FIR 濾波器,兩種觀點會平凡地重合:$p_k = h_k$、$d_0 = 1$,差分方程式就是 $y[n] = \sum_{k=0}^{M-1} h_k\,x[n-k]$。$h[n] = \{h_0, h_1, \dots, h_{M-1}\}$ の FIR フィルタでは、二つの見方は自明に一致する。$p_k = h_k$、$d_0 = 1$ であり、差分方程式は $y[n] = \sum_{k=0}^{M-1} h_k\,x[n-k]$ となる。
This slide ties three things together: the difference equation (LCCDE), the z-transform (ZT) and the frequency response. The bridge in the middle is the transfer function $H(z)$.這張投影片把三樣東西連起來:差分方程(LCCDE)、Z 轉換(ZT)、頻率響應。中間的橋樑就是轉移函數 $H(z)$。このスライドは三つのものをつなぐ。差分方程式(LCCDE)、z 変換(ZT)、周波数応答である。その間の橋渡しをするのが伝達関数 $H(z)$ である。
① LCCDE: the time-domain description. A linear constant-coefficient difference equation(線性常係數差分方程):① LCCDE:系統的時域描述。線性常係數差分方程(Linear Constant-Coefficient Difference Equation):① LCCDE:システムの時間領域での記述。線形定係数差分方程式(Linear Constant-Coefficient Difference Equation):
The left side is a weighted sum of the output and its past values, the right side a weighted sum of the input and its past values. Most practical digital filters can be written in this form.左邊是輸出及其過去值的加權和,右邊是輸入及其過去值的加權和。大多數實際的數位濾波器都能寫成這個形式。左辺は出力とその過去の値の重み付き和、右辺は入力とその過去の値の重み付き和である。実用的なディジタルフィルタのほとんどは、この形に書ける。
| System系統システム | Difference equation差分方程差分方程式 | Coefficients係數係数 |
|---|---|---|
| first-order IIR一階 IIR1 次 IIR | $y[n] - 0.5\,y[n-1] = x[n]$ | $d_0 = 1,\ d_1 = -0.5$; $p_0 = 1$ |
| MA filter ($M = 5$)MA 濾波器($M = 5$)MA フィルタ($M = 5$) | $y[n] = \frac{1}{5}\big(x[n] + \cdots + x[n-4]\big)$ | $d_0 = 1$; $p_0 = \cdots = p_4 = \frac{1}{5}$ |
When only $d_0$ is nonzero (no output feedback) the system is FIR; when some $d_k \ne 0$ with $k \ge 1$ it is IIR.只有 $d_0$ 非零時(輸出不回授),就是 FIR;有 $d_k \ne 0$($k \ge 1$)時,就是 IIR。$d_0$ だけが非零(出力の帰還がない)なら FIR であり、$k \ge 1$ で $d_k \ne 0$ となるものがあれば IIR である。
② Take the z-transform. The definition is② 取 Z 轉換。Z 轉換的定義:② z 変換をとる。z 変換の定義は次のとおり。
$$X(z) = \sum_{n=-\infty}^{\infty} x[n]\,z^{-n}$$and the key property is a delay of $k$ samples ↔ multiplication by $z^{-k}$:關鍵性質是延遲 $k$ 點 ↔ 乘上 $z^{-k}$:鍵となる性質は$k$ サンプルの遅延 ↔ $z^{-k}$ を掛けることである。
So every $y[n-k]$ in the difference equation becomes $z^{-k}Y(z)$ and every $x[n-k]$ becomes $z^{-k}X(z)$:所以差分方程中每個 $y[n-k]$ 變成 $z^{-k}Y(z)$、每個 $x[n-k]$ 變成 $z^{-k}X(z)$:したがって差分方程式の各 $y[n-k]$ は $z^{-k}Y(z)$ に、各 $x[n-k]$ は $z^{-k}X(z)$ になる。
$$\sum_k d_k\,z^{-k}\,Y(z) = \sum_k p_k\,z^{-k}\,X(z)$$The time-domain difference equation has become an algebraic equation in $z$.時域的差分方程,變成 $z$ 域的代數方程。時間領域の差分方程式が、$z$ 領域の代数方程式になった。
③ Solve for $Y(z)$. $Y(z)$ and $X(z)$ do not depend on $k$, so they come out of the sums:③ 解出 $Y(z)$。$Y(z)$ 和 $X(z)$ 與 $k$ 無關,可以提出求和:③ $Y(z)$ について解く。$Y(z)$ と $X(z)$ は $k$ に依存しないので、和の外へ出せる。
$$Y(z)\sum_k d_k\,z^{-k} = X(z)\sum_k p_k\,z^{-k} \quad\Longrightarrow\quad Y(z) = \underbrace{\frac{\sum_k p_k\,z^{-k}}{\sum_k d_k\,z^{-k}}}_{H(z)}\,X(z)$$These two slides are one complete IIR example: write $H(z)$ from the difference equation → factorize to find the poles and zeros → read the ROC and stability off the poles.這兩張投影片是一個完整的 IIR 例子:從差分方程寫出 $H(z)$ → 因式分解找出極零點 → 由極點判斷 ROC 與穩定性。この 2 枚のスライドは一つの完結した IIR の例である。差分方程式から $H(z)$ を書く → 因数分解して極と零点を求める → 極から ROC と安定性を判定する。
Slide 1, step 1 — move every $y$ to the left to get the standard LCCDE form. The $y$ coefficients change sign as they cross over:投影片 1,第 1 步 — 把所有 $y$ 移到左邊,整理成 LCCDE 的標準形式。移項時 $y$ 的係數變號:スライド 1、ステップ 1 — すべての $y$ を左辺へ移し、LCCDE の標準形に整える。移項すると $y$ の係数は符号が変わる。
$$y[n] - 1.3\,y[n-1] + 1.04\,y[n-2] - 0.222\,y[n-3] = x[n-1] - 1.2\,x[n-2] + x[n-3]$$| $k$ | 0 | 1 | 2 | 3 |
|---|---|---|---|---|
| $d_k$ ($y$ side)$d_k$($y$ 側)$d_k$($y$ 側) | 1 | −1.3 | 1.04 | −0.222 |
| $p_k$ ($x$ side)$p_k$($x$ 側)$p_k$($x$ 側) | 0 | 1 | −1.2 | 1 |
Step 2 — apply $H(z) = \sum p_k z^{-k} \big/ \sum d_k z^{-k}$:第 2 步 — 套用 $H(z) = \sum p_k z^{-k} \big/ \sum d_k z^{-k}$:ステップ 2 — $H(z) = \sum p_k z^{-k} \big/ \sum d_k z^{-k}$ を当てはめる。
$$H(z) = \frac{z^{-1} - 1.2\,z^{-2} + z^{-3}}{1 - 1.3\,z^{-1} + 1.04\,z^{-2} - 0.222\,z^{-3}}$$Why do the complex roots come in conjugate pairs? All the coefficients are real, so every complex root of the polynomial appears together with its conjugate. Multiplying out a conjugate pair gives a quadratic with real coefficients:為什麼複數根成對出現(共軛)?係數都是實數,多項式的複數根一定以共軛對出現。把一對共軛根乘開,得到實係數的二次式:なぜ複素根は共役の対で現れるのか。係数がすべて実数なので、多項式の複素根は必ず共役と対になって現れる。共役の対を掛け合わせると、実係数の 2 次式になる。
Check the numerator: $\text{Re}(\zeta_0) = 0.6$, $|\zeta_0|^2 = 0.36 + 0.64 = 1$:驗證分子:$\text{Re}(\zeta_0) = 0.6$,$|\zeta_0|^2 = 0.36 + 0.64 = 1$:分子の検算:$\text{Re}(\zeta_0) = 0.6$、$|\zeta_0|^2 = 0.36 + 0.64 = 1$。
$$z^{-1}\big(1 - 1.2\,z^{-1} + z^{-2}\big) = z^{-1} - 1.2\,z^{-2} + z^{-3} \quad\checkmark$$Check the denominator: $\text{Re}(\lambda_1) = 0.5$, $|\lambda_1|^2 = 0.25 + 0.49 = 0.74$:驗證分母:$\text{Re}(\lambda_1) = 0.5$,$|\lambda_1|^2 = 0.25 + 0.49 = 0.74$:分母の検算:$\text{Re}(\lambda_1) = 0.5$、$|\lambda_1|^2 = 0.25 + 0.49 = 0.74$。
$$\begin{aligned} (1 - 0.3\,z^{-1})(1 - z^{-1} + 0.74\,z^{-2}) &= 1 - z^{-1} + 0.74\,z^{-2} - 0.3\,z^{-1} + 0.3\,z^{-2} - 0.222\,z^{-3} \\ &= 1 - 1.3\,z^{-1} + 1.04\,z^{-2} - 0.222\,z^{-3} \quad\checkmark \end{aligned}$$Slide 2 — the pole-zero plot and the ROC.投影片 2 — 極零點圖與 ROC。スライド 2 — 極零配置図と ROC。
| Location位置位置 | Magnitude大小大きさ | Angle角度角度 | |
|---|---|---|---|
| zeros $\zeta_0, \zeta_0^*$ (○)零點 $\zeta_0, \zeta_0^*$(○)零点 $\zeta_0, \zeta_0^*$(○) | $0.6 \pm j0.8$ | 1 (on the unit circle)1(在單位圓上)1(単位円上) | $\pm 0.295\pi$ |
| pole $\lambda_0$ (×)極點 $\lambda_0$(×)極 $\lambda_0$(×) | $0.3$ | 0.3 | 0 |
| poles $\lambda_1, \lambda_1^*$ (×)極點 $\lambda_1, \lambda_1^*$(×)極 $\lambda_1, \lambda_1^*$(×) | $0.5 \pm j0.7$ | 0.860 | $\pm 0.303\pi$ |
About the $z^{-1}$ in the numerator. Multiply numerator and denominator by $z^3$ to rewrite in positive powers of $z$:關於分子的 $z^{-1}$:分子分母同乘 $z^3$,改寫成 $z$ 的正次方:分子の $z^{-1}$ について。分子と分母に $z^3$ を掛け、$z$ の正のべきで書き直す。
$$H(z) = \frac{(z - \zeta_0)(z - \zeta_0^*)}{(z - \lambda_0)(z - \lambda_1)(z - \lambda_1^*)}$$The numerator has degree 2 and the denominator degree 3, so there is no pole at the origin; the extra degree is a zero at $z = \infty$. The $z^{-1}$ is an overall one-sample delay: the right-hand side of the difference equation starts at $x[n-1]$, so $h[0] = 0$.分子 2 次、分母 3 次,所以原點沒有極點;多出的一個次數代表一個在 $z = \infty$ 的零點。$z^{-1}$ 的意義是整體延遲 1 個樣本(差分方程右邊從 $x[n-1]$ 開始,所以 $h[0] = 0$)。分子は 2 次、分母は 3 次なので、原点に極はない。余った 1 次は $z = \infty$ にある零点を表す。$z^{-1}$ は全体を 1 標本遅らせることを意味する(差分方程式の右辺は $x[n-1]$ から始まるので $h[0] = 0$)。
ROC. The system is causal, so the ROC lies outside the outermost pole:ROC:系統是因果的,ROC 在最外側極點之外:ROC。システムは因果的なので、ROC は最も外側の極の外側にある。
$$\text{ROC:}\ |z| > \max_i |\lambda_i| = |\lambda_1| \approx 0.860$$The orange region on the slide is this ROC: everything outside the circle of radius 0.86. The unit circle $|z| = 1$ lies inside it, so the system is stable. Equivalently, all the poles are inside the unit circle.投影片中的橘色區域就是 ROC:半徑 0.86 的圓外全部。單位圓 $|z| = 1$ 在其中,所以系統穩定。等價地說,所有極點都在單位圓內。スライドの橙色の領域がこの ROC、すなわち半径 0.86 の円の外側全体である。単位円 $|z| = 1$ はその中にあるので、システムは安定である。言い換えれば、すべての極が単位円の内側にある。
Predicting the frequency response from the poles and zeros.從極零點預測頻率響應。極と零点から周波数応答を予測する。
The interactive figure. Drag $\omega$ and the purple point moves around the unit circle; the middle panels show the frequency response and the bottom panel the impulse response computed recursively from the difference equation.互動圖。拖動 $\omega$,紫點沿單位圓移動;中間是頻率響應,下方是由差分方程遞迴算出的脈衝響應。インタラクティブな図。$\omega$ をドラッグすると紫の点が単位円上を動く。中央は周波数応答、下は差分方程式から再帰的に求めたインパルス応答である。
The two statements live at different levels:兩者講的是不同層次的事:両者は別の層の話である。
| Meaning意思意味 | Property of屬於何の性質か | |
|---|---|---|
| $z$ is in the ROC$z$ 在 ROC 內$z$ が ROC 内にある | at this value of $z$ the series $\sum h[n]z^{-n}$ converges在這個 $z$ 值,級數 $\sum h[n]z^{-n}$ 收斂この $z$ の値で級数 $\sum h[n]z^{-n}$ が収束する | one point某一個點ある一つの点 |
| the system is stable系統穩定システムが安定 | $\sum_n \lvert h[n]\rvert < \infty$ (bounded input → bounded output)$\sum_n \lvert h[n]\rvert < \infty$(有界輸入 → 有界輸出)$\sum_n \lvert h[n]\rvert < \infty$(有界入力 → 有界出力) | the whole system整個系統システム全体 |
There is only one link between them:兩者的連結只有一條:両者を結ぶのはただ一つである。
Other points of the ROC say nothing about stability.只看 ROC 裡的其他點,無法判斷穩定性。ROC の他の点を見ても、安定性は判断できない。
| Choice選擇選択 | $h[n]$ | ROC | Contains u.c.?包含單位圓?単位円を含む? | Stable?穩定?安定? | Causal?因果?因果的? |
|---|---|---|---|---|---|
| right-sided右邊序列右側系列 | $1.03^n\,u[n]$ | $\lvert z\rvert > 1.03$ | no否いいえ | unstable不穩定不安定 | yes是はい |
| left-sided左邊序列左側系列 | $-1.03^n\,u[-n-1]$ | $\lvert z\rvert < 1.03$ | yes是はい | stable穩定安定 | no否いいえ |
| Condition條件条件 | How to judge判斷方法判定方法 |
|---|---|
| a general system is stable一般系統穩定一般のシステムが安定 | does the ROC contain the unit circle?ROC 是否包含單位圓ROC が単位円を含むか |
| a causal system is stable因果系統穩定因果的システムが安定 | are all poles inside the unit circle? (the ROC is then outside the outermost pole and contains the unit circle automatically)所有極點是否都在單位圓內(此時 ROC 是最外側極點之外,自動包含單位圓)すべての極が単位円の内側にあるか(このとき ROC は最も外側の極より外であり、自動的に単位円を含む) |
| some $z$ is in the ROC某個 $z$ 在 ROC 內ある $z$ が ROC 内にある | only says $H(z)$ is defined at that point; unrelated to stability只表示 $H(z)$ 在那一點有定義,與穩定性無關$H(z)$ がその点で定義されることを示すだけで、安定性とは無関係 |
The DTFT is the z-transform evaluated on the unit circle(單位圓):DTFT 就是把 z 轉換在單位圓上求值:DTFT とは、z 変換を単位円上で評価したものである。
(This requires the ROC of $H(z)$ to contain the unit circle, i.e. the system to be stable.) Two checks with the systems of 10.2.1:(前提是 $H(z)$ 的收斂區包含單位圓,也就是系統穩定。)用 10.2.1 的兩個系統驗證:(ただし $H(z)$ の収束領域が単位円を含むこと、すなわちシステムが安定であることが前提である。)10.2.1 の二つのシステムで確かめる。
1. The definition feeds in a sinusoid. Complex exponentials are eigenfunctions of an LTI system (10.1.1): the input $e^{j\omega n}$ comes out as the same exponential times a complex number,1. 頻率響應的定義:輸入正弦波。複指數是 LTI 系統的特徵函數(10.1.1):輸入 $e^{j\omega n}$,輸出還是同一個頻率的複指數,只是乘上一個複數:1. 周波数応答の定義は正弦波を入れること。複素指数は LTI システムの固有関数である(10.1.1)。入力 $e^{j\omega n}$ は、同じ指数に複素数が掛かっただけのものとして出てくる。
$$x[n] = e^{j\omega n} \;\Rightarrow\; y[n] = \sum_k h[k]e^{j\omega(n-k)} = \underbrace{\Big(\sum_k h[k]e^{-j\omega k}\Big)}_{H(e^{j\omega})} e^{j\omega n}$$Its magnitude is the gain and its angle the phase shift. Compared with $H(z) = \sum_k h[k]z^{-k}$, this is $H(z)$ at $z = e^{j\omega}$ — and $|e^{j\omega}| = 1$ puts that point on the unit circle.這個複數 $H(e^{j\omega})$ 的大小是增益、角度是相移。和 $H(z) = \sum_k h[k]z^{-k}$ 比較,它就是把 $z = e^{j\omega}$ 代進去,而 $|e^{j\omega}| = 1$,正好落在單位圓上。その大きさが利得、角度が位相のずれである。$H(z) = \sum_k h[k]z^{-k}$ と比べると、これは $z = e^{j\omega}$ での $H(z)$ であり、$|e^{j\omega}| = 1$ だからその点は単位円上にある。
2. Off the circle, $z$ is not a pure frequency. Write a general $z = re^{j\omega}$; then $z^n = r^n e^{j\omega n}$:2. 單位圓以外的 $z$ 不是「純頻率」。一般的 $z$ 寫成 $z = re^{j\omega}$,那麼 $z^n = r^n e^{j\omega n}$:2. 単位円の外の $z$ は「純粋な周波数」ではない。一般の $z$ を $z = re^{j\omega}$ と書くと $z^n = r^n e^{j\omega n}$ となる。
| Where $z$ isz 的位置$z$ の位置 | What $z^n$ looks like$z^n$ 的樣子$z^n$ の形 | Meaning意義意味 |
|---|---|---|
| $|z| = 1$ (on the circle)$|z| = 1$(單位圓上)$|z| = 1$(単位円上) | sinusoid of constant amplitude振幅固定的正弦波振幅一定の正弦波 | pure frequency $\omega$純粹的頻率 $\omega$純粋な周波数 $\omega$ |
| $|z| < 1$ (inside)$|z| < 1$(圓內)$|z| < 1$(円の内側) | decaying sinusoid逐漸衰減的正弦波減衰する正弦波 | transient, damped oscillation暫態、阻尼振盪過渡、減衰振動 |
| $|z| > 1$ (outside)$|z| > 1$(圓外)$|z| > 1$(円の外側) | growing sinusoid逐漸增長的正弦波増大する正弦波 | diverging oscillation發散的振盪発散する振動 |
Only on the circle is the signal a steady, lasting sinusoid, so “how much is frequency $\omega$ amplified?” can only be read there. Put another way, the $z$-transform is a DTFT taken after weighting by $r^{-n}$:只有在單位圓上,訊號才是「穩定、持續的正弦波」;所以要問「頻率 $\omega$ 的正弦波被放大多少」,就只能看單位圓上的值。換個角度看,Z 轉換其實是「先乘上 $r^{-n}$ 再做 DTFT」:信号が安定して持続する正弦波になるのは単位円上だけなので、「周波数 $\omega$ の正弦波がどれだけ増幅されるか」はそこでしか読めない。言い換えれば、$z$ 変換は $r^{-n}$ で重み付けしてから DTFT をとったものである。
$$X(re^{j\omega}) = \sum_n \big(x[n]\,r^{-n}\big)e^{-j\omega n}$$At $r = 1$ the weight disappears and what remains is the spectrum (DTFT) of the original signal.$r = 1$ 時,乘上的權重消失,剩下的就是原始訊號的頻譜(DTFT)。$r = 1$ では重みが消え、残るのは元の信号のスペクトル(DTFT)である。
3. The continuous-time Bode plot does the same thing. It plots $H(s)$ on the imaginary axis $s = j\Omega$, for the same reason: $e^{j\Omega t}$ is a pure sinusoid, while $e^{\sigma t}e^{j\Omega t}$ decays or grows. After sampling, the two planes are linked by $z = e^{sT}$ ($T$ the sampling period):3. 和連續時間 Bode 圖的對應。連續時間的 Bode 圖畫的是 $H(s)$ 在虛軸 $s = j\Omega$ 上的值,理由完全一樣:$e^{j\Omega t}$ 是純正弦波,而 $e^{\sigma t}e^{j\Omega t}$ 會衰減或增長。取樣之後,s 平面和 z 平面的關係是 $z = e^{sT}$($T$ 是取樣週期):3. 連続時間のボード線図も同じことをしている。ボード線図は虚軸 $s = j\Omega$ 上の $H(s)$ を描く。理由は同じで、$e^{j\Omega t}$ は純粋な正弦波だが、$e^{\sigma t}e^{j\Omega t}$ は減衰または増大する。サンプリング後、二つの平面は $z = e^{sT}$($T$ はサンプリング周期)で結ばれる。
| $s$-planes 平面$s$ 平面 | $z$-planez 平面$z$ 平面 |
|---|---|
| imaginary axis $s = j\Omega$虛軸 $s = j\Omega$虚軸 $s = j\Omega$ | unit circle $z = e^{j\Omega T}$單位圓 $z = e^{j\Omega T}$単位円 $z = e^{j\Omega T}$ |
| left half-plane (stable)左半平面(穩定)左半平面(安定) | inside the unit circle單位圓內単位円の内側 |
| right half-plane (unstable)右半平面(不穩定)右半平面(不安定) | outside the unit circle單位圓外単位円の外側 |
| DC, $\Omega = 0$DC($\Omega = 0$)直流、$\Omega = 0$ | $z = 1$ |
| Nyquist, $\Omega = \pi/T$Nyquist 頻率($\Omega = \pi/T$)ナイキスト、$\Omega = \pi/T$ | $z = -1$ |
So the unit circle plays the role of the imaginary axis, with digital and analog frequency related by $\omega = \Omega T$.所以離散時間系統的單位圓,就扮演連續時間系統虛軸的角色;數位頻率和類比頻率的關係是 $\omega = \Omega T$。つまり離散時間の単位円は連続時間の虚軸の役割を果たし、ディジタル周波数とアナログ周波数は $\omega = \Omega T$ で結ばれる。
4. One lap of the circle sweeps every frequency.4. 沿單位圓走一圈等於掃過所有頻率。4. 単位円を一周するとすべての周波数を掃引する。
$\omega$ and $\omega + 2\pi$ are the same point, so the response is periodic and $0$ to $\pi$ (or $-\pi$ to $\pi$) is enough. That is why digital responses are drawn on a linear axis up to $\pi$, not on the logarithmic axis running to infinity of a continuous-time Bode plot.$\omega$ 和 $\omega + 2\pi$ 是同一個點,所以頻率響應是週期性的,只要看 $0$ 到 $\pi$(或 $-\pi$ 到 $\pi$)就夠了。這也是為什麼數位濾波器的頻率響應通常用線性頻率軸畫到 $\pi$,而不像連續時間 Bode 圖那樣用對數頻率軸一路延伸到無限大。$\omega$ と $\omega + 2\pi$ は同じ点なので応答は周期的であり、$0$ から $\pi$(または $-\pi$ から $\pi$)を見れば足りる。ディジタルフィルタの応答を、連続時間のボード線図のように無限大まで延びる対数軸ではなく、$\pi$ までの線形軸で描くのはこのためである。
5. Precondition: the unit circle must lie in the ROC. Only then does $H(e^{j\omega})$ exist — the system is stable and the DTFT converges (see the ROC view of the 3-D surface in 10.2.4). Reading the response along the unit circle therefore already assumes a stable system.5. 前提:單位圓必須在 ROC 內。只有當單位圓落在收斂區域(ROC)內,$H(e^{j\omega})$ 才存在,也就是系統穩定、DTFT 收斂(見 10.2.4 3D 曲面的 ROC 檢視)。所以「沿著單位圓看頻率變化」這件事,本身就隱含了系統是穩定的。5. 前提:単位円が ROC に含まれること。そのときに限り $H(e^{j\omega})$ が存在する。すなわちシステムが安定で DTFT が収束する(10.2.4 の 3 次元曲面の ROC 表示を参照)。したがって単位円に沿って応答を読むこと自体が、システムの安定を前提としている。
Factor $H(z)$ into first-order terms and evaluate at $z = e^{j\omega}$:把 $H(z)$ 因式分解成一階項,再代入 $z = e^{j\omega}$:$H(z)$ を 1 次の項に因数分解し、$z = e^{j\omega}$ を代入する。
Taking magnitude and argument term by term:逐項取振幅與輻角:項ごとに振幅と偏角をとる。
Each factor $(e^{j\omega} - \nu)$ is a vector drawn from the root $\nu$ (pole or zero) to the point $e^{j\omega}$ moving along the unit circle. The whole frequency response is a product/ratio of these vectors:每一個因式 $(e^{j\omega} - \nu)$,都是一個從根 $\nu$(極點或零點)指向沿單位圓移動的點 $e^{j\omega}$ 的向量。整個頻率響應,就是這些向量的乘積與比值:各因子 $(e^{j\omega} - \nu)$ は、根 $\nu$(極または零点)から、単位円上を動く点 $e^{j\omega}$ へ引いたベクトルである。周波数応答全体は、これらのベクトルの積と比にほかならない。
A worked case — the peaking filter of the lecture. A conjugate pole pair at radius 0.9 and a conjugate zero pair at radius 0.8, all at angles $\pm 0.3\pi$, so each zero sits just inside its pole on the same ray:實例:講義中的峰值濾波器。一對共軛極點在半徑 0.9、一對共軛零點在半徑 0.8,角度都是 $\pm 0.3\pi$,所以每個零點都落在同一條射線上、緊貼在對應極點的內側:実例:講義のピーキングフィルタ。共役な極の対が半径 0.9、共役な零点の対が半径 0.8 にあり、角度はどれも $\pm 0.3\pi$ である。したがって各零点は同じ半直線上で、対応する極のすぐ内側にある。
Why detour through the $z$ domain? $H(z)$ is a ratio of two polynomials, so it factors into zeros $\zeta_k$ (numerator is 0, $|H|$ drops to 0) and poles $\lambda_k$ (denominator is 0, $|H| \to \infty$). The frequency response is just the $|H(z)|$ seen while walking once round the unit circle: it shrinks near a zero and swells near a pole. So the pole and zero positions alone sketch the frequency response. The poles must lie inside the unit circle for stability — and the closer a pole is to the circle, the longer the transient lasts (10.1.3).為什麼要繞道 $z$ 域?$H(z)$ 是兩個多項式的比值,可以因式分解成零點 $\zeta_k$(分子為 0,$|H|$ 掉到 0)與極點 $\lambda_k$(分母為 0,$|H|$ 趨向無限大)。頻率響應就是沿著單位圓走一圈時看到的 $|H(z)|$:經過零點附近會變小,經過極點附近會變大。所以只要看零點和極點的位置,就能大致畫出頻率響應。極點必須在單位圓內,系統才穩定;極點越靠近單位圓,暫態越長(10.1.3)。なぜ $z$ 領域を経由するのか。$H(z)$ は二つの多項式の比なので、零点 $\zeta_k$(分子が 0、$|H|$ は 0 に落ちる)と極 $\lambda_k$(分母が 0、$|H| \to \infty$)に因数分解できる。周波数応答とは、単位円を一周しながら見た $|H(z)|$ にすぎない。零点の近くでは小さくなり、極の近くでは大きくなる。だから極と零点の位置だけで、周波数応答のおおよその形が描ける。安定であるためには極は単位円の内側になければならない。極が円に近いほど、過渡は長く続く(10.1.3)。
The factor $H(e^{j\omega}) = 1 - re^{j\theta}e^{-j\omega}$ (one zero at $z = re^{j\theta}$) is the elementary building block — every rational response is a product/ratio of such terms. Its log magnitude, phase and group delay all concentrate their activity around $\omega = \theta$:因式 $H(e^{j\omega}) = 1 - re^{j\theta}e^{-j\omega}$(一個位於 $z = re^{j\theta}$ 的零點)是最基本的元件——任何有理響應都是這種項的乘積或比值。它的對數振幅、相位與群延遲,活動全都集中在 $\omega = \theta$ 附近:因子 $H(e^{j\omega}) = 1 - re^{j\theta}e^{-j\omega}$($z = re^{j\theta}$ に零点がひとつ)は最も基本的な構成要素である。あらゆる有理応答はこうした項の積と比で書ける。その対数振幅・位相・群遅延は、いずれも $\omega = \theta$ の付近に働きが集中する。
Filtering means separating information in frequency with a constructed $H(e^{j\omega})$: given $x[n] = A\cos(\omega_1 n) + B\cos(\omega_2 n)$ where only the $\omega_1$ part is wanted, build $|H(e^{j\omega_1})| \approx 1$ and $|H(e^{j\omega_2})| \approx 0$, so that $y[n] \approx A\cos(\omega_1 n + \theta(\omega_1))$.所謂濾波,就是用一個刻意造出來的 $H(e^{j\omega})$ 在頻率上分離資訊:給定 $x[n] = A\cos(\omega_1 n) + B\cos(\omega_2 n)$,而我們只要 $\omega_1$ 那一部分,那就把 $|H(e^{j\omega_1})| \approx 1$、$|H(e^{j\omega_2})| \approx 0$ 做出來,於是 $y[n] \approx A\cos(\omega_1 n + \theta(\omega_1))$。フィルタリングとは、意図的に作った $H(e^{j\omega})$ によって周波数上で情報を分離することである。$x[n] = A\cos(\omega_1 n) + B\cos(\omega_2 n)$ のうち $\omega_1$ の成分だけがほしいなら、$|H(e^{j\omega_1})| \approx 1$ かつ $|H(e^{j\omega_2})| \approx 0$ となるように作ればよい。すると $y[n] \approx A\cos(\omega_1 n + \theta(\omega_1))$ となる。
Only two tools are needed here: factoring out $e^{-j\omega}$, and Euler's formula.這一步只用到兩個工具:把 $e^{-j\omega}$ 提出來,以及 Euler 公式。ここで使う道具は二つだけである。$e^{-j\omega}$ をくくり出すことと、オイラーの公式である。
1. Factor out $e^{-j\omega}$. The three exponents are $0$, $-\omega$ and $-2\omega$. Divide every term by $e^{-j\omega}$ and move that factor outside the bracket:1. 提出 $e^{-j\omega}$。原式三項的指數分別是 $0$、$-\omega$、$-2\omega$。把每一項都除以 $e^{-j\omega}$,再把它提到括號外:1. $e^{-j\omega}$ をくくり出す。三つの項の指数はそれぞれ $0$、$-\omega$、$-2\omega$ である。各項を $e^{-j\omega}$ で割り、その因子を括弧の外に出す。
$$\alpha + \beta e^{-j\omega} + \alpha e^{-j2\omega} = e^{-j\omega}\Big(\alpha\, e^{j\omega} + \beta + \alpha\, e^{-j\omega}\Big)$$Term by term:逐項檢查:項ごとに確かめると:
| Original term原本的項もとの項 | After dividing by $e^{-j\omega}$除以 $e^{-j\omega}$ 之後$e^{-j\omega}$ で割った後 |
|---|---|
| $\alpha$ | $\alpha e^{j\omega}$ |
| $\beta e^{-j\omega}$ | $\beta$ |
| $\alpha e^{-j2\omega}$ | $\alpha e^{-j\omega}$ |
Multiplying $e^{-j\omega}$ back into each term inside the bracket returns the original expression — worth checking once by hand.反過來把 $e^{-j\omega}$ 乘回括號內每一項,就會得到原式,可以自己驗證一次。逆に括弧の中の各項へ $e^{-j\omega}$ を掛け戻せばもとの式に戻る。一度は手で確かめておくとよい。
2. Merge with Euler's formula. Put the two $\alpha$ terms inside the bracket together:2. 用 Euler 公式合併。把括號內兩個 $\alpha$ 的項放在一起:2. オイラーの公式でまとめる。括弧の中の二つの $\alpha$ の項を並べる。
$$\alpha\big(e^{j\omega} + e^{-j\omega}\big)$$By Euler's formula:由 Euler 公式:オイラーの公式より:
$$e^{j\omega} = \cos\omega + j\sin\omega,\qquad e^{-j\omega} = \cos\omega - j\sin\omega$$Adding the two, the imaginary parts $\pm j\sin\omega$ cancel exactly:兩式相加,虛部 $\pm j\sin\omega$ 正好抵消:二つを足すと、虚部 $\pm j\sin\omega$ がちょうど打ち消し合う。
$$e^{j\omega} + e^{-j\omega} = 2\cos\omega$$So the bracket collapses to $\beta + 2\alpha\cos\omega$.所以括號內就變成 $\beta + 2\alpha\cos\omega$。したがって括弧の中は $\beta + 2\alpha\cos\omega$ になる。
Inside the bracket there is nothing but $\alpha$, $\beta$ and $\cos\omega$, all of them real — so the whole bracket is real, and the entire frequency dependence of the phase sits in the single factor $e^{-j\omega}$.括號內只有 $\alpha$、$\beta$ 和 $\cos\omega$,全部都是實數,所以整個括號是實數;相位對頻率的依賴,全部集中在 $e^{-j\omega}$ 這一個因子上。括弧の中にあるのは $\alpha$、$\beta$、$\cos\omega$ だけで、どれも実数である。よって括弧全体が実数であり、位相の周波数依存はすべて $e^{-j\omega}$ という一つの因子に収まっている。
Why factor out $e^{-j\omega}$ and not something else? Because $-\omega$ is exactly the middle of the three exponents $0$, $-\omega$, $-2\omega$ — that is, the centre of symmetry of $h[n]$, at $n = 1$. Pulling the centre out leaves the two outer terms as $e^{+j\omega}$ and $e^{-j\omega}$: same exponent magnitude, opposite sign. Since their coefficients are also equal (both $\alpha$), they pair up into $2\cos\omega$ and the imaginary parts vanish.為什麼要提出 $e^{-j\omega}$,而不是其他東西?因為 $-\omega$ 剛好是三個指數 $0$、$-\omega$、$-2\omega$ 的中間值,也就是 $h[n]$ 的對稱中心 $n = 1$。提出中心之後,左右兩邊的項變成 $e^{+j\omega}$ 和 $e^{-j\omega}$,指數大小相同、符號相反;再加上係數相同(都是 $\alpha$),就能配成 $2\cos\omega$,虛部完全抵消。なぜ $e^{-j\omega}$ をくくり出すのか。$-\omega$ が三つの指数 $0$、$-\omega$、$-2\omega$ のちょうど真ん中、すなわち $h[n]$ の対称中心 $n = 1$ に当たるからである。中心をくくり出すと、両端の項は $e^{+j\omega}$ と $e^{-j\omega}$ になる。指数の大きさが同じで符号が逆であり、係数も等しい(ともに $\alpha$)ので、$2\cos\omega$ にまとまり虚部は消える。
1. Impulse response. Only three samples are nonzero: $h[n] = \{\alpha, \beta, \alpha\} = \{6.76, -12.46, 6.76\}$. (The lecture plot puts them at $n = 1, 2, 3$ only because MATLAB indexes from 1; they are $n = 0, 1, 2$ here.)1. 脈衝響應。只有三個非零樣本:$h[n] = \{\alpha, \beta, \alpha\} = \{6.76, -12.46, 6.76\}$。(講義圖上畫在 $n = 1, 2, 3$,只是因為 MATLAB 從 1 開始編號,對應這裡的 $n = 0, 1, 2$。)1. インパルス応答。非零の標本は三つだけである:$h[n] = \{\alpha, \beta, \alpha\} = \{6.76, -12.46, 6.76\}$。(講義の図で $n = 1, 2, 3$ に描かれているのは MATLAB が 1 から番号を振るためで、ここでの $n = 0, 1, 2$ に当たる。)
| Frequency頻率周波数 | Sum計算計算 | $|H|$ | dB |
|---|---|---|---|
| $\omega = 0$ (DC)$\omega = 0$(DC)$\omega = 0$(直流) | $\alpha + \beta + \alpha = 6.76 - 12.46 + 6.76$ | $\approx 1.07$ | $\approx +0.6$ |
| $\omega = \pi$ | $\alpha - \beta + \alpha = 6.76 + 12.46 + 6.76$ | $\approx 25.98$ | $\approx +28.3$ |
How 0.4 is cancelled in the time domain. The filter computes在時域中,0.4 是怎麼被消掉的。濾波器的運算是:時間領域で 0.4 が消える仕組み。フィルタの計算は次のとおりである。
$$y[n] = 6.76\big(x[n] - 2\cos(0.4)\,x[n-1] + x[n-2]\big)$$and for a sinusoid at 0.4 the identity對頻率 0.4 的正弦波,三角恆等式周波数 0.4 の正弦波では、恒等式
$$\cos(0.4n) + \cos(0.4(n-2)) = 2\cos(0.4)\cos(0.4(n-1))$$makes the bracket exactly 0: the samples on either side add up to exactly $2\cos(0.4)$ times the middle one, and the three terms cancel. This is what the zeros on the unit circle at $e^{\pm j0.4}$ mean in the time domain.剛好讓括號變成 0:前後兩個樣本相加,恰好等於中間樣本的 $2\cos(0.4)$ 倍,三項完全抵消。這就是零點落在單位圓 $e^{\pm j0.4}$ 上,在時域中的意義。によって括弧の中がちょうど 0 になる。前後二つの標本の和が中央の標本のちょうど $2\cos(0.4)$ 倍になり、三つの項が打ち消し合う。これが、零点が単位円上の $e^{\pm j0.4}$ にあることの時間領域での意味である。
Large, alternating coefficients. The output is a small result of large numbers subtracted ($6.76 + 6.76 - 12.46 \approx 1$). Such heavy cancellation is sensitive to coefficient errors and input noise — the same fact as the high-frequency boost in the frequency domain.係數很大、正負交錯,代表什麼?輸出是「大數減大數」得到的小結果($6.76 + 6.76 - 12.46 \approx 1$)。這種大量互相抵消的結構,對係數誤差和輸入雜訊都很敏感,和頻域中高頻被放大是同一件事。大きく符号が交互する係数。出力は「大きな数から大きな数を引いた」小さな結果である($6.76 + 6.76 - 12.46 \approx 1$)。このように大きく打ち消し合う構造は係数誤差にも入力雑音にも敏感であり、周波数領域で高域が強められることと同じ事実である。
2. Frequency response (dB). At $\omega = 0.1$ it is about 0 dB ($|H| \approx 1$), so the tone is kept. At $\omega = 0.4$ the curve drops to $-\infty$ dB — the notch made by the zero on the unit circle. Beyond 0.4 it climbs all the way: about $+18$ dB at 1.2 rad (where the lecture plot stops) and $+28$ dB at $\pi$.2. 頻率響應(dB)。$\omega = 0.1$ 處約 0 dB($|H| \approx 1$),訊號保留。$\omega = 0.4$ 處曲線掉到 $-\infty$ dB,這就是單位圓上的零點造成的凹口(notch)。0.4 之後一路上升:到 1.2 rad(講義圖的橫軸只畫到這裡)已約 $+18$ dB,到 $\pi$ 時達 $+28$ dB。2. 周波数応答(dB)。$\omega = 0.1$ ではおよそ 0 dB($|H| \approx 1$)で、信号は残る。$\omega = 0.4$ で曲線は $-\infty$ dB まで落ちる。単位円上の零点によるノッチである。0.4 を過ぎると上がり続け、1.2 rad(講義の図はここまで)で約 $+18$ dB、$\pi$ では $+28$ dB に達する。
3. Input and output. The input is two sinusoids: a slow part with period $2\pi/0.1 \approx 63$ samples carrying a fast ripple with period $2\pi/0.4 \approx 16$ samples. In the output the ripple has vanished, leaving a smooth sinusoid of period about 63, amplitude about 1, delayed by one sample. The odd start at $n \approx 0$–$2$ is the FIR transient: the filter needs 3 samples to fill, and $x[n]$ is taken as 0 for $n < 0$, so the first two outputs are wrong. After that it is in steady state and matches the prediction exactly.3. 輸入輸出。輸入是兩個正弦波相加:慢的部分週期約 $2\pi/0.1 \approx 63$ 個樣本,上面疊加週期約 $2\pi/0.4 \approx 16$ 個樣本的快速抖動。輸出中快速抖動完全消失,只剩週期約 63、振幅約 1 的平滑正弦波,並延遲 1 個樣本。最開頭($n \approx 0$ 到 2)有一個突兀的起點,這是 FIR 的暫態:濾波器需要 3 個樣本才能「填滿」,$n < 0$ 時輸入被視為 0,所以前兩個輸出不正確。之後就是穩態,和預期完全一致。3. 入力と出力。入力は二つの正弦波の和で、周期約 $2\pi/0.1 \approx 63$ サンプルのゆっくりした成分に、周期約 $2\pi/0.4 \approx 16$ サンプルの速い揺れが乗っている。出力では速い揺れが完全に消え、周期約 63、振幅約 1 のなめらかな正弦波が 1 サンプル遅れて残る。最初($n \approx 0$~2)の不自然な立ち上がりは FIR の過渡である。フィルタが「満たされる」には 3 標本必要で、$n < 0$ では入力を 0 とみなすので、最初の二つの出力は正しくない。その後は定常状態で、予測と完全に一致する。
For a sinusoidal input $x[n] = \cos(\omega_0 n)$ we found $y[n] = |H(e^{j\omega_0})|\cos(\omega_0 n + \theta(\omega_0))$. The phase shift can be re-read as a time shift by factoring $\omega_0$ out:對弦波輸入 $x[n] = \cos(\omega_0 n)$,我們已經得到 $y[n] = |H(e^{j\omega_0})|\cos(\omega_0 n + \theta(\omega_0))$。把 $\omega_0$ 提出來,就能把相位平移改讀成時間平移:正弦波入力 $x[n] = \cos(\omega_0 n)$ に対しては $y[n] = |H(e^{j\omega_0})|\cos(\omega_0 n + \theta(\omega_0))$ が得られていた。$\omega_0$ をくくり出せば、位相のずれを時間のずれとして読み直せる。
The minus sign makes a positive $\tau_p$ mean a delay (causal behavior): a falling phase corresponds to the carrier arriving later.那個負號讓正的 $\tau_p$ 代表延遲(因果行為):相位往下掉,對應的就是載波較晚抵達。この負号のおかげで、正の $\tau_p$ が遅れ(因果的なふるまい)を意味する。位相が下がることは、搬送波が遅れて到着することに対応する。
These seven slides answer one question: after a signal passes through a system, how long has it actually been delayed?這七張投影片在回答一個問題:訊號通過系統之後,到底被延遲了多久?この 7 枚のスライドが答えているのは一つの問いである。信号はシステムを通ったあと、実際にどれだけ遅れるのか。
The answer depends on what you look at. Watch the crests of a single sinusoid and you get the phase delay; watch the envelope of a modulated signal — the part that actually carries the information — and you get the group delay. The two agree only when the phase is a straight line through the origin; in general they differ, and a group delay that changes with frequency causes phase distortion. The walk-through follows the slide order, split over §10.3.1–§10.3.3, with three interactive figures along the way.答案取決於你看的是什麼。如果看單一正弦波的波峰,得到的是相位延遲(phase delay);如果看調變訊號的包絡,也就是訊號實際攜帶的資訊,得到的是群延遲(group delay)。只有在相位是一條通過原點的直線時,兩者才相等;一般情況下它們不同,而群延遲隨頻率變化,就會造成相位失真。以下依投影片順序說明,分布在 §10.3.1–§10.3.3,中間穿插三個互動圖。答えは何を見るかによって変わる。単一の正弦波の山を見れば位相遅延(phase delay)が得られ、変調信号の包絡線、すなわち信号が実際に運んでいる情報を見れば群遅延(group delay)が得られる。両者が一致するのは位相が原点を通る直線のときだけであり、一般には異なる。そして群遅延が周波数によって変わると位相歪みが生じる。以下ではスライドの順に、§10.3.1–§10.3.3 にわたって説明し、途中に三つのインタラクティブな図を挟む。
Sinusoid in, sinusoid out. Start from the most basic fact: feed an LTI system a sinusoid at frequency $\omega_0$ and the steady-state output is still a sinusoid at the same frequency, with only its amplitude and phase changed:正弦波輸入,正弦波輸出。先回到最基本的事實:LTI 系統輸入一個頻率 $\omega_0$ 的正弦波,穩態輸出仍是同頻率的正弦波,只是振幅和相位改變:正弦波を入れれば正弦波が出る。最も基本的な事実に戻ろう。LTI システムに周波数 $\omega_0$ の正弦波を入れると、定常出力は同じ周波数の正弦波のままであり、変わるのは振幅と位相だけである。
Here $\theta(\omega) = \angle H(e^{j\omega})$. The derivation takes one step: split the cosine into two complex exponentials,其中 $\theta(\omega) = \angle H(e^{j\omega})$。推導只要一步:把餘弦拆成兩個複指數,ここで $\theta(\omega) = \angle H(e^{j\omega})$ である。導出は一歩で済む。余弦を二つの複素指数に分ける。
Complex exponentials are eigenfunctions of an LTI system, so each is multiplied by $H(e^{j\omega_0})$ and $H(e^{-j\omega_0})$ respectively. If $h[n]$ is real, $H(e^{-j\omega_0}) = H^*(e^{j\omega_0})$, the two terms are conjugates, and their sum is $|H|\cos(\omega_0 n + \theta)$.複指數是 LTI 系統的特徵函數,各自乘上 $H(e^{j\omega_0})$ 與 $H(e^{-j\omega_0})$。若 $h[n]$ 是實數,$H(e^{-j\omega_0}) = H^*(e^{j\omega_0})$,兩項互為共軛,相加後就是 $|H|\cos(\omega_0 n + \theta)$。複素指数は LTI システムの固有関数なので、それぞれ $H(e^{j\omega_0})$ と $H(e^{-j\omega_0})$ を掛けられる。$h[n]$ が実数なら $H(e^{-j\omega_0}) = H^*(e^{j\omega_0})$ であり、二つの項は互いに共役で、足し合わせると $|H|\cos(\omega_0 n + \theta)$ になる。
The slide marks two things here: the gain in blue, and the phase shift or time shift in red. The second is the point of this slide.投影片在這裡標出兩件事:藍色的 gain(增益)與紅色的 phase shift or time shift(相移,或者說時間平移)。後者是本頁重點。スライドはここで二つを示している。青の gain(利得)と、赤の phase shift or time shift(位相のずれ、すなわち時間のずれ)である。このページの要点は後者である。
Rewriting the phase shift as a time shift. Factor out $\omega_0$:把相移改寫成時間平移。把 $\omega_0$ 提出來:位相のずれを時間のずれに書き直す。$\omega_0$ をくくり出す。
So the output is the input shifted along the time axis by $\theta(\omega_0)/\omega_0$ samples. To make the “delay” a positive number, the slide defines所以輸出就是輸入在時間軸上平移了 $\theta(\omega_0)/\omega_0$ 個樣本。為了讓「延遲」是正數,投影片定義:したがって出力は、入力を時間軸上で $\theta(\omega_0)/\omega_0$ サンプルずらしたものである。「遅延」を正の数にするため、スライドは次のように定義する。
Why the minus sign? The output of a causal system always comes later than its input, so the phase usually lags, i.e. $\theta \lt 0$. With the minus sign, $\tau_p \gt 0$ means a delay — which is what the green note on the slide says: “subtraction so positive $\tau_p$ means delay (causal)”.為什麼要加負號?因果系統的輸出一定比輸入晚,相位通常是「落後」的,也就是 $\theta \lt 0$。加上負號後,$\tau_p \gt 0$ 就代表延遲,這是投影片綠字的意思:「subtraction so positive $\tau_p$ means delay (causal)」。なぜ負号を付けるのか。因果的なシステムの出力は必ず入力より遅れるので、位相は通常「遅れ」、つまり $\theta \lt 0$ である。負号を付ければ $\tau_p \gt 0$ が遅延を表す。スライドの緑字「subtraction so positive $\tau_p$ means delay (causal)」はこのことを言っている。
The simplest check is the pure delay $h[n] = \delta[n - n_0]$:最簡單的驗證是純延遲系統 $h[n] = \delta[n - n_0]$:最も簡単な確認は純粋な遅延 $h[n] = \delta[n - n_0]$ である。
— exactly what intuition says.和直覺完全一致。直観とまったく一致する。
Two intrinsic limits of phase delay.相位延遲的兩個本質限制。位相遅延の二つの本質的な限界。
First, it need not be an integer. $\tau_p = 2.7$ means: shift the continuous sinusoid that lies behind the input samples by 2.7 samples, sample it again, and you have the output. A discrete-time signal has no “instant 2.7”, but shifting the sinusoidal curve does make sense.第一,它可以不是整數。$\tau_p = 2.7$ 的意思是:把輸入樣本背後那條連續的正弦曲線平移 2.7 個樣本,再重新取樣,就是輸出。離散時間訊號本身沒有「2.7 這個時刻」,但正弦曲線的平移是有意義的。第一に、整数でなくてもよい。$\tau_p = 2.7$ とは、入力サンプルの背後にある連続な正弦曲線を 2.7 サンプルずらし、それを標本化し直したものが出力だという意味である。離散時間信号そのものには「時刻 2.7」はないが、正弦曲線をずらすことには意味がある。
Second, it is inherently ambiguous. A sinusoid is periodic: shifting it by one period $2\pi/\omega_0$ looks exactly like not shifting it at all. In terms of phase, $\theta$ is fixed only up to an integer multiple of $2\pi$. So from a single steady-state sinusoid you cannot tell a delay of 1 sample from a delay of 1 sample plus one period.第二,它本質上有歧義。正弦波是週期性的,平移一個週期 $2\pi/\omega_0$ 和完全沒平移看起來一模一樣。對應到相位上,$\theta$ 只能確定到 $2\pi$ 的整數倍。因此單看一個穩態正弦波,你無法分辨它延遲了 1 個樣本,還是延遲了 1 加上一個週期。第二に、本質的にあいまいである。正弦波は周期的なので、一周期 $2\pi/\omega_0$ ずらしたものはまったくずらさないものと見分けがつかない。位相で言えば、$\theta$ は $2\pi$ の整数倍の不定性を残してしか決まらない。したがって定常正弦波を一つ見ただけでは、1 サンプル遅れたのか、1 サンプルと一周期遅れたのかを区別できない。
The usual remedy is the unwrapped phase: set $\theta(0) = 0$ and let $\theta(\omega)$ vary continuously with $\omega$, with no $2\pi$ jumps allowed. Only then does $\tau_p$ have a consistent meaning.一般的處理方式是使用展開相位(unwrapped phase):規定 $\theta(0) = 0$,然後讓 $\theta(\omega)$ 隨 $\omega$ 連續變化,不允許 $2\pi$ 的跳躍。這樣算出的 $\tau_p$ 才有一致的意義。普通の対処はアンラップ位相(unwrapped phase)を使うことである。$\theta(0) = 0$ と決め、$\theta(\omega)$ を $\omega$ とともに連続に変化させ、$2\pi$ の跳びを許さない。こうして求めた $\tau_p$ だけが一貫した意味をもつ。
The deeper problem, though, is that a sinusoid that goes on forever, with no start and no end, carries no time marker. You cannot tell which output crest belongs to which input crest. That is exactly why the modulated carrier is brought in later: its envelope supplies the time marker.但更根本的問題是:一個永遠持續、沒有起點也沒有終點的正弦波,沒有任何「時間標記」。你看不出哪一個波峰對應輸入的哪一個波峰。這正是後面要引入「調變載波」的原因:包絡提供了時間標記。しかし、より根本的な問題は、始まりも終わりもなく永遠に続く正弦波には「時刻の目印」がまったくないことである。出力のどの山が入力のどの山に対応するのかがわからない。後で「変調された搬送波」を持ち込むのはまさにこのためであり、包絡線が時刻の目印を与える。
We already derived (Example 10-D)前面已經推過(範例 10-D):すでに導いたとおり(例題 10-D)、
Where $A(\omega) \gt 0$, $\theta(\omega) = -\omega$, so當 $A(\omega) \gt 0$ 時,$\theta(\omega) = -\omega$,所以$A(\omega) \gt 0$ のところでは $\theta(\omega) = -\omega$ なので、
Every frequency is delayed by 1 sample, just as the input/output plots showed earlier (“as observed”).所有頻率都延遲 1 個樣本,和之前輸入輸出圖中看到的完全一致(as observed)。すべての周波数が 1 サンプル遅れる。これは先に入出力の図で見たこと(as observed)とまったく一致する。
A detail the slide skips: the region where $A(\omega) \lt 0$. The slide says “at all frequencies”; strictly, that holds only where $A(\omega) \gt 0$. In our example $A(\omega) = -12.46 + 13.52\cos\omega$, which turns negative for $\omega \gt 0.4$. There投影片沒說的細節:$A(\omega) \lt 0$ 的區域。投影片寫「at all frequencies」,嚴格來說只對 $A(\omega) \gt 0$ 的頻率成立。在我們的例子中,$A(\omega) = -12.46 + 13.52\cos\omega$,在 $\omega \gt 0.4$ 時變成負數,此時スライドが触れていない細部:$A(\omega) \lt 0$ の領域。スライドは「at all frequencies」と書いているが、厳密には $A(\omega) \gt 0$ の周波数でしか成り立たない。この例では $A(\omega) = -12.46 + 13.52\cos\omega$ であり、$\omega \gt 0.4$ で負になる。そこでは
At $\omega = 1.5$, for instance, $\tau_p \approx -1.09$ — it looks like an advance.例如 $\omega = 1.5$ 時,$\tau_p \approx -1.09$,看起來像是「提前」了。たとえば $\omega = 1.5$ では $\tau_p \approx -1.09$ となり、「進んだ」ように見える。
Of course nothing has really come early. The output is “delayed by 1 sample and then turned upside down”:這當然不是真的提前。輸出其實是「延遲 1 個樣本再上下翻轉」:もちろん本当に進んだわけではない。出力は実は「1 サンプル遅らせてから上下を反転したもの」である。
Flipping is multiplying by $-1 = e^{j\pi}$, and for a single sinusoid a flip cannot be told apart from a shift of half a period. The phase-delay formula counts the flip as a time shift too, and so produces a strange number.翻轉等於乘上 $-1 = e^{j\pi}$,而對單一正弦波來說,翻轉和平移半個週期無法區分。相位延遲的公式把翻轉也算成了時間平移,所以得到一個奇怪的數字。反転は $-1 = e^{j\pi}$ を掛けることであり、単一の正弦波にとって反転と半周期のずれは区別できない。位相遅延の式は反転まで時間のずれとして数えてしまうので、奇妙な数が出てくる。
The group delay, by contrast, looks only at the slope of the phase; the constant $\pi$ disappears on differentiation, so the group delay is still 1 for $\omega \gt 0.4$. This is the first example of the group delay being more reliable than the phase delay.相比之下,群延遲只看相位的斜率,常數 $\pi$ 微分後消失,所以群延遲在 $\omega \gt 0.4$ 時仍然是 1。這是群延遲比相位延遲更可靠的第一個例子。これに対し群遅延は位相の傾きしか見ないので、定数 $\pi$ は微分で消え、$\omega \gt 0.4$ でも群遅延は 1 のままである。これが、群遅延のほうが位相遅延より信頼できることの最初の例である。
The interactive figure below shows this. On the left is the phase $\theta(\omega)$; the blue segment runs from the origin to the current point, and its slope equals $-\tau_p$. On the right are the input (grey) and the output (blue, divided by the gain so the positions are easy to compare).下面的互動圖可以觀察這一點。左邊是相位 $\theta(\omega)$,藍色線段從原點連到目前的點,它的斜率等於 $-\tau_p$;右邊是輸入(灰色)和輸出(藍色,已除以增益,方便比較位置)。下のインタラクティブな図でこれを確かめられる。左は位相 $\theta(\omega)$ で、青の線分は原点から現在の点へ引いたものであり、その傾きが $-\tau_p$ に等しい。右は入力(灰色)と出力(青。位置を比べやすいよう利得で割ってある)である。
Things to try:建議這樣操作:次のように操作してみるとよい。
The phase-delay formula cannot tell these two situations apart, but the group delay in the readout stays at 1 throughout.相位延遲的公式無法分辨這兩件事,但讀數中的群延遲始終是 1。位相遅延の式はこの二つを区別できないが、表示の群遅延は終始 1 のままである。
Phase delay describes a single sinusoid. Real signals carry information in envelopes. Consider a modulated carrier(調變載波):相位延遲描述的是單一弦波。真實訊號則把資訊帶在包絡裡。考慮一個調變載波:位相遅延が記述するのは単一の正弦波である。現実の信号は情報を包絡線に載せて運ぶ。変調された搬送波を考えよう。
The signal occupies two close frequencies $\omega_c \pm \omega_m$. Assume $|H(e^{j\omega})| \approx 1$ around $\omega_c$, and call the phases at the two side frequencies $\theta(\omega_c - \omega_m) = \theta_l$ and $\theta(\omega_c + \omega_m) = \theta_u$. The output recombines into這個訊號佔據兩個相近的頻率 $\omega_c \pm \omega_m$。假設 $\omega_c$ 附近有 $|H(e^{j\omega})| \approx 1$,並把兩個邊頻處的相位記為 $\theta(\omega_c - \omega_m) = \theta_l$ 與 $\theta(\omega_c + \omega_m) = \theta_u$。輸出重新合成為この信号は近接した二つの周波数 $\omega_c \pm \omega_m$ を占める。$\omega_c$ の周辺で $|H(e^{j\omega})| \approx 1$ と仮定し、二つの側波周波数での位相を $\theta(\omega_c - \omega_m) = \theta_l$、$\theta(\omega_c + \omega_m) = \theta_u$ と書く。出力は次のようにまとまる。
The average phase shifts the carrier; the difference shifts the envelope. If $\theta(\omega)$ is locally linear around $\omega_c$, i.e. $\theta(\omega_c + \Delta\omega) \approx \theta(\omega_c) + S\,\Delta\omega$ with slope $S = \frac{d\theta}{d\omega}\big|_{\omega_c}$, then:相位的平均移動載波,相位的差移動包絡。若 $\theta(\omega)$ 在 $\omega_c$ 附近是局部線性的,也就是 $\theta(\omega_c + \Delta\omega) \approx \theta(\omega_c) + S\,\Delta\omega$,斜率為 $S = \frac{d\theta}{d\omega}\big|_{\omega_c}$,那麼:位相の平均が搬送波を、位相の差が包絡線をずらす。$\theta(\omega)$ が $\omega_c$ の近傍で局所的に線形、すなわち傾きを $S = \frac{d\theta}{d\omega}\big|_{\omega_c}$ として $\theta(\omega_c + \Delta\omega) \approx \theta(\omega_c) + S\,\Delta\omega$ であれば、
| Phase delay $\tau_p(\omega)$相位延遲 $\tau_p(\omega)$位相遅延 $\tau_p(\omega)$ | Group delay $\tau_g(\omega)$群延遲 $\tau_g(\omega)$群遅延 $\tau_g(\omega)$ | |
|---|---|---|
| Definition定義定義 | $-\theta(\omega)/\omega$ | $-d\theta(\omega)/d\omega$ |
| Geometric meaning幾何意義幾何学的な意味 | (negative) slope of the chord from origin to $(\omega, \theta(\omega))$從原點到 $(\omega, \theta(\omega))$ 這條弦的(負)斜率原点から $(\omega, \theta(\omega))$ へ引いた弦の(負の)傾き | (negative) local tangent slope of $\theta(\omega)$$\theta(\omega)$ 的(負)局部切線斜率$\theta(\omega)$ の局所的な接線の(負の)傾き |
| What it delays延遲的對象何を遅らせるか | the carrier / a single sinusoid(載波)載波/單一弦波(載波)搬送波/単一の正弦波 | the envelope / information(包絡、訊息)包絡/資訊(包絡、訊息)包絡線/情報そのもの |
| Equal when兩者相等的條件等しくなる条件 | $\theta(\omega) = -\tau\omega$ exactly linear through the origin → $\tau_p = \tau_g = \tau$$\theta(\omega) = -\tau\omega$ 恰好是通過原點的直線 → $\tau_p = \tau_g = \tau$$\theta(\omega) = -\tau\omega$ が原点を通るちょうど線形のとき → $\tau_p = \tau_g = \tau$ | |
In discrete time $\tau_g$ is measured in samples: it is how long a small group of frequencies near $\omega$ is held up by the system. The AM pair above is one case of a general narrowband signal, a slowly varying envelope on a fast carrier:離散時間系統中 $\tau_g$ 的單位是樣本數。直覺上,它表示頻率在 $\omega$ 附近的一小群訊號,通過系統後被延遲了多久。上面的 AM 兩頻訊號只是一般窄頻訊號的特例:一個高頻載波乘上一個變化緩慢的包絡:離散時間では $\tau_g$ の単位はサンプル数である。直観的には、$\omega$ 付近の周波数の小さな群がシステムを通ってどれだけ遅れるかを表す。上の AM の 2 周波信号は、速い搬送波にゆっくり変わる包絡線を掛けた、一般の狭帯域信号の一例である。
$$x[n] = s[n]\cos(\omega_0 n)$$Its spectrum occupies only a narrow band around $\omega_0$ — a “group” of frequencies. Across that band the phase can be expanded to first order,它的頻譜只集中在 $\omega_0$ 附近一小段,也就是一「群」頻率。在這一小段範圍內,可以把相位做一階泰勒展開:そのスペクトルは $\omega_0$ のまわりの狭い帯域、すなわち周波数の「群」だけを占める。その帯域内では位相を 1 次まで展開できる。
$$\angle H(e^{j\omega}) \approx \angle H(e^{j\omega_0}) - \tau_g(\omega_0)\,(\omega-\omega_0)$$and substituting gives代進去之後,輸出大約是:これを代入すると、出力はおよそ次のようになる。
The group delay usually matters more, because what a signal “looks like” is set mainly by its envelope. The magnitude response says how much each frequency is amplified or attenuated; the group delay says how long each frequency is delayed.一般更關心群延遲,因為訊號「長什麼樣子」主要由包絡決定。大小響應告訴你「各頻率被放大或衰減多少」,群延遲告訴你「各頻率被延遲多久」。ふつう重要なのは群遅延のほうである。信号の「見た目」は主に包絡線で決まるからである。振幅応答は各周波数がどれだけ増幅・減衰されるかを、群遅延は各周波数がどれだけ遅れるかを教える。
Why check whether $\tau_g$ is constant.為什麼要看群延遲是不是常數。なぜ $\tau_g$ が一定かどうかを見るのか。
Link to the lecture — group delays add. Cascaded phases add, so cascaded group delays add too:和講義的關係 — 群延遲可以逐個因式相加。串接系統的相位相加,因此群延遲也直接相加:講義とのつながり — 群遅延は足し合わせられる。縦続接続した系の位相は足し合わされるので、群遅延もそのまま足し合わされる。
$$H = H_1 H_2 \cdots \;\Rightarrow\; \tau_g = \tau_{g,1} + \tau_{g,2} + \cdots$$That is why the lecture analyses only the single zero $1 - re^{j\theta}z^{-1}$ (Fig. 10-4): the group delay of any rational system is the sum of the contributions of its poles and zeros. For one zero,這就是講義只分析單一零點 $1 - re^{j\theta}z^{-1}$(圖 10-4)的原因:任何有理系統的群延遲,都是各個極點、零點的貢獻加總。單一零點的群延遲為:講義がひとつの零点 $1 - re^{j\theta}z^{-1}$(図 10-4)だけを解析するのはこのためである。どんな有理システムの群遅延も、各極・零点の寄与の和になる。ひとつの零点では次のとおりである。
With $r = 0.9$ this is $-9$ samples — the downward spike in Examples I and II. Geometrically, as $e^{j\omega}$ passes beside the zero, the zero vector's angle turns very fast; the steeper the phase change, the larger $|\tau_g|$.$r = 0.9$ 時是 $-9$,這就是 Example I、II 圖上往下的尖刺。從幾何上看,$e^{j\omega}$ 經過零點旁邊時,零點向量的角度轉得非常快;相位變化越急,群延遲的絕對值越大。$r = 0.9$ では $-9$ サンプルとなり、これが例 I・II の下向きの尖りである。幾何的には、$e^{j\omega}$ が零点のそばを通るとき零点ベクトルの角度が非常に速く回る。位相の変化が急なほど $|\tau_g|$ は大きい。
Why a modulated carrier. The conclusion of §10.3.1 was that a single sinusoid has no time marker, so its phase delay is ambiguous. To really measure “how long the signal is delayed” we need a signal with a feature we can follow.為什麼要用調變載波。§10.3.1 的結論是:單一正弦波沒有時間標記,所以相位延遲有歧義。要真正量測「訊號延遲多久」,需要一個有特徵、可以追蹤的訊號。なぜ変調された搬送波を使うのか。§10.3.1 の結論は、単一の正弦波には時刻の目印がないので位相遅延があいまいになる、というものだった。「信号がどれだけ遅れるか」を本当に測るには、追跡できる特徴をもつ信号が必要である。
The simplest choice is amplitude modulation (AM):最簡單的選擇是振幅調變(AM):最も簡単な選択は振幅変調(AM)である。
In the figure of slide 34 the fast blue oscillation is the carrier and the slowly varying red dashed curve is the envelope. The envelope's period $2\pi/\omega_m$ is much longer than the carrier's, so its crest is a clear time marker: see when the output envelope's crest turns up and you know how long the signal was delayed.第 34 頁的圖中,藍色是快速振盪的載波(carrier),紅色虛線是緩慢變化的包絡(envelope)。包絡的週期是 $2\pi/\omega_m$,比載波長很多,所以它的波峰就是一個清楚的時間標記:看輸出包絡的波峰出現在什麼時候,就知道訊號延遲了多久。スライド 34 の図では、速く振動する青が搬送波(carrier)、ゆっくり変化する赤の破線が包絡線(envelope)である。包絡線の周期 $2\pi/\omega_m$ は搬送波よりずっと長いので、その山ははっきりした時刻の目印になる。出力の包絡線の山がいつ現れるかを見れば、信号がどれだけ遅れたかがわかる。
This is not just a teaching example. Real communication and speech signals almost all have this form: a high-frequency carrier, or a band, carrying slowly varying information. The information is in the envelope, not in each crest of the carrier.這不只是教學用的例子。實際的通訊和語音訊號幾乎都是這種形式:一個高頻載波,或一段頻帶,上面承載著緩慢變化的資訊。資訊在包絡裡,不在載波的每個波峰裡。これは単なる教材用の例ではない。実際の通信信号や音声信号はほとんどがこの形をしている。高周波の搬送波、あるいは一つの帯域の上に、ゆっくり変化する情報が乗っている。情報は包絡線の中にあり、搬送波の一つひとつの山にあるのではない。
The spectrum: two lines. Slide 35 splits the signal with the product-to-sum identity:頻譜:兩條線。第 35 頁用積化和差把訊號拆開:スペクトル:2 本の線。スライド 35 は積和公式で信号を分解する。
In the frequency domain these are two spectral lines on either side of $\omega_c$, $\omega_m$ away from it (the small plot at the top right of the slide): the lower and the upper sideband.在頻域中,這是 $\omega_c$ 兩側、相距 $\omega_m$ 的兩條譜線(投影片右上角的小圖),分別稱為下邊帶和上邊帶。周波数領域では、$\omega_c$ の両側に $\omega_m$ 離れて立つ 2 本のスペクトル線(スライド右上の小さな図)であり、それぞれ下側波帯、上側波帯と呼ばれる。
This step is the key: the modulated signal is just a sum of two sinusoids. And we already know what a system does to a single sinusoid: multiply by the gain, add the phase shift.這一步很關鍵:調變訊號只是兩個正弦波的和。而我們已經知道系統對單一正弦波做什麼:乘上增益、加上相移。この一歩が鍵である。変調信号は二つの正弦波の和にすぎない。そしてシステムが単一の正弦波に何をするかはもうわかっている。利得を掛け、位相のずれを加えるだけである。
Through the system. Each sinusoid passes through on its own:通過系統。兩個正弦波各自通過系統:システムを通す。二つの正弦波はそれぞれ別々にシステムを通る。
with $\theta_l = \theta(\omega_c - \omega_m)$ and $\theta_u = \theta(\omega_c + \omega_m)$. The slide's $H(\cdot)\cos(\cdot)$ is shorthand for “multiply by the magnitude, add the phase”.其中 $\theta_l = \theta(\omega_c - \omega_m)$,$\theta_u = \theta(\omega_c + \omega_m)$。投影片寫成 $H(\cdot)\cos(\cdot)$ 是簡寫,意思就是「乘上大小、加上相位」。ここで $\theta_l = \theta(\omega_c - \omega_m)$、$\theta_u = \theta(\omega_c + \omega_m)$ である。スライドの $H(\cdot)\cos(\cdot)$ は略記であり、「大きさを掛け、位相を加える」という意味である。
The slide then assumes $|H| \approx 1$. The assumption is there to separate the effect of the magnitude from that of the phase: if the two lines got different gains the depth of the envelope would change, and that is amplitude distortion, not the subject here. Under the assumption only the phase is left:投影片接著假設 $|H| \approx 1$。這個假設是為了把大小的影響和相位的影響分開:如果兩條譜線的增益不同,包絡的深度會改變,那是振幅失真,不是這裡要討論的。在這個假設下,只剩相位:スライドは続いて $|H| \approx 1$ と仮定する。この仮定は大きさの影響と位相の影響を切り分けるためのものである。2 本の線の利得が違えば包絡線の深さが変わるが、それは振幅歪みであり、ここでの話題ではない。この仮定のもとでは位相だけが残る。
Recombining: sum-to-product. Slide 36 turns the two terms back into a product with重新組合:和差化積。第 36 頁用和差化積公式把兩項合回乘積形式:組み立て直す:和積公式。スライド 36 は和積公式で二つの項を積の形に戻す。
Let $a = (\omega_c + \omega_m)n + \theta_u$ and $b = (\omega_c - \omega_m)n + \theta_l$:令 $a = (\omega_c + \omega_m)n + \theta_u$,$b = (\omega_c - \omega_m)n + \theta_l$:$a = (\omega_c + \omega_m)n + \theta_u$、$b = (\omega_c - \omega_m)n + \theta_l$ とおく。
The leading $\frac{A}{2}$ times the 2 of the identity is exactly $A$:前面的 $\frac{A}{2}$ 乘上公式中的 2,剛好是 $A$:前の $\frac{A}{2}$ に公式の 2 を掛けると、ちょうど $A$ になる。
The output is still “carrier times envelope”, the same shape as the input, but the two parts get different phase shifts:輸出仍然是「載波乘上包絡」,形狀和輸入一樣,但兩者各自得到不同的相移:出力はやはり「搬送波 × 包絡線」で入力と同じ形をしているが、両者はそれぞれ異なる位相のずれを受ける。
What the result means. This is the most important observation of the whole topic: the envelope cares only about the difference in phase between the two frequencies, not about the size of the phase itself.這個結果的意義。這是整個主題中最重要的觀察:包絡只在乎兩個頻率之間的相位「差」,不在乎相位本身的大小。この結果の意味。これがこの話題全体で最も重要な観察である。包絡線が気にするのは二つの周波数の間の位相の「差」だけであり、位相そのものの大きさは関係ない。
If the system adds the same phase to both frequencies — say $\pi/2$ to each — then $\theta_u - \theta_l$ is unchanged, the envelope is untouched and only the carrier moves. Conversely, as soon as the phases at the two frequencies differ, the envelope is shifted.如果系統讓兩個頻率都多了同一個相位,例如都加上 $\pi/2$,那麼 $\theta_u - \theta_l$ 不變,包絡完全不受影響,只有載波被移動。反過來說,只要兩個頻率的相位有差異,包絡就會被平移。システムが二つの周波数に同じ位相、たとえばどちらにも $\pi/2$ を加えるなら、$\theta_u - \theta_l$ は変わらず、包絡線はまったく影響を受けず、動くのは搬送波だけである。逆に、二つの周波数の位相に差があれば、包絡線はずらされる。
This echoes exactly what was said before: the group delay is the rate at which the steady-state phase shift differs from one frequency to the next. Slide 36 is the full derivation of that sentence.這正好呼應我們先前討論過的:群延遲是穩態相移在不同頻率之間的差異率。第 36 頁就是這句話的完整推導。これは前に述べたことにちょうど対応する。群遅延とは、定常状態の位相のずれが周波数ごとに異なる割合である。スライド 36 はこの一文の完全な導出である。
Taylor expansion. If $\omega_m$ is small, both lines sit close to $\omega_c$, and over such a short range $\theta(\omega)$ can be approximated by a straight line:泰勒展開。如果 $\omega_m$ 很小,兩條譜線都很靠近 $\omega_c$,在這麼小的範圍內,$\theta(\omega)$ 可以近似成直線:テイラー展開。$\omega_m$ が小さければ 2 本の線はどちらも $\omega_c$ のすぐ近くにあり、その狭い範囲では $\theta(\omega)$ を直線で近似できる。
Substituting $\Delta\omega = \pm\omega_m$:代入 $\Delta\omega = \pm\omega_m$:$\Delta\omega = \pm\omega_m$ を代入する。
The carrier gives the phase delay.載波:得到相位延遲。搬送波からは位相遅延が得られる。
The $S\omega_m$ terms cancel, the carrier's phase shift is $\theta(\omega_c)$ itself, and the corresponding delay is precisely the phase delay of slide 32:$S\omega_m$ 正負抵消,載波的相移就是 $\theta(\omega_c)$ 本身,對應的延遲正是第 32 頁的相位延遲:$S\omega_m$ は正負で打ち消し合い、搬送波の位相のずれは $\theta(\omega_c)$ そのものになる。対応する遅延はまさにスライド 32 の位相遅延である。
The envelope gives the group delay.包絡:得到群延遲。包絡線からは群遅延が得られる。
The envelope becomes包絡變成:包絡線は次のようになる。
so the envelope is delayed by所以包絡被延遲:したがって包絡線の遅れは
— which is the definition of the group delay.這就是群延遲的定義。であり、これが群遅延の定義である。
An easily missed point: $\omega_m$ has disappeared. Note that $\tau_g = -S$ contains no $\omega_m$. As long as local linearity holds, the delay is the same no matter how fast or slowly the envelope varies.一個容易忽略的重點:$\omega_m$ 消失了。注意 $\tau_g = -S$ 裡沒有 $\omega_m$。只要局部線性成立,不論包絡變化多快或多慢,延遲都一樣。見落としやすい点:$\omega_m$ が消えた。$\tau_g = -S$ には $\omega_m$ が含まれていないことに注意しよう。局所的な線形性が成り立つかぎり、包絡線がどれほど速く、あるいは遅く変化しても、遅延は同じである。
This matters. A real envelope is not a single cosine; it is made of many different $\omega_m$ (its Fourier decomposition). Since every component is delayed by the same $\tau_g$, the whole envelope is shifted intact and keeps its shape. So the group delay is a property of the system near $\omega_c$, not a property of a particular signal.這很重要。真實的包絡不只是一個餘弦,而是由許多不同的 $\omega_m$ 組成(傅立葉分解)。既然每個成分都被延遲同樣的 $\tau_g$,整個包絡就被完整地平移,形狀不變。所以群延遲是系統在 $\omega_c$ 附近的性質,不是特定訊號的性質。これは重要である。現実の包絡線は一つの余弦ではなく、多くの異なる $\omega_m$ からなる(フーリエ分解)。どの成分も同じ $\tau_g$ だけ遅れるのだから、包絡線全体はそのまま平行移動し、形は変わらない。したがって群遅延は$\omega_c$ 付近におけるシステムの性質であり、特定の信号の性質ではない。
Conversely, if $\theta(\omega)$ is not a straight line across the signal's bandwidth, different $\omega_m$ get different delays and the shape of the envelope is destroyed. This sets up the last sentence of slide 38.反過來說,如果 $\theta(\omega)$ 在訊號頻寬內不是直線,不同的 $\omega_m$ 就會得到不同的延遲,包絡形狀就會被破壞。這是第 38 頁最後一句話的伏筆。逆に、$\theta(\omega)$ が信号の帯域幅の中で直線でなければ、異なる $\omega_m$ は異なる遅延を受け、包絡線の形は壊れる。これがスライド 38 の最後の一文への伏線である。
Interactive: a modulated carrier through an allpass filter. To match the slide's $|H| \approx 1$ assumption exactly, the figure below uses a second-order allpass filter: poles at $\rho e^{\pm j1.2}$, zeros at the conjugate-reciprocal positions. Its $|H|$ is exactly 1 at every frequency; only the phase changes, so everything you see comes from the phase alone.互動圖:調變載波通過全通濾波器。為了完全符合投影片 $|H| \approx 1$ 的假設,下面用二階全通濾波器(allpass):極點在 $\rho e^{\pm j1.2}$,零點在共軛倒數位置。它的 $|H|$ 在所有頻率都精確等於 1,只改變相位,所以看到的效果純粹來自相位。インタラクティブ:全域通過フィルタを通る変調された搬送波。スライドの $|H| \approx 1$ の仮定に正確に合わせるため、下の図では2 次の全域通過フィルタ(allpass)を使う。極は $\rho e^{\pm j1.2}$、零点は共役逆数の位置にある。その $|H|$ はすべての周波数で厳密に 1 であり、変わるのは位相だけなので、見える効果は純粋に位相から来ている。
Near the resonance at 1.2 the phase changes fast, so the group delay is large; the phase itself is not large there, so the phase delay is small. This is a case where the two differ clearly.在共振頻率 1.2 附近,相位變化很快,群延遲很大;而相位本身的大小不大,所以相位延遲較小。這正是兩者明顯不同的情況。共振周波数 1.2 の付近では位相が急に変わるので群遅延は大きいが、位相そのものはそれほど大きくないので位相遅延は小さい。両者がはっきり異なる場面である。
Try these in order:建議依序嘗試:次の順に試してみるとよい。
Chord and tangent. The lower figure of slide 38 is the most intuitive way to understand the two delays:割線與切線。第 38 頁下方的圖是理解兩種延遲最直觀的方式:弦と接線。スライド 38 の下の図は、二つの遅延を理解するいちばん直観的な方法である。
In the slide's figure the red tangent is steeper than the blue chord, so $\tau_g \gt \tau_p$, matching the upper time-domain figure where the red arrow (envelope delay) is longer than the blue one (carrier delay). It is exactly what the interactive figure above shows at $\omega_c = 1.2$.投影片的圖中,紅色切線比藍色割線陡,所以 $\tau_g \gt \tau_p$,對應上方時域圖中紅色箭頭(包絡延遲)比藍色箭頭(載波延遲)長。這正是前一個互動圖在 $\omega_c = 1.2$ 時看到的情況。スライドの図では赤の接線が青の弦より急なので $\tau_g \gt \tau_p$ であり、上の時間領域の図で赤の矢印(包絡線の遅延)が青の矢印(搬送波の遅延)より長いことに対応する。上のインタラクティブな図で $\omega_c = 1.2$ のときに見えたのがまさにこれである。
How the two are related. From $\theta(\omega) = -\omega\,\tau_p(\omega)$, differentiate with respect to $\omega$:兩者的數學關係。由 $\theta(\omega) = -\omega\,\tau_p(\omega)$,對 $\omega$ 微分:両者の関係。$\theta(\omega) = -\omega\,\tau_p(\omega)$ を $\omega$ で微分する。
So $\tau_g = \tau_p$ if and only if $\tau_p$ does not change with frequency, i.e. the phase is a straight line through the origin.所以 $\tau_g = \tau_p$ 若且唯若 $\tau_p$ 不隨頻率改變,也就是相位是一條通過原點的直線。したがって $\tau_g = \tau_p$ となるのは $\tau_p$ が周波数によって変わらないとき、かつそのときに限る。つまり位相が原点を通る直線のときである。
Three cases.三種情況。三つの場合。
| Form of the phase相位形式位相の形 | Examples例子例 | τp | τg | Output輸出出力 |
|---|---|---|---|---|
| $\theta = -\omega n_0$ (a line through the origin)$\theta = -\omega n_0$(過原點的直線)$\theta = -\omega n_0$(原点を通る直線) | pure delay; symmetric FIR where $A \gt 0$純延遲、對稱 FIR 在 $A \gt 0$ 區域純粋な遅延、$A \gt 0$ の領域の対称 FIR | $n_0$ | $n_0$ | the whole waveform is shifted, completely undistorted整個波形平移,完全不變形波形全体が平行移動し、まったく歪まない |
| $\theta = -\omega n_0 + c$ (a line not through the origin)$\theta = -\omega n_0 + c$(不過原點的直線)$\theta = -\omega n_0 + c$(原点を通らない直線) | the 3-point filter for $\omega \gt 0.4$ ($c = \pi$); the Hilbert transformer ($c = -\pi/2$)3 點濾波器在 $\omega \gt 0.4$($c = \pi$)、Hilbert 轉換器($c = -\pi/2$)$\omega \gt 0.4$ の 3 点フィルタ($c = \pi$)、ヒルベルト変換器($c = -\pi/2$) | $n_0 - c/\omega$ | $n_0$ | the envelope is shifted intact; the carrier gets an extra fixed phase relative to the envelope包絡完整平移;載波相對包絡多了固定相移包絡線はそのまま平行移動し、搬送波は包絡線に対して一定の位相のずれを余分に受ける |
| $\theta$ not a straight lineθ 不是直線θ が直線でない | IIR filters, allpass filtersIIR 濾波器、全通濾波器IIR フィルタ、全域通過フィルタ | varies with ω隨 ω 變化ω とともに変わる | varies with ω隨 ω 變化ω とともに変わる | the envelope may be deformed (phase distortion)包絡可能變形(相位失真)包絡線が変形しうる(位相歪み) |
The second case deserves attention. The envelope is not distorted at all, because the group delay is constant, but the relative position of carrier and envelope changes. For signals whose information is in the envelope, such as AM or speech, that has no effect whatsoever, so the second case is treated as linear phase just like the first — strictly, generalized linear phase. This also explains the earlier question of why the 3-point filter still counts as linear phase although its phase jumps by $\pi$ at $\omega = 0.4$.第二種情況特別值得注意。包絡完全不失真,因為群延遲是常數,但載波和包絡的相對位置改變了。對 AM 或語音這類「資訊在包絡」的訊號,這完全沒有影響,所以第二種情況和第一種一樣,都被視為線性相位,嚴格稱為廣義線性相位。這也解釋了先前的問題:3 點濾波器在 $\omega \gt 0.4$ 處相位跳 $\pi$,為什麼仍然算線性相位。二番目の場合は特に注目に値する。群遅延が一定なので包絡線はまったく歪まないが、搬送波と包絡線の相対的な位置は変わる。AM や音声のように「情報が包絡線にある」信号にとってはまったく影響がないので、二番目の場合も一番目と同じく線形位相とみなされ、厳密には一般化線形位相と呼ばれる。3 点フィルタの位相が $\omega = 0.4$ で $\pi$ 跳ぶのに、なぜ線形位相とみなされるのかという先の問いも、これで説明がつく。
If $\theta(\omega)$ is not linear across the band occupied by $A[n]$, different envelope frequencies suffer different delays and the envelope is smeared — phase distortion(相位失真). The lecture closes with a striking example, worked through in full at the end of this subsection: a signal made of three pulse bursts at different carrier frequencies passes through a system whose group delay differs by some 140 samples between two of those bands — at the output the bursts come out in the opposite order, and nothing in the magnitude response hints that they would.若 $\theta(\omega)$ 在 $A[n]$ 所佔的頻帶上不是線性的,不同的包絡頻率就會承受不同的延遲,包絡因而被塗抹開來——這就是相位失真。投影片以一個很鮮明的例子收尾,本小節最後會完整推演一遍:一個由三段不同載波頻率的脈衝所組成的訊號,通過一個在其中兩個頻帶之間群延遲相差約 140 個樣本的系統——輸出端的脈衝以相反的順序出場,而大小響應完全看不出這件事會發生。$\theta(\omega)$ が $A[n]$ の占める帯域にわたって線形でない場合、包絡線の周波数ごとに遅れが異なり、包絡線はにじんでしまう。これが位相歪みである。講義は印象的な例で締めくくられる。本小節の最後で最後まで追う例である。搬送波周波数の異なる三つのパルスからなる信号を、そのうち二つの帯域のあいだで群遅延がおよそ 140 サンプルも違うシステムに通す。出力ではパルスが逆の順序で出てきて、そうなることは振幅応答からはまったく読み取れない。
From two lines to a real signal. The last sentence of slide 38: “If θ(ω) is not linear around ωc, A[n] suffers phase distortion → correction…”從兩條譜線到真實訊號。第 38 頁最後一句:「If θ(ω) is not linear around ωc, A[n] suffers phase distortion → correction…」2 本の線から現実の信号へ。スライド 38 の最後の一文:「If θ(ω) is not linear around ωc, A[n] suffers phase distortion → correction…」
In the two-line analysis of §10.3.2 the envelope had a single frequency $\omega_m$; whatever the system, the output envelope was still a cosine, at most with a different amount of delay. A real envelope is far more complicated — a pulse, a syllable of speech, a digital symbol. Its spectrum occupies a bandwidth $B$ around $\omega_c$ and can be thought of as many pairs of lines with different $\omega_m$ stacked together.§10.3.2 的兩譜線分析中,包絡只有一個頻率 $\omega_m$,不論系統如何,輸出包絡都還是一個餘弦,最多只是延遲量不同。真實的包絡則複雜得多,例如一個脈衝、一個語音音節、一個數位符號。它的頻譜在 $\omega_c$ 附近佔據一段頻寬 $B$,可以想成很多對不同 $\omega_m$ 的譜線疊加。§10.3.2 の 2 本線の解析では包絡線の周波数は $\omega_m$ 一つだけで、システムが何であれ出力の包絡線は余弦のままであり、違うのはせいぜい遅延量だけだった。現実の包絡線はずっと複雑で、たとえば一つのパルス、音声の一音節、一つのディジタル符号である。そのスペクトルは $\omega_c$ のまわりに帯域幅 $B$ を占め、異なる $\omega_m$ をもつ多くの線の対が重なったものと考えられる。
Most noteworthy: even if $|H| = 1$ at every frequency, the signal can still be badly deformed. A flat magnitude response does not mean the waveform is preserved.最值得注意的是:即使 $|H|$ 在所有頻率都等於 1,訊號仍然可能嚴重變形。大小響應平坦,不代表波形會被保留。最も注目すべきは、すべての周波数で $|H| = 1$ であっても、信号はなお大きく変形しうるということである。振幅応答が平坦だからといって、波形が保たれるわけではない。
Judging how bad the distortion is. A practical criterion compares two quantities:失真程度的判斷。一個實用的判斷標準是比較兩個量:歪みの程度の見積もり。実用的な目安は二つの量を比べることである。
If $\Delta\tau_g$ is much smaller than the envelope's width, the misalignment between components is negligible and the envelope hardly changes shape. If $\Delta\tau_g$ is comparable to the width or larger, the distortion is obvious.如果 $\Delta\tau_g$ 遠小於包絡寬度,不同成分之間的錯開微不足道,包絡幾乎不變形。如果 $\Delta\tau_g$ 和包絡寬度差不多或更大,失真就很明顯。$\Delta\tau_g$ が包絡線の幅よりずっと小さければ、成分どうしのずれは無視でき、包絡線はほとんど変形しない。$\Delta\tau_g$ が包絡線の幅と同程度かそれ以上なら、歪みははっきり現れる。
This brings out an important trade-off: the narrower (shorter) the envelope, the wider its bandwidth, the more $\tau_g$ variation it covers, and the more easily it is distorted.這帶出一個重要的取捨:包絡越窄(越短促),頻寬越寬,涵蓋的 $\tau_g$ 變化越多,越容易失真。ここから重要なトレードオフが出てくる。包絡線が狭い(短い)ほど帯域幅は広くなり、覆う $\tau_g$ の変化が大きくなり、歪みやすくなる。
Interactive: a Gaussian pulse through the allpass filter. The input is a Gaussian envelope times a carrier,互動圖:高斯脈衝通過全通濾波器。輸入是高斯包絡乘上載波:インタラクティブ:全域通過フィルタを通るガウスパルス。入力はガウス形の包絡線に搬送波を掛けたものである。
$\sigma$ is the envelope width, and the bandwidth is roughly $1/\sigma$. The system is still the allpass of §10.3.2, but this time the output is computed sample by sample from the difference equation, with no approximation at all:$\sigma$ 是包絡寬度,頻寬約為 $1/\sigma$。系統仍是 §10.3.2 的全通濾波器,這次直接用差分方程逐點計算輸出,不使用任何近似:$\sigma$ は包絡線の幅で、帯域幅はおよそ $1/\sigma$ である。システムは §10.3.2 と同じ全域通過フィルタだが、今回は差分方程式から出力を一点ずつ計算し、近似は一切使わない。
The red dashed line is the “ideal result”: the input envelope shifted by $\tau_g(\omega_c)$. If the output hugs the red dashed line, the envelope has survived intact; if it does not, that is phase distortion.紅色虛線是「理想結果」:把輸入包絡平移 $\tau_g(\omega_c)$。如果輸出貼合紅色虛線,代表包絡完整保留;如果不貼合,就是相位失真。赤の破線は「理想的な結果」、すなわち入力の包絡線を $\tau_g(\omega_c)$ だけずらしたものである。出力が赤の破線にぴったり沿えば包絡線はそのまま保たれており、沿わなければそれが位相歪みである。
Try the three buttons in order:建議依序嘗試三個按鈕:三つのボタンを順に試してみるとよい。
Once more: this filter's $|H|$ is exactly 1 at every frequency. No frequency component of the signal is amplified or attenuated, yet the waveform is wrecked — only because the components are delayed by different amounts. You can also fix $\sigma = 12$ and slowly pull $\rho$ from 0.5 to 0.97: the larger $\rho$, the sharper the resonance and the higher and narrower the $\tau_g$ peak, so the same pulse is distorted more easily.再強調一次:這個濾波器的 $|H|$ 在所有頻率都精確等於 1。訊號的每個頻率成分都沒有被放大或衰減,但波形還是被破壞了,只因為各成分的延遲不同。你也可以固定 $\sigma = 12$,然後慢慢把 $\rho$ 從 0.5 拉到 0.97。$\rho$ 越大,共振越尖銳,$\tau_g$ 峰值越高也越窄,同樣的脈衝就越容易被扭曲。もう一度強調しておく。このフィルタの $|H|$ はすべての周波数で厳密に 1 である。信号のどの周波数成分も増幅も減衰もされていないのに、波形は壊れる。ただ各成分の遅延が違うというだけの理由で。$\sigma = 12$ に固定して、$\rho$ を 0.5 から 0.97 へゆっくり上げてみてもよい。$\rho$ が大きいほど共振は鋭くなり、$\tau_g$ のピークは高く狭くなり、同じパルスがより歪みやすくなる。
Oppenheim & Schafer close §5.1 with this example, which takes everything above — “group delay = envelope delay”, applied band by band — and runs it on one real signal. The whole point fits in a sentence: the system delays the envelopes at different frequencies by different amounts, and as a result the three bursts change places.Oppenheim & Schafer 在 §5.1 的最後用這個例子收尾。它把前面「群延遲 = 包絡延遲」的觀念逐頻帶套用在一個實際的訊號上。重點只有一句話:系統對不同頻率的包絡造成不同的延遲,結果三個脈衝的出場順序被改變了。Oppenheim & Schafer は §5.1 をこの例で締めくくる。これまでの「群遅延 = 包絡線の遅延」という考えを、帯域ごとに一つの現実の信号へ当てはめたものである。要点は一文に収まる。システムは周波数ごとに包絡線を違う量だけ遅らせ、その結果、三つのパルスの出てくる順序が入れ替わる。
The input: three narrowband bursts. $x[n]$ is three bursts laid end to end, each about 61 samples long, each a Hamming envelope (a bell, large in the middle and small at the edges) filled with a carrier of its own. The rate of oscillation tells them apart:輸入訊號:三個窄頻脈衝。$x[n]$ 由三段脈衝串接而成,每段長約 61 個樣本,包絡是 Hamming 窗(中間大、兩邊小的鐘形),裡面是各自不同頻率的載波。從振盪的快慢可以分辨:入力信号:三つの狭帯域パルス。$x[n]$ は三つのパルスを端から端へつないだもので、各区間はおよそ 61 サンプル、包絡線は Hamming 窓(中央が大きく両端が小さい釣鐘形)、その中にそれぞれ別の搬送波が入っている。振動の速さで見分けられる。
| Burst脈衝パルス | Samples樣本範圍サンプル範囲 | Centre中心中心 | Carrier載波頻率搬送波周波数 | Period週期周期 |
|---|---|---|---|---|
| 1st第 1 個1 番目 | 0 – 60 | $n = 30$ | $0.8\pi$ | 2.5 samples (nearly alternating sign)2.5 個樣本(幾乎正負交錯)2.5 サンプル(ほぼ正負が交互) |
| 2nd第 2 個2 番目 | 61 – 121 | $n = 91$ | $0.2\pi$ | 10 samples10 個樣本10 サンプル |
| 3rd第 3 個3 番目 | 122 – 182 | $n = 152$ | $0.4\pi$ | 5 samples5 個樣本5 サンプル |
$|X(e^{j\omega})|$ confirms it: six peaks, at $\pm 0.2\pi$, $\pm 0.4\pi$ and $\pm 0.8\pi$, and essentially nothing anywhere else. Each peak is narrow — a Hamming window of length 61 has a main lobe of half-width $4\pi/60 \approx 0.067\pi$. The longer the envelope lasts, the narrower its band. All three bursts are narrowband signals, so the result derived above applies to each one directly: after passing through the system, a burst's amplitude is multiplied by $|H(e^{j\omega_c})|$ and its envelope is delayed by $\tau_g(\omega_c)$.$|X(e^{j\omega})|$ 驗證了這一點:只有在 $\pm 0.2\pi$、$\pm 0.4\pi$、$\pm 0.8\pi$ 有六個尖峰,其他頻率幾乎是 0。每個尖峰都很窄——長度 61 的 Hamming 窗,主瓣半寬約 $4\pi/60 \approx 0.067\pi$。包絡持續的時間越長,頻寬越窄。三個脈衝都是窄頻訊號,所以前面推導的結論可以逐一直接套用:每個脈衝通過系統後,振幅乘上 $|H(e^{j\omega_c})|$,包絡延遲 $\tau_g(\omega_c)$。$|X(e^{j\omega})|$ がそれを裏づける。$\pm 0.2\pi$、$\pm 0.4\pi$、$\pm 0.8\pi$ に六つの峰があるだけで、ほかの周波数はほぼ 0 である。どの峰も狭い。長さ 61 の Hamming 窓のメインローブは半値幅でおよそ $4\pi/60 \approx 0.067\pi$ である。包絡線が長く続くほど帯域は狭くなる。三つのパルスはいずれも狭帯域信号なので、上で導いた結果をそれぞれにそのまま当てはめられる。システムを通ると、振幅は $|H(e^{j\omega_c})|$ 倍になり、包絡線は $\tau_g(\omega_c)$ だけ遅れる。
The system. It is best read as two parts in cascade: one shapes the magnitude, the other does nothing but pile up delay near $0.2\pi$.系統。這個系統最好看成兩部分的串接:一部分負責塑造大小,另一部分只負責在 $0.2\pi$ 附近製造大量延遲。システム。二つの部分の縦続接続と読むのがいちばんわかりやすい。一方は振幅を形づくり、もう一方は $0.2\pi$ の近くに遅延を積み上げるだけの働きをする。
with $r_z = 0.98$ (a pair of zeros just inside the circle at $\pm 0.8\pi$), $r_p = 0.80$ (a pair of poles at $\pm 0.4\pi$), $G = 0.3283$, and the eight allpass pole pairs at radius $\rho = 0.95$, their angles $\phi_k$ spread evenly from $0.164\pi$ to $0.236\pi$.其中 $r_z = 0.98$(一對零點,在 $\pm 0.8\pi$ 方向、剛好落在單位圓內),$r_p = 0.80$(一對極點在 $\pm 0.4\pi$),$G = 0.3283$,而八組全通極點的半徑 $\rho = 0.95$,角度 $\phi_k$ 從 $0.164\pi$ 到 $0.236\pi$ 均勻分佈。ここで $r_z = 0.98$($\pm 0.8\pi$ の方向で単位円のすぐ内側にある一対の零点)、$r_p = 0.80$($\pm 0.4\pi$ にある一対の極)、$G = 0.3283$ であり、八組の全域通過の極は半径 $\rho = 0.95$、角度 $\phi_k$ は $0.164\pi$ から $0.236\pi$ まで等間隔に並ぶ。
| $\omega$ | $|H(e^{j\omega})|$ | $\tau_g(\omega)$ | Why原因理由 |
|---|---|---|---|
| $0.2\pi$ | 1.21 | 150 | the allpass cluster — a tall, narrow spike with no feature at all in the magnitude全通群集——又高又窄的尖峰,而大小響應上完全沒有對應的特徵全域通過の集まり — 高く狭い尖りで、振幅側には対応する特徴がまったくない |
| $0.4\pi$ | 2.10 | 7 | poles near this direction form the peak; delay is only a few samples極點靠近這個方向,形成峰值;延遲只有幾個樣本この方向に近い極が峰をつくる。遅延は数サンプルだけ |
| $0.8\pi$ | 0.004 | −49 | zeros at $r_z = 0.98$ almost on the circle: the band is removed, and the delay dives negative$r_z = 0.98$ 的零點幾乎貼在圓上:這個頻帶被濾掉,延遲同時深深轉負$r_z = 0.98$ の零点が円のすぐ内側にある。この帯域は取り除かれ、同時に遅延は大きく負へ沈む |
| $0$, 、、$\pi$ | 1.02, 、、0.06 | 5, 、、0.5 | far from everything遠離所有極零點極にも零点にも遠い |
Predicting the output. Handle the three bursts separately; each needs only two numbers, $|H(e^{j\omega_c})|$ and $\tau_g(\omega_c)$.預測輸出。把三個脈衝分開處理,每一個只需要兩個數字:$|H(e^{j\omega_c})|$ 和 $\tau_g(\omega_c)$。出力を予測する。三つのパルスを別々に扱えばよい。必要な数字はそれぞれ二つ、$|H(e^{j\omega_c})|$ と $\tau_g(\omega_c)$ だけである。
| Burst脈衝パルス | Input centre輸入中心入力の中心 | $|H|$ | $\tau_g$ | Predicted centre預測輸出中心予測される出力の中心 | Measured centre實測輸出中心実測の出力の中心 | Measured peak實測振幅実測の振幅 |
|---|---|---|---|---|---|---|
| $0.8\pi$ | 30 | ≈ 0 | −49 (meaningless)−49(無意義)−49(意味をもたない) | removed被濾除取り除かれる | — | 0.02 |
| $0.2\pi$ | 91 | 1.21 | 150 | $91 + 150 = \mathbf{241}$ | 232 | 1.00 |
| $0.4\pi$ | 152 | 2.10 | 7 | $152 + 7 = \mathbf{159}$ | 159 | 1.98 |
The $0.4\pi$ burst lands exactly where predicted and comes out $2.1\times$ taller. The $0.2\pi$ burst arrives 9 samples early and 0.2 short of its predicted height — §“Three details” below explains why, and it is not an error in the prediction so much as a limit of it.$0.4\pi$ 的脈衝精準落在預測位置,而且振幅放大了 $2.1$ 倍。$0.2\pi$ 的脈衝則早到 9 個樣本、振幅比預測少 0.2——下面「三個值得注意的細節」會解釋原因,那與其說是預測錯了,不如說是預測本身的適用界限。$0.4\pi$ のパルスは予測どおりの位置にぴたりと現れ、振幅は $2.1$ 倍になる。$0.2\pi$ のパルスは 9 サンプル早く着き、振幅は予測より 0.2 低い。その理由は下の「注目すべき三つの細部」で説明する。これは予測が誤っているというより、予測が通用する範囲の限界である。
What actually comes out. Compare the output panel against the prediction:實際的輸出。把輸出面板和預測對照:実際に出てくるもの。出力のパネルを予測と照らし合わせる。
The order is reversed.脈衝的順序被顛倒了。パルスの順序が入れ替わる。
| First先出現先に現れる | Second後出現あとに現れる | |
|---|---|---|
| Input輸入入力 | $0.2\pi$ ($n \approx 91$) | $0.4\pi$ ($n \approx 152$) |
| Output輸出出力 | $0.4\pi$ ($n \approx 159$) | $0.2\pi$ ($n \approx 232$) |
The $0.2\pi$ burst entered the system 61 samples before the $0.4\pi$ burst and came out 73 samples after it. This is phase distortion in its most extreme form. Each burst on its own still looks intact, because each is narrowband and the phase is nearly linear within its own band. But the signal as a whole is made of components at different frequencies, their delays differ by some 140 samples, and the waveform they combine into is a completely different one.$0.2\pi$ 的脈衝比 $0.4\pi$ 早 61 個樣本進入系統,卻晚了 73 個樣本才出來。這是相位失真最極端的形式。每個脈衝單獨看起來都還算完整,因為各自頻寬很窄,在自己的頻帶內相位近似線性。但整個訊號是由不同頻率的成分組成的,各成分的延遲差了約 140 個樣本,組合起來的波形就完全不一樣了。$0.2\pi$ のパルスは $0.4\pi$ のパルスより 61 サンプル早くシステムに入り、73 サンプル遅れて出てきた。これは位相歪みのもっとも極端な形である。パルスは一つずつ見ればまだ無事に見える。どれも狭帯域で、それぞれの帯域の中では位相がほぼ線形だからである。しかし信号全体は異なる周波数の成分からできており、その遅延はおよそ 140 サンプルも違う。組み合わさってできる波形はまったく別のものになる。
How this differs from the Gaussian-pulse figure above. There, one pulse had a bandwidth that straddled a region where the group delay was changing, and the pulse itself was stretched and deformed. Here, different pulses sit at frequencies with different group delays; each keeps its shape more or less intact, but their relative timing is thrown into disarray. Both are phase distortion, happening at different scales.和上面高斯脈衝互動圖的差別。在那裡,是同一個脈衝的頻寬跨越了群延遲變化的區域,結果脈衝本身被拉長、變形。在這個例子中,是不同的脈衝落在群延遲不同的頻率,每個脈衝大致保持形狀,但彼此的相對時間被打亂了。兩者都是相位失真,只是發生在不同的尺度上。上のガウスパルスの図との違い。あちらは一つのパルスの帯域幅が群遅延の変化する領域にまたがり、パルスそのものが引き伸ばされて変形した。この例では別々のパルスが群遅延の異なる周波数に置かれており、各パルスはおおむね形を保つものの、互いの時間関係が崩される。どちらも位相歪みであり、起きている尺度が違うだけである。
Three details worth noticing.三個值得注意的細節。注目すべき三つの細部。
The figure below reproduces the textbook example in full. The system is a cascade of two parts:下面用課本的系統把這個例子完整重現。系統由兩部分串接而成:下の図は教科書の例をそのまま再現したものである。システムは二つの部分の縦続接続でできている。
The output is computed sample by sample from the difference equations, with no approximation.輸出是用差分方程逐點算出來的,沒有使用任何近似。出力は差分方程式から一点ずつ計算しており、近似は一切使っていない。
How to read the figure.怎麼讀這張圖。図の読み方。
Try these in order:建議依序嘗試:次の順に試してみるとよい。
| burst脈衝パルス | input centre輸入中心入力の中心 | |H(ωc)| | τg(ωc) | predicted centre預測中心予測中心 | measured centre實測中心実測中心 | output peak輸出峰值出力のピーク |
|---|
Why the blue burst's measured centre differs from its predicted centre. The predicted centre uses $\tau_g$ at the centre frequency $0.2\pi$ — the very top of the group-delay spike. But the blue burst's bandwidth (about $\pm 0.066\pi$, see the band) is about as wide as the spike itself, so the components at the edges of the band are delayed by only a few tens of samples, much less than the centre.為什麼藍色的「實測中心」和「預測中心」有差距。預測中心用的是中心頻率 $0.2\pi$ 的 $\tau_g$,也就是群延遲尖峰的頂端。但藍色脈衝的頻寬(約 $\pm 0.066\pi$,見色帶)和群延遲尖峰的寬度差不多,所以頻帶邊緣的成分只被延遲了幾十個樣本,比中心少很多。青の「実測中心」と「予測中心」がずれる理由。予測中心は中心周波数 $0.2\pi$ での $\tau_g$、すなわち群遅延の尖りの頂上を使っている。しかし青のパルスの帯域幅(約 $\pm 0.066\pi$、色の帯を参照)は尖りの幅とほぼ同じなので、帯域の端の成分は数十サンプルしか遅れず、中心よりずっと少ない。
The measured centre is the energy-weighted mean position of the output, roughly “the $\tau_g$ of the components in the band, averaged with their energies as weights”, so it comes out somewhat below the peak value. The size of this gap is how far the local-linear approximation has failed — the dispersion discussed above:實測中心是輸出能量的加權平均位置,大致等於「頻帶內各成分的 $\tau_g$ 依能量加權平均」,所以會比尖峰值小一些。這個差距就是局部線性近似失效的程度,也就是前面說的色散:実測中心は出力のエネルギーで重み付けした平均位置であり、おおよそ「帯域内の各成分の $\tau_g$ をエネルギーで重み付けして平均したもの」なので、ピークの値より少し小さくなる。このずれの大きさが局所線形近似が破れている程度、すなわち上で述べた分散である。
| Observation觀察観察 | Underlying principle對應的原理対応する原理 |
|---|---|
| the $0.8\pi$ burst disappears$0.8\pi$ 脈衝消失$0.8\pi$ のパルスが消える | zeros close to the unit circle, $|H| \approx 0$零點在單位圓附近,$|H| \approx 0$零点が単位円の近くにあり、$|H| \approx 0$ |
| the $0.4\pi$ burst is amplified about 2×$0.4\pi$ 脈衝被放大約 2 倍$0.4\pi$ のパルスが約 2 倍に増幅される | poles near the $0.4\pi$ direction極點靠近 $0.4\pi$ 方向極が $0.4\pi$ の方向の近くにある |
| the $0.2\pi$ burst is about 150 samples late$0.2\pi$ 脈衝晚了約 150 個樣本$0.2\pi$ のパルスが約 150 サンプル遅れる | the allpass poles cluster near $0.2\pi$ and create a large group delay全通極點聚集在 $0.2\pi$ 附近,產生大量群延遲全域通過の極が $0.2\pi$ 付近に集まり、大きな群遅延を作る |
| the order of the bursts is reversed脈衝順序顛倒パルスの順序が入れ替わる | the group delay differs strongly between frequencies — phase distortion不同頻率的群延遲差異很大,這是相位失真周波数によって群遅延が大きく違う。これが位相歪みである |
| removing the allpass leaves $|H|$ unchanged, yet the order is restored移除全通後 $|H|$ 不變,順序卻恢復全域通過部を外しても $|H|$ は変わらないのに、順序は元に戻る | the information about time structure is hidden entirely in the phase時間結構的資訊完全藏在相位中時間構造の情報はすべて位相の中に隠れている |
| the blue burst's prediction and measurement differ, and its shape is distorted藍色脈衝預測與實測有落差,形狀被扭曲青のパルスの予測と実測がずれ、形が歪む | $\tau_g$ is not constant across the band; local linearity fails (dispersion)頻帶內 $\tau_g$ 不是常數,局部線性失效(色散)帯域内で $\tau_g$ が一定でなく、局所線形性が破れる(分散) |
The slide ends with “→ correction…”, a pointer to the next topic. Here is the direction in brief.投影片最後寫「→ correction…」,這是下一個主題的伏筆,這裡先簡單說明方向。スライドは「→ correction…」で終わっており、次の話題への伏線になっている。ここでは方向だけを簡単に述べる。
Group delays add. Frequency responses of cascaded systems multiply, their phases add, and so their group delays add directly:群延遲可以相加。串接系統的頻率響應相乘,相位相加,所以群延遲也直接相加:群遅延は足し合わせられる。縦続接続したシステムの周波数応答は掛け合わされ、位相は足し合わされるので、群遅延もそのまま足し合わされる。
That gives a natural way to correct: if the filter's $\tau_g$ has a hump somewhere in the passband, cascade an allpass equalizer whose $\tau_g$ is low there and high elsewhere, so the sum is close to constant. Because an allpass has $|H| = 1$, the magnitude response is untouched; only the phase is repaired.這給了一個很自然的校正方法:如果濾波器的 $\tau_g$ 在通帶中某處有隆起,就串接一個全通等化器(allpass equalizer),讓它的 $\tau_g$ 在那裡低、在其他地方高,兩者相加後變成接近常數。因為全通濾波器的 $|H| = 1$,大小響應完全不受影響,只修正相位。ここから自然な補正の方法が得られる。フィルタの $\tau_g$ が通過域のどこかで盛り上がっていれば、そこで $\tau_g$ が低く、ほかで高い全域通過イコライザ(allpass equalizer)を縦続接続し、和がほぼ一定になるようにする。全域通過フィルタは $|H| = 1$ なので、振幅応答はまったく変わらず、位相だけが直される。
Other ways to avoid phase distortion.其他避免相位失真的方法。位相歪みを避けるほかの方法。
filtfilt works this way. Because it needs “future” samples it can only be used for non-real-time processing.零相位濾波(離線處理):把訊號先順向通過濾波器,再把結果反轉後再通過一次。兩次的相位互相抵消,總相位為 0,群延遲為 0。MATLAB 的 filtfilt 就是這樣做的。因為需要用到「未來」的樣本,所以只能用於非即時處理。ゼロ位相フィルタリング(オフライン処理):信号をまず順方向にフィルタに通し、結果を反転してもう一度通す。二回の位相が打ち消し合い、全体の位相は 0、群遅延も 0 になる。MATLAB の filtfilt はこの方法である。「未来」のサンプルを使うので、非実時間の処理にしか使えない。Where it matters.在哪些地方重要。どこで重要になるか。
| Field領域分野 | Problem caused by a non-flat group delay群延遲不平坦造成的問題平坦でない群遅延が引き起こす問題 |
|---|---|
| Digital communications數位通訊ディジタル通信 | symbol envelopes are stretched into neighbouring symbols, causing intersymbol interference (ISI)符號的包絡被拉長,侵入相鄰符號,造成符號間干擾(ISI)符号の包絡線が引き伸ばされて隣の符号に入り込み、符号間干渉(ISI)を起こす |
| Audio音訊オーディオ | transients such as drum hits and plucked strings become blurred; the phase design of loudspeaker crossovers鼓聲、撥弦等暫態聲音變得模糊;喇叭分音器的相位設計ドラムや撥弦などの過渡的な音がぼやける。スピーカーのクロスオーバーの位相設計 |
| Measurement, biomedical signals量測、生醫訊號計測、生体信号 | for waveforms such as the ECG the shape itself is the information; distortion leads to misreadings心電圖等波形的形狀本身就是資訊,失真會造成誤判心電図などは波形の形そのものが情報であり、歪むと誤判定につながる |
| Radar雷達レーダー | dispersion is used deliberately: transmit a stretched chirp and compress it back into a narrow pulse with the opposite dispersion at the receiver (pulse compression)反過來刻意利用色散:發射被拉長的 chirp,接收端用相反的色散壓縮回窄脈衝(pulse compression)逆に分散を意図的に利用する。引き伸ばしたチャープを送信し、受信側で逆の分散によって狭いパルスに圧縮し直す(パルス圧縮) |
The last row is the interesting one: dispersion is not necessarily bad. As long as you know its form, you can exploit it, or compensate for it exactly.最後一項很有意思:色散不一定是壞事,只要知道它的形式,就能利用它,或者精確地補償回來。最後の行は興味深い。分散は必ずしも悪者ではない。その形がわかっていれば、利用することも、正確に補償して元に戻すこともできる。
“Phase delay and group delay are two different times.” Both are steady-state properties, but they look at different things. The phase delay looks at the carrier — the crests of a single frequency; the group delay looks at the envelope — the information. For carrying information, the group delay is what really matters.「相位延遲和群延遲是兩種不同的時間」:它們都是穩態性質,但看的對象不同。相位延遲看載波,也就是單一頻率的波峰;群延遲看包絡,也就是資訊。對傳遞資訊而言,群延遲才是真正重要的。「位相遅延と群遅延は二種類の違う時間である」:どちらも定常状態の性質だが、見ている対象が違う。位相遅延は搬送波、つまり単一周波数の山を見ており、群遅延は包絡線、つまり情報を見ている。情報を伝えるうえで本当に重要なのは群遅延である。
“The phase delay can be measured directly.” A single steady-state sinusoid has no time marker; its phase delay is ambiguous by integer multiples of $2\pi/\omega$ and can even count a sign flip as a time shift (the $\omega \gt 0.4$ case in the first interactive figure, §10.3.1). The group delay depends only on the slope of the phase, is insensitive to $2\pi$ wrapping and to constant phase offsets, and can actually be measured from the envelope.「相位延遲可以直接量測」:單一穩態正弦波沒有時間標記,相位延遲有 $2\pi/\omega$ 整數倍的歧義,還會把反相誤算成時間平移(§10.3.1 互動圖一的 $\omega \gt 0.4$ 情況)。群延遲只依賴相位的斜率,對 $2\pi$ 折返和常數相位偏移都不敏感,可以透過包絡實際量測。「位相遅延は直接測れる」:定常正弦波一つには時刻の目印がなく、位相遅延には $2\pi/\omega$ の整数倍のあいまいさがあり、符号の反転まで時間のずれと数えてしまう(§10.3.1 の一つ目の図の $\omega \gt 0.4$ の場合)。群遅延は位相の傾きだけに依存し、$2\pi$ の折り返しにも一定の位相オフセットにも影響されず、包絡線から実際に測ることができる。
“If $|H|$ is flat there is no distortion.” Wrong. The allpass in the Gaussian-pulse figure above has $|H|$ exactly 1, yet a wideband pulse is badly distorted. Distortionless transmission needs two conditions: $|H|$ constant across the signal band, and $\tau_g$ constant across the signal band.「|H| 平坦就不會失真」:錯。上面高斯脈衝互動圖中的全通濾波器 $|H|$ 精確等於 1,但寬頻脈衝仍然被嚴重扭曲。波形不失真需要兩個條件:$|H|$ 在訊號頻帶內是常數,而且 $\tau_g$ 在訊號頻帶內是常數。「|H| が平坦なら歪まない」:誤り。上のガウスパルスの図の全域通過フィルタは $|H|$ が厳密に 1 だが、広帯域パルスはやはり大きく歪む。波形を歪ませないには二つの条件が要る。信号帯域内で $|H|$ が一定であること、かつ信号帯域内で $\tau_g$ が一定であることである。
“A negative group delay violates causality.” As discussed in §10.3.2: a negative group delay appears near a zero inside the unit circle, where $|H|$ is small and the signal heavily attenuated, and it holds only under the narrowband approximation. The signal's energy never comes out earlier than the input.「負的群延遲違反因果律」:先前(§10.3.2)討論過:負群延遲出現在單位圓內的零點附近,那裡 $|H|$ 很小,訊號被大幅衰減;而且它只在窄頻近似下成立。訊號的能量不會比輸入更早出現。「負の群遅延は因果律に反する」:§10.3.2 で述べたとおり、負の群遅延は単位円内の零点の近くに現れ、そこでは $|H|$ が小さく信号は大きく減衰している。しかも狭帯域近似のもとでしか成り立たない。信号のエネルギーが入力より早く出てくることはない。
| Item項目項目 | Phase delay τp相位延遲 τp位相遅延 τp | Group delay τg群延遲 τg群遅延 τg |
|---|---|---|
| Definition定義定義 | $-\dfrac{\theta(\omega)}{\omega}$ | $-\dfrac{d\theta}{d\omega}$ |
| Geometric meaning幾何意義幾何学的な意味 | slope of the chord from the origin to the point on the curve (negated)原點到曲線上點的割線斜率(取負)原点から曲線上の点への弦の傾き(符号反転) | slope of the tangent to the curve (negated)曲線的切線斜率(取負)曲線の接線の傾き(符号反転) |
| Describes描述對象表すもの | the carrier (where the crests of a single frequency are)載波(單一頻率的波峰位置)搬送波(単一周波数の山の位置) | the envelope (when the information and energy arrive)包絡(資訊、能量到達的時間)包絡線(情報やエネルギーが届く時刻) |
| Derived from推導來源導出のもと | the average of the two sideband phases $\frac{\theta_u+\theta_l}{2}$兩邊帶相位的平均 $\frac{\theta_u+\theta_l}{2}$二つの側波帯の位相の平均 $\frac{\theta_u+\theta_l}{2}$ | half the difference of the two sideband phases $\frac{\theta_u-\theta_l}{2}$兩邊帶相位的半差 $\frac{\theta_u-\theta_l}{2}$二つの側波帯の位相の差の半分 $\frac{\theta_u-\theta_l}{2}$ |
| $2\pi$ wrapping or constant phase offset對 $2\pi$ 折返或常數相位偏移$2\pi$ の折り返しや一定の位相オフセット | sensitive, ambiguous敏感,有歧義影響を受け、あいまい | unaffected不受影響影響を受けない |
| 3-point FIR filter3 點 FIR 濾波器3 点 FIR フィルタ | 1 for $\omega \lt 0.4$, $1 - \pi/\omega$ for $\omega \gt 0.4$$\omega \lt 0.4$ 時為 1,$\omega \gt 0.4$ 時為 $1 - \pi/\omega$$\omega \lt 0.4$ で 1、$\omega \gt 0.4$ で $1 - \pi/\omega$ | 1 everywhere (except at the zero)處處為 1(零點處除外)至るところ 1(零点を除く) |
| Allpass filter at its resonance全通濾波器在共振頻率共振周波数での全域通過フィルタ | small (about 3)小(約 3)小さい(約 3) | large (about 19)大(約 19)大きい(約 19) |
(Supplementary to the slides.) Three classical families of systems are defined directly by their pole-zero geometry, and together they organize everything in 10.2–10.3: all-pass systems carry only phase, minimum-phase systems carry the least possible phase for a given magnitude, and linear-phase FIR systems carry a pure delay phase.(投影片之外的補充。)有三個經典的系統家族,直接由極零幾何定義,而它們合起來正好把 10.2–10.3 的內容整理了一遍:全通系統只帶相位,最小相位系統在給定振幅下帶最少可能的相位,線性相位 FIR 系統則帶一個純延遲式的相位。(スライドへの補足。)極零配置によって直接定義される古典的なシステムの族が三つあり、それらを合わせると 10.2–10.3 の内容がきれいに整理される。全域通過システムは位相だけを担い、最小位相システムは与えられた振幅に対してとりうる最小の位相を担い、線形位相 FIR システムは純粋な遅延としての位相を担う。
An all-pass system(全通系統) has $|H_{ap}(e^{j\omega})| = $ constant for all $\omega$. The construction: pair every pole $a$ with a zero at its conjugate reciprocal $1/a^{*}$:全通系統是指對所有 $\omega$ 都有 $|H_{ap}(e^{j\omega})| = $ 常數的系統。做法是:把每個極點 $a$ 與一個位於其共軛倒數 $1/a^{*}$ 的零點配成一對:全域通過システムとは、すべての $\omega$ について $|H_{ap}(e^{j\omega})| = $ 一定となるシステムである。作り方はこうである。各々の極 $a$ に対し、その共役逆数 $1/a^{*}$ に零点を置いて対にする。
On the unit circle the vector to the zero at $1/a^{*}$ is always proportional in length to the vector to the pole at $a$ — the ratio of lengths is constant, so the magnitude never moves. Only the phase (and hence the group delay) varies.在單位圓上,指向 $1/a^{*}$ 那個零點的向量,長度永遠與指向極點 $a$ 的向量成固定比例——長度之比是常數,所以振幅動也不動。只有相位(連帶群延遲)會變。単位円上では、$1/a^{*}$ の零点へ向かうベクトルの長さは、極 $a$ へ向かうベクトルの長さにつねに比例する。長さの比が一定なので、振幅はまったく動かない。変わるのは位相(したがって群遅延)だけである。
A causal stable system is minimum-phase(最小相位) if all its zeros (as well as poles) are inside the unit circle. Then the inverse system $1/H(z)$ — whose poles are the zeros of $H(z)$ — is also causal and stable: minimum-phase systems are exactly the ones that can be undone by a realizable filter (e.g. channel equalization).一個因果穩定的系統,若它的零點(連同極點)全都在單位圓內,就稱為最小相位系統。此時反系統 $1/H(z)$ —— 它的極點就是 $H(z)$ 的零點 —— 也是因果且穩定的:最小相位系統正好就是那些能被可實現濾波器還原回去的系統(例如通道等化)。因果的かつ安定なシステムは、その零点が(極と同様に)すべて単位円の内側にあるとき最小位相であるという。このとき逆システム $1/H(z)$ —— その極は $H(z)$ の零点である —— もまた因果的かつ安定になる。最小位相システムとは、実現可能なフィルタによって元に戻せるシステムにほかならない(たとえば通信路の等化)。
Reflecting a zero $\zeta$ to $1/\zeta^{*}$ multiplies the response by an all-pass factor, which leaves $|H(e^{j\omega})|$ unchanged. Hence many systems share one magnitude, and:把零點 $\zeta$ 反射到 $1/\zeta^{*}$,等於把響應乘上一個全通因子,而這並不改變 $|H(e^{j\omega})|$。因此許多系統會共用同一個振幅,於是:零点 $\zeta$ を $1/\zeta^{*}$ へ反射させることは、応答に全域通過の因子を掛けることにあたり、$|H(e^{j\omega})|$ は変わらない。したがって多くのシステムがひとつの振幅を共有することになり、次が成り立つ。
Every rational $H(z)$ is a minimum-phase system (carrying all of the magnitude) cascaded with an all-pass (carrying the excess phase). Among all systems with the same $|H(e^{j\omega})|$, the minimum-phase one has:任何有理的 $H(z)$,都可以寫成一個最小相位系統(承擔全部振幅)串接一個全通系統(承擔多餘的相位)。在所有具有相同 $|H(e^{j\omega})|$ 的系統之中,最小相位的那一個具有:任意の有理な $H(z)$ は、最小位相システム(振幅のすべてを担う)と全域通過システム(余分な位相を担う)の縦続接続として書ける。同じ $|H(e^{j\omega})|$ をもつすべてのシステムのうち、最小位相のものは次の性質をもつ。
A system has (generalized) linear phase(線性相位) if $H(e^{j\omega}) = A(\omega)\,e^{-j(\alpha\omega - \beta)}$ with $A(\omega)$ real. Then $\tau_g(\omega) = \alpha$ is constant — no phase distortion, every envelope is delayed by the same $\alpha$ samples. For causal FIR filters of length $M$ this happens exactly when the impulse response is symmetric or antisymmetric:若一個系統可以寫成 $H(e^{j\omega}) = A(\omega)\,e^{-j(\alpha\omega - \beta)}$,其中 $A(\omega)$ 為實數,就說它具有(廣義)線性相位。此時 $\tau_g(\omega) = \alpha$ 為常數——沒有相位失真,每個包絡都被延遲同樣的 $\alpha$ 個樣本。對長度為 $M$ 的因果 FIR 濾波器,這件事發生的充要條件就是脈衝響應對稱或反對稱:システムが実数の $A(\omega)$ を用いて $H(e^{j\omega}) = A(\omega)\,e^{-j(\alpha\omega - \beta)}$ と書けるとき、そのシステムは(一般化)線形位相をもつという。このとき $\tau_g(\omega) = \alpha$ は一定であり、位相歪みはなく、どの包絡線も同じ $\alpha$ 標本だけ遅れる。長さ $M$ の因果的 FIR フィルタでは、これが成り立つのはインパルス応答が対称または反対称であるとき、かつそのときにかぎる。
The symmetry constrains the zero locations: $z^{-(M-1)}H(1/z) = \pm H(z)$, so if $\zeta$ is a zero, then $1/\zeta$, $\zeta^{*}$, $1/\zeta^{*}$ are all zeros. The zero constellation therefore consists of:這個對稱性限制了零點的位置:$z^{-(M-1)}H(1/z) = \pm H(z)$,所以只要 $\zeta$ 是零點,$1/\zeta$、$\zeta^{*}$、$1/\zeta^{*}$ 也都是零點。因此零點的分布必然由下列幾種組合構成:この対称性は零点の位置を縛る。$z^{-(M-1)}H(1/z) = \pm H(z)$ なので、$\zeta$ が零点ならば $1/\zeta$、$\zeta^{*}$、$1/\zeta^{*}$ もすべて零点である。したがって零点の配置は次のものから成る。
filtfilt).線性相位 FIR 不可能是最小相位的(它的零點成對出現為 $\zeta, 1/\zeta$,所以除了退化情形之外,總有一些落在圓上或圓外)——零相位失真的代價,就是固定的延遲 $(M-1)/2$。更一般地說,因果性本身就逼出非零的相位響應:除非 $h[n] = c\,\delta[n]$,否則因果系統不可能有 $\theta(\omega) \equiv 0$。零相位濾波只能離線做(非因果處理,例如 MATLAB 的 filtfilt)。線形位相 FIR は最小位相にはなれない(零点が $\zeta, 1/\zeta$ の対で現れるため、退化した場合を除いていくつかは円上か円の外に出る)。位相歪みをゼロにする代価が、固定遅延 $(M-1)/2$ なのである。より一般に、因果性そのものが非零の位相応答を強いる。$h[n] = c\,\delta[n]$ でないかぎり、因果的なシステムが $\theta(\omega) \equiv 0$ をもつことはない。零位相フィルタリングはオフラインでしかできない(非因果的な処理、たとえば MATLAB の filtfilt)。
Problem: A 3-point moving-average filter $h[n] = \{\tfrac13, \tfrac13, \tfrac13\}$ (for $n = 0,1,2$) is driven by $x[n] = \cos\!\big(\tfrac{\pi}{4}n + \tfrac{\pi}{6}\big)$. Find the steady-state output $y[n]$ exactly, and interpret the phase shift as a time delay.題目:一個 3 點移動平均濾波器 $h[n] = \{\tfrac13, \tfrac13, \tfrac13\}$($n = 0,1,2$)由 $x[n] = \cos\!\big(\tfrac{\pi}{4}n + \tfrac{\pi}{6}\big)$ 驅動。請精確求出穩態輸出 $y[n]$,並把相位平移解讀為時間延遲。問題:3 点移動平均フィルタ $h[n] = \{\tfrac13, \tfrac13, \tfrac13\}$($n = 0,1,2$)を $x[n] = \cos\!\big(\tfrac{\pi}{4}n + \tfrac{\pi}{6}\big)$ で駆動する。定常状態の出力 $y[n]$ を厳密に求め、位相のずれを時間遅れとして解釈せよ。
Frequency response (Example 10-A with $M = 3$, or directly):頻率響應(取範例 10-A 的 $M = 3$,或直接計算):周波数応答(例題 10-A で $M = 3$ とするか、直接求める)。
$$H(e^{j\omega}) = \frac13\big(1 + e^{-j\omega} + e^{-2j\omega}\big) = \frac13 e^{-j\omega}\big(1 + 2\cos\omega\big)$$At $\omega_0 = \pi/4$: $\cos(\pi/4) = \tfrac{\sqrt2}{2}$, so在 $\omega_0 = \pi/4$ 處:$\cos(\pi/4) = \tfrac{\sqrt2}{2}$,所以$\omega_0 = \pi/4$ では $\cos(\pi/4) = \tfrac{\sqrt2}{2}$ なので、
$$H(e^{j\pi/4}) = \frac{1+\sqrt{2}}{3}\,e^{-j\pi/4} \approx 0.805\,e^{-j\pi/4}$$Real input, real $h[n]$ → apply the boxed result of 10.1.2:實數輸入、實數 $h[n]$ → 套用 10.1.2 的方框結論:実入力で $h[n]$ も実数 → 10.1.2 の囲みの結論を適用する。
$$y[n] = 0.805\cos\Big(\frac{\pi}{4}n + \frac{\pi}{6} - \frac{\pi}{4}\Big) = 0.805\cos\Big(\frac{\pi}{4}(n-1) + \frac{\pi}{6}\Big)$$The phase shift $-\pi/4$ at frequency $\pi/4$ is exactly a one-sample delay ($\tau_p = \frac{\pi/4}{\pi/4} = 1$), consistent with linear phase $\theta(\omega) = -\omega$ and the symmetry center of $h[n]$ at $n = 1$.在頻率 $\pi/4$ 處的相位平移 $-\pi/4$,恰好就是一個樣本的延遲($\tau_p = \frac{\pi/4}{\pi/4} = 1$),這與線性相位 $\theta(\omega) = -\omega$ 以及 $h[n]$ 的對稱中心落在 $n = 1$ 完全一致。周波数 $\pi/4$ における位相のずれ $-\pi/4$ は、ちょうど1 標本の遅れ($\tau_p = \frac{\pi/4}{\pi/4} = 1$)である。これは線形位相 $\theta(\omega) = -\omega$ とも、$h[n]$ の対称中心が $n = 1$ にあることとも整合している。
Problem: (a) Derive the closed-form frequency response of the $M$-point moving average $y[n] = \frac1M\sum_{\ell=0}^{M-1}x[n-\ell]$. (b) For $M = 4$, list all frequencies in $(0, \pi]$ at which the response is exactly zero, and explain them with the zero locations of $H(z)$. (c) Which input sinusoids are completely blocked?題目:(a) 推導 $M$ 點移動平均 $y[n] = \frac1M\sum_{\ell=0}^{M-1}x[n-\ell]$ 的閉合形式頻率響應。(b) 取 $M = 4$,列出 $(0, \pi]$ 中所有使響應恰為零的頻率,並用 $H(z)$ 的零點位置解釋它們。(c) 哪些輸入弦波會被完全擋掉?問題:(a) $M$ 点移動平均 $y[n] = \frac1M\sum_{\ell=0}^{M-1}x[n-\ell]$ の周波数応答を閉じた形で導け。(b) $M = 4$ のとき、$(0, \pi]$ のうち応答がちょうど 0 になる周波数をすべて挙げ、$H(z)$ の零点の位置から説明せよ。(c) 完全に阻止される入力正弦波はどれか。
(a) Geometric sum:(a) 等比級數和:(a) 等比級数和を使う。
$$H(e^{j\omega}) = \frac1M\,\frac{1-e^{-j\omega M}}{1-e^{-j\omega}} = \frac1M\,e^{-j\omega\frac{M-1}{2}}\,\frac{\sin(M\omega/2)}{\sin(\omega/2)}$$(b) Zeros where $\sin(M\omega/2) = 0$ but $\sin(\omega/2) \ne 0$: $\omega = 2\pi k/M$, $k = 1, \dots, M-1$. For $M = 4$:(b) 零點出現在 $\sin(M\omega/2) = 0$ 但 $\sin(\omega/2) \ne 0$ 之處:$\omega = 2\pi k/M$,$k = 1, \dots, M-1$。取 $M = 4$:(b) 零点は $\sin(M\omega/2) = 0$ かつ $\sin(\omega/2) \ne 0$ の点、すなわち $\omega = 2\pi k/M$、$k = 1, \dots, M-1$ である。$M = 4$ なら、
$$\omega = \frac{\pi}{2}, \ \pi \quad (k = 1, 2; \ k=3 \text{ gives } \tfrac{3\pi}{2} \notin (0,\pi])$$In the z-plane, $H(z) = \frac{z^4 - 1}{4z^3(z-1)}$: the 4th roots of unity $\{1, j, -1, -j\}$ minus the cancelled root at $z=1$ leave zeros at $z = j, -1, -j$, i.e. on the unit circle at angles $\pm\pi/2$ and $\pi$ — exactly the magnitude nulls (the vector from an on-circle zero to $e^{j\omega}$ has length zero there).在 z 平面上 $H(z) = \frac{z^4 - 1}{4z^3(z-1)}$:4 次單位根 $\{1, j, -1, -j\}$ 扣掉在 $z=1$ 被對消的那個根,剩下 $z = j, -1, -j$ 這三個零點,也就是單位圓上角度 $\pm\pi/2$ 與 $\pi$ 的位置——正好就是振幅零點(從圓上零點指向 $e^{j\omega}$ 的向量在那裡長度為零)。z 平面では $H(z) = \frac{z^4 - 1}{4z^3(z-1)}$ である。1 の 4 乗根 $\{1, j, -1, -j\}$ から $z=1$ で相殺された根を除くと、零点は $z = j, -1, -j$ に残る。すなわち単位円上の角度 $\pm\pi/2$ と $\pi$ の位置であり、これが振幅の零点そのものである(円上の零点から $e^{j\omega}$ へのベクトルの長さがそこで 0 になる)。
(c) Any sinusoid $A\cos(\omega n + \phi)$ with $\omega = \pi/2$ or $\omega = \pi$ produces zero steady-state output — e.g. $\{\dots, 1, 0, -1, 0, \dots\}$ and $\{\dots, 1, -1, 1, -1, \dots\}$: every 4 consecutive samples average to 0.(c) 任何 $\omega = \pi/2$ 或 $\omega = \pi$ 的弦波 $A\cos(\omega n + \phi)$ 的穩態輸出都是零——例如 $\{\dots, 1, 0, -1, 0, \dots\}$ 與 $\{\dots, 1, -1, 1, -1, \dots\}$:任意連續 4 個樣本的平均都是 0。(c) $\omega = \pi/2$ または $\omega = \pi$ の正弦波 $A\cos(\omega n + \phi)$ は、どれも定常状態の出力が 0 になる。たとえば $\{\dots, 1, 0, -1, 0, \dots\}$ や $\{\dots, 1, -1, 1, -1, \dots\}$ であり、連続する 4 標本の平均がつねに 0 になるからである。
Problem: A causal system satisfies題目:某因果系統滿足問題:ある因果的なシステムが次を満たすとする。 $$y[n] = x[n] - 0.5x[n-1] + 1.1y[n-1] - 0.3y[n-2].$$ (a) Find $H(z)$ and factor it. (b) Locate poles and zeros. (c) Give the ROC and decide stability. (d) Could a different (non-causal) system with the same $H(z)$ expression be stable?(a) 求出 $H(z)$ 並因式分解。(b) 找出極點與零點。(c) 給出收斂區域並判斷穩定性。(d) 具有相同 $H(z)$ 表示式的另一個(非因果)系統有可能穩定嗎?(a) $H(z)$ を求めて因数分解せよ。(b) 極と零点の位置を求めよ。(c) 収束領域を示し、安定性を判定せよ。(d) 同じ $H(z)$ の式をもつ別の(非因果的な)システムが安定になりうるか。
(a) Transforming with $Y$ terms collected:(a) 把 $Y$ 的各項集中之後取轉換:(a) $Y$ の項をまとめてから変換する。
$$H(z) = \frac{1 - 0.5z^{-1}}{1 - 1.1z^{-1} + 0.3z^{-2}} = \frac{1 - 0.5z^{-1}}{(1 - 0.5z^{-1})(1 - 0.6z^{-1})} = \frac{1}{1 - 0.6z^{-1}}$$The denominator factors via $\lambda^2 - 1.1\lambda + 0.3 = (\lambda - 0.5)(\lambda - 0.6)$.分母可以分解為 $\lambda^2 - 1.1\lambda + 0.3 = (\lambda - 0.5)(\lambda - 0.6)$。分母は $\lambda^2 - 1.1\lambda + 0.3 = (\lambda - 0.5)(\lambda - 0.6)$ と因数分解できる。
(b) Zero at $z = 0.5$, poles at $z = 0.5$ and $z = 0.6$. The pole at $0.5$ is cancelled by the zero (a pole-zero cancellation — "near poles and zeros cancel", here exactly). Effective system: single pole at $z = 0.6$, $h[n] = (0.6)^n u[n]$.(b) 零點在 $z = 0.5$,極點在 $z = 0.5$ 與 $z = 0.6$。位於 $0.5$ 的極點被零點對消掉了(極零對消——「靠得很近的極零對會互相抵消」,這裡則是完全抵消)。等效的系統是:單一極點位於 $z = 0.6$,$h[n] = (0.6)^n u[n]$。(b) 零点は $z = 0.5$、極は $z = 0.5$ と $z = 0.6$ にある。$0.5$ の極は零点によって相殺される(極零相殺であり、「近接した極と零点は打ち消し合う」が、ここでは厳密に打ち消し合う)。実効的なシステムは $z = 0.6$ に極をひとつもつもので、$h[n] = (0.6)^n u[n]$ となる。
(c) Causal → ROC $|z| > 0.6$. This includes the unit circle → stable.(c) 因果 → 收斂區域 $|z| > 0.6$。它包含單位圓 → 穩定。(c) 因果的 → 収束領域は $|z| > 0.6$。これは単位円を含む → 安定である。
(d) The other ROC choice $|z| < 0.6$ (left-sided $h[n] = -(0.6)^n u[-n-1]$) does not include the unit circle → unstable. So for this $H(z)$ only the causal version is stable. (In general a stable non-causal choice exists only when some pole lies outside the unit circle and the annular ROC between pole radii still contains $|z| = 1$.)(d) 另一種收斂區域的選擇 $|z| < 0.6$(左側序列 $h[n] = -(0.6)^n u[-n-1]$)不包含單位圓 → 不穩定。所以對這個 $H(z)$ 而言,只有因果的版本是穩定的。(一般來說,只有當某個極點落在單位圓外、而極點半徑之間的環狀收斂區域仍然包含 $|z| = 1$ 時,才會存在穩定的非因果選擇。)(d) もう一方の収束領域 $|z| < 0.6$(左側系列 $h[n] = -(0.6)^n u[-n-1]$)は単位円を含まない → 不安定である。つまりこの $H(z)$ では、安定なのは因果的な版だけである。(一般に、安定な非因果的の選択が存在するのは、いずれかの極が単位円の外にあり、極の半径のあいだの円環状の収束領域がなお $|z| = 1$ を含む場合にかぎる。)
Problem: Consider the two-pole resonator題目:考慮下面這個雙極點共振器問題:次の 2 極共振器を考える。 $$H(z) = \frac{1}{(1 - 0.9e^{j\pi/3}z^{-1})(1 - 0.9e^{-j\pi/3}z^{-1})}.$$ Using only pole-vector geometry (no full evaluation): (a) near which frequency does $|H(e^{j\omega})|$ peak, and why? (b) Estimate $|H|$ at $\omega = \pi/3$. (c) What happens to the peak as the pole radius $r \to 1$? (d) Sketch the expected phase behavior near the peak.只用極點向量的幾何(不做完整計算):(a) $|H(e^{j\omega})|$ 會在哪個頻率附近出現尖峰?為什麼?(b) 估計 $\omega = \pi/3$ 處的 $|H|$。(c) 當極點半徑 $r \to 1$ 時,尖峰會怎麼變?(d) 畫出尖峰附近預期的相位行為。極ベクトルの幾何だけを使って(完全な評価はせずに)答えよ。(a) $|H(e^{j\omega})|$ はどの周波数の近くでピークをもつか。またそれはなぜか。(b) $\omega = \pi/3$ での $|H|$ を見積もれ。(c) 極の半径が $r \to 1$ となるとピークはどうなるか。(d) ピーク付近で予想される位相のふるまいを図示せよ。
(a) $|H(e^{j\omega})| = \dfrac{1}{|e^{j\omega} - 0.9e^{j\pi/3}|\;|e^{j\omega} - 0.9e^{-j\pi/3}|}$. The first pole-vector length is minimized when $e^{j\omega}$ passes closest to the pole, i.e. $\omega \approx \pi/3$ — the magnitude peaks there (the true maximum is a hair below $\pi/3$ because the conjugate pole's vector also shrinks slightly as $\omega$ decreases, but for $r = 0.9$ the shift is tiny).(a) $|H(e^{j\omega})| = \dfrac{1}{|e^{j\omega} - 0.9e^{j\pi/3}|\;|e^{j\omega} - 0.9e^{-j\pi/3}|}$。第一個極點向量的長度在 $e^{j\omega}$ 離該極點最近時達到最小,也就是 $\omega \approx \pi/3$——振幅就在那裡出現尖峰(真正的極大值會比 $\pi/3$ 低一點點,因為 $\omega$ 減小時共軛極點的向量也會稍微變短,但在 $r = 0.9$ 的情況下這個偏移小得可以忽略)。(a) $|H(e^{j\omega})| = \dfrac{1}{|e^{j\omega} - 0.9e^{j\pi/3}|\;|e^{j\omega} - 0.9e^{-j\pi/3}|}$ である。一つめの極ベクトルの長さは、$e^{j\omega}$ がその極に最も近づくとき、すなわち $\omega \approx \pi/3$ で最小になる。振幅はそこでピークをもつ(真の最大はごくわずかに $\pi/3$ より下である。$\omega$ が小さくなると共役な極のベクトルもわずかに縮むからだが、$r = 0.9$ ではこのずれはごく小さい)。
(b) At $\omega = \pi/3$:(b) 在 $\omega = \pi/3$ 處:(b) $\omega = \pi/3$ では次のようになる。
(c) As $r \to 1$ the near-pole distance $\to 0$, so the peak height $\to \infty$ (pole on the unit circle = oscillator) and the peak bandwidth shrinks $\propto (1-r)$.(c) 當 $r \to 1$ 時,到近極點的距離 $\to 0$,所以尖峰高度 $\to \infty$(極點落在單位圓上就是振盪器),而尖峰的頻寬則以 $\propto (1-r)$ 的方式收窄。(c) $r \to 1$ とすると、近い極までの距離が $\to 0$ となるので、ピークの高さは $\to \infty$ になる(極が単位円上に乗れば発振器である)。ピークの帯域幅は $\propto (1-r)$ で狭くなる。
(d) Phase = $-$(sum of pole-vector angles)$+2\omega$. As $e^{j\omega}$ sweeps past the near pole, that pole's vector angle rotates rapidly by nearly $\pi$, so $\theta(\omega)$ drops steeply through the resonance — and correspondingly $\tau_g = -d\theta/d\omega$ has a tall positive peak at $\omega \approx \pi/3$.(d) 相位 = $-$(極點向量角度之和)$+2\omega$。當 $e^{j\omega}$ 掃過近極點時,該極點的向量角度會快速轉過接近 $\pi$,所以 $\theta(\omega)$ 在共振處陡降——相對應地,$\tau_g = -d\theta/d\omega$ 在 $\omega \approx \pi/3$ 處出現一個高聳的正峰。(d) 位相 = $-$(極ベクトルの角度の和)$+2\omega$ である。$e^{j\omega}$ が近い極のそばを通り過ぎるとき、その極のベクトルの角度は $\pi$ 近くまで一気に回る。したがって $\theta(\omega)$ は共振を抜けるところで急降下し、対応して $\tau_g = -d\theta/d\omega$ は $\omega \approx \pi/3$ で高い正のピークをもつ。
Problem: For the symmetric 3-pt filter $h[n] = \{\alpha, \beta, \alpha\}$ with $H(e^{j\omega}) = e^{-j\omega}(\beta + 2\alpha\cos\omega)$: (a) Show that wherever $\beta + 2\alpha\cos\omega > 0$, $\tau_p(\omega) = \tau_g(\omega) = 1$ sample, and explain why this filter causes no phase distortion. (b) What happens to $\theta(\omega)$, $\tau_p$, and $\tau_g$ at a frequency $\omega_z$ where $\beta + 2\alpha\cos\omega$ changes sign? (c) An AM signal $A[n]\cos(\omega_c n)$ with narrowband $A[n]$ passes through a system with $\theta(\omega_c) = -1.2$ rad, $\frac{d\theta}{d\omega}|_{\omega_c} = -8$, $\omega_c = 0.3$. By how many samples are the carrier and the envelope delayed?題目:對於對稱的 3 點濾波器 $h[n] = \{\alpha, \beta, \alpha\}$,其 $H(e^{j\omega}) = e^{-j\omega}(\beta + 2\alpha\cos\omega)$:(a) 證明只要 $\beta + 2\alpha\cos\omega > 0$,就有 $\tau_p(\omega) = \tau_g(\omega) = 1$ 個樣本,並說明為何這個濾波器不會造成相位失真。(b) 在 $\beta + 2\alpha\cos\omega$ 變號的頻率 $\omega_z$ 處,$\theta(\omega)$、$\tau_p$ 與 $\tau_g$ 會發生什麼事?(c) 一個 AM 訊號 $A[n]\cos(\omega_c n)$($A[n]$ 為窄頻)通過一個系統,該系統有 $\theta(\omega_c) = -1.2$ rad、$\frac{d\theta}{d\omega}|_{\omega_c} = -8$、$\omega_c = 0.3$。載波與包絡各被延遲幾個樣本?問題:対称な 3 点フィルタ $h[n] = \{\alpha, \beta, \alpha\}$ は $H(e^{j\omega}) = e^{-j\omega}(\beta + 2\alpha\cos\omega)$ をもつ。(a) $\beta + 2\alpha\cos\omega > 0$ であるかぎり $\tau_p(\omega) = \tau_g(\omega) = 1$ 標本となることを示し、なぜこのフィルタが位相歪みを生じないのかを説明せよ。(b) $\beta + 2\alpha\cos\omega$ が符号を変える周波数 $\omega_z$ で、$\theta(\omega)$、$\tau_p$、$\tau_g$ はどうなるか。(c) 狭帯域の $A[n]$ をもつ AM 信号 $A[n]\cos(\omega_c n)$ が、$\theta(\omega_c) = -1.2$ rad、$\frac{d\theta}{d\omega}|_{\omega_c} = -8$、$\omega_c = 0.3$ のシステムを通る。搬送波と包絡線はそれぞれ何標本遅れるか。
(a) Where the real factor is positive it contributes no phase, so $\theta(\omega) = -\omega$:(a) 在實數因子為正的地方,它不貢獻任何相位,所以 $\theta(\omega) = -\omega$:(a) 実数因子が正のところでは位相をまったく寄与しないので、$\theta(\omega) = -\omega$ である。
$$\tau_p = -\frac{-\omega}{\omega} = 1, \qquad \tau_g = -\frac{d(-\omega)}{d\omega} = 1$$Both delays are constant and equal: every sinusoid and every envelope is delayed exactly 1 sample, so the output waveform shape is preserved (only amplitude-shaped by $|\beta + 2\alpha\cos\omega|$). Constant $\tau_g$ ⇔ linear phase ⇔ no phase distortion.兩個延遲都是常數而且相等:每個弦波、每個包絡都剛好延遲 1 個樣本,因此輸出波形的形狀被保留下來(只是振幅被 $|\beta + 2\alpha\cos\omega|$ 塑形)。$\tau_g$ 為常數 ⇔ 線性相位 ⇔ 沒有相位失真。二つの遅延はともに一定で、しかも等しい。どの正弦波もどの包絡線もちょうど 1 標本だけ遅れるので、出力波形の形は保たれる(振幅が $|\beta + 2\alpha\cos\omega|$ で整形されるだけである)。$\tau_g$ が一定 ⇔ 線形位相 ⇔ 位相歪みなし。
(b) At a sign change the phase jumps by $\pi$ (the magnitude $|\beta + 2\alpha\cos\omega|$ has a corner touching 0). $\tau_p$ becomes $\frac{\omega + \pi \cdot(\text{accumulated jumps})}{\omega}$-like — ill-defined/discontinuous at $\omega_z$ — and $\tau_g$, the derivative of a step, has an impulse there: group delay is simply not meaningful at frequencies where the response is zero (nothing gets through to be delayed).(b) 變號處相位跳 $\pi$(振幅 $|\beta + 2\alpha\cos\omega|$ 在那裡有一個觸及 0 的尖角)。$\tau_p$ 變成 $\frac{\omega + \pi \cdot(\text{accumulated jumps})}{\omega}$ 那樣的東西——在 $\omega_z$ 處定義不清、不連續——而 $\tau_g$ 作為階梯函數的導數,在那裡會出現一個脈衝:在響應為零的頻率上,群延遲根本沒有意義(根本沒有東西通過,自然談不上延遲)。(b) 符号が変わるところで位相は $\pi$ だけ跳ぶ(振幅 $|\beta + 2\alpha\cos\omega|$ はそこで 0 に触れる角をもつ)。$\tau_p$ は $\frac{\omega + \pi \cdot(\text{accumulated jumps})}{\omega}$ のような形になり、$\omega_z$ では定義が定まらず不連続になる。また $\tau_g$ は階段関数の微分として、そこにインパルスをもつ。応答がゼロになる周波数では、群遅延はそもそも意味をもたない(遅れるべきものが何も通っていないからである)。
(c) Carrier delay $\tau_p(\omega_c) = -\theta(\omega_c)/\omega_c = 1.2/0.3 = 4$ samples. Envelope delay $\tau_g(\omega_c) = -d\theta/d\omega = 8$ samples. The carrier slips 4 samples while the information envelope arrives 8 samples late.(c) 載波延遲 $\tau_p(\omega_c) = -\theta(\omega_c)/\omega_c = 1.2/0.3 = 4$ 個樣本。包絡延遲 $\tau_g(\omega_c) = -d\theta/d\omega = 8$ 個樣本。載波滑移 4 個樣本,而承載資訊的包絡則晚到 8 個樣本。(c) 搬送波の遅れは $\tau_p(\omega_c) = -\theta(\omega_c)/\omega_c = 1.2/0.3 = 4$ 標本。包絡線の遅れは $\tau_g(\omega_c) = -d\theta/d\omega = 8$ 標本である。搬送波は 4 標本ずれ、情報を担う包絡線は 8 標本遅れて到着する。
Problem: For the lecture system題目:對投影片中的系統問題:講義のシステム $$H(z) = \frac{z^{-1} - 1.2z^{-2} + z^{-3}}{1 - 1.3z^{-1} + 1.04z^{-2} - 0.222z^{-3}},$$ write MATLAB code to (a) plot the log-magnitude, principal-value phase, and group delay over $0 \le \omega \le \pi$, and (b) draw the pole-zero plot. (c) From the plots, explain the deep magnitude notch and the adjacent peak near $\omega \approx 0.3\pi$, and the sign of the group delay around the notch.請寫出 MATLAB 程式碼:(a) 在 $0 \le \omega \le \pi$ 上畫出對數振幅、主值相位與群延遲;(b) 畫出極零圖。(c) 由這些圖說明 $\omega \approx 0.3\pi$ 附近那個很深的振幅凹口與緊鄰的尖峰,以及凹口周圍群延遲的正負號。について、MATLAB のコードを書け。(a) $0 \le \omega \le \pi$ で対数振幅・主値位相・群遅延を描くこと。(b) 極零図を描くこと。(c) 得られた図から、$\omega \approx 0.3\pi$ 付近の深い振幅のノッチと隣り合うピーク、そしてノッチ周辺での群遅延の符号を説明せよ。
b = [0 1 -1.2 1]; % numerator z^-1 - 1.2 z^-2 + z^-3
a = [1 -1.3 1.04 -0.222]; % denominator
[H, w] = freqz(b, a, 1024); % frequency response on [0, pi)
[gd, wg] = grpdelay(b, a, 1024);
figure;
subplot(3,1,1); plot(w/pi, 20*log10(abs(H)));
ylabel('dB'); title('Log magnitude'); grid on;
subplot(3,1,2); plot(w/pi, angle(H));
ylabel('radians'); title('Phase (principal value)'); grid on;
subplot(3,1,3); plot(wg/pi, gd);
ylabel('samples'); xlabel('\omega/\pi'); title('Group delay'); grid on;
figure; zplane(b, a); % poles x, zeros o
% roots(b(2:end)) -> zeros 0.6 +/- 0.8i (|z| = 1, angle 0.2952*pi)
% roots(a) -> poles 0.5 +/- 0.7i, 0.3 (|z| = 0.86, angle 0.3026*pi)
Regenerated result (numpy/scipy equivalent of freqz/grpdelay):重新產生的結果(freqz/grpdelay 的 numpy/scipy 對應版本):再生成した結果(freqz/grpdelay の numpy/scipy 版)。
(c) The zeros $0.6 \pm j0.8$ sit exactly on the unit circle at angle $0.295\pi$ → $|H| = 0$ there: $20\log_{10}|H| \to -\infty$, the deep notch. The complex poles $0.5 \pm j0.7$ sit at radius $0.86$ at almost the same angle ($0.303\pi$) → a short pole vector boosts the magnitude immediately around the notch, producing the adjacent peak. Around the on-circle zero the phase jumps by $\pi$ (sign change of the zero factor), so the group delay shows a large negative spike right at the notch — group delay is not meaningful at a frequency where nothing is transmitted; just outside the notch the pole dominates and $\tau_g$ is large and positive.(c) 零點 $0.6 \pm j0.8$ 恰好落在單位圓上、角度為 $0.295\pi$ → 那裡 $|H| = 0$:$20\log_{10}|H| \to -\infty$,形成那個很深的凹口。複極點 $0.5 \pm j0.7$ 位於半徑 $0.86$、角度幾乎相同($0.303\pi$)→ 很短的極點向量把凹口緊鄰處的振幅推高,於是出現相鄰的尖峰。在圓上零點附近,相位跳了 $\pi$(零點因子變號),所以群延遲在凹口正上方出現一根很大的負向尖刺——在什麼都傳不過去的頻率上,群延遲本就沒有意義;而一離開凹口,極點便佔上風,$\tau_g$ 又大又正。(c) 零点 $0.6 \pm j0.8$ はちょうど単位円上、角度 $0.295\pi$ にある → そこで $|H| = 0$ となり、$20\log_{10}|H| \to -\infty$、これが深いノッチである。複素極 $0.5 \pm j0.7$ は半径 $0.86$、ほぼ同じ角度($0.303\pi$)にある → 短い極ベクトルがノッチのすぐ周りで振幅を押し上げ、隣り合うピークを作る。円上の零点のまわりでは位相が $\pi$ だけ跳ぶので(零点因子の符号変化による)、群遅延はノッチのちょうどその場所で大きな負のスパイクを見せる。何も伝わらない周波数で群遅延は意味をもたないのである。ノッチを少し外れれば極が支配的になり、$\tau_g$ は大きな正の値になる。