112 年 國立臺灣大學電信工程學研究所乙組《通信原理》

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

第 1 題

Suppose that y[n]−0.5y[n−1]=x[n]y[n] - 0.5y[n-1] = x[n].
(a) Please determine the causal sequence h[n]h[n] such that y[n]=x[n]∗h[n]y[n] = x[n] * h[n] where ∗* means the convolution.
(b) Also suppose that 3f[n]−f[n−1]=y[n]−y[n−1]3f[n] - f[n-1] = y[n] - y[n-1]. Determine the causal sequence g[n]g[n] such that f[n]=x[n]∗g[n]f[n] = x[n] * g[n].

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

這一題的完整詳解

這題主要考查離散時間系統的脈衝響應(h[n]h[n])與系統的串接(g[n]g[n])的計算。

(a) 求系統 h[n]h[n]

題目給定的差分方程式為 y[n]−0.5y[n−1]=x[n]y[n] - 0.5y[n-1] = x[n]。
我們知道,對於 LTI 系統,若輸入為 x[n]x[n],輸出為 y[n]y[n],且 y[n]=x[n]∗h[n]y[n] = x[n] * h[n],其中 h[n]h[n] 是系統的脈衝響應。
我們可以對差分方程式進行 Z 轉換。
Z{y[n]}−0.5Z{y[n−1]}=Z{x[n]}Z\{y[n]\} - 0.5Z\{y[n-1]\} = Z\{x[n]\}
Y(z)−0.5z−1Y(z)=X(z)Y(z) - 0.5z^{-1}Y(z) = X(z)
Y(z)(1−0.5z−1)=X(z)Y(z)(1 - 0.5z^{-1}) = X(z)
系統的傳遞函數 H(z)H(z) 定義為 Y(z)X(z)\frac{Y(z)}{X(z)}:
H(z)=Y(z)X(z)=11−0.5z−1H(z) = \frac{Y(z)}{X(z)} = \frac{1}{1 - 0.5z^{-1}}
為了求取脈衝響應 h[n]h[n],我們需要對 H(z)H(z) 進行反 Z 轉換。
H(z)=11−0.5z−1H(z) = \frac{1}{1 - 0.5z^{-1}}
這是一個幾何級數的 Z 轉換形式 ∑n=0∞anz−n=11−az−1\sum_{n=0}^{\infty} a^n z^{-n} = \frac{1}{1-az^{-1}}。
因此,對於因果系統(causal system),其脈衝響應為 h[n]=(0.5)nu[n]h[n] = (0.5)^n u[n]。

【答案】h[n]=(0.5)nu[n]h[n] = (0.5)^n u[n]

(b) 求系統 g[n]g[n]

題目給定兩個方程式:

  1. y[n]−0.5y[n−1]=x[n]y[n] - 0.5y[n-1] = x[n]
  2. 3f[n]−f[n−1]=y[n]−y[n−1]3f[n] - f[n-1] = y[n] - y[n-1]
    我們要求 f[n]=x[n]∗g[n]f[n] = x[n] * g[n] 中的 g[n]g[n],其中 f[n]f[n] 是輸出,x[n]x[n] 是輸入。

首先,我們利用 Z 轉換來分析第二個方程式。
Z{3f[n]−f[n−1]}=Z{y[n]−y[n−1]}Z\{3f[n] - f[n-1]\} = Z\{y[n] - y[n-1]\}
3F(z)−z−1F(z)=Y(z)−z−1Y(z)3F(z) - z^{-1}F(z) = Y(z) - z^{-1}Y(z)
F(z)(3−z−1)=Y(z)(1−z−1)F(z)(3 - z^{-1}) = Y(z)(1 - z^{-1})

我們已經知道從第一個方程式 y[n]−0.5y[n−1]=x[n]y[n] - 0.5y[n-1] = x[n] 得到 Y(z)=H(z)X(z)Y(z) = H(z)X(z),其中 H(z)=11−0.5z−1H(z) = \frac{1}{1 - 0.5z^{-1}}。
將 Y(z)Y(z) 代入上面的方程式:
F(z)(3−z−1)=H(z)X(z)(1−z−1)F(z)(3 - z^{-1}) = H(z)X(z)(1 - z^{-1})
F(z)=H(z)(1−z−1)3−z−1X(z)F(z) = \frac{H(z)(1 - z^{-1})}{3 - z^{-1}} X(z)

因為我們要求 f[n]=x[n]∗g[n]f[n] = x[n] * g[n],所以 F(z)=G(z)X(z)F(z) = G(z)X(z)。
因此,系統 g[n]g[n] 的傳遞函數 G(z)G(z) 為:

🔒

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

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

免費註冊

第 2 題

Determine the following convolution results (* means the convolution).
(a) (6%)(6\%) [3,3,2,2,1,1]∗[−1,2,−1][3, 3, 2, 2, 1, 1] * [-1, 2, -1];
(b) (6%)(6\%) sinc(t)∗sinc(3t)∗sinc(5t)∗δ(t−6)\text{sinc}(t) * \text{sinc}(3t) * \text{sinc}(5t) * \delta(t -6);
(c) (6%)(6\%) sinc2(3t/2)∗(1+sin⁡(2πt)+cos⁡(4πt))\text{sinc}^2(3t/2) * (1+\sin(2\pi t) + \cos(4\pi t)).

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

這一題的完整詳解

這題主要考查離散時間和連續時間的卷積運算。

(a) 離散時間卷積

給定兩個離散序列:x[n]=[3,3,2,2,1,1]x[n] = [3, 3, 2, 2, 1, 1] 和 h[n]=[−1,2,−1]h[n] = [-1, 2, -1]。
序列 x[n]x[n] 的長度為 N=6N=6,序列 h[n]h[n] 的長度為 M=3M=3。
卷積結果 y[n]=x[n]∗h[n]y[n] = x[n] * h[n] 的長度為 N+M−1=6+3−1=8N+M-1 = 6+3-1 = 8。
卷積的定義為 y[n]=∑k=−∞∞x[k]h[n−k]y[n] = \sum_{k=-\infty}^{\infty} x[k]h[n-k]。

我們可以逐項計算:
y[0]=x[0]h[0]=3×(−1)=−3y[0] = x[0]h[0] = 3 \times (-1) = -3
y[1]=x[0]h[1]+x[1]h[0]=3×2+3×(−1)=6−3=3y[1] = x[0]h[1] + x[1]h[0] = 3 \times 2 + 3 \times (-1) = 6 - 3 = 3
y[2]=x[0]h[2]+x[1]h[1]+x[2]h[0]=3×(−1)+3×2+2×(−1)=−3+6−2=1y[2] = x[0]h[2] + x[1]h[1] + x[2]h[0] = 3 \times (-1) + 3 \times 2 + 2 \times (-1) = -3 + 6 - 2 = 1
y[3]=x[1]h[2]+x[2]h[1]+x[3]h[0]=3×(−1)+2×2+2×(−1)=−3+4−2=−1y[3] = x[1]h[2] + x[2]h[1] + x[3]h[0] = 3 \times (-1) + 2 \times 2 + 2 \times (-1) = -3 + 4 - 2 = -1
y[4]=x[2]h[2]+x[3]h[1]+x[4]h[0]=2×(−1)+2×2+1×(−1)=−2+4−1=1y[4] = x[2]h[2] + x[3]h[1] + x[4]h[0] = 2 \times (-1) + 2 \times 2 + 1 \times (-1) = -2 + 4 - 1 = 1
y[5]=x[3]h[2]+x[4]h[1]+x[5]h[0]=2×(−1)+1×2+1×(−1)=−2+2−1=−1y[5] = x[3]h[2] + x[4]h[1] + x[5]h[0] = 2 \times (-1) + 1 \times 2 + 1 \times (-1) = -2 + 2 - 1 = -1
y[6]=x[4]h[2]+x[5]h[1]=1×(−1)+1×2=−1+2=1y[6] = x[4]h[2] + x[5]h[1] = 1 \times (-1) + 1 \times 2 = -1 + 2 = 1
y[7]=x[5]h[2]=1×(−1)=−1y[7] = x[5]h[2] = 1 \times (-1) = -1

所以,卷積結果為 [−3,3,1,−1,1,−1,1,−1][-3, 3, 1, -1, 1, -1, 1, -1]。

【答案】[−3,3,1,−1,1,−1,1,−1][-3, 3, 1, -1, 1, -1, 1, -1]

(b) 連續時間卷積

題目要求計算 sinc(t)∗sinc(3t)∗sinc(5t)∗δ(t−6)\text{sinc}(t) * \text{sinc}(3t) * \text{sinc}(5t) * \delta(t -6)。
我們知道卷積的性質:

  1. sinc(at)∗sinc(bt)=1max⁡(a,b)sinc(max⁡(a,b)t)\text{sinc}(at) * \text{sinc}(bt) = \frac{1}{\max(a,b)} \text{sinc}(\max(a,b)t)。
  2. x(t)∗δ(t−t0)=x(t−t0)x(t) * \delta(t-t_0) = x(t-t_0)。

首先,我們計算 sinc(t)∗sinc(3t)\text{sinc}(t) * \text{sinc}(3t)。
這裡 a=1,b=3a=1, b=3。
sinc(t)∗sinc(3t)=1max⁡(1,3)sinc(max⁡(1,3)t)=13sinc(3t)\text{sinc}(t) * \text{sinc}(3t) = \frac{1}{\max(1,3)} \text{sinc}(\max(1,3)t) = \frac{1}{3} \text{sinc}(3t)。

接著,計算 (13sinc(3t))∗sinc(5t)(\frac{1}{3} \text{sinc}(3t)) * \text{sinc}(5t)。
這裡我們可以將常數提出:13(sinc(3t)∗sinc(5t))\frac{1}{3} (\text{sinc}(3t) * \text{sinc}(5t))。
對於 sinc(3t)∗sinc(5t)\text{sinc}(3t) * \text{sinc}(5t),這裡 a=3,b=5a=3, b=5。
sinc(3t)∗sinc(5t)=1max⁡(3,5)sinc(max⁡(3,5)t)=15sinc(5t)\text{sinc}(3t) * \text{sinc}(5t) = \frac{1}{\max(3,5)} \text{sinc}(\max(3,5)t) = \frac{1}{5} \text{sinc}(5t)。
所以,(13sinc(3t))∗sinc(5t)=13⋅15sinc(5t)=115sinc(5t)(\frac{1}{3} \text{sinc}(3t)) * \text{sinc}(5t) = \frac{1}{3} \cdot \frac{1}{5} \text{sinc}(5t) = \frac{1}{15} \text{sinc}(5t)。

最後,計算 (115sinc(5t))∗δ(t−6)(\frac{1}{15} \text{sinc}(5t)) * \delta(t -6)。
利用卷積與衝激響應的性質:
(115sinc(5t))∗δ(t−6)=115sinc(5(t−6))(\frac{1}{15} \text{sinc}(5t)) * \delta(t -6) = \frac{1}{15} \text{sinc}(5(t-6))。

【答案】115sinc(5(t−6))\frac{1}{15} \text{sinc}(5(t-6))

(c) 連續時間卷積

題目要求計算 sinc2(3t/2)∗(1+sin⁡(2πt)+cos⁡(4πt))\text{sinc}^2(3t/2) * (1+\sin(2\pi t) + \cos(4\pi t))。
我們將卷積看成三個部分的卷積之和:

  1. sinc2(3t/2)∗1\text{sinc}^2(3t/2) * 1
  2. sinc2(3t/2)∗sin⁡(2πt)\text{sinc}^2(3t/2) * \sin(2\pi t)
  3. sinc2(3t/2)∗cos⁡(4πt)\text{sinc}^2(3t/2) * \cos(4\pi t)

首先,計算 sinc2(3t/2)∗1\text{sinc}^2(3t/2) * 1。
我們知道 x(t)∗δ(t)=x(t)x(t) * \delta(t) = x(t)。這裡的 11 指的是 δ(t)\delta(t)。
所以,sinc2(3t/2)∗1=sinc2(3t/2)\text{sinc}^2(3t/2) * 1 = \text{sinc}^2(3t/2)。

接著,我們需要計算 sinc2(3t/2)\text{sinc}^2(3t/2) 與 sin⁡(2πt)\sin(2\pi t) 和 cos⁡(4πt)\cos(4\pi t) 的卷積。
我們需要知道 sinc2(at)\text{sinc}^2(at) 的傅立葉轉換。
sinc(t)=sin⁡(πt)πt\text{sinc}(t) = \frac{\sin(\pi t)}{\pi t}。
sinc(at)=sin⁡(πat)πat\text{sinc}(at) = \frac{\sin(\pi at)}{\pi at}。
sinc2(at)=(sin⁡(πat)πat)2\text{sinc}^2(at) = (\frac{\sin(\pi at)}{\pi at})^2。

我們知道 sinc(t)\text{sinc}(t) 的傅立葉轉換是 1πrect(ω2π)\frac{1}{\pi} \text{rect}(\frac{\omega}{2\pi})。
sinc(at)\text{sinc}(at) 的傅立葉轉換是 1a1πrect(ω2πa)=1πarect(ω2πa)\frac{1}{a} \frac{1}{\pi} \text{rect}(\frac{\omega}{2\pi a}) = \frac{1}{\pi a} \text{rect}(\frac{\omega}{2\pi a})。
sinc2(at)\text{sinc}^2(at) 的傅立葉轉換是 sinc(at)∗sinc(at)\text{sinc}(at) * \text{sinc}(at) 的傅立葉轉換,這會得到一個三角形函數。

🔒

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

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

免費註冊

第 3 題

Determine x[n]x[n] and y[n]y[n] if their Z transform are as follows. Only consider the case where x[n]x[n] and y[n]y[n] are causal.
(a) X(z)=cos⁡(z3)X(z) = \cos(z^3);
(b) Y(z)=z2−z−1z3−12z2−z−2Y(z) = \frac{z^2 - z - 1}{z^3 - 12z^2 - z - 2}

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

這一題的完整詳解

核心觀念

本題考查:

  1. 雙邊 Z 轉換定義
    X(z)=∑n=−∞∞x[n]z−nX(z)=\sum_{n=-\infty}^{\infty}x[n]z^{-n}

  2. 因果序列的判斷:若 x[n]x[n] 為因果序列,則 x[n]=0x[n]=0(n<0n<0),其 Z 轉換只能包含 z−nz^{-n} 的非負次方。

  3. 有理式 Z 轉換的部分分式分解,以及基本反轉換公式

11−pz−1↔Z−1pnu[n]\frac{1}{1-pz^{-1}} \xleftrightarrow{\mathcal{Z}^{-1}} p^n u[n]

(a) X(z)=cos⁡(z3)X(z)=\cos(z^3)

解題方法

將 cosine 展開成冪級數:

cos⁡(z3)=∑k=0∞(−1)k(z3)2k(2k)!=∑k=0∞(−1)kz6k(2k)!\cos(z^3) =\sum_{k=0}^{\infty}\frac{(-1)^k(z^3)^{2k}}{(2k)!} =\sum_{k=0}^{\infty}\frac{(-1)^kz^{6k}}{(2k)!}

根據 Z 轉換定義,z6kz^{6k} 可寫成

z6k=z−(−6k)z^{6k}=z^{-(-6k)}

因此可直接對照得到

x[−6k]=(−1)k(2k)!,k=0,1,2,…x[-6k]=\frac{(-1)^k}{(2k)!}, \qquad k=0,1,2,\ldots

其餘位置皆為零。若不考慮因果限制,序列為

x[n]={(−1)k(2k)!,n=−6k,k=0,1,2,…0,其他x[n]= \begin{cases} \dfrac{(-1)^k}{(2k)!}, & n=-6k,\quad k=0,1,2,\ldots\\[6pt] 0, & \text{其他} \end{cases}

也可寫成

x[n]=∑k=0∞(−1)k(2k)!δ[n+6k]x[n]=\sum_{k=0}^{\infty} \frac{(-1)^k}{(2k)!}\delta[n+6k]

其中 k=0k=0 對應 x[0]=1x[0]=1,k=1k=1 對應 x[−6]=−1/2!x[-6]=-1/2!,k=2k=2 對應 x[−12]=1/4!x[-12]=1/4!。

因果性判斷

因果序列必須滿足

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

但上述 Z 轉換明確包含 z6,z12,z18,…z^6,z^{12},z^{18},\ldots 等正次方項,這些項分別對應到 x[−6],x[−12],x[−18],…x[-6],x[-12],x[-18],\ldots,因此必然存在負時間索引的非零樣本。

所以,X(z)=cos⁡(z3)X(z)=\cos(z^3) 不存在同時符合題目所要求的因果序列 x[n]x[n]。

解題技巧

對因果序列而言,Z 轉換應表示成 z−1z^{-1} 的冪級數:

X(z)=x[0]+x[1]z−1+x[2]z−2+⋯X(z)=x[0]+x[1]z^{-1}+x[2]z^{-2}+\cdots

只要看到轉換式含有非零的正次方 z,z2,…z,z^2,\ldots,即可立即判斷其對應序列含有負時間樣本,不可能是因果序列。


(b)

Y(z)=z2−z−1z3−12z2−z−2Y(z)=\frac{z^2-z-1}{z^3-12z^2-z-2}

化為 z−1z^{-1} 形式

分子、分母同除以 z3z^3:

Y(z)=z−1−z−2−z−31−12z−1−z−2−2z−3Y(z) = \frac{z^{-1}-z^{-2}-z^{-3}} {1-12z^{-1}-z^{-2}-2z^{-3}}

因為題目指定 y[n]y[n] 為因果序列,所以其 ROC 位於最外側極點之外,並採用右側序列的反轉換形式。


方法一:部分分式分解

令分母的三個根為 p1,p2,p3p_1,p_2,p_3:

z3−12z2−z−2=(z−p1)(z−p2)(z−p3)z^3-12z^2-z-2 =(z-p_1)(z-p_2)(z-p_3)

三根約為

p1≈12.09634p_1\approx 12.09634 p2≈−0.04817+j0.40376,p3≈−0.04817−j0.40376p_2\approx -0.04817+j0.40376, \qquad p_3\approx -0.04817-j0.40376

將 Y(z)Y(z) 寫成

Y(z)=C11−p1z−1+C21−p2z−1+C31−p3z−1Y(z)= \frac{C_1}{1-p_1z^{-1}} +\frac{C_2}{1-p_2z^{-1}} +\frac{C_3}{1-p_3z^{-1}}

其中係數可由留數公式求得:

Ci=pi2−pi−1pi(3pi2−24pi−1)C_i= \frac{p_i^2-p_i-1} {p_i\left(3p_i^2-24p_i-1\right)}
🔒

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

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

免費註冊

第 4 題

Suppose that X(jω)X(j\omega) is the continuous-time Fourier transform of x(t)x(t) and X(jω)≠0X(j\omega) \neq 0 for ∣ω∣<6000π|\omega| < 6000\pi, X(jω)=0X(j\omega) = 0 for ∣ω∣>6000π|\omega| > 6000\pi. Determine what is the lower bound of the sampling interval if we want to sample y(t)y(t) and z(t)z(t) as follows without the aliasing effect. In (a) and (b), ∗* means the convolution and ′' means the derivative.
(a) y(t)=x(t)∗cos⁡(4000πt)y(t) = x(t) * \cos(4000\pi t)
(b) z(t)=x(t−1)x(2t+2)∗x(4t)z(t) = x(t-1)x(2t+2) * x(4t)

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

這一題的完整詳解

這題主要考查奈奎斯特採樣定理在不同信號形式下的應用。核心是確定每個信號的頻寬,然後計算採樣間隔的下限。

奈奎斯特採樣定理
一個帶寬為 BB Hz 的帶限信號,若要無失真地恢復,則採樣頻率 fsf_s 必須大於或等於 2B2B(即 fs≥2Bf_s \ge 2B)。對應的採樣間隔 Ts=1/fsT_s = 1/f_s,則 Ts≤12BT_s \le \frac{1}{2B}。
在頻率 ω\omega (弧度/秒) 的單位下,帶寬為 W=B⋅2πW = B \cdot 2\pi 弧度/秒。採樣角頻率 ωs=fs⋅2π\omega_s = f_s \cdot 2\pi。則 ωs≥2W\omega_s \ge 2W。採樣間隔 TsT_s 滿足 Ts≤πWT_s \le \frac{\pi}{W}。

題目給定 X(jω)X(j\omega) 的頻寬為 ∣ω∣<6000π|\omega| < 6000\pi。這意味著 x(t)x(t) 的最高角頻率是 Wx=6000πW_x = 6000\pi rad/s。
根據奈奎斯特準則,為了無失真地採樣 x(t)x(t),採樣角頻率 ωs\omega_s 必須滿足 ωs>2Wx=12000π\omega_s > 2 W_x = 12000\pi rad/s。
對應的採樣間隔 TsT_s 的上限為 Ts≤πWx=π6000π=16000T_s \le \frac{\pi}{W_x} = \frac{\pi}{6000\pi} = \frac{1}{6000} 秒。

我們需要分別計算 y(t)y(t) 和 z(t)z(t) 的頻寬,然後確定各自的採樣間隔下限。

(a) y(t)=x(t)∗cos⁡(4000πt)y(t) = x(t) * \cos(4000\pi t)

卷積在時域對應於傅立葉轉換域的乘積。
Y(jω)=X(jω)⋅F{cos⁡(4000πt)}Y(j\omega) = X(j\omega) \cdot \mathcal{F}\{\cos(4000\pi t)\}。
cos⁡(4000πt)=12(ej4000πt+e−j4000πt)\cos(4000\pi t) = \frac{1}{2} (e^{j4000\pi t} + e^{-j4000\pi t})。
F{cos⁡(4000πt)}=π(δ(ω−4000π)+δ(ω+4000π))\mathcal{F}\{\cos(4000\pi t)\} = \pi (\delta(\omega - 4000\pi) + \delta(\omega + 4000\pi))。

Y(jω)=X(jω)⋅π(δ(ω−4000π)+δ(ω+4000π))Y(j\omega) = X(j\omega) \cdot \pi (\delta(\omega - 4000\pi) + \delta(\omega + 4000\pi))。
Y(jω)=π[X(j(4000π))δ(ω−4000π)+X(j(−4000π))δ(ω+4000π)]Y(j\omega) = \pi [X(j(4000\pi)) \delta(\omega - 4000\pi) + X(j(-4000\pi)) \delta(\omega + 4000\pi)]。

我們知道 X(jω)X(j\omega) 的頻寬是 ∣ω∣<6000π|\omega| < 6000\pi。
因此,X(j4000π)X(j4000\pi) 和 X(j(−4000π))X(j(-4000\pi)) 都是非零的(因為 ∣±4000π∣<6000π| \pm 4000\pi | < 6000\pi)。
Y(jω)Y(j\omega) 的頻譜由兩個衝激函數組成,位於 ω=4000π\omega = 4000\pi 和 ω=−4000π\omega = -4000\pi。
這意味著 y(t)y(t) 的頻寬是 Wy=4000πW_y = 4000\pi rad/s。
注意:嚴格來說,如果 X(jω)X(j\omega) 在 ∣ω∣≥4000π|\omega| \ge 4000\pi 的地方有值,卷積的頻寬會擴展。但是,這裡的 cos⁡\cos 函數在頻域只包含兩個離散頻率點 ±4000π\pm 4000\pi。
Y(jω)Y(j\omega) 的頻譜由 X(jω)X(j\omega) 在 ω=4000π\omega = 4000\pi 和 ω=−4000π\omega = -4000\pi 的值決定。
如果 X(jω)X(j\omega) 在 ∣ω∣<6000π|\omega| < 6000\pi 範圍內是非零的,那麼 Y(jω)Y(j\omega) 將包含在 ω=±4000π\omega = \pm 4000\pi 的衝激,並且其值由 X(j(±4000π))X(j(\pm 4000\pi)) 決定。
Y(jω)Y(j\omega) 的非零部分僅存在於 ω=4000π\omega = 4000\pi 和 ω=−4000π\omega = -4000\pi。
因此,y(t)y(t) 的有效頻寬為 Wy=4000πW_y = 4000\pi rad/s。

根據奈奎斯特準則,採樣角頻率 ωs\omega_s 必須大於 2Wy=2×4000π=8000π2 W_y = 2 \times 4000\pi = 8000\pi rad/s。
採樣間隔 TsT_s 的下限為 Ts≤πWy=π4000π=14000T_s \le \frac{\pi}{W_y} = \frac{\pi}{4000\pi} = \frac{1}{4000} 秒。

🔒

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

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

免費註冊

第 5 題

Let X(t)=rsin⁡(ωt+θ)X(t) = r \sin(\omega t + \theta), where rr and ω\omega are non-zero constants and θ\theta is uniformly distributed between 00 and 2π2\pi.
(a) (4%) please find function R(t,s)=E[X(t)X(s)]R(t,s) = E[X(t)X(s)].
(b) (5%) Assume that R(0,3)=2R(0,3)=2, please find R(0.5,3.5)R(0.5,3.5).

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

這一題的完整詳解

這題考查隨機信號的自相關函數的計算,以及利用自相關函數的性質求解未知參數。

核心概念:

  1. 隨機變數的期望值: 對於一個連續型隨機變數 θ\theta 均勻分佈在 [a,b][a, b] 上,其概率密度函數 (PDF) 為 f(θ)=1b−af(\theta) = \frac{1}{b-a}。期望值 E[θ]=∫abθf(θ)dθE[\theta] = \int_a^b \theta f(\theta) d\theta。
  2. 自相關函數: 對於一個隨機過程 X(t)X(t),其自相關函數定義為 RXX(t,s)=E[X(t)X(s)]R_{XX}(t, s) = E[X(t)X(s)]。
  3. 統計獨立性: 如果 X(t)X(t) 和 X(s)X(s) 的統計特性(例如 E[X(t)]E[X(t)] 和 E[X(s)]E[X(s)])是確定的,並且 X(t)X(t) 和 X(s)X(s) 的隨機性來源於同一個隨機變數,則它們的乘積的期望值需要對該隨機變數進行平均。

給定條件:
X(t)=rsin⁡(ωt+θ)X(t) = r \sin(\omega t + \theta)
r,ωr, \omega 是常數。
θ\theta 是均勻分佈在 [0,2π][0, 2\pi] 上的隨機變數。
θ\theta 的 PDF 為 f(θ)=12πf(\theta) = \frac{1}{2\pi},對於 0≤θ≤2π0 \le \theta \le 2\pi,否則為 0。

(a) 求 R(t,s)=E[X(t)X(s)]R(t,s) = E[X(t)X(s)]

X(t)X(s)=(rsin⁡(ωt+θ))(rsin⁡(ωs+θ))X(t)X(s) = (r \sin(\omega t + \theta))(r \sin(\omega s + \theta))
X(t)X(s)=r2sin⁡(ωt+θ)sin⁡(ωs+θ)X(t)X(s) = r^2 \sin(\omega t + \theta) \sin(\omega s + \theta)

我們需要計算這個乘積的期望值,其中隨機性來自 θ\theta。
R(t,s)=E[X(t)X(s)]=E[r2sin⁡(ωt+θ)sin⁡(ωs+θ)]R(t,s) = E[X(t)X(s)] = E[r^2 \sin(\omega t + \theta) \sin(\omega s + \theta)]
由於 r2r^2 是常數,可以提出期望值運算:
R(t,s)=r2E[sin⁡(ωt+θ)sin⁡(ωs+θ)]R(t,s) = r^2 E[\sin(\omega t + \theta) \sin(\omega s + \theta)]

使用三角恆等式 2sin⁡Asin⁡B=cos⁡(A−B)−cos⁡(A+B)2\sin A \sin B = \cos(A-B) - \cos(A+B):
令 A=ωt+θA = \omega t + \theta,B=ωs+θB = \omega s + \theta。
A−B=(ωt+θ)−(ωs+θ)=ωt−ωs=ω(t−s)A-B = (\omega t + \theta) - (\omega s + \theta) = \omega t - \omega s = \omega(t-s)。
A+B=(ωt+θ)+(ωs+θ)=ωt+ωs+2θ=ω(t+s)+2θA+B = (\omega t + \theta) + (\omega s + \theta) = \omega t + \omega s + 2\theta = \omega(t+s) + 2\theta。

所以,sin⁡(ωt+θ)sin⁡(ωs+θ)=12[cos⁡(ω(t−s))−cos⁡(ω(t+s)+2θ)]\sin(\omega t + \theta) \sin(\omega s + \theta) = \frac{1}{2} [\cos(\omega(t-s)) - \cos(\omega(t+s) + 2\theta)]。

現在計算期望值:
E[sin⁡(ωt+θ)sin⁡(ωs+θ)]=E[12(cos⁡(ω(t−s))−cos⁡(ω(t+s)+2θ))]E[\sin(\omega t + \theta) \sin(\omega s + \theta)] = E[\frac{1}{2} (\cos(\omega(t-s)) - \cos(\omega(t+s) + 2\theta))]
=12E[cos⁡(ω(t−s))]−12E[cos⁡(ω(t+s)+2θ)]= \frac{1}{2} E[\cos(\omega(t-s))] - \frac{1}{2} E[\cos(\omega(t+s) + 2\theta)]

因為 ω(t−s)\omega(t-s) 不依賴於 θ\theta,所以 E[cos⁡(ω(t−s))]=cos⁡(ω(t−s))E[\cos(\omega(t-s))] = \cos(\omega(t-s))。

接下來計算 E[cos⁡(ω(t+s)+2θ)]E[\cos(\omega(t+s) + 2\theta)]:

🔒

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

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

免費註冊

第 6 題

Consider a binary communication system where a random bit bb is transmitted. Suppose the 2×12 \times 1 transmitted signal is
s=[(−1)bdmin/20]s = \begin{bmatrix} (-1)^{b d_{min}/2} \\ 0 \end{bmatrix}
where b∈{0,1}b \in \{0,1\} and dmin>0d_{min} > 0. The received signal is r=s+nr = s + n, where noise
n=[n1n2]n = \begin{bmatrix} n_1 \\ n_2 \end{bmatrix}
is a white Gaussian random vector that is independent of ss. T is the vector transpose.
(a) (4%) Assume that nin_i has mean μi\mu_i and variance σi2\sigma_i^2, i=1,2i=1, 2. Please find the mean of noise vector nn.
(b) (2%) What is the definition of the covariance matrix of a random vector?
(c) (4%) Please find the covariance matrix of nn using your answer in (b).
(d) (10%) Suppose that received
r=[r1r2]r = \begin{bmatrix} r_1 \\ r_2 \end{bmatrix}
Show that given r1r_1, r2r_2 is irrelevant; that is, r2r_2 is independent of ss.

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

這一題的完整詳解

這題考查隨機向量的期望、協方差矩陣的計算,以及獨立性的概念,特別是在通信系統的背景下。

核心概念:

  1. 隨機向量的期望: 對於一個 N×1N \times 1 的隨機向量 v=[v1,v2,…,vN]Tv = [v_1, v_2, \dots, v_N]^T,其期望向量為 E[v]=[E[v1],E[v2],…,E[vN]]TE[v] = [E[v_1], E[v_2], \dots, E[v_N]]^T。
  2. 協方差矩陣: 對於一個 N×1N \times 1 的隨機向量 vv,其協方差矩陣 CvC_v 是一個 N×NN \times N 的矩陣,定義為 Cv=E[(v−E[v])(v−E[v])T]C_v = E[(v - E[v])(v - E[v])^T]。其對角線元素是各分量的方差,非對角線元素是各分量之間的協方差。
  3. 白雜訊 (White Noise): 指的是一個隨機過程,其不同時間點的取值是統計獨立的,並且具有相同的概率分佈。對於離散時間的白雜訊向量,通常意味著各分量是統計獨立的,並且具有相同的均值(通常為 0)和方差。題目中提到「white Gaussian random vector」,意味著其分量是獨立的,且服從高斯分佈。
  4. 統計獨立性: 兩個隨機變數 XX 和 YY 是統計獨立的,如果 P(X∈A,Y∈B)=P(X∈A)P(Y∈B)P(X \in A, Y \in B) = P(X \in A)P(Y \in B) 對於所有事件 A,BA, B 成立。對於隨機向量,如果其所有分量之間都是統計獨立的,則稱向量是統計獨立的。
  5. 條件概率與獨立性: 如果 r2r_2 和 ss 是獨立的,那麼知道 r1r_1 的值不應該改變 r2r_2 的概率分佈。

給定條件:
s=[(−1)bdmin/20]s = \begin{bmatrix} (-1)^{b d_{min}/2} \\ 0 \end{bmatrix}, b∈{0,1}b \in \{0,1\}, dmin>0d_{min} > 0.
n=[n1n2]n = \begin{bmatrix} n_1 \\ n_2 \end{bmatrix}, n1,n2n_1, n_2 是隨機變數。
r=s+n=[r1r2]=[s1+n1s2+n2]r = s + n = \begin{bmatrix} r_1 \\ r_2 \end{bmatrix} = \begin{bmatrix} s_1 + n_1 \\ s_2 + n_2 \end{bmatrix}。
nn 是白高斯隨機向量,獨立於 ss。
E[n1]=μ1E[n_1] = \mu_1, Var(n1)=σ12Var(n_1) = \sigma_1^2.
E[n2]=μ2E[n_2] = \mu_2, Var(n2)=σ22Var(n_2) = \sigma_2^2.

(a) 求雜訊向量 nn 的期望

根據期望的線性性質,雜訊向量 nn 的期望是其分量期望組成的向量。
E[n]=E[[n1n2]]=[E[n1]E[n2]]E[n] = E\left[\begin{bmatrix} n_1 \\ n_2 \end{bmatrix}\right] = \begin{bmatrix} E[n_1] \\ E[n_2] \end{bmatrix}。
已知 E[n1]=μ1E[n_1] = \mu_1 和 E[n2]=μ2E[n_2] = \mu_2。
所以,E[n]=[μ1μ2]E[n] = \begin{bmatrix} \mu_1 \\ \mu_2 \end{bmatrix}。

【答案】E[n]=[μ1μ2]E[n] = \begin{bmatrix} \mu_1 \\ \mu_2 \end{bmatrix}。

(b) 隨機向量協方差矩陣的定義

對於一個 N×1N \times 1 的隨機向量 vv,其協方差矩陣 CvC_v 定義為:
Cv=E[(v−E[v])(v−E[v])T]C_v = E[(v - E[v])(v - E[v])^T]。
如果 v=[v1,v2,…,vN]Tv = [v_1, v_2, \dots, v_N]^T,則 E[v]=[μ1,μ2,…,μN]TE[v] = [\mu_1, \mu_2, \dots, \mu_N]^T。
v−E[v]=[v1−μ1,v2−μ2,…,vN−μN]Tv - E[v] = [v_1-\mu_1, v_2-\mu_2, \dots, v_N-\mu_N]^T。
CvC_v 是一個 N×NN \times N 的矩陣,其元素 Cv,ijC_{v,ij} 為:
Cv,ij=Cov(vi,vj)=E[(vi−E[vi])(vj−E[vj])]C_{v,ij} = \text{Cov}(v_i, v_j) = E[(v_i - E[v_i])(v_j - E[v_j])]。
特別地,對角線元素是方差:Cv,ii=Var(vi)=E[(vi−E[vi])2]C_{v,ii} = \text{Var}(v_i) = E[(v_i - E[v_i])^2]。

【答案】Cv=E[(v−E[v])(v−E[v])T]C_v = E[(v - E[v])(v - E[v])^T]。

(c) 求雜訊向量 nn 的協方差矩陣

n=[n1n2]n = \begin{bmatrix} n_1 \\ n_2 \end{bmatrix}。
E[n]=[μ1μ2]E[n] = \begin{bmatrix} \mu_1 \\ \mu_2 \end{bmatrix}。
n−E[n]=[n1−μ1n2−μ2]n - E[n] = \begin{bmatrix} n_1 - \mu_1 \\ n_2 - \mu_2 \end{bmatrix}。
(n−E[n])(n−E[n])T=[n1−μ1n2−μ2][n1−μ1n2−μ2](n - E[n])(n - E[n])^T = \begin{bmatrix} n_1 - \mu_1 \\ n_2 - \mu_2 \end{bmatrix} \begin{bmatrix} n_1 - \mu_1 & n_2 - \mu_2 \end{bmatrix}
=[(n1−μ1)2(n1−μ1)(n2−μ2)(n2−μ2)(n1−μ1)(n2−μ2)2]= \begin{bmatrix} (n_1 - \mu_1)^2 & (n_1 - \mu_1)(n_2 - \mu_2) \\ (n_2 - \mu_2)(n_1 - \mu_1) & (n_2 - \mu_2)^2 \end{bmatrix}。

Cn=E[(n−E[n])(n−E[n])T]=[E[(n1−μ1)2]E[(n1−μ1)(n2−μ2)]E[(n2−μ2)(n1−μ1)]E[(n2−μ2)2]]C_n = E[(n - E[n])(n - E[n])^T] = \begin{bmatrix} E[(n_1 - \mu_1)^2] & E[(n_1 - \mu_1)(n_2 - \mu_2)] \\ E[(n_2 - \mu_2)(n_1 - \mu_1)] & E[(n_2 - \mu_2)^2] \end{bmatrix}。

已知 Var(n1)=E[(n1−μ1)2]=σ12Var(n_1) = E[(n_1 - \mu_1)^2] = \sigma_1^2。
已知 Var(n2)=E[(n2−μ2)2]=σ22Var(n_2) = E[(n_2 - \mu_2)^2] = \sigma_2^2。

題目中提到 nn 是「white Gaussian random vector」。這意味著 n1n_1 和 n2n_2 是統計獨立的。
因此,E[(n1−μ1)(n2−μ2)]=E[n1−μ1]E[n2−μ2]E[(n_1 - \mu_1)(n_2 - \mu_2)] = E[n_1 - \mu_1] E[n_2 - \mu_2]。
由於 E[n1−μ1]=E[n1]−μ1=μ1−μ1=0E[n_1 - \mu_1] = E[n_1] - \mu_1 = \mu_1 - \mu_1 = 0。
同理,E[n2−μ2]=0E[n_2 - \mu_2] = 0。
所以,E[(n1−μ1)(n2−μ2)]=0×0=0E[(n_1 - \mu_1)(n_2 - \mu_2)] = 0 \times 0 = 0。
這也表示 n1n_1 和 n2n_2 的協方差為 0。

因此,協方差矩陣為:
Cn=[σ1200σ22]C_n = \begin{bmatrix} \sigma_1^2 & 0 \\ 0 & \sigma_2^2 \end{bmatrix}。

【答案】Cn=[σ1200σ22]C_n = \begin{bmatrix} \sigma_1^2 & 0 \\ 0 & \sigma_2^2 \end{bmatrix}。

(d) 證明 r2r_2 與 ss 無關(獨立)

r=s+nr = s + n。
r1=s1+n1r_1 = s_1 + n_1
r2=s2+n2r_2 = s_2 + n_2

其中 s1=(−1)bdmin/2s_1 = (-1)^{b d_{min}/2}, s2=0s_2 = 0。
n1,n2n_1, n_2 是獨立的高斯隨機變數,且獨立於 ss。

我們需要證明,在已知 r1r_1 的情況下,r2r_2 與 ss 是獨立的。
如果 r2r_2 與 ss 獨立,那麼知道 r1r_1 的值不應該影響 r2r_2 的概率分佈,特別是 r2r_2 的條件期望和條件方差。

考慮 r2r_2 的條件期望 E[r2∣r1]E[r_2 | r_1]。
E[r2∣r1]=E[s2+n2∣r1]E[r_2 | r_1] = E[s_2 + n_2 | r_1]。
由於 s2=0s_2 = 0 是一個常數,且 n2n_2 獨立於 ss(因此也獨立於 s1s_1, s2s_2),所以 n2n_2 也獨立於 r1=s1+n1r_1 = s_1 + n_1。
因此,E[n2∣r1]=E[n2]=μ2E[n_2 | r_1] = E[n_2] = \mu_2。
E[r2∣r1]=E[0∣r1]+E[n2∣r1]=0+μ2=μ2E[r_2 | r_1] = E[0 | r_1] + E[n_2 | r_1] = 0 + \mu_2 = \mu_2。
由於 E[r2∣r1]E[r_2 | r_1] 不依賴於 r1r_1,這暗示了 r2r_2 的分佈可能與 r1r_1 無關(或至少其均值與 r1r_1 無關)。

更嚴格地說,我們需要證明 P(r2≤y∣r1)=P(r2≤y)P(r_2 \le y | r_1) = P(r_2 \le y) 對於所有 yy 成立。
r2=s2+n2=0+n2=n2r_2 = s_2 + n_2 = 0 + n_2 = n_2。
所以,r2r_2 的分佈就是 n2n_2 的分佈。
r1=s1+n1r_1 = s_1 + n_1。
s1=(−1)bdmin/2s_1 = (-1)^{b d_{min}/2},其中 b∈{0,1}b \in \{0,1\}。
如果 dmind_{min} 是偶數,則 dmin/2d_{min}/2 是整數。
如果 dmin/2d_{min}/2 是整數,則 s1s_1 只能取值 11 (當 b=0b=0) 或 −1-1 (當 b=1b=1)。
如果 dmind_{min} 是奇數,則 dmin/2d_{min}/2 不是整數。但題目要求 dmin>0d_{min}>0。
通常在通信系統中,dmind_{min} 指的是漢明距離,是整數。所以 dmin/2d_{min}/2 可能是整數或半整數。
如果 bdmin/2b d_{min}/2 始終是整數,則 s1s_1 只能是 11 或 −1-1。
如果 bdmin/2b d_{min}/2 不是整數(例如 dmin=1d_{min}=1),則 (−1)bdmin/2(-1)^{b d_{min}/2} 的定義需要更仔細。假設 dmind_{min} 是偶數,以簡化問題。
假設 dmind_{min} 是偶數,則 s1∈{1,−1}s_1 \in \{1, -1\}。
r1=s1+n1r_1 = s_1 + n_1。
r2=n2r_2 = n_2。

由於 n1n_1 和 n2n_2 是獨立的,且 s1s_1 是與 n1,n2n_1, n_2 無關的隨機變數(它取決於傳輸的比特 bb),
那麼 r1=s1+n1r_1 = s_1 + n_1 的值取決於 s1s_1 和 n1n_1。
r2=n2r_2 = n_2 的值僅僅取決於 n2n_2。

我們需要證明 P(r2∈A∣r1=r1,value)=P(r2∈A)P(r_2 \in A | r_1 = r_{1,value}) = P(r_2 \in A)。
P(r2∈A∣r1=r1,value)=P(n2∈A∣s1+n1=r1,value)P(r_2 \in A | r_1 = r_{1,value}) = P(n_2 \in A | s_1 + n_1 = r_{1,value})。

🔒

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

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

免費註冊

第 7 題

Let Q(x)Q(x) be the probability that a zero-mean unit-variance Gaussian random variable is no less than xx.
(a) (5%) The complementary error function is defined by
erfc(x)=2π∫x∞e−z2dz\text{erfc}(x) = \frac{2}{\sqrt{\pi}} \int_{x}^{\infty} e^{-z^2} dz
Please express Q(x)Q(x) using erfc(x)\text{erfc}(x). In other words, show how to transform function erfc(x)\text{erfc}(x) to Q(x)Q(x)?
(b) (4%) Please prove your answer in (a).
(c) (12%) For 1>p0>01 > p_0 > 0, find the value of rr which minimizes
p0Q(1−r)+(1−p0)Q(1+r)p_0 Q(1-r) + (1-p_0) Q(1+r)
Please express your answer in p0p_0.

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

這一題的完整詳解

核心觀念

本題考查標準常態分布的尾機率函數 Q(x)Q(x) 與互補誤差函數 erfc⁡(x)\operatorname{erfc}(x) 之間的轉換,以及利用微分求加權尾機率和的最小值。

令 Z∼N(0,1)Z\sim\mathcal{N}(0,1),其機率密度函數為

ϕ(t)=12πe−t2/2.\phi(t)=\frac{1}{\sqrt{2\pi}}e^{-t^2/2}.

因此

Q(x)=Pr⁡(Z≥x)=∫x∞ϕ(t) dt.Q(x)=\Pr(Z\ge x)=\int_x^\infty \phi(t)\,dt.

題目定義的互補誤差函數為

erfc⁡(x)=2π∫x∞e−z2 dz.\operatorname{erfc}(x)=\frac{2}{\sqrt{\pi}}\int_x^\infty e^{-z^2}\,dz.

(a)Q(x)Q(x) 與 erfc⁡(x)\operatorname{erfc}(x) 的關係

將 Q(x)Q(x) 中常態密度的指數 e−t2/2e^{-t^2/2},透過變數縮放轉成 e−z2e^{-z^2},可得

Q(x)=12erfc⁡(x2).Q(x)=\frac{1}{2}\operatorname{erfc}\left(\frac{x}{\sqrt{2}}\right).

(b)證明

從 Q(x)Q(x) 的積分定義出發,令

z=t2,t=2z,dt=2 dz.z=\frac{t}{\sqrt{2}}, \qquad t=\sqrt{2}z, \qquad dt=\sqrt{2}\,dz.

當 t=xt=x 時,z=x/2z=x/\sqrt{2};當 t→∞t\to\infty 時,z→∞z\to\infty。代入積分:

🔒

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

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

免費註冊

其他考古題