109 年 國立臺灣聯合大學系統(清華、政治、陽明交通、中央四校聯招)研究所電機類《訊號與系統》

📄 試題原卷 免費註冊後即可對照原始考卷 PDF免費註冊

第 一、(10%) 題10 分

The output y[n]y[n] of the three-point moving average system is related to the input x[n]x[n] according to the
formula
y[n]=13∑k=02x[n−k]y[n] = \frac{1}{3} \sum_{k=0}^{2} x[n-k]
(一) (5%) Please determine the impulse response of this system.
(二) (5%) Determine the output of the system when the input is the rectangular function defined as x[n]=u[n]−u[n−5]x[n] = u[n] - u[n-5].

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗離散時間系統的系統響應,包含求取系統的脈衝響應與特定輸入下的輸出。

(一) 脈衝響應 h[n]h[n] 的定義是當輸入為單位脈衝 δ[n]\delta[n] 時的輸出。
由題目給定的差分方程式 y[n]=13∑k=02x[n−k]y[n] = \frac{1}{3} \sum_{k=0}^{2} x[n-k],我們可以寫成
y[n]=13(x[n]+x[n−1]+x[n−2])y[n] = \frac{1}{3} (x[n] + x[n-1] + x[n-2])。
若輸入為 x[n]=δ[n]x[n] = \delta[n],則輸出為
y[n]=13(δ[n]+δ[n−1]+δ[n−2])y[n] = \frac{1}{3} (\delta[n] + \delta[n-1] + \delta[n-2])。
因此,此系統的脈衝響應為 h[n]=13(δ[n]+δ[n−1]+δ[n−2])h[n] = \frac{1}{3} (\delta[n] + \delta[n-1] + \delta[n-2])。

(二) 當輸入為 x[n]=u[n]−u[n−5]x[n] = u[n] - u[n-5] 時,輸出的計算可以透過卷積(convolution)進行,即 y[n]=x[n]∗h[n]y[n] = x[n] * h[n]。
x[n]x[n] 的定義為:
x[n]=1x[n] = 1 for 0≤n≤40 \le n \le 4
x[n]=0x[n] = 0 otherwise.

卷積的定義為 y[n]=∑k=−∞∞x[k]h[n−k]y[n] = \sum_{k=-\infty}^{\infty} x[k] h[n-k]。
將 h[n]=13(δ[n]+δ[n−1]+δ[n−2])h[n] = \frac{1}{3} (\delta[n] + \delta[n-1] + \delta[n-2]) 代入,
y[n]=∑k=−∞∞x[k]13(δ[n−k]+δ[n−k−1]+δ[n−k−2])y[n] = \sum_{k=-\infty}^{\infty} x[k] \frac{1}{3} (\delta[n-k] + \delta[n-k-1] + \delta[n-k-2])
y[n]=13∑k=−∞∞x[k](δ[n−k]+δ[n−(k+1)]+δ[n−(k+2)])y[n] = \frac{1}{3} \sum_{k=-\infty}^{\infty} x[k] (\delta[n-k] + \delta[n-(k+1)] + \delta[n-(k+2)])
y[n]=13(∑k=−∞∞x[k]δ[n−k]+∑k=−∞∞x[k]δ[n−(k+1)]+∑k=−∞∞x[k]δ[n−(k+2)])y[n] = \frac{1}{3} \left( \sum_{k=-\infty}^{\infty} x[k]\delta[n-k] + \sum_{k=-\infty}^{\infty} x[k]\delta[n-(k+1)] + \sum_{k=-\infty}^{\infty} x[k]\delta[n-(k+2)] \right)
y[n]=13(x[n]+x[n−1]+x[n−2])y[n] = \frac{1}{3} (x[n] + x[n-1] + x[n-2])。

現在我們將 x[n]x[n] 的值代入:
x[n]=1x[n] = 1 for 0≤n≤40 \le n \le 4
x[n−1]=1x[n-1] = 1 for 0≤n−1≤4  ⟹  1≤n≤50 \le n-1 \le 4 \implies 1 \le n \le 5

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 一、(15%) 題15 分

Suppose that an LTI system has the following output y(t)y(t) when the input is the unit step x(t)=u(t)x(t) = u(t):
y(t)=e−tu(t)+u(−1−t)y(t) = e^{-t}u(t) + u(-1-t)
Determine the response of the system to the input x(t)x(t) shown in Figure 1.

🖼️【此處有附圖,請對照原卷】
Figure 1: x(t)x(t) waveform. The waveform is a rectangular pulse that starts at t=0t=0, ends at t=2t=2, and has an amplitude of 1.

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 1 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗 LTI 系統的特性,特別是利用已知輸入與輸出關係來求取系統對其他輸入的響應。

首先,我們需要找出系統的脈衝響應 h(t)h(t)。
已知當輸入為單位階躍信號 x(t)=u(t)x(t) = u(t) 時,輸出為 y(t)=e−tu(t)+u(−1−t)y(t) = e^{-t}u(t) + u(-1-t)。
對於 LTI 系統,輸入 x(t)x(t) 與輸出 y(t)y(t) 的關係為 y(t)=x(t)∗h(t)y(t) = x(t) * h(t)。
若 x(t)=u(t)x(t) = u(t),則 y(t)=u(t)∗h(t)y(t) = u(t) * h(t)。
我們知道 u(t)∗h(t)=∫−∞th(τ)dτu(t) * h(t) = \int_{-\infty}^{t} h(\tau) d\tau。
所以,我們有 ∫−∞th(τ)dτ=e−tu(t)+u(−1−t)\int_{-\infty}^{t} h(\tau) d\tau = e^{-t}u(t) + u(-1-t)。

為了求得 h(t)h(t),我們對等號兩邊進行微分。
ddt(∫−∞th(τ)dτ)=h(t)\frac{d}{dt} \left( \int_{-\infty}^{t} h(\tau) d\tau \right) = h(t)。
對右式微分:
ddt(e−tu(t))=−e−tu(t)+e−tδ(t)\frac{d}{dt} (e^{-t}u(t)) = -e^{-t}u(t) + e^{-t}\delta(t)。
ddt(u(−1−t))\frac{d}{dt} (u(-1-t))。令 τ=−1−t\tau = -1-t,則 dτdt=−1\frac{d\tau}{dt} = -1。
ddt(u(−1−t))=du(τ)dτdτdt=δ(τ)×(−1)=−δ(−1−t)\frac{d}{dt} (u(-1-t)) = \frac{du(\tau)}{d\tau} \frac{d\tau}{dt} = \delta(\tau) \times (-1) = -\delta(-1-t)。
由於 δ(−t)=δ(t)\delta(-t) = \delta(t),所以 −δ(−1−t)=−δ(t+1)-\delta(-1-t) = -\delta(t+1)。
因此,
h(t)=−e−tu(t)+e−tδ(t)−δ(t+1)h(t) = -e^{-t}u(t) + e^{-t}\delta(t) - \delta(t+1)。

然而,單位脈衝 δ(t)\delta(t) 的定義是 e−tδ(t)=δ(t)e^{-t}\delta(t) = \delta(t) (因為 δ(t)\delta(t) 只在 t=0t=0 時有值,且 e0=1e^0=1)。
所以,h(t)=−e−tu(t)+δ(t)−δ(t+1)h(t) = -e^{-t}u(t) + \delta(t) - \delta(t+1)。

現在我們需要求取當輸入為圖一所示的 x(t)x(t) 時的輸出。
圖一的 x(t)x(t) 可以表示為 x(t)=u(t)−u(t−2)x(t) = u(t) - u(t-2)。
輸出 ynew(t)y_{new}(t) 為 ynew(t)=x(t)∗h(t)y_{new}(t) = x(t) * h(t)。
ynew(t)=(u(t)−u(t−2))∗h(t)y_{new}(t) = (u(t) - u(t-2)) * h(t)
ynew(t)=(u(t)∗h(t))−(u(t−2)∗h(t))y_{new}(t) = (u(t) * h(t)) - (u(t-2) * h(t))。

我們已經知道 u(t)∗h(t)=y(t)=e−tu(t)+u(−1−t)u(t) * h(t) = y(t) = e^{-t}u(t) + u(-1-t)。
現在計算 u(t−2)∗h(t)u(t-2) * h(t)。
利用時移性質,若 f(t)∗h(t)=g(t)f(t) * h(t) = g(t),則 f(t−t0)∗h(t)=g(t−t0)f(t-t_0) * h(t) = g(t-t_0)。
這裡 f(t)=u(t)f(t) = u(t),所以 u(t−2)∗h(t)=y(t−2)u(t-2) * h(t) = y(t-2)。
y(t−2)=e−(t−2)u(t−2)+u(−1−(t−2))y(t-2) = e^{-(t-2)}u(t-2) + u(-1-(t-2))

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 一、(15%) 題15 分

For a discrete time LTI system:
y[n]−13y[n−1]+14y[n−2]=x[n−1]+12x[n−2]y[n] - \frac{1}{3}y[n-1] + \frac{1}{4}y[n-2] = x[n-1] + \frac{1}{2}x[n-2]
(1) (4%) Determine the frequency response of system.
(2) (4%) Determine the time domain impulse response of system.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查離散時間 LTI 系統的:

  1. 頻率響應:令輸入為複指數訊號 x[n]=ejωnx[n]=e^{j\omega n},輸出為 y[n]=H(ejω)ejωny[n]=H(e^{j\omega})e^{j\omega n},即可求得頻率響應。
  2. 單位脈衝響應:令輸入 x[n]=δ[n]x[n]=\delta[n],求出輸出 y[n]=h[n]y[n]=h[n]。
  3. 差分方程與系統函數:利用延遲運算子 z−1z^{-1} 將差分方程轉為 H(z)H(z),再反轉換得到 h[n]h[n]。

以下採用系統初始靜止且為因果系統的標準假設。


解題方法一:求頻率響應

題目給定差分方程:

y[n]−13y[n−1]+14y[n−2]=x[n−1]+12x[n−2]y[n]-\frac{1}{3}y[n-1]+\frac{1}{4}y[n-2] = x[n-1]+\frac{1}{2}x[n-2]

對 LTI 系統令:

x[n]=ejωnx[n]=e^{j\omega n}

則輸出可寫為:

y[n]=H(ejω)ejωny[n]=H(e^{j\omega})e^{j\omega n}

利用延遲性質:

x[n−k]=ejω(n−k)=e−jωkejωnx[n-k]=e^{j\omega(n-k)} =e^{-j\omega k}e^{j\omega n} y[n−k]=H(ejω)e−jωkejωny[n-k] = H(e^{j\omega})e^{-j\omega k}e^{j\omega n}

代入原差分方程:

H(ejω)ejωn−13H(ejω)e−jωejωn+14H(ejω)e−j2ωejωnH(e^{j\omega})e^{j\omega n} -\frac{1}{3}H(e^{j\omega})e^{-j\omega}e^{j\omega n} +\frac{1}{4}H(e^{j\omega})e^{-j2\omega}e^{j\omega n} =e−jωejωn+12e−j2ωejωn= e^{-j\omega}e^{j\omega n} +\frac{1}{2}e^{-j2\omega}e^{j\omega n}

消去 ejωne^{j\omega n},得:

H(ejω)(1−13e−jω+14e−j2ω)=e−jω+12e−j2ωH(e^{j\omega}) \left( 1-\frac{1}{3}e^{-j\omega} +\frac{1}{4}e^{-j2\omega} \right) = e^{-j\omega} +\frac{1}{2}e^{-j2\omega}

因此頻率響應為:

H(ejω)=e−jω+12e−j2ω1−13e−jω+14e−j2ω\boxed{ H(e^{j\omega}) = \frac{e^{-j\omega}+\frac{1}{2}e^{-j2\omega}} {1-\frac{1}{3}e^{-j\omega}+\frac{1}{4}e^{-j2\omega}} }

解題方法二:求系統函數

在初始靜止條件下,對差分方程取 ZZ 轉換:

Y(z)−13z−1Y(z)+14z−2Y(z)=z−1X(z)+12z−2X(z)Y(z)-\frac{1}{3}z^{-1}Y(z)+\frac{1}{4}z^{-2}Y(z) = z^{-1}X(z)+\frac{1}{2}z^{-2}X(z)

整理得:

Y(z)(1−13z−1+14z−2)=X(z)(z−1+12z−2)Y(z) \left( 1-\frac{1}{3}z^{-1} +\frac{1}{4}z^{-2} \right) = X(z) \left( z^{-1}+\frac{1}{2}z^{-2} \right)

所以系統函數為:

H(z)=Y(z)X(z)=z−1+12z−21−13z−1+14z−2H(z)=\frac{Y(z)}{X(z)} = \frac{z^{-1}+\frac{1}{2}z^{-2}} {1-\frac{1}{3}z^{-1}+\frac{1}{4}z^{-2}}

將 z=ejωz=e^{j\omega} 代入,即可得到上一節的頻率響應。


求時間域單位脈衝響應

令:

x[n]=δ[n],y[n]=h[n]x[n]=\delta[n],\qquad y[n]=h[n]

原差分方程變為:

h[n]−13h[n−1]+14h[n−2]=δ[n−1]+12δ[n−2]h[n]-\frac{1}{3}h[n-1]+\frac{1}{4}h[n-2] = \delta[n-1]+\frac{1}{2}\delta[n-2]

因系統為因果且初始靜止:

h[n]=0,n<0h[n]=0,\qquad n<0

逐點計算初始值。

n=0n=0

h[0]−13h[−1]+14h[−2]=0h[0]-\frac{1}{3}h[-1]+\frac{1}{4}h[-2]=0

因此:

h[0]=0h[0]=0

n=1n=1

h[1]−13h[0]+14h[−1]=1h[1]-\frac{1}{3}h[0]+\frac{1}{4}h[-1]=1
🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 三、(15%) 題15 分

(一) (8%) Determine the Fourier series coefficients of x(t)=(sin⁡(t))(cos⁡(3t))x(t) = (\sin(t))(\cos(3t)).
(二) (7%) Given the discrete-time Fourier series coefficients X[k]=cos⁡(3π16k)+jsin⁡(5π16k)X[k] = \cos\left(\frac{3\pi}{16}k\right) + j\sin\left(\frac{5\pi}{16}k\right), determine the discrete-time periodic signal x[n]x[n].

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗連續時間與離散時間的傅立葉級數。

(一) 求取 x(t)=(sin⁡(t))(cos⁡(3t))x(t) = (\sin(t))(\cos(3t)) 的傅立葉級數係數。
我們可以使用三角恆等式來簡化 x(t)x(t)。
sin⁡(A)cos⁡(B)=12[sin⁡(A+B)+sin⁡(A−B)]\sin(A)\cos(B) = \frac{1}{2}[\sin(A+B) + \sin(A-B)]。
令 A=tA=t, B=3tB=3t。
x(t)=12[sin⁡(t+3t)+sin⁡(t−3t)]x(t) = \frac{1}{2}[\sin(t+3t) + \sin(t-3t)]
x(t)=12[sin⁡(4t)+sin⁡(−2t)]x(t) = \frac{1}{2}[\sin(4t) + \sin(-2t)]
x(t)=12[sin⁡(4t)−sin⁡(2t)]x(t) = \frac{1}{2}[\sin(4t) - \sin(2t)]。

傅立葉級數的定義為 x(t)=∑k=−∞∞ckejkω0tx(t) = \sum_{k=-\infty}^{\infty} c_k e^{jk\omega_0 t},其中 ω0\omega_0 是基頻。
對於 x(t)=12sin⁡(4t)−12sin⁡(2t)x(t) = \frac{1}{2}\sin(4t) - \frac{1}{2}\sin(2t),我們可以將 sin⁡\sin 函數用複指數形式表示:
sin⁡(ωt)=ejωt−e−jωt2j\sin(\omega t) = \frac{e^{j\omega t} - e^{-j\omega t}}{2j}。
所以,
x(t)=12(ej4t−e−j4t2j)−12(ej2t−e−j2t2j)x(t) = \frac{1}{2} \left( \frac{e^{j4t} - e^{-j4t}}{2j} \right) - \frac{1}{2} \left( \frac{e^{j2t} - e^{-j2t}}{2j} \right)
x(t)=14j(ej4t−e−j4t)−14j(ej2t−e−j2t)x(t) = \frac{1}{4j} (e^{j4t} - e^{-j4t}) - \frac{1}{4j} (e^{j2t} - e^{-j2t})
x(t)=14jej4t−14je−j4t−14jej2t+14je−j2tx(t) = \frac{1}{4j} e^{j4t} - \frac{1}{4j} e^{-j4t} - \frac{1}{4j} e^{j2t} + \frac{1}{4j} e^{-j2t}。

比較 x(t)=∑k=−∞∞ckejkω0tx(t) = \sum_{k=-\infty}^{\infty} c_k e^{jk\omega_0 t},我們可以看到主要的頻率是 22 和 44。
假設基頻 ω0\omega_0 是 22 (或者 22 的因數)。如果 ω0=2\omega_0 = 2,則 kk 可以是 11 或 22。
如果 ω0=2\omega_0 = 2,則 x(t)=14jej2⋅2t−14jej(−2)⋅2t−14jej1⋅2t+14jej(−1)⋅2tx(t) = \frac{1}{4j} e^{j2 \cdot 2 t} - \frac{1}{4j} e^{j(-2) \cdot 2 t} - \frac{1}{4j} e^{j1 \cdot 2 t} + \frac{1}{4j} e^{j(-1) \cdot 2 t}。
那麼,傅立葉級數係數 ckc_k 為:
c2=14jc_2 = \frac{1}{4j}
c−2=−14jc_{-2} = -\frac{1}{4j}
c1=−14jc_1 = -\frac{1}{4j}
c−1=14jc_{-1} = \frac{1}{4j}
所有其他的 ckc_k 係數都為 0。

如果我們取 ω0=1\omega_0 = 1,那麼 jω0tj\omega_0 t 的形式是 jt,j2t,j3t,…jt, j2t, j3t, \dots。
x(t)=14jej4t−14je−j4t−14jej2t+14je−j2tx(t) = \frac{1}{4j} e^{j4t} - \frac{1}{4j} e^{-j4t} - \frac{1}{4j} e^{j2t} + \frac{1}{4j} e^{-j2t}
x(t)=14jej4ω0t−14jej(−4)ω0t−14jej2ω0t+14jej(−2)ω0tx(t) = \frac{1}{4j} e^{j4\omega_0 t} - \frac{1}{4j} e^{j(-4)\omega_0 t} - \frac{1}{4j} e^{j2\omega_0 t} + \frac{1}{4j} e^{j(-2)\omega_0 t}
所以,
c4=14j=−j4c_4 = \frac{1}{4j} = -\frac{j}{4}
c−4=−14j=j4c_{-4} = -\frac{1}{4j} = \frac{j}{4}
c2=−14j=j4c_2 = -\frac{1}{4j} = \frac{j}{4}
c−2=14j=−j4c_{-2} = \frac{1}{4j} = -\frac{j}{4}
所有其他的 ckc_k 係數都為 0。

通常,傅立葉級數的週期 TT 決定了基頻 ω0=2πT\omega_0 = \frac{2\pi}{T}。
這裡的函數是 sin⁡(2t)\sin(2t) 和 sin⁡(4t)\sin(4t)。
sin⁡(2t)\sin(2t) 的週期是 T1=2π2=πT_1 = \frac{2\pi}{2} = \pi。
sin⁡(4t)\sin(4t) 的週期是 T2=2π4=π2T_2 = \frac{2\pi}{4} = \frac{\pi}{2}。
這兩個函數的最小公倍週期是 T=πT = \pi。
所以,基頻 ω0=2ππ=2\omega_0 = \frac{2\pi}{\pi} = 2。
因此,我們應該使用 ω0=2\omega_0 = 2 來表示。
x(t)=c1ej1⋅2t+c−1ej(−1)⋅2t+c2ej2⋅2t+c−2ej(−2)⋅2tx(t) = c_1 e^{j1 \cdot 2t} + c_{-1} e^{j(-1) \cdot 2t} + c_2 e^{j2 \cdot 2t} + c_{-2} e^{j(-2) \cdot 2t}。
x(t)=−14jej2t+14je−j2t+14jej4t−14je−j4tx(t) = -\frac{1}{4j} e^{j2t} + \frac{1}{4j} e^{-j2t} + \frac{1}{4j} e^{j4t} - \frac{1}{4j} e^{-j4t}。
所以,
c1=−14j=j4c_1 = -\frac{1}{4j} = \frac{j}{4}
c−1=14j=−j4c_{-1} = \frac{1}{4j} = -\frac{j}{4}
c2=14j=−j4c_2 = \frac{1}{4j} = -\frac{j}{4}
c−2=−14j=j4c_{-2} = -\frac{1}{4j} = \frac{j}{4}
所有其他的 ckc_k 係數都為 0。

(二) 給定離散時間傅立葉級數係數 X[k]=cos⁡(3π16k)+jsin⁡(5π16k)X[k] = \cos\left(\frac{3\pi}{16}k\right) + j\sin\left(\frac{5\pi}{16}k\right),求取離散時間

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 三、(15%) 題15 分

(一) (5%) Prove the Parseval relationship:
∫−∞∞∣x(t)∣2dt=12π∫−∞∞∣X(jω)∣2dω\int_{-\infty}^{\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega)|^2 d\omega
(二) (5%) Determine the Fourier transform of x(t)=e−2tu(t)x(t) = e^{-2t}u(t).
(三) (5%) Let X(ejΩ)X(e^{j\Omega}) be the discrete-time Fourier transform of a real signal x[n]x[n]. Given
Y(ejΩ)=X(ejΩ)ej2Ω+X(e−jΩ)e−j2ΩY(e^{j\Omega}) = X(e^{j\Omega})e^{j2\Omega} + X(e^{-j\Omega})e^{-j2\Omega}
determine the discrete-time signal y[n]y[n], expressed by x[n]x[n].

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗連續時間傅立葉轉換的帕斯伐定理、基本函數的傅立葉轉換,以及離散時間傅立葉轉換的性質。

(一) 證明帕斯伐定理 (Parseval's Theorem)。
帕斯伐定理是傅立葉轉換的能量守恆定理。
連續時間傅立葉轉換的定義是 X(jω)=∫−∞∞x(t)e−jωtdtX(j\omega) = \int_{-\infty}^{\infty} x(t) e^{-j\omega t} dt。
其逆轉換是 x(t)=12π∫−∞∞X(jω)ejωtdωx(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} X(j\omega) e^{j\omega t} d\omega。
我們需要證明 ∫−∞∞∣x(t)∣2dt=12π∫−∞∞∣X(jω)∣2dω\int_{-\infty}^{\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega)|^2 d\omega。

考慮左式:∫−∞∞∣x(t)∣2dt=∫−∞∞x(t)x∗(t)dt\int_{-\infty}^{\infty} |x(t)|^2 dt = \int_{-\infty}^{\infty} x(t) x^*(t) dt。
將 x∗(t)x^*(t) 用逆轉換表示:
x∗(t)=(12π∫−∞∞X(jω)ejωtdω)∗x^*(t) = \left( \frac{1}{2\pi} \int_{-\infty}^{\infty} X(j\omega) e^{j\omega t} d\omega \right)^*
x∗(t)=12π∫−∞∞X∗(jω)e−jωtdωx^*(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega) e^{-j\omega t} d\omega。

代入左式:
∫−∞∞x(t)(12π∫−∞∞X∗(jω′)e−jω′tdω′)dt\int_{-\infty}^{\infty} x(t) \left( \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega') e^{-j\omega' t} d\omega' \right) dt
=12π∫−∞∞X∗(jω′)(∫−∞∞x(t)e−jω′tdt)dω′= \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega') \left( \int_{-\infty}^{\infty} x(t) e^{-j\omega' t} dt \right) d\omega'。
括號內的積分就是 X(jω′)X(j\omega') 的定義。
=12π∫−∞∞X∗(jω′)X(jω′)dω′= \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega') X(j\omega') d\omega'
=12π∫−∞∞∣X(jω′)∣2dω′= \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega')|^2 d\omega'。
將積分變數改為 ω\omega,即得證:
∫−∞∞∣x(t)∣2dt=12π∫−∞∞∣X(jω)∣2dω\int_{-\infty}^{\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega)|^2 d\omega。

(二) 求取 x(t)=e−2tu(t)x(t) = e^{-2t}u(t) 的傅立葉轉換。
傅立葉轉換的定義是 X(jω)=∫−∞∞x(t)e−jωtdtX(j\omega) = \int_{-\infty}^{\infty} x(t) e^{-j\omega t} dt。
X(jω)=∫−∞∞e−2tu(t)e−jωtdtX(j\omega) = \int_{-\infty}^{\infty} e^{-2t}u(t) e^{-j\omega t} dt。
由於 u(t)u(t) 的存在,積分範圍變為從 00 到 ∞\infty。
X(jω)=∫0∞e−2te−jωtdtX(j\omega) = \int_{0}^{\infty} e^{-2t} e^{-j\omega t} dt
X(jω)=∫0∞e−(2+jω)tdtX(j\omega) = \int_{0}^{\infty} e^{-(2+j\omega)t} dt。
此積分收斂的條件是 Re(2+jω)>0Re(2+j\omega) > 0,即 2>02 > 0,這總是成立的。
X(jω)=[e−(2+jω)t−(2+jω)]0∞X(j\omega) = \left[ \frac{e^{-(2+j\omega)t}}{-(2+j\omega)} \right]_{0}^{\infty}

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 四、(15%) 題15 分

(一) (5%) Prove the Parseval relationship:
∫−∞∞∣x(t)∣2dt=12π∫−∞∞∣X(jω)∣2dω\int_{-\infty}^{\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega)|^2 d\omega
(二) (5%) Determine the Fourier transform of x(t)=e−2tu(t)x(t) = e^{-2t}u(t).
(三) (5%) Let X(ejΩ)X(e^{j\Omega}) be the discrete-time Fourier transform of a real signal x[n]x[n]. Given
Y(ejΩ)=X(ejΩ)ej2Ω+X(e−jΩ)e−j2ΩY(e^{j\Omega}) = X(e^{j\Omega})e^{j2\Omega} + X(e^{-j\Omega})e^{-j2\Omega}
determine the discrete-time signal y[n]y[n], expressed by x[n]x[n].

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗連續時間傅立葉轉換的帕斯伐定理、基本函數的傅立葉轉換,以及離散時間傅立葉轉換的性質。

(一) 證明帕斯伐定理 (Parseval's Theorem)。
帕斯伐定理是傅立葉轉換的能量守恆定理。
連續時間傅立葉轉換的定義是 X(jω)=∫−∞∞x(t)e−jωtdtX(j\omega) = \int_{-\infty}^{\infty} x(t) e^{-j\omega t} dt。
其逆轉換是 x(t)=12π∫−∞∞X(jω)ejωtdωx(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} X(j\omega) e^{j\omega t} d\omega。
我們需要證明 ∫−∞∞∣x(t)∣2dt=12π∫−∞∞∣X(jω)∣2dω\int_{-\infty}^{\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega)|^2 d\omega。

考慮左式:∫−∞∞∣x(t)∣2dt=∫−∞∞x(t)x∗(t)dt\int_{-\infty}^{\infty} |x(t)|^2 dt = \int_{-\infty}^{\infty} x(t) x^*(t) dt。
將 x∗(t)x^*(t) 用逆轉換表示:
x∗(t)=(12π∫−∞∞X(jω)ejωtdω)∗x^*(t) = \left( \frac{1}{2\pi} \int_{-\infty}^{\infty} X(j\omega) e^{j\omega t} d\omega \right)^*
x∗(t)=12π∫−∞∞X∗(jω)e−jωtdωx^*(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega) e^{-j\omega t} d\omega。

代入左式:
∫−∞∞x(t)(12π∫−∞∞X∗(jω′)e−jω′tdω′)dt\int_{-\infty}^{\infty} x(t) \left( \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega') e^{-j\omega' t} d\omega' \right) dt
=12π∫−∞∞X∗(jω′)(∫−∞∞x(t)e−jω′tdt)dω′= \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega') \left( \int_{-\infty}^{\infty} x(t) e^{-j\omega' t} dt \right) d\omega'。
括號內的積分就是 X(jω′)X(j\omega') 的定義。
=12π∫−∞∞X∗(jω′)X(jω′)dω′= \frac{1}{2\pi} \int_{-\infty}^{\infty} X^*(j\omega') X(j\omega') d\omega'
=12π∫−∞∞∣X(jω′)∣2dω′= \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega')|^2 d\omega'。
將積分變數改為 ω\omega,即得證:
∫−∞∞∣x(t)∣2dt=12π∫−∞∞∣X(jω)∣2dω\int_{-\infty}^{\infty} |x(t)|^2 dt = \frac{1}{2\pi} \int_{-\infty}^{\infty} |X(j\omega)|^2 d\omega。

(二) 求取 x(t)=e−2tu(t)x(t) = e^{-2t}u(t) 的傅立葉轉換。
傅立葉轉換的定義是 X(jω)=∫−∞∞x(t)e−jωtdtX(j\omega) = \int_{-\infty}^{\infty} x(t) e^{-j\omega t} dt。
X(jω)=∫−∞∞e−2tu(t)e−jωtdtX(j\omega) = \int_{-\infty}^{\infty} e^{-2t}u(t) e^{-j\omega t} dt。
由於 u(t)u(t) 的存在,積分範圍變為從 00 到 ∞\infty。
X(jω)=∫0∞e−2te−jωtdtX(j\omega) = \int_{0}^{\infty} e^{-2t} e^{-j\omega t} dt
X(jω)=∫0∞e−(2+jω)tdtX(j\omega) = \int_{0}^{\infty} e^{-(2+j\omega)t} dt。
此積分收斂的條件是 Re(2+jω)>0Re(2+j\omega) > 0,即 2>02 > 0,這總是成立的。
X(jω)=[e−(2+jω)t−(2+jω)]0∞X(j\omega) = \left[ \frac{e^{-(2+j\omega)t}}{-(2+j\omega)} \right]_{0}^{\infty}

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 五、(5%) 題5 分

Given two continuous-time periodic signals a(t)a(t) and b(t)b(t), is the following signal: a(t)eb(t)a(t)e^{b(t)} always a periodic signal? Justify your answer.

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

連續時間訊號 x(t)x(t) 為週期訊號的定義,是存在某個有限且正的 TT,使得對所有 tt 均滿足

x(t+T)=x(t).x(t+T)=x(t).

若 a(t)a(t) 的週期為 TaT_a,b(t)b(t) 的週期為 TbT_b,則:

  • eb(t)e^{b(t)} 與 b(t)b(t) 具有相同的週期;
  • a(t)eb(t)a(t)e^{b(t)} 至少需要存在一個共同週期,使兩個因子同時重複;
  • 兩個週期訊號不一定具有共同週期。

當兩個基本週期的比值為有理數時,通常可找到共同週期;若週期比值為無理數,則不一定存在共同週期。


解題方法

先取兩個週期比值為無理數的訊號:

a(t)=ecos⁡t,b(t)=cos⁡(2t).a(t)=e^{\cos t}, \qquad b(t)=\cos(\sqrt{2}t).

兩者皆為連續時間週期訊號,其中

Ta=2π,Tb=2π2=2π.T_a=2\pi, \qquad T_b=\frac{2\pi}{\sqrt{2}}=\sqrt{2}\pi.

其週期比值為

TaTb=2π2π=2,\frac{T_a}{T_b} = \frac{2\pi}{\sqrt{2}\pi} = \sqrt{2},

而 2\sqrt{2} 為無理數,因此兩者沒有有限的共同基本週期。

題目所給訊號為

x(t)=a(t)eb(t)=ecos⁡tecos⁡(2t)=ecos⁡t+cos⁡(2t).x(t)=a(t)e^{b(t)} =e^{\cos t}e^{\cos(\sqrt{2}t)} =e^{\cos t+\cos(\sqrt{2}t)}.

由於 x(t)>0x(t)>0,若 x(t)x(t) 是週期訊號,則其自然對數也必須是週期訊號:

ln⁡x(t)=cos⁡t+cos⁡(2t).\ln x(t)=\cos t+\cos(\sqrt{2}t).

訊號 cos⁡t\cos t 的頻率為 11,訊號 cos⁡(2t)\cos(\sqrt{2}t) 的頻率為 2\sqrt{2}。若它們的和具有週期 TT,則兩個頻率都必須是某個基本角頻率 ω0\omega_0 的整數倍,因此應存在整數 m,nm,n 使得

1=mω0,2=nω0.1=m\omega_0, \qquad \sqrt{2}=n\omega_0.

兩式相除可得

2=nm.\sqrt{2}=\frac{n}{m}.

右側為有理數,左側為無理數,矛盾。因此 cos⁡t+cos⁡(2t)\cos t+\cos(\sqrt{2}t) 不具週期性,進而

x(t)=a(t)eb(t)x(t)=a(t)e^{b(t)}

也不具週期性。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 六、(15%) 題15 分

Given a real signal x(t)x(t) which has a spectrum in Figure 2(a) with ω1>ω2−ω1\omega_1 > \omega_2 - \omega_1, consider the sampling system in Figure 2(b), where δ(t)\delta(t) denotes the continuous-time unit impulse function and H(jω)H(j\omega) indicates a filter to be determined by you. Find the maximum value of TT such that xs(t)=x(t)x_s(t) = x(t). Does your answer violate the sampling theorem? Note that answers without justification get no credit.

🖼️【此處有附圖,請對照原卷】
Figure 2(a): Spectrum of x(t)x(t). The spectrum is non-zero between −ω2-\omega_2 and −ω1-\omega_1, and between ω1\omega_1 and ω2\omega_2. The shape is a rectangle of height 1.
Figure 2(b): Sampling system. Input x(t)x(t) is multiplied by p(t)=∑k=−∞∞δ(t−kT)p(t) = \sum_{k=-\infty}^{\infty} \delta(t-kT), then filtered by H(jω)H(j\omega) to produce xs(t)x_s(t).

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 2 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  • 週期脈衝取樣造成頻譜以取樣角頻率 ωs\omega_s 為間隔週期性複製。
  • 帶通取樣定理(bandpass sampling theorem)。
  • 只要各頻譜複本不互相重疊,即可用濾波器 H(jω)H(j\omega) 選出原頻譜,並補償取樣造成的比例因子。

圖中原頻譜為

X(jω)≠0,ω1≤∣ω∣≤ω2X(j\omega)\neq 0,\qquad \omega_1\leq |\omega|\leq \omega_2

令頻帶寬度為

B=ω2−ω1.B=\omega_2-\omega_1.

取樣脈衝列為

p(t)=∑k=−∞∞δ(t−kT),p(t)=\sum_{k=-\infty}^{\infty}\delta(t-kT),

其取樣角頻率為

ωs=2πT.\omega_s=\frac{2\pi}{T}.

解題方法

時域相乘對應到頻域複製,因此

Xp(jω)=1T∑k=−∞∞X(j(ω−kωs)).X_p(j\omega) = \frac{1}{T} \sum_{k=-\infty}^{\infty} X\bigl(j(\omega-k\omega_s)\bigr).

每一個原頻帶寬度都是 BB,取樣後會以 ωs\omega_s 為間隔重複出現。為了讓濾波器能夠完整取出原本的 X(jω)X(j\omega),相鄰頻譜複本不可重疊。

帶通取樣條件可寫成

2ω2n≤ωs≤2ω1n−1,n=1,2,3,…\frac{2\omega_2}{n} \leq \omega_s \leq \frac{2\omega_1}{n-1}, \qquad n=1,2,3,\ldots

其中 n=1n=1 時右側上限視為無限大。

此條件有可行解的必要條件為

2ω2n≤2ω1n−1.\frac{2\omega_2}{n} \leq \frac{2\omega_1}{n-1}.

整理得

ω2(n−1)≤ω1n,\omega_2(n-1)\leq \omega_1n,

因此

n(ω2−ω1)≤ω1,n(\omega_2-\omega_1)\leq \omega_1,

也就是

nB≤ω1.nB\leq \omega_1.

為使 T=2π/ωsT=2\pi/\omega_s 最大,必須使可行的 nn 最大。令

nmax⁡=⌊ω1ω2−ω1⌋.n_{\max} = \left\lfloor\frac{\omega_1}{\omega_2-\omega_1}\right\rfloor.

此時可採用的最小取樣角頻率為

ωs,min⁡=2ω2nmax⁡.\omega_{s,\min} = \frac{2\omega_2}{n_{\max}}.

因此

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 七、(10%) 題10 分

Consider a continuous-time stable and causal system with impulse response h(t)h(t) and system function H(s)H(s). Suppose H(s)H(s) is rational, has a pole at s=−4+4js = -4 + 4j, and does not have a zero at the origin (at the origin, s=0+0js = 0 + 0j). The location of all the other poles and zeros is unknown. For each of the following statements, please determine it is "true" or "false" and justify your decision.
(一) (4%) F{h(t)e4t}F\{h(t)e^{4t}\} converges (F{}F\{\} represents Fourier transform)
(二) (3%) ∫−∞∞h(t)dt=0\int_{-\infty}^{\infty} h(t)dt = 0
(三) (3%) H(s)=H(−s)H(s) = H(-s)

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗對 LTI 系統的穩定性、因果性、系統函數、傅立葉轉換以及拉氏轉換的理解。

已知條件:

  1. 系統是連續時間、穩定且因果的。
  2. 系統函數 H(s)H(s) 是有理函數。
  3. H(s)H(s) 在 s=−4+4js = -4 + 4j 處有一個極點。
  4. H(s)H(s) 在 s=0s = 0 處沒有零點。
  5. 其他極點和零點的位置未知。

(一) F{h(t)e4t}F\{h(t)e^{4t}\} converges (F{}F\{\} 代表傅立葉轉換)。
傅立葉轉換 F{g(t)}F\{g(t)\} 收斂,若且唯若 g(t)g(t) 是絕對可積的,即 ∫−∞∞∣g(t)∣dt<∞\int_{-\infty}^{\infty} |g(t)| dt < \infty。
令 g(t)=h(t)e4tg(t) = h(t)e^{4t}。
我們知道系統是穩定的,這意味著 h(t)h(t) 的所有極點(在 ss-平面上)的實部都必須是負的。
我們也知道系統是因果的,這意味著 h(t)=0h(t)=0 對於 t<0t<0。
由於 H(s)H(s) 是有理函數且系統是因果穩定的,其拉氏轉換 H(s)H(s) 的收斂域 (ROC) 必定包含右半平面,並且包含虛軸。
即,對於所有 Re(s)>σmin\text{Re}(s) > \sigma_{min},其中 σmin<0\sigma_{min} < 0。
由於系統是穩定的,所有極點的實部都必須小於 0。
給定的極點是 s=−4+4js = -4 + 4j。它的實部是 −4-4,這是一個負數,符合穩定性要求。
h(t)h(t) 可以表示為一系列指數項的和(或包含振盪項),形式為 ckepktu(t)c_k e^{p_k t} u(t),其中 pkp_k 是 H(s)H(s) 的極點,且 Re(pk)<0\text{Re}(p_k) < 0。
那麼,g(t)=h(t)e4t=∑kckepktu(t)e4t=∑kcke(pk+4)tu(t)g(t) = h(t)e^{4t} = \sum_k c_k e^{p_k t} u(t) e^{4t} = \sum_k c_k e^{(p_k+4)t} u(t)。
為了使 g(t)g(t) 絕對可積,我們需要 Re(pk+4)<0\text{Re}(p_k+4) < 0 對於所有的 kk。
也就是說,Re(pk)<−4\text{Re}(p_k) < -4 對於所有的 kk。

然而,我們只知道其中一個極點是 p1=−4+4jp_1 = -4 + 4j,其 Re(p1)=−4\text{Re}(p_1) = -4。
那麼,對於這個極點對應的項,我們有 c1e(−4+4j)tu(t)e4t=c1e4jtu(t)c_1 e^{(-4+4j)t} u(t) e^{4t} = c_1 e^{4j t} u(t)。
這個函數是振盪的,並且是絕對不可積的。
例如,如果 h(t)h(t) 包含一個項 e(−4+4j)tu(t)e^{(-4+4j)t} u(t),那麼 h(t)e4th(t)e^{4t} 包含一個項 e4jtu(t)e^{4jt} u(t)。
∫−∞∞∣e4jtu(t)∣dt=∫0∞∣e4jt∣dt=∫0∞1dt=∞\int_{-\infty}^{\infty} |e^{4jt} u(t)| dt = \int_{0}^{\infty} |e^{4jt}| dt = \int_{0}^{\infty} 1 dt = \infty。
因此,F{h(t)e4t}F\{h(t)e^{4t}\} 不收斂。

要使 F{h(t)e4t}F\{h(t)e^{4t}\} 收斂,需要 h(t)e4th(t)e^{4t} 是絕對可積的。
這意味著對於 h(t)h(t) 的所有極點 pkp_k,必須滿足 Re(pk)<−4\text{Re}(p_k) < -4。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 八、(5%) 題5 分

Consider a discrete-time noncausal seven-point weighted moving average system as follows.
y[n]=∑k=−33(a3)∣k∣x[n−k]y[n] = \sum_{k=-3}^{3} \left(\frac{a}{3}\right)^{|k|} x[n-k]
where x[n]x[n] is the input, y[n]y[n] is the output, ∣k∣|k| represents the absolute value of the integer kk, and aa is real.
Find the system function of its causal version. Note that the causal seven-point weighted moving average
system should have the same magnitude response as the noncasual one has, and should have minimum
group delay when being compared with all the other causal versions.

登入後即可作答並保存紀錄。

這一題的完整詳解

此題考驗離散時間系統的系統函數、因果性、非因果性以及群延遲。

題目給定一個非因果的七點移動平均系統:
y[n]=∑k=−33(a3)∣k∣x[n−k]y[n] = \sum_{k=-3}^{3} \left(\frac{a}{3}\right)^{|k|} x[n-k]
其中 aa 是實數。
我們可以將這個方程式寫成:
y[n]=(a3)∣−3∣x[n−(−3)]+(a3)∣−2∣x[n−(−2)]+⋯+(a3)∣0∣x[n−0]+⋯+(a3)∣3∣x[n−3]y[n] = (\frac{a}{3})^{|-3|} x[n-(-3)] + (\frac{a}{3})^{|-2|} x[n-(-2)] + \dots + (\frac{a}{3})^{|0|} x[n-0] + \dots + (\frac{a}{3})^{|3|} x[n-3]
y[n]=(a3)3x[n+3]+(a3)2x[n+2]+(a3)1x[n+1]+(a3)0x[n]+(a3)1x[n−1]+(a3)2x[n−2]+(a3)3x[n−3]y[n] = (\frac{a}{3})^3 x[n+3] + (\frac{a}{3})^2 x[n+2] + (\frac{a}{3})^1 x[n+1] + (\frac{a}{3})^0 x[n] + (\frac{a}{3})^1 x[n-1] + (\frac{a}{3})^2 x[n-2] + (\frac{a}{3})^3 x[n-3]
y[n]=(a3)3x[n+3]+(a3)2x[n+2]+a3x[n+1]+x[n]+a3x[n−1]+(a3)2x[n−2]+(a3)3x[n−3]y[n] = \left(\frac{a}{3}\right)^3 x[n+3] + \left(\frac{a}{3}\right)^2 x[n+2] + \frac{a}{3} x[n+1] + x[n] + \frac{a}{3} x[n-1] + \left(\frac{a}{3}\right)^2 x[n-2] + \left(\frac{a}{3}\right)^3 x[n-3]

這個系統是非因果的,因為它包含了未來輸入 x[n+3],x[n+2],x[n+1]x[n+3], x[n+2], x[n+1]。
我們需要找到其因果版本。
系統函數 H(z)H(z) 是對差分方程式進行 Z 轉換得到的。
H(z)=∑k=−33(a3)∣k∣z−kH(z) = \sum_{k=-3}^{3} \left(\frac{a}{3}\right)^{|k|} z^{-k}
H(z)=(a3)3z−(−3)+(a3)2z−(−2)+(a3)1z−(−1)+(a3)0z0+(a3)1z−1+(a3)2z−2+(a3)3z−3H(z) = \left(\frac{a}{3}\right)^3 z^{-(-3)} + \left(\frac{a}{3}\right)^2 z^{-(-2)} + \left(\frac{a}{3}\right)^1 z^{-(-1)} + \left(\frac{a}{3}\right)^0 z^{0} + \left(\frac{a}{3}\right)^1 z^{-1} + \left(\frac{a}{3}\right)^2 z^{-2} + \left(\frac{a}{3}\right)^3 z^{-3}
H(z)=(a3)3z3+(a3)2z2+a3z+1+a3z−1+(a3)2z−2+(a3)3z−3H(z) = \left(\frac{a}{3}\right)^3 z^3 + \left(\frac{a}{3}\right)^2 z^2 + \frac{a}{3} z + 1 + \frac{a}{3} z^{-1} + \left(\frac{a}{3}\right)^2 z^{-2} + \left(\frac{a}{3}\right)^3 z^{-3}

為了使系統因果,我們需要將所有 zkz^k 的項(其中 k>0k>0)移到分母。
非因果系統的系統函數 H(z)H(z) 通常可以寫成 H(z)=B(z)A(z)H(z) = \frac{B(z)}{A(z)},其中 B(z)B(z) 和 A(z)A(z) 是 zz 的多項式。
在這個例子中,H(z)H(z) 是一個雙邊 Z 轉換的形式。

我們需要找到一個因果系統 Hcausal(z)H_{causal}(z),它具有相同的幅度響應 ∣H(ejΩ)∣=∣Hcausal(ejΩ)∣|H(e^{j\Omega})| = |H_{causal}(e^{j\Omega})|,並且具有最小群延遲。
幅度響應 ∣H(ejΩ)∣|H(e^{j\Omega})| 是 H(z)H(z) 在單位圓上的值。
H(z)H(z) 的幅度響應是 ∣H(ejΩ)∣=∣∑k=−33(a3)∣k∣e−jkΩ∣|H(e^{j\Omega})| = |\sum_{k=-3}^{3} (\frac{a}{3})^{|k|} e^{-jk\Omega}|。
注意到 H(z)H(z) 的係數是關於 k=0k=0 對稱的(因為 ∣k∣=∣−k∣|k|=|-k|)。
H(z)=1+∑k=13(a3)k(zk+z−k)H(z) = 1 + \sum_{k=1}^{3} (\frac{a}{3})^k (z^k + z^{-k})。
這表明 H(z)H(z) 是一個對稱的(或稱為回文)多項式。
對於這種對稱的系統,其頻率響應 H(ejΩ)H(e^{j\Omega}) 是實數的(如果 aa 是實數),這意味著相位響應是線性的,群延遲是常數。

一個非因果的 FIR 濾波器 h[n]h[n],如果其係數是對稱的(h[n]=h[N−n]h[n] = h[N-n]),則其頻率響應的相位是線性的。
在這個例子中,中心是 n=0n=0,係數關於 k=0k=0 對稱。
h[n]=(a3)∣n∣h[n] = (\frac{a}{3})^{|n|} for ∣n∣≤3|n| \le 3,且 h[n]=0h[n]=0 for ∣n∣>3|n|>3。
h[3]=h[−3]h[3] = h[-3], h[2]=h[−2]h[2] = h[-2], h[1]=h[−1]h[1] = h[-1]。
中心點是 h[0]h[0]。
這個系統的脈衝響應是 h[n]=(a3)∣n∣h[n] = (\frac{a}{3})^{|n|} for n∈{−3,−2,−1,0,1,2,3}n \in \{-3, -2, -1, 0, 1, 2, 3\}。
這是一個非因果的 FIR 系統。
為了得到一個因果的 FIR 系統,我們需要將脈衝響應向右平移,使其第一個非零點在 n=0n=0 或之後。
由於非因果系統的脈衝響應是 h[n]h[n],對於 n=−3,−2,−1,0,1,2,3n=-3, -2, -1, 0, 1, 2, 3,我們需要將它平移 3 個單位,得到 hcausal[n]=h[n−3]h_{causal}[n] = h[n-3]。
那麼,hcausal[n]=(a3)∣n−3∣h_{causal}[n] = (\frac{a}{3})^{|n-3|} for n−3∈{−3,−2,−1,0,1,2,3}n-3 \in \{-3, -2, -1, 0, 1, 2, 3\}。
hcausal[n]=(a3)∣n−3∣h_{causal}[n] = (\frac{a}{3})^{|n-3|} for n∈{0,1,2,3,4,5,6}n \in \{0, 1, 2, 3, 4, 5, 6\}。
這個因果系統的脈衝響應是 hcausal[n]=(a3)∣n−3∣h_{causal}[n] = (\frac{a}{3})^{|n-3|} for 0≤n≤60 \le n \le 6。
其系統函數是 Hcausal(z)=∑n=06(a3)∣n−3∣z−nH_{causal}(z) = \sum_{n=0}^{6} (\frac{a}{3})^{|n-3|} z^{-n}。

這個因果系統 hcausal[n]h_{causal}[n] 的幅度響應 ∣Hcausal(ejΩ)∣|H_{causal}(e^{j\Omega})| 會和原來的非因果系統的幅度響應 ∣H(ejΩ)∣|H(e^{j\Omega})| 相同嗎?
對於 FIR 濾波器,如果係數是對稱的,那麼相位響應是線性的。
原來的脈衝響應 h[n]h[n] 是關於 n=0n=0 對稱的。
平移後的脈衝響應 hcausal[n]h_{causal}[n] 是關於 n=3n=3 對稱的。
hcausal[n]=(a3)∣n−3∣h_{causal}[n] = (\frac{a}{3})^{|n-3|} for 0≤n≤60 \le n \le 6。
hcausal[0]=(a3)3h_{causal}[0] = (\frac{a}{3})^3, hcausal[1]=(a3)2h_{causal}[1] = (\frac{a}{3})^2, hcausal[2]=(a3)1h_{causal}[2] = (\frac{a}{3})^1, hcausal[3]=(a3)0=1h_{causal}[3] = (\frac{a}{3})^0 = 1, hcausal[4]=(a3)1h_{causal}[4] = (\frac{a}{3})^1, hcausal[5]=(a3)2h_{causal}[5] = (\frac{a}{3})^2, hcausal[6]=(a3)3h_{causal}[6] = (\frac{a}{3})^3。
這些係數是對稱的。
對於 FIR 濾波器,如果係數是對稱的,則其頻率響應 H(ejΩ)H(e^{j\Omega}) 的相位響應是線性的,群延遲是常數。
而且,幅度響應 ∣H(ejΩ)∣|H(e^{j\Omega})| 與係數的對稱性有關。
如果 h[n]h[n] 是關於 n=0n=0 對稱的,那麼 H(ejΩ)H(e^{j\Omega}) 的相位是線性的(或為零)。
如果 hcausal[n]h_{causal}[n] 是關於 n=N/2n=N/2 對稱的(其中 NN 是濾波器的階數),那麼其頻率響應的相位也是線性的。

要求「相同的幅度響應」。
非因果系統的系統函數是 H(z)=∑k=−33(a3)∣k∣z−kH(z) = \sum_{k=-3}^{3} (\frac{a}{3})^{|k|} z^{-k}。
因果系統的系統函數是 Hcausal(z)=∑n=06(a3)∣n−3∣z−nH_{causal}(z) = \sum_{n=0}^{6} (\frac{a}{3})^{|n-3|} z^{-n}。
我們需要 ∣H(ejΩ)∣=∣Hcausal(ejΩ)∣|H(e^{j\Omega})| = |H_{causal}(e^{j\Omega})|。
這通常要求 H(z)H(z) 和 Hcausal(z)H_{causal}(z) 之間有特定的關係。

考慮一個 FIR 濾波器 h[n]h[n]。
如果 h[n]h[n] 的係數是對稱的,h[n]=h[N−n]h[n] = h[N-n],則其頻率響應為 H(ejΩ)=e−jΩN/2∑n=0Nh[n]e−jΩ(n−N/2)H(e^{j\Omega}) = e^{-j\Omega N/2} \sum_{n=0}^{N} h[n] e^{-j\Omega (n-N/2)}。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 九、(10%) 題10 分

The input x[n]x[n] and output y[n]y[n] of a causal linear time invariant system are related through the block-diagram representation shown in the following figure.

🖼️【此處有附圖,請對照原卷】
Figure: A block diagram of a causal LTI system. The input x[n]x[n] goes through a summing junction. One input to the summer is x[n]x[n]. Another input comes from a path with a gain of −1/3-1/3 and a z−1z^{-1} block, which feeds into another summing junction. The output of the first summing junction goes through a z−1z^{-1} block, and also through a path with a gain of 2/92/9 and a z−1z^{-1} block. The output of the z−1z^{-1} block from the first summing junction and the output of the 2/92/9 gain block are fed into the second summing junction. The output of the second summing junction goes through a gain of −5/8-5/8 and then into the final summing junction, which produces the output y[n]y[n]. The input x[n]x[n] also directly feeds into the final summing junction.

(一) (5%) Determine a difference equation relating y[n]y[n] and x[n]x[n].
(二) (5%) Is this system stable?

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 3 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  • 方塊圖轉換為差分方程。
  • 延遲單元 z−1z^{-1} 的時域意義:
    z−1v[n]=v[n−1]z^{-1}v[n]=v[n-1]
  • 利用 ZZ 轉換消去內部訊號。
  • 因果 LTI 系統的 BIBO 穩定條件:所有極點必須位於單位圓內。

解題方法

令第一個加法器輸出為 v[n]v[n]。由圖中的回授路徑可得,第一個加法器接收:

  • 輸入 x[n]x[n];
  • 經過 z−1z^{-1} 與增益 −13-\frac{1}{3} 的訊號;
  • 經過兩個 z−1z^{-1} 與增益 29\frac{2}{9} 的訊號。

因此:

v[n]=x[n]−13v[n−1]+29v[n−2]v[n]=x[n]-\frac{1}{3}v[n-1]+\frac{2}{9}v[n-2]

整理為:

v[n]+13v[n−1]−29v[n−2]=x[n]v[n]+\frac{1}{3}v[n-1]-\frac{2}{9}v[n-2]=x[n]

最後一個加法器將 v[n]v[n] 與經過 z−1z^{-1}、增益 −58-\frac{5}{8} 的訊號相加,因此:

y[n]=v[n]−58v[n−1]y[n]=v[n]-\frac{5}{8}v[n-1]

取 ZZ 轉換:

V(z)(1+13z−1−29z−2)=X(z)V(z)\left(1+\frac{1}{3}z^{-1}-\frac{2}{9}z^{-2}\right)=X(z)

且

Y(z)=(1−58z−1)V(z)Y(z)=\left(1-\frac{5}{8}z^{-1}\right)V(z)

消去 V(z)V(z),得到:

(1+13z−1−29z−2)Y(z)=(1−58z−1)X(z)\left(1+\frac{1}{3}z^{-1}-\frac{2}{9}z^{-2}\right)Y(z) = \left(1-\frac{5}{8}z^{-1}\right)X(z)

轉回時域:

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

其他考古題