113 年 國立中山大學通訊工程研究所甲組《通訊理論》

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

第 1 題20 分

  1. (20%) We are given the complex baseband signal x2(t)=sinc(t−1)+j2sinc(t)x_2(t) = \text{sinc}(t - 1) + j2\text{sinc}(t).
    (A) (5%) Calculate the real and imaginary parts of the Fourier transform X2(f)X_2(f). To express the transforms, please use the function rect(t)\text{rect}(t), which is defined as a rectangle of unit height and spanning the interval [-1/2,1/2].
    (B) (5%) Plot the real and imaginary parts of X2(f)X_2(f).
    (C) (10%) Plot the real and imaginary parts of the Fourier transform of the bandpass signal obtained by upconverting x2(t)x_2(t) to the carrier frequency of 10 Hz.

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

這一題的完整詳解

本題主要考驗學生對複數基頻訊號的傅立葉轉換及其應用,包含基本訊號轉換、頻域分析與升頻轉換的理解。

首先,我們需要知道 sinc(t)\text{sinc}(t) 的傅立葉轉換是 rect(f)\text{rect}(f),即
F{sinc(t)}=∫−∞∞sinc(t)e−j2πftdt=rect(f)\mathcal{F}\{\text{sinc}(t)\} = \int_{-\infty}^{\infty} \text{sinc}(t) e^{-j2\pi ft} dt = \text{rect}(f)
其中 rect(f)\text{rect}(f) 的定義為
rect(f)={1if ∣f∣≤1/20if ∣f∣>1/2\text{rect}(f) = \begin{cases} 1 & \text{if } |f| \le 1/2 \\ 0 & \text{if } |f| > 1/2 \end{cases}
而題目給定的 rect(t)\text{rect}(t) 定義為單位高度且範圍為 [−1/2,1/2][-1/2, 1/2],這與標準的 rect(f)\text{rect}(f) 定義略有不同,標準的 rect(f)\text{rect}(f) 是單位高度且頻寬為 1。題目要求使用 rect(t)\text{rect}(t) 來表示傅立葉轉換,這裡我們假設題目是指標準的 rect(f)\text{rect}(f) 函數,只是其參數寫成 tt。若題目嚴格要求使用 rect(t)\text{rect}(t) 來表示,則 rect(t)\text{rect}(t) 的傅立葉轉換為 sinc(f)\text{sinc}(f)。然而,根據題目提供的訊號形式,更符合工程上的習慣是 sinc(t)\text{sinc}(t) 轉換為 rect(f)\text{rect}(f)。我們將採用後者。

若題目指的 rect(t)\text{rect}(t) 是單位高度且範圍為 [−1/2,1/2][-1/2, 1/2],則其傅立葉轉換為 sinc(f)\text{sinc}(f)。
F{rect(t)}=∫−1/21/21⋅e−j2πftdt=[e−j2πft−j2πf]−1/21/2=e−jπf−ejπf−j2πf=−2jsin⁡(πf)−j2πf=sin⁡(πf)πf=sinc(f)\mathcal{F}\{\text{rect}(t)\} = \int_{-1/2}^{1/2} 1 \cdot e^{-j2\pi ft} dt = \left[ \frac{e^{-j2\pi ft}}{-j2\pi f} \right]_{-1/2}^{1/2} = \frac{e^{-j\pi f} - e^{j\pi f}}{-j2\pi f} = \frac{-2j \sin(\pi f)}{-j2\pi f} = \frac{\sin(\pi f)}{\pi f} = \text{sinc}(f)
這與標準的 sinc(t)\text{sinc}(t) 轉換為 rect(f)\text{rect}(f) 是對偶的。

題目中的 sinc(t)\text{sinc}(t) 通常定義為 sin⁡(πt)πt\frac{\sin(\pi t)}{\pi t}。其傅立葉轉換為 rect(f)\text{rect}(f)。
F{sin⁡(πt)πt}=rect(f)\mathcal{F}\left\{\frac{\sin(\pi t)}{\pi t}\right\} = \text{rect}(f)
若 sinc(t)\text{sinc}(t) 定義為 sin⁡(t)t\frac{\sin(t)}{t},則其傅立葉轉換為 π⋅rect(f2π)\pi \cdot \text{rect}(\frac{f}{2\pi})。
根據題目敘述,使用 rect(t)\text{rect}(t) 定義為單位高度且跨度為 [−1/2,1/2][-1/2, 1/2],這暗示了 sinc(t)\text{sinc}(t) 的頻譜是 rect(f)\text{rect}(f)。我們假設 sinc(t)=sin⁡(πt)πt\text{sinc}(t) = \frac{\sin(\pi t)}{\pi t}。

(A) 計算 X2(f)X_2(f) 的實部與虛部。
我們有 x2(t)=sinc(t−1)+j2sinc(t)x_2(t) = \text{sinc}(t - 1) + j2\text{sinc}(t)。
利用傅立葉轉換的時域位移性質:F{f(t−t0)}=e−j2πft0F(f)\mathcal{F}\{f(t-t_0)\} = e^{-j2\pi ft_0} F(f)。
則 F{sinc(t−1)}=e−j2πf(1)F{sinc(t)}=e−j2πfrect(f)\mathcal{F}\{\text{sinc}(t - 1)\} = e^{-j2\pi f(1)} \mathcal{F}\{\text{sinc}(t)\} = e^{-j2\pi f} \text{rect}(f)。
而 F{2sinc(t)}=2F{sinc(t)}=2rect(f)\mathcal{F}\{2\text{sinc}(t)\} = 2\mathcal{F}\{\text{sinc}(t)\} = 2\text{rect}(f)。

所以,X2(f)=F{x2(t)}=F{sinc(t−1)}+jF{2sinc(t)}X_2(f) = \mathcal{F}\{x_2(t)\} = \mathcal{F}\{\text{sinc}(t - 1)\} + j\mathcal{F}\{2\text{sinc}(t)\}
X2(f)=e−j2πfrect(f)+j2rect(f)X_2(f) = e^{-j2\pi f} \text{rect}(f) + j2\text{rect}(f)
X2(f)=(e−j2πf+j2)rect(f)X_2(f) = (e^{-j2\pi f} + j2) \text{rect}(f)

其中 rect(f)={1if ∣f∣≤1/20if ∣f∣>1/2\text{rect}(f) = \begin{cases} 1 & \text{if } |f| \le 1/2 \\ 0 & \text{if } |f| > 1/2 \end{cases}。
e−j2πf=cos⁡(2πf)−jsin⁡(2πf)e^{-j2\pi f} = \cos(2\pi f) - j\sin(2\pi f)。

所以,X2(f)=(cos⁡(2πf)−jsin⁡(2πf)+j2)rect(f)X_2(f) = (\cos(2\pi f) - j\sin(2\pi f) + j2) \text{rect}(f)
X2(f)=(cos⁡(2πf)+j(2−sin⁡(2πf)))rect(f)X_2(f) = (\cos(2\pi f) + j(2 - \sin(2\pi f))) \text{rect}(f)

當 ∣f∣≤1/2|f| \le 1/2 時,rect(f)=1\text{rect}(f) = 1。
X2(f)=cos⁡(2πf)+j(2−sin⁡(2πf))X_2(f) = \cos(2\pi f) + j(2 - \sin(2\pi f)),當 ∣f∣≤1/2|f| \le 1/2。
X2(f)=0X_2(f) = 0,當 ∣f∣>1/2|f| > 1/2。

X2(f)X_2(f) 的實部為 Re{X2(f)}=cos⁡(2πf)⋅rect(f)\text{Re}\{X_2(f)\} = \cos(2\pi f) \cdot \text{rect}(f)。
X2(f)X_2(f) 的虛部為 Im{X2(f)}=(2−sin⁡(2πf))⋅rect(f)\text{Im}\{X_2(f)\} = (2 - \sin(2\pi f)) \cdot \text{rect}(f)。

【答案】
X2(f)X_2(f) 的實部為 cos⁡(2πf)⋅rect(f)\cos(2\pi f) \cdot \text{rect}(f)。
X2(f)X_2(f) 的虛部為 (2−sin⁡(2πf))⋅rect(f)(2 - \sin(2\pi f)) \cdot \text{rect}(f)。
其中 rect(f)=1\text{rect}(f) = 1 for ∣f∣≤1/2|f| \le 1/2 and 00 otherwise.

(B) 繪製 X2(f)X_2(f) 的實部與虛部。
實部:Re{X2(f)}=cos⁡(2πf)\text{Re}\{X_2(f)\} = \cos(2\pi f),當 ∣f∣≤1/2|f| \le 1/2。
虛部:Im{X2(f)}=2−sin⁡(2πf)\text{Im}\{X_2(f)\} = 2 - \sin(2\pi f),當 ∣f∣≤1/2|f| \le 1/2。

繪製圖形:
頻率範圍是 [−1/2,1/2][-1/2, 1/2]。
對於實部 cos⁡(2πf)\cos(2\pi f):
當 f=0f = 0, cos⁡(0)=1\cos(0) = 1.
當 f=1/2f = 1/2, cos⁡(π)=−1\cos(\pi) = -1.
當 f=−1/2f = -1/2, cos⁡(−π)=−1\cos(-\pi) = -1.
這是一個餘弦函數在 [−1/2,1/2][-1/2, 1/2] 區間的圖形,從 f=−1/2f=-1/2 的 −1-1 開始,上升到 f=0f=0 的 11,再下降到 f=1/2f=1/2 的 −1-1。

對於虛部 2−sin⁡(2πf)2 - \sin(2\pi f):
當 f=0f = 0, 2−sin⁡(0)=22 - \sin(0) = 2.
當 f=1/2f = 1/2, 2−sin⁡(π)=22 - \sin(\pi) = 2.
當 f=−1/2f = -1/2, 2−sin⁡(−π)=22 - \sin(-\pi) = 2.
這是一個常數函數 22 在 [−1/2,1/2][-1/2, 1/2] 區間。

【答案】
繪製圖形如下:
實部 Re{X2(f)}\text{Re}\{X_2(f)\}:在 f∈[−1/2,1/2]f \in [-1/2, 1/2] 區間,為 cos⁡(2πf)\cos(2\pi f) 的曲線,峰值為 1 在 f=0f=0,兩端點為 -1 在 f=±1/2f=\pm 1/2。在區間外為 0。
虛部 Im{X2(f)}\text{Im}\{X_2(f)\}:在 f∈[−1/2,1/2]f \in [-1/2, 1/2] 區間,為常數 2。在區間外為 0。

🔒

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

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

免費註冊

第 2 題20 分

  1. (20%) Let x(t)=m(t)+cos⁡(ωct)x(t) = m(t) + \cos(\omega_c t). Let WW be the bandwidth of m(t)m(t). Assume that the average value of m(t)m(t) is zero and that the maximum value of ∣m(t)∣|m(t)| is MM. Also assume that the square-law device is defined by y(t)=4x(t)+2x2(t)y(t) = 4x(t) + 2x^2(t).
    (A) (5%) Write the equation for y(t)y(t).
    (B) (10%) Describe the filter with input signal y(t)y(t) that produces an AM signal for g(t)g(t), where g(t)g(t) represents the output of the filter.
    (C) (5%) Specify the requirement of MM to ensure no distortion when using envelope demodulation.

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

這一題的完整詳解

核心觀念

本題考查平方律元件產生 AM(Amplitude Modulation,振幅調變)訊號的原理。

給定輸入

x(t)=m(t)+cos⁡(ωct)x(t)=m(t)+\cos(\omega_c t)

其中 m(t)m(t) 為訊息訊號,頻寬為 WW;平方律元件輸出為

y(t)=4x(t)+2x2(t)y(t)=4x(t)+2x^2(t)

平方律元件會產生:

  • 基頻訊息成分;
  • 載波成分;
  • 訊息與載波相乘所形成的 ωc\omega_c 附近側帶;
  • 二倍載波頻率 2ωc2\omega_c 成分。

利用帶通濾波器選出 ωc\omega_c 附近的頻率成分,即可得到 AM 訊號。


解題方法

將 x(t)=m(t)+cos⁡(ωct)x(t)=m(t)+\cos(\omega_c t) 代入平方律元件:

y(t)=4[m(t)+cos⁡(ωct)]+2[m(t)+cos⁡(ωct)]2=4m(t)+4cos⁡(ωct)+2m2(t)+4m(t)cos⁡(ωct)+2cos⁡2(ωct)\begin{aligned} y(t) &=4\left[m(t)+\cos(\omega_c t)\right] +2\left[m(t)+\cos(\omega_c t)\right]^2\\ &=4m(t)+4\cos(\omega_c t) +2m^2(t)+4m(t)\cos(\omega_c t) +2\cos^2(\omega_c t) \end{aligned}

再利用恆等式

2cos⁡2(ωct)=1+cos⁡(2ωct)2\cos^2(\omega_c t)=1+\cos(2\omega_c t)

即可整理出各頻率區域的成分。


(A) 求 y(t)y(t)

因此

y(t)=1+4m(t)+2m2(t)+4[1+m(t)]cos⁡(ωct)+cos⁡(2ωct)\boxed{ y(t)=1+4m(t)+2m^2(t) +4\left[1+m(t)\right]\cos(\omega_c t) +\cos(2\omega_c t) }

也可寫成分群形式:

y(t)=1+4m(t)+2m2(t)⏟基頻及直流成分+4[1+m(t)]cos⁡(ωct)⏟ωc 附近的 AM 成分+cos⁡(2ωct)⏟2ωc 成分\boxed{ y(t)= \underbrace{1+4m(t)+2m^2(t)}_{\text{基頻及直流成分}} + \underbrace{4[1+m(t)]\cos(\omega_c t)}_{\omega_c\text{ 附近的 AM 成分}} + \underbrace{\cos(2\omega_c t)}_{2\omega_c\text{ 成分}} }

其中 4m(t)cos⁡(ωct)4m(t)\cos(\omega_c t) 是訊息搬移至載波頻率附近所形成的雙邊帶成分,而 4cos⁡(ωct)4\cos(\omega_c t) 則是載波成分。


(B) 濾波器設計與輸出 g(t)g(t)

1. 應保留的頻率成分

由於 m(t)m(t) 的頻寬為 WW,乘上 cos⁡(ωct)\cos(\omega_c t) 後,其頻譜會搬移至 ωc\omega_c 附近,範圍為

ωc−W≤∣ω∣≤ωc+W\omega_c-W\leq |\omega|\leq \omega_c+W

因此應使用以 ωc\omega_c 為中心頻率、通帶寬度為 2W2W 的理想帶通濾波器。

其正頻率通帶可表示為

ωc−W≤ω≤ωc+W\omega_c-W\leq \omega\leq \omega_c+W

並應同時保留對稱的負頻率通帶。濾波器增益取 11,且抑制直流、基頻以及 2ωc2\omega_c 附近的成分。

2. 濾波後輸出

通過帶通濾波器的部分為

🔒

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

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

免費註冊

第 3 題10 分

  1. (10%) A transmitter uses a carrier frequency of 1000 Hz, with the unmodulated carrier represented as Accos⁡(2πfct)A_c \cos(2\pi f_c t). Determine both the phase and frequency deviation for each of the following transmitter outputs:
    (A) (5%) x(t)=cos⁡[2π(1000)t+40sin⁡(5t2)]x(t) = \cos[2\pi(1000)t + 40 \sin(5t^2)]
    (B) (5%) x(t)=cos⁡[2π(600)t]x(t) = \cos[2\pi(600)t]

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

這一題的完整詳解

核心觀念

將輸出表示為

x(t)=Accos⁡[2πfct+ϕd(t)]x(t)=A_c\cos\left[2\pi f_ct+\phi_d(t)\right]

其中:

  • ϕd(t)\phi_d(t):相對於未調變載波的瞬時相位偏移,單位為 rad。
  • 瞬時頻率為
fi(t)=fc+Δf(t)f_i(t)=f_c+\Delta f(t)
  • 瞬時頻率偏移為
Δf(t)=12πdϕd(t)dt\Delta f(t)=\frac{1}{2\pi}\frac{d\phi_d(t)}{dt}

本題載波頻率為 fc=1000f_c=1000 Hz,因此須將各輸出相位與 2π(1000)t2\pi(1000)t 比較。


(A) x(t)=cos⁡[2π(1000)t+40sin⁡(5t2)]x(t)=\cos[2\pi(1000)t+40\sin(5t^2)]

相位偏移

與標準形式比較:

ϕd(t)=40sin⁡(5t2)\phi_d(t)=40\sin(5t^2)

因此瞬時相位偏移為

ϕd(t)=40sin⁡(5t2) rad\boxed{\phi_d(t)=40\sin(5t^2)\ \text{rad}}

由於 −1≤sin⁡(5t2)≤1-1\leq \sin(5t^2)\leq 1,其最大相位偏移量為

Δϕmax⁡=40 rad\boxed{\Delta\phi_{\max}=40\ \text{rad}}

頻率偏移

對相位偏移微分:

dϕd(t)dt=40ddtsin⁡(5t2)\frac{d\phi_d(t)}{dt} =40\frac{d}{dt}\sin(5t^2) =40cos⁡(5t2)⋅10t=400tcos⁡(5t2)=40\cos(5t^2)\cdot 10t =400t\cos(5t^2)

因此瞬時頻率偏移為

Δf(t)=12π⋅400tcos⁡(5t2)\Delta f(t) =\frac{1}{2\pi}\cdot 400t\cos(5t^2) Δf(t)=200πtcos⁡(5t2) Hz\boxed{\Delta f(t)=\frac{200}{\pi}t\cos(5t^2)\ \text{Hz}}

瞬時頻率則為

fi(t)=1000+200πtcos⁡(5t2) Hz\boxed{f_i(t)=1000+\frac{200}{\pi}t\cos(5t^2)\ \text{Hz}}

由於題目未限定時間範圍,tcos⁡(5t2)t\cos(5t^2) 的大小沒有有限上限,所以頻率偏移的最大值不存在有限數值;其幅度會隨 ∣t∣|t| 增大而增大。


(B) x(t)=cos⁡[2π(600)t]x(t)=\cos[2\pi(600)t]

相位偏移

以 10001000 Hz 載波為基準,改寫輸出相位:

🔒

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

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

免費註冊

第 4 題10 分

  1. (10%) Consider the following three signals:
    [Image of three signals S1(t), S2(t), S3(t)]
    (A) (6%) Use Gram-Schmidt procedure to find the set of basis functions from the three signals and determine the dimensionality of the set.
    (B) (2%) Find the signal-space representation of the three signals based on the basis functions obtained in (a). (Represent the signals in terms of vectors)
    (C) (2%) Determine the minimum distance between any pair of waveforms.
🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 2 頁原卷第 3 頁

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

這一題的完整詳解

核心觀念

本題考 Gram–Schmidt 正交化、訊號空間表示,以及正交基底下的歐氏距離。訊號內積與能量定義為

⟨x,y⟩=∫−∞∞x(t)y(t) dt,∥x∥=⟨x,x⟩.\langle x,y\rangle=\int_{-\infty}^{\infty}x(t)y(t)\,dt, \qquad \|x\|=\sqrt{\langle x,x\rangle}.

由圖可讀出三個訊號皆只在 0≤t≤40\le t\le 4 非零:

s1(t)=1,s2(t)=54−34t,s3(t)=−74+34t.s_1(t)=1,\qquad s_2(t)=\frac54-\frac34t,\qquad s_3(t)=-\frac74+\frac34t.

(A) Gram–Schmidt 正交化與維度

先以 s1(t)s_1(t) 建立第一個正交單位基底。其能量為

∥s1∥2=∫041 dt=4,\|s_1\|^2=\int_0^4 1\,dt=4,

因此

ϕ1(t)=s1(t)∥s1∥={12,0≤t≤4,0,其他.\phi_1(t)=\frac{s_1(t)}{\|s_1\|} = \begin{cases} \frac12,&0\le t\le4,\\ 0,&\text{其他}. \end{cases}

接著將 s2(t)s_2(t) 對 ϕ1(t)\phi_1(t) 投影:

⟨s2,ϕ1⟩=∫04(54−34t)12 dt=−12.\langle s_2,\phi_1\rangle =\int_0^4\left(\frac54-\frac34t\right)\frac12\,dt =-\frac12.

扣除投影後的殘差為

g2(t)=s2(t)−⟨s2,ϕ1⟩ϕ1(t)=s2(t)+12ϕ1(t)=32−34t,0≤t≤4.g_2(t)=s_2(t)-\langle s_2,\phi_1\rangle\phi_1(t) =s_2(t)+\frac12\phi_1(t) =\frac32-\frac34t,\quad 0\le t\le4.

其能量為

∥g2∥2=∫04(32−34t)2dt=3,\|g_2\|^2 =\int_0^4\left(\frac32-\frac34t\right)^2dt =3,

所以第二個正交單位基底為

ϕ2(t)=g2(t)3={3/2−3t/43,0≤t≤4,0,其他.\phi_2(t)=\frac{g_2(t)}{\sqrt3} = \begin{cases} \frac{3/2-3t/4}{\sqrt3},&0\le t\le4,\\ 0,&\text{其他}. \end{cases}

由圖中的兩條斜線可得

s2(t)+s3(t)=−12s1(t),s_2(t)+s_3(t)=-\frac12s_1(t),

故 s3s_3 不會產生新的獨立方向。此訊號集合的維度為 22,正交單位基底為 {ϕ1,ϕ2}\{\phi_1,\phi_2\}。

(B) 訊號空間向量表示

🔒

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

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

免費註冊

第 5 題20 分

  1. (20%) Let x(t)x(t) denote a real valued WSS random process with an autocorrelation function Rx(τ)R_x(\tau) and y(t)=x(t)cos⁡(2πfct+θ)y(t) = x(t) \cos(2\pi f_c t + \theta), θ∼U(0,2π)\theta \sim U(0,2\pi).
    (A) (5%) Find Ry(τ)R_y(\tau).
    (B) (5%) y(t)y(t) is passed through a low-pass filter with a frequency-domain response H(f)H(f)
    H(f)={1if ∣f∣<fc0if ∣f∣>fcH(f) = \begin{cases} 1 & \text{if } |f| < f_c \\ 0 & \text{if } |f| > f_c \end{cases}
    and y~(t)\tilde{y}(t) denotes the filter output. Find E[y~(t)2]E[\tilde{y}(t)^2].
    (C) (10%) If we let π4=π4\frac{\pi}{4} = \frac{\pi}{4} and r(t)=s(t)cos⁡(2πfct)+y(t)r(t) = s(t) \cos(2\pi f_c t) + y(t), please show how to demodulate s(t)s(t) based on r(t)r(t) in detail.

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

這一題的完整詳解

核心觀念

  1. 廣義平穩隨機程序(WSS Random Process)之載波調變:
    • 設實數隨機程序 x(t)x(t) 為 WSS,其與均勻分佈相位常態隨機變數 θ∼U(0,2π)\theta \sim U(0, 2\pi) 彼此獨立。
    • 經載波調變 y(t)=x(t)cos⁡(2πfct+θ)y(t) = x(t)\cos(2\pi f_c t + \theta) 後,y(t)y(t) 亦為 WSS,其自相關函數為:
      Ry(τ)=12Rx(τ)cos⁡(2πfcτ)R_y(\tau) = \frac{1}{2} R_x(\tau) \cos(2\pi f_c \tau)
  2. 隨機程序通過 LTI 系統與平均功率計算:
    • 根據 Wiener-Khinchin 定理,功率頻譜密度(PSD)Sy(f)=F{Ry(τ)}=14[Sx(f−fc)+Sx(f+fc)]S_y(f) = \mathcal{F}\{R_y(\tau)\} = \frac{1}{4}[S_x(f-f_c) + S_x(f+f_c)]。
    • 輸出程序 y~(t)\tilde{y}(t) 的平均功率為其 PSD 在低通濾波器(LPF)通帶內的積分:
      E[y~(t)2]=Ry~(0)=∫−∞∞Sy(f)∣H(f)∣2dfE[\tilde{y}(t)^2] = R_{\tilde{y}}(0) = \int_{-\infty}^{\infty} S_y(f)|H(f)|^2 df
  3. 同頻干擾下的正交載波同步解調(Orthogonal Coherent Demodulation):
    • 當接收訊號包含目標訊號 s(t)cos⁡(2πfct)s(t)\cos(2\pi f_c t) 與具有固定相位差的干擾訊號 x(t)cos⁡(2πfct+θ)x(t)\cos(2\pi f_c t + \theta) 時,可利用正交載波投影或 I/Q 雙通道解調消除法,在基頻將干擾項完全抵消並精確還原 s(t)s(t)。

各小題詳解

(A) 求 Ry(τ)R_y(\tau)

根據自相關函數定義:
Ry(t+τ,t)=E[y(t+τ)y(t)]=E[x(t+τ)cos⁡(2πfc(t+τ)+θ)⋅x(t)cos⁡(2πfct+θ)]R_y(t+\tau, t) = E[y(t+\tau)y(t)] = E\big[ x(t+\tau)\cos(2\pi f_c(t+\tau) + \theta) \cdot x(t)\cos(2\pi f_c t + \theta) \big]

因為 x(t)x(t) 與 θ\theta 彼此獨立,期望值可分開計算:
Ry(t+τ,t)=E[x(t+τ)x(t)]⋅E[cos⁡(2πfc(t+τ)+θ)cos⁡(2πfct+θ)]R_y(t+\tau, t) = E[x(t+\tau)x(t)] \cdot E\big[ \cos(2\pi f_c(t+\tau) + \theta)\cos(2\pi f_c t + \theta) \big]

已知 x(t)x(t) 為 WSS,故 E[x(t+τ)x(t)]=Rx(τ)E[x(t+\tau)x(t)] = R_x(\tau)。利用三角函數積化和差公式 cos⁡Acos⁡B=12[cos⁡(A−B)+cos⁡(A+B)]\cos A \cos B = \frac{1}{2}[\cos(A-B) + \cos(A+B)]:
cos⁡(2πfc(t+τ)+θ)cos⁡(2πfct+θ)=12cos⁡(2πfcτ)+12cos⁡(4πfct+2πfcτ+2θ)\cos(2\pi f_c(t+\tau) + \theta)\cos(2\pi f_c t + \theta) = \frac{1}{2}\cos(2\pi f_c \tau) + \frac{1}{2}\cos(4\pi f_c t + 2\pi f_c \tau + 2\theta)

對第二項取關於 θ∼U(0,2π)\theta \sim U(0, 2\pi) 的期望值:
E[cos⁡(4πfct+2πfcτ+2θ)]=12π∫02πcos⁡(4πfct+2πfcτ+2θ)dθ=0E[\cos(4\pi f_c t + 2\pi f_c \tau + 2\theta)] = \frac{1}{2\pi}\int_{0}^{2\pi} \cos(4\pi f_c t + 2\pi f_c \tau + 2\theta) d\theta = 0

因此,Ry(t+τ,t)R_y(t+\tau, t) 僅與時間差 τ\tau 有關:
Ry(τ)=12Rx(τ)cos⁡(2πfcτ)R_y(\tau) = \frac{1}{2} R_x(\tau) \cos(2\pi f_c \tau)


(B) 求 E[y~(t)2]E[\tilde{y}(t)^2]

由 (A) 的結果,取傅立葉轉換求 y(t)y(t) 的功率頻譜密度(PSD):
Sy(f)=F{Ry(τ)}=14[Sx(f−fc)+Sx(f+fc)]S_y(f) = \mathcal{F}\{R_y(\tau)\} = \frac{1}{4} \big[ S_x(f - f_c) + S_x(f + f_c) \big]

低通濾波器之頻率響應為:
H(f)={1,∣f∣<fc0,∣f∣>fcH(f) = \begin{cases} 1, & |f| < f_c \\ 0, & |f| > f_c \end{cases}

輸出訊號 y~(t)\tilde{y}(t) 的平均功率即為輸出 PSD 在全頻域的積分:
E[y~(t)2]=∫−∞∞Sy~(f)df=∫−∞∞Sy(f)∣H(f)∣2df=∫−fcfcSy(f)dfE[\tilde{y}(t)^2] = \int_{-\infty}^{\infty} S_{\tilde{y}}(f) df = \int_{-\infty}^{\infty} S_y(f) |H(f)|^2 df = \int_{-f_c}^{f_c} S_y(f) df

代入 Sy(f)S_y(f):
E[y~(t)2]=14[∫−fcfcSx(f−fc)df+∫−fcfcSx(f+fc)df]E[\tilde{y}(t)^2] = \frac{1}{4} \left[ \int_{-f_c}^{f_c} S_x(f - f_c) df + \int_{-f_c}^{f_c} S_x(f + f_c) df \right]

進行變數代換:

  • 第一項令 u=f−fc  ⟹  ∫−fcfcSx(f−fc)df=∫−2fc0Sx(u)duu = f - f_c \implies \int_{-f_c}^{f_c} S_x(f - f_c) df = \int_{-2f_c}^{0} S_x(u) du
  • 第二項令 v=f+fc  ⟹  ∫−fcfcSx(f+fc)df=∫02fcSx(v)dvv = f + f_c \implies \int_{-f_c}^{f_c} S_x(f + f_c) df = \int_{0}^{2f_c} S_x(v) dv

將兩積分區間合併:
E[y~(t)2]=14∫−2fc2fcSx(f)dfE[\tilde{y}(t)^2] = \frac{1}{4} \int_{-2f_c}^{2f_c} S_x(f) df

標準通訊系統條件:
一般情況下,隨機程序 x(t)x(t) 為基頻訊號,其最高頻寬 W≤2fcW \le 2f_c(通常滿足 W<fc≪2fcW < f_c \ll 2f_c),此時 Sx(f)S_x(f) 的所有功率皆落在 [−2fc,2fc][-2f_c, 2f_c] 內:
∫−2fc2fcSx(f)df=∫−∞∞Sx(f)df=Rx(0)=E[x(t)2]\int_{-2f_c}^{2f_c} S_x(f) df = \int_{-\infty}^{\infty} S_x(f) df = R_x(0) = E[x(t)^2]

因此:
E[y~(t)2]=14Rx(0)=14E[x(t)2]E[\tilde{y}(t)^2] = \frac{1}{4} R_x(0) = \frac{1}{4} E[x(t)^2]


(C) 詳細說明如何由 r(t)r(t) 解調出 s(t)s(t)

題目勘誤說明:題幹「If we let π4=π4\frac{\pi}{4} = \frac{\pi}{4}」為原考卷印刷筆誤,實為指定相位常數 θ=π4\theta = \frac{\pi}{4}。

🔒

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

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

免費註冊

第 6 題10 分

  1. (10%) Let x(t)=g(t)+w(t)x(t) = g(t) + w(t), 0≤t≤T0 \le t \le T, be the received noisy signal, where g(t)g(t) denotes the transmitted pulse that represents a binary symbol 0 or 1 and w(t)w(t) denotes an additive white noise process with zero mean and power spectral density (PSD) N02\frac{N_0}{2}. Since the filter is linear, the result output can be express as y(t)=go(t)+n(t)y(t) = g_o(t) + n(t), where go(t)g_o(t) denotes the response to g(t)g(t) and n(t)n(t) denotes the response to w(t)w(t). We know that the peak pulse signal to noise ratio of the match filter is
    η=∣go(T)∣2E[n2(t)]\eta = \frac{|g_o(T)|^2}{E[n^2(t)]}
    Please show that η≤1N0∫−∞∞∣G(f)∣2df\eta \le \frac{1}{N_0} \int_{-\infty}^{\infty} |G(f)|^2 df.
🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 3 頁

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

這一題的完整詳解

核心觀念

本題考查:

  • 線性時不變系統的輸出表示;
  • 白雜訊通過濾波器後的輸出功率;
  • Cauchy–Schwarz 不等式;
  • Parseval 定理;
  • 匹配濾波器在取樣時刻達到最大訊雜比。

題目圖中原式為

η≤2N0∫−∞∞∣G(f)∣2 df.\eta\leq \frac{2}{N_0}\int_{-\infty}^{\infty}|G(f)|^2\,df.

解題方法

令濾波器脈衝響應為 h(t)h(t),並在 t=Tt=T 取樣。

由線性時不變系統,訊號部分輸出為

go(T)=∫−∞∞g(τ)h(T−τ) dτ.g_o(T)=\int_{-\infty}^{\infty}g(\tau)h(T-\tau)\,d\tau.

雜訊 w(t)w(t) 的雙邊功率譜密度為 N0/2N_0/2,因此通過濾波器後,在任意時刻的雜訊變異數為

E[n2(t)]=N02∫−∞∞∣h(t)∣2 dt.E[n^2(t)] =\frac{N_0}{2}\int_{-\infty}^{\infty}|h(t)|^2\,dt.

令

hT(τ)=h(T−τ),h_T(\tau)=h(T-\tau),

則時間反轉與平移不會改變能量,因此

∫−∞∞∣hT(τ)∣2 dτ=∫−∞∞∣h(t)∣2 dt.\int_{-\infty}^{\infty}|h_T(\tau)|^2\,d\tau = \int_{-\infty}^{\infty}|h(t)|^2\,dt.

利用 Cauchy–Schwarz 不等式,

∣go(T)∣2=∣∫−∞∞g(τ)hT(τ) dτ∣2|g_o(T)|^2 = \left|\int_{-\infty}^{\infty}g(\tau)h_T(\tau)\,d\tau\right|^2 ≤(∫−∞∞∣g(τ)∣2 dτ)(∫−∞∞∣hT(τ)∣2 dτ).\leq \left(\int_{-\infty}^{\infty}|g(\tau)|^2\,d\tau\right) \left(\int_{-\infty}^{\infty}|h_T(\tau)|^2\,d\tau\right).

因此

∣go(T)∣2≤(∫−∞∞∣g(t)∣2 dt)(∫−∞∞∣h(t)∣2 dt).|g_o(T)|^2 \leq \left(\int_{-\infty}^{\infty}|g(t)|^2\,dt\right) \left(\int_{-\infty}^{\infty}|h(t)|^2\,dt\right).

代入峰值訊雜比

🔒

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

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

免費註冊

第 7 題10 分

  1. (10%) Consider two discrete random variables X and Y with the joint distribution:
    P(x, y) | X = -1 | X = 0 | X = 1
    ------- | ------ | ------ | ------
    Y = 2 | 0.1 | 0.15 | 0.15
    Y = 4 | 0.05 | 0.2 | 0.15
    Y = 6 | 0.05 | 0.05 | 0.1

(A) (2%) Find the entropy H(X, Y).
(B) (2%) Find the entropy H(X).
(C) (4%) Find the entropy H(X|Y).
(D) (2%) Find the mutual information I(X; Y).

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

這一題的完整詳解

本題主要考驗學生對資訊理論基本概念的理解,包含聯合熵、邊際熵、條件熵和互資訊的計算。

首先,我們需要根據給定的聯合機率分佈表,計算出所有邊際機率和條件機率。
聯合機率分佈表:

P(x, y)X = -1X = 0X = 1P(Y=y)
Y = 20.10.150.150.4
Y = 40.050.20.150.4
Y = 60.050.050.10.2
P(X=x)0.20.40.41.0 (總和)

計算邊際機率:
P(Y=2) = 0.1 + 0.15 + 0.15 = 0.4
P(Y=4) = 0.05 + 0.2 + 0.15 = 0.4
P(Y=6) = 0.05 + 0.05 + 0.1 = 0.2
總和 = 0.4 + 0.4 + 0.2 = 1.0

P(X=-1) = 0.1 + 0.05 + 0.05 = 0.2
P(X=0) = 0.15 + 0.2 + 0.05 = 0.4
P(X=1) = 0.15 + 0.15 + 0.1 = 0.4
總和 = 0.2 + 0.4 + 0.4 = 1.0

(A) 計算聯合熵 H(X, Y)。
聯合熵的定義為 H(X,Y)=−∑x,yP(x,y)log⁡2P(x,y)H(X, Y) = -\sum_{x,y} P(x,y) \log_2 P(x,y)。
我們需要計算所有非零的 P(x,y)P(x,y) 項。
H(X,Y)=−[P(−1,2)log⁡2P(−1,2)+P(0,2)log⁡2P(0,2)+P(1,2)log⁡2P(1,2)+H(X, Y) = - [ P(-1,2)\log_2 P(-1,2) + P(0,2)\log_2 P(0,2) + P(1,2)\log_2 P(1,2) +
P(−1,4)log⁡2P(−1,4)+P(0,4)log⁡2P(0,4)+P(1,4)log⁡2P(1,4)+P(-1,4)\log_2 P(-1,4) + P(0,4)\log_2 P(0,4) + P(1,4)\log_2 P(1,4) +
P(−1,6)log⁡2P(−1,6)+P(0,6)log⁡2P(0,6)+P(1,6)log⁡2P(1,6)]P(-1,6)\log_2 P(-1,6) + P(0,6)\log_2 P(0,6) + P(1,6)\log_2 P(1,6) ]

H(X,Y)=−[0.1log⁡20.1+0.15log⁡20.15+0.15log⁡20.15+H(X, Y) = - [ 0.1 \log_2 0.1 + 0.15 \log_2 0.15 + 0.15 \log_2 0.15 +
0.05log⁡20.05+0.2log⁡20.2+0.15log⁡20.15+0.05 \log_2 0.05 + 0.2 \log_2 0.2 + 0.15 \log_2 0.15 +
0.05log⁡20.05+0.05log⁡20.05+0.1log⁡20.1]0.05 \log_2 0.05 + 0.05 \log_2 0.05 + 0.1 \log_2 0.1 ]

計算 log⁡2\log_2 值:
log⁡20.1≈−3.3219\log_2 0.1 \approx -3.3219
log⁡20.15≈−2.7370\log_2 0.15 \approx -2.7370
log⁡20.05≈−4.3219\log_2 0.05 \approx -4.3219
log⁡20.2≈−2.3219\log_2 0.2 \approx -2.3219

H(X,Y)=−[0.1(−3.3219)+0.15(−2.7370)+0.15(−2.7370)+H(X, Y) = - [ 0.1(-3.3219) + 0.15(-2.7370) + 0.15(-2.7370) +
0.05(−4.3219)+0.2(−2.3219)+0.15(−2.7370)+0.05(-4.3219) + 0.2(-2.3219) + 0.15(-2.7370) +
0.05(−4.3219)+0.05(−4.3219)+0.1(−3.3219)]0.05(-4.3219) + 0.05(-4.3219) + 0.1(-3.3219) ]

H(X,Y)=−[−0.33219−0.41055−0.41055−0.216095−0.46438−0.41055−0.216095−0.216095−0.33219]H(X, Y) = - [ -0.33219 - 0.41055 - 0.41055 - 0.216095 - 0.46438 - 0.41055 - 0.216095 - 0.216095 - 0.33219 ]
H(X,Y)=−[−2.99865]≈3.0H(X, Y) = - [ -2.99865 ] \approx 3.0 bits。

精確計算:
0.1log⁡20.1≈0.332190.1 \log_2 0.1 \approx 0.33219
0.15log⁡20.15≈0.410550.15 \log_2 0.15 \approx 0.41055
0.05log⁡20.05≈0.2160950.05 \log_2 0.05 \approx 0.216095
0.2log⁡20.2≈0.464380.2 \log_2 0.2 \approx 0.46438
0.1log⁡20.1≈0.332190.1 \log_2 0.1 \approx 0.33219

H(X,Y)=−(0.33219+2×0.41055+0.216095+0.46438+0.41055+2×0.216095+0.33219)H(X, Y) = - (0.33219 + 2 \times 0.41055 + 0.216095 + 0.46438 + 0.41055 + 2 \times 0.216095 + 0.33219)
H(X,Y)=−(0.33219+0.8211+0.216095+0.46438+0.41055+0.43219+0.33219)H(X, Y) = - (0.33219 + 0.8211 + 0.216095 + 0.46438 + 0.41055 + 0.43219 + 0.33219)
H(X,Y)=−(3.008695)H(X, Y) = - (3.008695)
取近似值:H(X,Y)≈3.009H(X,Y) \approx 3.009 bits。

【答案】
H(X,Y)≈3.009H(X, Y) \approx 3.009 bits。

(B) 計算邊際熵 H(X)。
H(X)=−∑xP(X=x)log⁡2P(X=x)H(X) = -\sum_x P(X=x) \log_2 P(X=x)。
P(X=−1)=0.2P(X=-1) = 0.2
P(X=0)=0.4P(X=0) = 0.4
P(X=1)=0.4P(X=1) = 0.4

H(X)=−[0.2log⁡20.2+0.4log⁡20.4+0.4log⁡20.4]H(X) = - [ 0.2 \log_2 0.2 + 0.4 \log_2 0.4 + 0.4 \log_2 0.4 ]
H(X)=−[0.2log⁡20.2+0.8log⁡20.4]H(X) = - [ 0.2 \log_2 0.2 + 0.8 \log_2 0.4 ]
log⁡20.2≈−2.3219\log_2 0.2 \approx -2.3219
log⁡20.4≈−1.3219\log_2 0.4 \approx -1.3219

H(X)=−[0.2(−2.3219)+0.8(−1.3219)]H(X) = - [ 0.2 (-2.3219) + 0.8 (-1.3219) ]
H(X)=−[−0.46438−1.05752]H(X) = - [ -0.46438 - 1.05752 ]
H(X)=−[−1.5219]=1.5219H(X) = - [ -1.5219 ] = 1.5219 bits。

精確計算:
0.2log⁡20.2≈0.464380.2 \log_2 0.2 \approx 0.46438
0.4log⁡20.4≈0.528770.4 \log_2 0.4 \approx 0.52877
H(X)=−(0.46438+2×0.52877)=−(0.46438+1.05754)=−1.52192H(X) = - (0.46438 + 2 \times 0.52877) = - (0.46438 + 1.05754) = -1.52192.
H(X)≈1.522H(X) \approx 1.522 bits。

【答案】
H(X)≈1.522H(X) \approx 1.522 bits。

(C) 計算條件熵 H(X|Y)。
條件熵的定義為 H(X∣Y)=∑yP(Y=y)H(X∣Y=y)H(X|Y) = \sum_y P(Y=y) H(X|Y=y)。
H(X∣Y=y)=−∑xP(X=x∣Y=y)log⁡2P(X=x∣Y=y)H(X|Y=y) = -\sum_x P(X=x|Y=y) \log_2 P(X=x|Y=y)。
P(X=x∣Y=y)=P(x,y)P(Y=y)P(X=x|Y=y) = \frac{P(x,y)}{P(Y=y)}。

計算 H(X∣Y=2)H(X|Y=2): P(Y=2)=0.4P(Y=2) = 0.4.
P(X=−1∣Y=2)=0.1/0.4=0.25P(X=-1|Y=2) = 0.1/0.4 = 0.25
P(X=0∣Y=2)=0.15/0.4=0.375P(X=0|Y=2) = 0.15/0.4 = 0.375
P(X=1∣Y=2)=0.15/0.4=0.375P(X=1|Y=2) = 0.15/0.4 = 0.375

🔒

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

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

免費註冊

其他考古題