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

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

第 1 題10 分

Consider a signal xs(t)=x(t)s(t)x_s(t) = x(t)s(t), where x(t)x(t) is a bandlimited signal with bandwidth WBW_B Hz and s(t)s(t) is a periodic signal, expressed as s(t)=∑k=−∞∞Akδ(t−kT0)s(t) = \sum_{k=-\infty}^{\infty} A_k \delta(t - k T_0), where
Ak={1,if k is even0.8,if k is oddA_k = \begin{cases} 1, & \text{if } k \text{ is even} \\ 0.8, & \text{if } k \text{ is odd} \end{cases}
(a) Find the spectrum of s(t)s(t).
(b) What is the maximum value of T0T_0 such that x(t)x(t) can be fully recovered from xs(t)x_s(t) by passing xs(t)x_s(t) through a linear time-invariant filter?
(c) What is the impulse response of the linear filter in (b)?

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

這一題的完整詳解

核心觀念

本題考查:

  1. 週期脈衝串的傅立葉轉換。
  2. 時域相乘對應頻域摺積。
  3. 取樣頻譜的頻譜重複與避免混疊條件。
  4. 理想低通濾波器的設計與脈衝響應。

假設 x(t)x(t) 為低通限頻訊號,其頻譜滿足

X(f)=0,∣f∣>WBX(f)=0,\qquad |f|>W_B

因此原始訊號頻譜範圍為 [−WB,WB][-W_B,W_B]。


(a)求 s(t)s(t) 的頻譜

題目給定

s(t)=∑k=−∞∞Akδ(t−kT0)s(t)=\sum_{k=-\infty}^{\infty}A_k\delta(t-kT_0)

其中

Ak={1,k 為偶數0.8,k 為奇數A_k= \begin{cases} 1, & k\text{ 為偶數}\\ 0.8, & k\text{ 為奇數} \end{cases}

將 AkA_k 改寫成

Ak=0.9+0.1(−1)kA_k=0.9+0.1(-1)^k

因為:

  • kk 為偶數時,Ak=0.9+0.1=1A_k=0.9+0.1=1
  • kk 為奇數時,Ak=0.9−0.1=0.8A_k=0.9-0.1=0.8

所以

s(t)=0.9∑k=−∞∞δ(t−kT0)+0.1∑k=−∞∞(−1)kδ(t−kT0)s(t) = 0.9\sum_{k=-\infty}^{\infty}\delta(t-kT_0) + 0.1\sum_{k=-\infty}^{\infty}(-1)^k\delta(t-kT_0)

利用週期脈衝串的傅立葉轉換:

∑k=−∞∞δ(t−kT0)⟷1T0∑n=−∞∞δ(f−nT0)\sum_{k=-\infty}^{\infty}\delta(t-kT_0) \longleftrightarrow \frac{1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n}{T_0}\right)

而

(−1)k=ejπk(-1)^k=e^{j\pi k}

會使頻譜向右平移 1/(2T0)1/(2T_0),因此

∑k=−∞∞(−1)kδ(t−kT0)⟷1T0∑n=−∞∞δ(f−n+12T0)\sum_{k=-\infty}^{\infty}(-1)^k\delta(t-kT_0) \longleftrightarrow \frac{1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n+\frac12}{T_0}\right)

故

S(f)=0.9T0∑n=−∞∞δ(f−nT0)+0.1T0∑n=−∞∞δ(f−n+12T0)\boxed{ S(f)= \frac{0.9}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n}{T_0}\right) + \frac{0.1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n+\frac12}{T_0}\right) }

也可將頻率位置寫成

S(f)=0.9T0∑n=−∞∞δ(f−nT0)+0.1T0∑n=−∞∞δ(f−2n+12T0)\boxed{ S(f)= \frac{0.9}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n}{T_0}\right) + \frac{0.1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{2n+1}{2T_0}\right) }

頻譜結構

s(t)s(t) 的週期為 2T02T_0,所以頻譜基本間隔為

12T0\frac{1}{2T_0}

其中:

  • 偶數次諧波位於 f=n/T0f=n/T_0,權重為 0.9/T00.9/T_0
  • 奇數次諧波位於 f=(2n+1)/(2T0)f=(2n+1)/(2T_0),權重為 0.1/T00.1/T_0

(b)T0T_0 的最大值

由時域相乘定理:

xs(t)=x(t)s(t)⟺Xs(f)=X(f)∗S(f)x_s(t)=x(t)s(t) \quad\Longleftrightarrow\quad X_s(f)=X(f)*S(f)

因此

Xs(f)=0.9T0∑n=−∞∞X(f−nT0)+0.1T0∑n=−∞∞X(f−n+12T0)X_s(f) = \frac{0.9}{T_0} \sum_{n=-\infty}^{\infty} X\left(f-\frac{n}{T_0}\right) + \frac{0.1}{T_0} \sum_{n=-\infty}^{\infty} X\left(f-\frac{n+\frac12}{T_0}\right)

原始頻譜 X(f)X(f) 位於 [−WB,WB][-W_B,W_B]。由於 S(f)S(f) 中最近的非零頻譜線位於

f=±12T0f=\pm\frac{1}{2T_0}
🔒

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

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

免費註冊

第 2 題10 分

Assume the DSB demodulation process is modeled as below, where the received signal is expressed as xr(t)=m(t)cos⁡(2πfct)+n(t)x_r(t) = m(t) \cos(2\pi f_c t) + n(t) with the noise n(t)n(t) being modeled as white, Gaussian noise with two-sided power spectral density of N0/2=1×10−10N_0/2 = 1 \times 10^{-10} W/Hz.

🖼️【此處有附圖,請對照原卷】
The figure shows a block diagram of a demodulator. The received signal xr(t)x_r(t) is first passed through a bandpass filter H1(f)H_1(f) centered at fcf_c. The output of the bandpass filter is then multiplied by cos⁡(2πfdt)\cos(2\pi f_d t). The result is then passed through a lowpass filter H2(f)H_2(f) to obtain the output y(t)y(t). The magnitude responses of H1(f)H_1(f) and H2(f)H_2(f) are given in the figure.

Assume the demodulation output y(t)y(t) can be decomposed into two parts, the signal part ys(t)y_s(t) and the noise part yn(t)y_n(t).

Case 1: fc=fd=200,000f_c = f_d = 200,000 Hz.
(a) Find ys(t)y_s(t).
(b) Find the total power of yn(t)y_n(t).

Case 2: fc=200,000f_c = 200,000 Hz but fd=200,100f_d = 200, 100 Hz.
(c) Plot the spectrum of ys(t)y_s(t).
(d) Find the auto-correlation of yn(t)y_n(t).

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

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

這一題的完整詳解

核心觀念

本題考查:

  1. DSB-SC 訊號的頻譜平移:
    xs(t)=m(t)cos⁡(2πfct)x_s(t)=m(t)\cos(2\pi f_ct)
    Xs(f)=12[M(f−fc)+M(f+fc)]X_s(f)=\frac{1}{2}\left[M(f-f_c)+M(f+f_c)\right]

  2. 乘上載波 cos⁡(2πfdt)\cos(2\pi f_dt) 所造成的頻譜平移:

F{v(t)cos⁡(2πfdt)}=12[V(f−fd)+V(f+fd)]\mathcal{F}\{v(t)\cos(2\pi f_dt)\} =\frac{1}{2}\left[V(f-f_d)+V(f+f_d)\right]
  1. 白色高斯雜訊通過濾波器後,其功率頻譜密度為
    Sn(f)=N02∣H(f)∣2S_n(f)=\frac{N_0}{2}|H(f)|^2

圖中資訊為:

  • M(f)M(f) 為三角形頻譜,支撐於 ∣f∣≤500|f|\leq 500 Hz,且 M(0)=1M(0)=1。
  • H1(f)H_1(f) 在 ±fc\pm f_c 附近通過,頻寬各為 10001000 Hz,增益為 22。
  • H2(f)H_2(f) 通過 ∣f∣≤500|f|\leq 500 Hz,增益為 11。
  • 雜訊雙邊 PSD:
    N02=10−10 W/Hz\frac{N_0}{2}=10^{-10}\ \text{W/Hz}

解題方法

先利用 H1(f)H_1(f) 保留 DSB 訊號的兩個旁帶,再乘上 cos⁡(2πfdt)\cos(2\pi f_dt) 將頻譜搬移至低頻,最後由 H2(f)H_2(f) 保留基頻部分。

對訊號而言,通過 H1(f)H_1(f) 後,由於 H1(f)H_1(f) 增益為 22,恰好抵消 DSB 調變頻譜中的係數 1/21/2:

X1(f)=M(f−fc)+M(f+fc)X_1(f)=M(f-f_c)+M(f+f_c)

再乘上 cos⁡(2πfdt)\cos(2\pi f_dt):

Ys(f)=12[X1(f−fd)+X1(f+fd)]Y_s(f)=\frac{1}{2}\left[X_1(f-f_d)+X_1(f+f_d)\right]


Case 1:fc=fd=200000f_c=f_d=200000 Hz

(a) 求 ys(t)y_s(t)

將 fc=fdf_c=f_d 代入:

Ys(f)=12[M(f)+M(f−2fc)+M(f+2fc)+M(f)]Y_s(f)=\frac{1}{2} \left[ M(f)+M(f-2f_c)+M(f+2f_c)+M(f) \right]

經過低通濾波器 H2(f)H_2(f) 後,只保留基頻部分:

Ys(f)=12[M(f)+M(f)]=M(f)Y_s(f)=\frac{1}{2}[M(f)+M(f)]=M(f)

因此

ys(t)=m(t)\boxed{y_s(t)=m(t)}

此時為理想同步解調,載波頻率完全匹配,所以原始訊號可完整恢復。


(b) 求 yn(t)y_n(t) 的總功率

H1(f)H_1(f) 的增益為 22,因此通過 H1(f)H_1(f) 後:

Sn1(f)=10−10×∣2∣2=4×10−10 W/HzS_{n1}(f)=10^{-10}\times |2|^2 =4\times10^{-10}\ \text{W/Hz}

在乘法器之後,每一個頻譜平移分量乘上 1/21/2,所以其 PSD 乘上 1/41/4。在低通濾波器通過的 ∣f∣≤500|f|\leq500 Hz 範圍內,正、負頻帶均搬移至基頻,因此共有兩個雜訊分量:

Syn(f)=2(14)(4×10−10)=2×10−10 W/HzS_{yn}(f) =2\left(\frac14\right) \left(4\times10^{-10}\right) =2\times10^{-10}\ \text{W/Hz}

低通頻寬為 10001000 Hz,因此雜訊總功率為

Pyn=∫−500500Syn(f) dfP_{yn} =\int_{-500}^{500}S_{yn}(f)\,df Pyn=(2×10−10)(1000)=2×10−7 WP_{yn} =(2\times10^{-10})(1000) =2\times10^{-7}\ \text{W}

因此

Pyn=2×10−7 W\boxed{P_{yn}=2\times10^{-7}\ \text{W}}


Case 2:fc=200000f_c=200000 Hz,fd=200100f_d=200100 Hz

此時解調頻率比載波頻率高 100100 Hz:

fd−fc=100 Hzf_d-f_c=100\ \text{Hz}

因此解調後的訊號頻譜會偏移 ±100\pm100 Hz。

(c) 繪出 ys(t)y_s(t) 的頻譜

低通範圍內的訊號頻譜為

Ys(f)=12[M(f−100)+M(f+100)],∣f∣≤500Y_s(f) =\frac12\left[M(f-100)+M(f+100)\right], \qquad |f|\leq500

其中

🔒

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

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

免費註冊

第 3 題13 分

Assume X(t) and Y(t) are two wide-sense stationary random processes with the autocorrelation functions RX(τ)=10e−2∣τ∣+2R_X(\tau) = 10e^{-2|\tau|} + 2 and RY(τ)=5e−10∣τ∣+5R_Y(\tau) = 5e^{-10|\tau|} + 5, respectively.

(a) Find the mean functions E[X(t)]E[X(t)] and E[Y(t)]E[Y(t)].
(b) Find the variance functions Var[X(t)]Var[X(t)] and Var[Y(t)]Var[Y(t)].
(c) Which random process has a wider bandwidth? Why?
(d) Which random process has a larger power? Why?
(e) Which random process is more like a Gaussian random process? Why?

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

這一題的完整詳解

核心觀念

本題考查寬平穩隨機程序(WSS)的自相關函數、平均功率、變異數、功率譜密度與頻寬,以及高斯隨機程序的判定。

對寬平穩隨機程序:

RX(τ)=E[X(t)X(t+τ)]R_X(\tau)=E[X(t)X(t+\tau)]

若自相關函數可分解為

RX(τ)=CX(τ)+mX2R_X(\tau)=C_X(\tau)+m_X^2

其中 mX=E[X(t)]m_X=E[X(t)],CX(τ)C_X(\tau) 為自共變異函數,則

CX(τ)=RX(τ)−mX2C_X(\tau)=R_X(\tau)-m_X^2

且

Var⁡[X(t)]=CX(0)=RX(0)−mX2\operatorname{Var}[X(t)]=C_X(0)=R_X(0)-m_X^2

平均功率為

PX=E[X2(t)]=RX(0)P_X=E[X^2(t)]=R_X(0)

解題方法

題目中的指數項會隨著 ∣τ∣|\tau| 增大而趨近於零,因此常數項代表平均值平方:

mX2=lim⁡∣τ∣→∞RX(τ)m_X^2=\lim_{|\tau|\to\infty}R_X(\tau) mY2=lim⁡∣τ∣→∞RY(τ)m_Y^2=\lim_{|\tau|\to\infty}R_Y(\tau)

採用隨機程序平均值為非負的慣例,取正平方根。

(a) 平均函數

對 X(t)X(t):

RX(τ)=10e−2∣τ∣+2R_X(\tau)=10e^{-2|\tau|}+2

因此

mX2=lim⁡∣τ∣→∞RX(τ)=2m_X^2=\lim_{|\tau|\to\infty}R_X(\tau)=2 E[X(t)]=mX=2E[X(t)]=m_X=\sqrt{2}

對 Y(t)Y(t):

RY(τ)=5e−10∣τ∣+5R_Y(\tau)=5e^{-10|\tau|}+5

因此

mY2=lim⁡∣τ∣→∞RY(τ)=5m_Y^2=\lim_{|\tau|\to\infty}R_Y(\tau)=5 E[Y(t)]=mY=5E[Y(t)]=m_Y=\sqrt{5}

嚴格而言,自相關函數只能決定平均值的平方,因此平均值也可能是負值;本題依一般教科書慣例取非負平均值。


(b) 變異數函數

由於兩個程序皆為寬平穩,平均值與自相關函數均與時間 tt 無關,所以變異數也是常數。

對 X(t)X(t):

RX(0)=10e0+2=12R_X(0)=10e^0+2=12 Var⁡[X(t)]=RX(0)−mX2=12−2=10\operatorname{Var}[X(t)] =R_X(0)-m_X^2 =12-2 =10

對 Y(t)Y(t):

RY(0)=5e0+5=10R_Y(0)=5e^0+5=10 Var⁡[Y(t)]=RY(0)−mY2=10−5=5\operatorname{Var}[Y(t)] =R_Y(0)-m_Y^2 =10-5 =5

因此

Var⁡[X(t)]=10\boxed{\operatorname{Var}[X(t)]=10} Var⁡[Y(t)]=5\boxed{\operatorname{Var}[Y(t)]=5}

(c) 哪一個隨機程序具有較寬的頻寬?

先扣除平均值所造成的常數項,得到自共變異函數:

CX(τ)=10e−2∣τ∣C_X(\tau)=10e^{-2|\tau|} CY(τ)=5e−10∣τ∣C_Y(\tau)=5e^{-10|\tau|}

利用傅立葉轉換公式:

F{e−a∣τ∣}=2aa2+ω2\mathcal{F}\{e^{-a|\tau|}\} =\frac{2a}{a^2+\omega^2}

可得連續部分的功率譜密度:

SX(ω)=10⋅44+ω2=404+ω2S_X(\omega)=10\cdot\frac{4}{4+\omega^2} =\frac{40}{4+\omega^2}
🔒

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

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

免費註冊

第 4 題10 分

Consider a bandlimited communication system with the overall response of the transmitter, channel, and receiver filters being represented by the pulse shaping waveform p(t)p(t). The Nyquist pulse shaping criterion
∑n=−∞∞P(f+nT)=T\sum_{n=-\infty}^{\infty} P\left(f + \frac{n}{T}\right) = T
provides a rule to design p(t)p(t) to avoid inter-symbol interference (ISI) in the bandlimited system, where TT is the sampling interval and P(f)P(f) is the Fourier transform of the overall pulse p(t)p(t).

(a) A commonly used pulse shaping waveform is the raised cosine pulse defined by
P(f)={a,∣f∣≤W−βb+ccos⁡(π(∣f∣−W+β)2β),W−β≤∣f∣≤W+β0,∣f∣≥W+βP(f) = \begin{cases} a, & |f| \le W - \beta \\ b + c \cos\left(\frac{\pi( |f|-W+\beta)}{2\beta}\right), & W - \beta \le |f| \le W + \beta \\ 0, & |f| \ge W + \beta \end{cases}
where W=12TW = \frac{1}{2T} and β\beta is the excessive bandwidth with 0≤β≤W0 \le \beta \le W. Find the constants a,b,ca, b, c and dd such that the above raised cosine spectrum satisfies the Nyquist pulse shaping criterion.

(b) Consider another pulse shaping spectrum P(f)P(f) shown below, where the vertical axis of the plot corresponds to the value of 1TP(f)\frac{1}{T}P(f) and the labeled frequencies correspond to the corner points of the spectrum. Find one signaling interval TT such that the bandlimited system with the pulse p(t)p(t) results in zero ISI.

🖼️【此處有附圖,請對照原卷】

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

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

這一題的完整詳解

核心觀念

零 ISI 的頻域奈奎斯特條件是

∑n=−∞∞P(f+nT)=T.\sum_{n=-\infty}^{\infty}P\left(f+\frac{n}{T}\right)=T.

也就是將頻譜每隔 1/T1/T 平移後相加,總和必須對所有 ff 都等於 TT。本題 (a) 利用相鄰頻譜在重疊區互補;(b) 則找出圖形平移後能互補的頻率間隔。

(a) 升餘弦頻譜常數

依原卷圖中的式子,轉換帶寫成

P(f)=b+T2cos⁡(π(∣f∣+c)d),W−β≤∣f∣≤W+β.P(f)=b+\frac{T}{2}\cos\left(\frac{\pi(|f|+c)}{d}\right), \qquad W-\beta\le |f|\le W+\beta.

在通帶中央,奈奎斯特條件要求頻譜值為 TT,因此

a=T.a=T.

考慮正頻率一側的轉換帶,令 f∈[W−β,W+β]f\in[W-\beta,W+\beta]。此時與 ff 配對、相隔 1/T=2W1/T=2W 的頻率為 f−2Wf-2W,其絕對值為 2W−f2W-f。兩者的頻譜相加必須等於 TT:

2b+T2[cos⁡(π(f+c)d)+cos⁡(π(2W−f+c)d)]=T.2b+\frac{T}{2}\left[ \cos\left(\frac{\pi(f+c)}{d}\right) +\cos\left(\frac{\pi(2W-f+c)}{d}\right) \right]=T.

兩個餘弦項必須互相抵消,因此兩個角度相加應為 π\pi,且兩個常數項相加為 TT。由此得到

🔒

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

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

免費註冊

第 5 題24 分

Consider a communication system with the vector representation (all with dimension N×1N \times 1):
r=s+n\mathbf{r} = \mathbf{s} + \mathbf{n}
where r\mathbf{r} is the received signal, s∈{S1,S2,…,SM}\mathbf{s} \in \{S_1, S_2, \dots, S_M\} is the transmitted symbol and the noise n\mathbf{n} has i.i.d. Gaussian components each with zero mean and variance N0/2N_0/2. The transmitted signal s\mathbf{s} and the noise n\mathbf{n} are assumed to be independent. All the possible transmitted symbols S1,S2,…,SMS_1, S_2, \dots, S_M have equal prior probability, but may not have equal energy. The signaling interval for each symbol s\mathbf{s} is TT seconds.

Maximum-likelihood (ML) decision rule is often employed at the receiver for symbol decoding and can be represented by
S^ML=arg⁡max⁡smp(r∣s=sm)\hat{S}_{ML} = \arg \max_{\mathbf{s}_m} p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m)
where p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) is the likelihood function associated with sm\mathbf{s}_m.

(a) What is the average signal energy of s\mathbf{s}? What is the transmission bit rate of this system?
(b) Find the likelihood function p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m). (Note: You can express the likelihood function using vector 2-norm ∣∣⋅∣∣2||\cdot||_2, or define appropriate notations by yourself. The probability density function of a N(0,σ2)N(0, \sigma^2) Gaussian random variable ZZ is pZ(z)=12πσ2e−z22σ2p_Z(z) = \frac{1}{\sqrt{2\pi\sigma^2}} e^{-\frac{z^2}{2\sigma^2}}.)
(c) A decoding rule based on correlation of the received signal r\mathbf{r} and each sm\mathbf{s}_m for m=1,2,…,Mm=1, 2, \dots, M can be devised as follows:
"If rTsm≥rTsn\mathbf{r}^T \mathbf{s}_m \ge \mathbf{r}^T \mathbf{s}_n for all n≠mn \neq m, then the decoder decides that sm\mathbf{s}_m is the transmitted symbol."
Is this decoder optimum in terms of symbol error probability? Please explain your answers.
(d) One vector representation of QPSK signaling is designed as S1=[2,0]TS_1 = [\sqrt{2}, 0]^T, S2=[0,2]TS_2 = [0, \sqrt{2}]^T, S3=[2,22]TS_3 = [\sqrt{2}, 2\sqrt{2}]^T and S4=[22,2]TS_4 = [2\sqrt{2}, \sqrt{2}]^T with allocation of two information bits to the four symbols being arranged as follows: (00)→S1(00) \rightarrow S_1, (01)→S2(01) \rightarrow S_2, (10)→S3(10) \rightarrow S_3 and (11)→S4(11) \rightarrow S_4.
i. Given that the transmitted symbol is S1S_1. Find the probability that the ML decoder erroneously gives the decision S^ML=S4\hat{S}_{ML} = S_4 when N0=2N_0 = 2 (Express your answer either in terms of the Gaussian Q-function or the erfc function).
ii. What would be wrong with this design? What would you do to improve any of the system performance in a reasonable way? Please explain.

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

這一題的完整詳解

此題考驗對數位通訊系統中最大概似接收器(ML receiver)、符號錯誤機率、最佳解碼準則、QPSK 調變、以及高斯雜訊下的錯誤機率計算。

核心觀念:

  1. 最大概似接收器 (ML Receiver) 的原理。
  2. 高斯雜訊下的似然函數 (Likelihood Function)。
  3. Correlation Receiver 與 ML Receiver 的關係。
  4. Symbol Error Probability (SEP) 的計算。
  5. QPSK 調變的符號表示與能量。
  6. 高斯 Q-function 和 erfc function。

解題步驟:

系統模型:r=s+n\mathbf{r} = \mathbf{s} + \mathbf{n}
s∈{S1,…,SM}\mathbf{s} \in \{S_1, \dots, S_M\} 是傳送符號,具有相等先驗機率 P(Si)=1/MP(S_i) = 1/M。
n\mathbf{n} 是零均值、對角協方差矩陣為 (N0/2)IN(N_0/2)I_N 的高斯雜訊。
N0/2N_0/2 是單邊功率譜密度。

(a) 平均信號能量與傳輸位元率

平均信號能量 EsE_s 定義為所有可能符號能量的平均值:
Es=E[∣∣s∣∣22]=1M∑m=1M∣∣Sm∣∣22E_s = E[||\mathbf{s}||_2^2] = \frac{1}{M} \sum_{m=1}^M ||S_m||_2^2
題目提到符號能量可能不相等。

傳輸位元率 RbR_b:
假設每個符號 SmS_m 代表 kk 個位元。則 M=2kM = 2^k。
傳輸速率(符號速率)是 Rs=1/TR_s = 1/T (符號/秒)。
傳輸位元率 Rb=k×Rs=k/TR_b = k \times R_s = k/T (位元/秒)。
題目沒有給出 MM 或 kk,也沒有給出 TT。
但是,題目有提到 TT 秒是信號間隔。
所以 Rs=1/TR_s = 1/T。
如果 MM 個符號代表 kk 個位元,則 M=2kM = 2^k。
位元率 Rb=k×(1/T)R_b = k \times (1/T)。
如果 MM 和 kk 未知,則無法計算。
重新閱讀題目: 題目沒有給出 MM 或 kk。所以這部分可能需要用符號表示。
平均信號能量: Es=1M∑m=1M∣∣Sm∣∣22E_s = \frac{1}{M} \sum_{m=1}^M ||S_m||_2^2。
傳輸位元率: Rb=log⁡2MTR_b = \frac{\log_2 M}{T} bits/sec。

答案 (a):
平均信號能量:Es=1M∑m=1M∣∣Sm∣∣22E_s = \frac{1}{M} \sum_{m=1}^M ||S_m||_2^2。
傳輸位元率:Rb=log⁡2MTR_b = \frac{\log_2 M}{T}。

(b) 尋找似然函數 p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m)

似然函數是給定傳送符號 sm\mathbf{s}_m 時,接收到信號 r\mathbf{r} 的條件機率密度函數。
我們有 r=sm+n\mathbf{r} = \mathbf{s}_m + \mathbf{n}。
雜訊 n\mathbf{n} 是 NN 維高斯向量,其各分量獨立,均值為 0,方差為 N0/2N_0/2。
所以 n∼N(0,N02IN)\mathbf{n} \sim N(\mathbf{0}, \frac{N_0}{2} I_N)。
則 r−sm=n\mathbf{r} - \mathbf{s}_m = \mathbf{n}。
r−sm\mathbf{r} - \mathbf{s}_m 是一個均值為 sm\mathbf{s}_m、協方差為 N02IN\frac{N_0}{2} I_N 的高斯向量。
一個 NN 維高斯隨機向量 x\mathbf{x},其均值為 μ\boldsymbol{\mu},協方差矩陣為 Σ\Sigma,其 PDF 為:
p(x)=1(2π)Ndet⁡(Σ)exp⁡(−12(x−μ)TΣ−1(x−μ))p(\mathbf{x}) = \frac{1}{\sqrt{(2\pi)^N \det(\Sigma)}} \exp\left(-\frac{1}{2}(\mathbf{x}-\boldsymbol{\mu})^T \Sigma^{-1} (\mathbf{x}-\boldsymbol{\mu})\right)
在這裡,x=r\mathbf{x} = \mathbf{r},μ=sm\boldsymbol{\mu} = \mathbf{s}_m,Σ=N02IN\Sigma = \frac{N_0}{2} I_N。
det⁡(Σ)=(N02)N\det(\Sigma) = (\frac{N_0}{2})^N。
Σ−1=2N0IN\Sigma^{-1} = \frac{2}{N_0} I_N。
所以,似然函數為:
p(r∣s=sm)=1(2π)N(N02)Nexp⁡(−12(r−sm)T(2N0IN)(r−sm))p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{\sqrt{(2\pi)^N (\frac{N_0}{2})^N}} \exp\left(-\frac{1}{2}(\mathbf{r}-\mathbf{s}_m)^T \left(\frac{2}{N_0} I_N\right) (\mathbf{r}-\mathbf{s}_m)\right)
p(r∣s=sm)=1(N02)N/2(2π)Nexp⁡(−1N0(r−sm)T(r−sm))p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\frac{N_0}{2})^{N/2} \sqrt{(2\pi)^N}} \exp\left(-\frac{1}{N_0}(\mathbf{r}-\mathbf{s}_m)^T (\mathbf{r}-\mathbf{s}_m)\right)
p(r∣s=sm)=1(N02)N/2(2π)N/2exp⁡(−1N0∣∣r−sm∣∣22)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\frac{N_0}{2})^{N/2} (2\pi)^{N/2}} \exp\left(-\frac{1}{N_0} ||\mathbf{r}-\mathbf{s}_m||_2^2\right)
p(r∣s=sm)=1(πN0)N/2exp⁡(−1N0∣∣r−sm∣∣22)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\pi N_0)^{N/2}} \exp\left(-\frac{1}{N_0} ||\mathbf{r}-\mathbf{s}_m||_2^2\right)
其中 ∣∣r−sm∣∣22=(r−sm)T(r−sm)=∑i=1N(ri−sm,i)2||\mathbf{r}-\mathbf{s}_m||_2^2 = (\mathbf{r}-\mathbf{s}_m)^T (\mathbf{r}-\mathbf{s}_m) = \sum_{i=1}^N (r_i - s_{m,i})^2。

ML 解碼器是最大化這個似然函數。由於指數函數是單調遞增的,最大化似然函數等價於最小化指數的指數項。
最大化 p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) 等價於最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
s^ML=arg⁡min⁡sm∣∣r−sm∣∣22\hat{\mathbf{s}}_{ML} = \arg \min_{\mathbf{s}_m} ||\mathbf{r}-\mathbf{s}_m||_2^2
這也稱為最小距離準則 (Minimum Distance Criterion)。

答案 (b):
p(r∣s=sm)=1(πN0)N/2exp⁡(−1N0∣∣r−sm∣∣22)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\pi N_0)^{N/2}} \exp\left(-\frac{1}{N_0} ||\mathbf{r}-\mathbf{s}_m||_2^2\right)。

(c) 相關性解碼器是否最佳(在符號錯誤機率方面)

相關性解碼器規則:
若 rTsm≥rTsn\mathbf{r}^T \mathbf{s}_m \ge \mathbf{r}^T \mathbf{s}_n 對所有 n≠mn \neq m,則選擇 sm\mathbf{s}_m。
我們知道 ML 解碼器是最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
∣∣r−sm∣∣22=(r−sm)T(r−sm)=rTr−2rTsm+smTsm||\mathbf{r}-\mathbf{s}_m||_2^2 = (\mathbf{r}-\mathbf{s}_m)^T (\mathbf{r}-\mathbf{s}_m) = \mathbf{r}^T \mathbf{r} - 2\mathbf{r}^T \mathbf{s}_m + \mathbf{s}_m^T \mathbf{s}_m
∣∣r−sm∣∣22=∣∣r∣∣22−2rTsm+∣∣sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 = ||\mathbf{r}||_2^2 - 2\mathbf{r}^T \mathbf{s}_m + ||\mathbf{s}_m||_2^2
最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 等價於最大化 2rTsm−∣∣sm∣∣222\mathbf{r}^T \mathbf{s}_m - ||\mathbf{s}_m||_2^2。
或者,如果所有符號的能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 都相等,則最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 等價於最大化 rTsm\mathbf{r}^T \mathbf{s}_m。
此時,ML 解碼器就變成了相關性解碼器。

然而,題目明確指出「符號能量可能不相等」。
如果符號能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 不相等,那麼最大化 rTsm\mathbf{r}^T \mathbf{s}_m 並不一定等同於最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
ML 解碼器是最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
相關性解碼器是最大化 rTsm\mathbf{r}^T \mathbf{s}_m。

考慮以下情況:
令 r=[1,0]T\mathbf{r} = [1, 0]^T。
令 S1=[1,0]TS_1 = [1, 0]^T,則 ∣∣s1∣∣22=1||\mathbf{s}_1||_2^2 = 1。rTs1=1\mathbf{r}^T \mathbf{s}_1 = 1。∣∣r−s1∣∣22=0||\mathbf{r}-\mathbf{s}_1||_2^2 = 0。
令 S2=[2,0]TS_2 = [2, 0]^T,則 ∣∣s2∣∣22=4||\mathbf{s}_2||_2^2 = 4。rTs2=2\mathbf{r}^T \mathbf{s}_2 = 2。∣∣r−s2∣∣22=∣∣[1,0]T−[2,0]T∣∣22=∣∣[−1,0]T∣∣22=1||\mathbf{r}-\mathbf{s}_2||_2^2 = ||[1,0]^T - [2,0]^T||_2^2 = ||[-1,0]^T||_2^2 = 1。
ML 解碼器會選擇 S1S_1 (因為 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 最小,為 0)。
相關性解碼器會選擇 S2S_2 (因為 rTsm\mathbf{r}^T \mathbf{s}_m 最大,2 > 1)。
在這種情況下,相關性解碼器不是最佳的。

解釋:
ML 解碼器選擇 sm\mathbf{s}_m 使得 p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) 最大。
這等價於最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
∣∣r−sm∣∣22=∣∣r∣∣22−2rTsm+∣∣sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 = ||\mathbf{r}||_2^2 - 2\mathbf{r}^T \mathbf{s}_m + ||\mathbf{s}_m||_2^2。
ML 解碼器選擇 sm\mathbf{s}_m 使得 rTsm−12∣∣sm∣∣22\mathbf{r}^T \mathbf{s}_m - \frac{1}{2}||\mathbf{s}_m||_2^2 最大。
相關性解碼器選擇 sm\mathbf{s}_m 使得 rTsm\mathbf{r}^T \mathbf{s}_m 最大。

只有當所有符號能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 都相等時,最大化 rTsm\mathbf{r}^T \mathbf{s}_m 才等價於最大化 rTsm−12∣∣sm∣∣22\mathbf{r}^T \mathbf{s}_m - \frac{1}{2}||\mathbf{s}_m||_2^2。
由於題目明確指出符號能量可能不相等,所以相關性解碼器不是最佳的。

答案 (c):
否。這個解碼器不是最佳的。
解釋:最大概似解碼器是最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2,這等價於最大化 rTsm−12∣∣sm∣∣22\mathbf{r}^T \mathbf{s}_m - \frac{1}{2}||\mathbf{s}_m||_2^2。相關性解碼器則是最大化 rTsm\mathbf{r}^T \mathbf{s}_m。只有當所有符號的能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 都相等時,這兩者才等價。由於題目說明符號能量可能不相等,因此相關性解碼器不保證是最佳的(即能獲得最小的符號錯誤機率)。

(d) QPSK 符號表示與錯誤機率計算

QPSK 符號:
S1=[2,0]TS_1 = [\sqrt{2}, 0]^T (00)
S2=[0,2]TS_2 = [0, \sqrt{2}]^T (01)
S3=[2,22]TS_3 = [\sqrt{2}, 2\sqrt{2}]^T (10)
S4=[22,2]TS_4 = [2\sqrt{2}, \sqrt{2}]^T (11)
N=2N=2 (向量維度)。
N0=2N_0 = 2。所以 N0/2=1N_0/2 = 1 (雜訊方差)。

(d) i. P(S^ML=S4∣S1 transmitted)P(\hat{S}_{ML} = S_4 | S_1 \text{ transmitted})

ML 解碼器選擇 sm\mathbf{s}_m 使得 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 最小。
我們知道 r=s1+n\mathbf{r} = \mathbf{s}_1 + \mathbf{n}。
ML 解碼器錯誤地選擇 S4S_4 而不是 S1S_1 當:
∣∣r−s4∣∣22<∣∣r−s1∣∣22||\mathbf{r}-\mathbf{s}_4||_2^2 < ||\mathbf{r}-\mathbf{s}_1||_2^2
代入 r=s1+n\mathbf{r} = \mathbf{s}_1 + \mathbf{n}:
∣∣s1+n−s4∣∣22<∣∣s1+n−s1∣∣22||\mathbf{s}_1 + \mathbf{n} - \mathbf{s}_4||_2^2 < ||\mathbf{s}_1 + \mathbf{n} - \mathbf{s}_1||_2^2
∣∣n−(s4−s1)∣∣22<∣∣n∣∣22||\mathbf{n} - (\mathbf{s}_4 - \mathbf{s}_1)||_2^2 < ||\mathbf{n}||_2^2

計算向量差:
s4−s1=[22,2]T−[2,0]T=[2,2]T\mathbf{s}_4 - \mathbf{s}_1 = [2\sqrt{2}, \sqrt{2}]^T - [\sqrt{2}, 0]^T = [\sqrt{2}, \sqrt{2}]^T。
令 Δ41=s4−s1=[2,2]T\Delta_{41} = \mathbf{s}_4 - \mathbf{s}_1 = [\sqrt{2}, \sqrt{2}]^T。
不等式變為:
∣∣n−Δ41∣∣22<∣∣n∣∣22||\mathbf{n} - \Delta_{41}||_2^2 < ||\mathbf{n}||_2^2
(n−Δ41)T(n−Δ41)<nTn(\mathbf{n} - \Delta_{41})^T (\mathbf{n} - \Delta_{41}) < \mathbf{n}^T \mathbf{n}
nTn−2nTΔ41+Δ41TΔ41<nTn\mathbf{n}^T \mathbf{n} - 2 \mathbf{n}^T \Delta_{41} + \Delta_{41}^T \Delta_{41} < \mathbf{n}^T \mathbf{n}
−2nTΔ41+∣∣Δ41∣∣22<0- 2 \mathbf{n}^T \Delta_{41} + ||\Delta_{41}||_2^2 < 0
∣∣Δ41∣∣22<2nTΔ41||\Delta_{41}||_2^2 < 2 \mathbf{n}^T \Delta_{41}

計算 ∣∣Δ41∣∣22||\Delta_{41}||_2^2:
∣∣Δ41∣∣22=(2)2+(2)2=2+2=4||\Delta_{41}||_2^2 = (\sqrt{2})^2 + (\sqrt{2})^2 = 2 + 2 = 4。
所以,4<2nTΔ414 < 2 \mathbf{n}^T \Delta_{41}。
2<nTΔ412 < \mathbf{n}^T \Delta_{41}。
nTΔ41=n1(2)+n2(2)=2(n1+n2)\mathbf{n}^T \Delta_{41} = n_1 (\sqrt{2}) + n_2 (\sqrt{2}) = \sqrt{2} (n_1 + n_2)。
所以,2<2(n1+n2)2 < \sqrt{2} (n_1 + n_2)。
2<n1+n2\sqrt{2} < n_1 + n_2
令 Z=n1+n2Z = n_1 + n_2。

🔒

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

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

免費註冊

第 6 題33 分

Consider the binary symmetric channel (BSC) and the binary erasure channel (BEC) shown below. In the BSC, the transmitted bit X (a zero or a one) will be received as Y with a crossover (error) probability of p. In the BEC, the transmitted bit will be lost (erased) with a probability p.

🖼️【此處有附圖,請對照原卷】
The figure shows two channel models: BSC and BEC.
For BSC: P(Y=0∣X=1)=pP(Y=0|X=1) = p, P(Y=1∣X=0)=pP(Y=1|X=0) = p.
For BEC: P(Y=Erasure∣X=0)=pP(Y=\text{Erasure}|X=0) = p, P(Y=Erasure∣X=1)=pP(Y=\text{Erasure}|X=1) = p.

(a) The channel capacity is the maximized mutual information I(X;Y)=H(Y)−H(Y∣X)I(X; Y) = H(Y) - H(Y|X). With known p and p, explain when the maximization can be achieved.
(b) What is the channel capacity of the BSC and BEC respectively? You need to show your derivation.
(Note: the conditional entropy is defined by H(Y∣X)=−∑i∑jP(xi,yj)log⁡2P(yj∣xi)H(Y|X) = -\sum_{i}\sum_{j} P(x_i, y_j) \log_2 P(y_j|x_i).)
(c) If 0<p=p<0.50 < p = p < 0.5, which channel has the higher capacity? Why?
(d) Consider the composite channel shown below. The channel W's are two independent identical BEC considered above. The pair (X1,X2)(X_1, X_2) is encoded from the source via X1=XOR(U1,U2)X_1 = \text{XOR}(U_1, U_2) and X2=U2X_2 = U_2. Please write down the generator (encoding) matrix.
🖼️【此處有附圖,請對照原卷】
The figure shows a composite channel with two inputs U1,U2U_1, U_2 and two outputs Y1,Y2Y_1, Y_2. The intermediate variables are X1,X2X_1, X_2. The channels W are BECs.

(e) Following (d), we can deduce the source (U1,U2)(U_1, U_2) by decoding the received (Y1,Y2)(Y_1, Y_2). As a matter of fact, we can regard the channel from U1U_1 to (Y1,Y2)(Y_1, Y_2) as a BEC. What is its erasure probability?
(f) Following (e), if the correct U1U_1 is already known by the receiver, then the channel from U2U_2 to (Y1,Y2)(Y_1, Y_2) can also be regarded as a BEC. What is its erasure probability?

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

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

這一題的完整詳解

核心觀念

本題考查:

  • 通道容量定義:C=max⁡P(X)I(X;Y)C=\max_{P(X)} I(X;Y)。

  • 二元熵函數:

    h2(q)=−qlog⁡2q−(1−q)log⁡2(1−q)h_2(q)=-q\log_2q-(1-q)\log_2(1-q)
  • BSC 的交叉錯誤機率為 pp。

  • BEC 的抹除機率依圖記為 ρ\rho。

  • 線性區塊碼的生成矩陣,以及串聯通道的等效抹除機率。


(a) 容量最大化的條件

令

P(X=0)=q,P(X=1)=1−q.P(X=0)=q,\qquad P(X=1)=1-q.

BSC

BSC 的輸出為

P(Y=1)=qp+(1−q)(1−p).P(Y=1)=q p+(1-q)(1-p).

因此

I(X;Y)=H(Y)−H(Y∣X).I(X;Y)=H(Y)-H(Y|X).

由於 BSC 的條件熵為

H(Y∣X)=h2(p),H(Y|X)=h_2(p),

且當輸出 YY 為均勻分布時,H(Y)H(Y) 最大。輸出均勻的條件是

P(Y=0)=P(Y=1)=12.P(Y=0)=P(Y=1)=\frac12.

對稱通道下,這要求輸入也採用均勻分布:

P(X=0)=P(X=1)=12.P(X=0)=P(X=1)=\frac12.

BEC

BEC 的輸出符號為 0,1,Erasure0,1,\text{Erasure}。其抹除機率與輸入無關,因此

P(Y=Erasure)=ρ.P(Y=\text{Erasure})=\rho.

非抹除的總機率為 1−ρ1-\rho。要使輸出資訊量最大,非抹除的 00、11 必須等機率出現,因此仍須令

P(X=0)=P(X=1)=12.P(X=0)=P(X=1)=\frac12.

所以,BSC 與 BEC 的容量最大化皆在輸入位元均勻分布時達成。


(b) BSC 與 BEC 的通道容量

BSC 容量

在均勻輸入下,BSC 的輸出也是均勻分布,因此

H(Y)=1.H(Y)=1.

給定輸入 XX 後,輸出錯誤的機率為 pp,故

H(Y∣X)=h2(p).H(Y|X)=h_2(p).

因此 BSC 容量為

CBSC=1−h2(p)\boxed{C_{\mathrm{BSC}}=1-h_2(p)}

亦即

CBSC=1+plog⁡2p+(1−p)log⁡2(1−p)\boxed{ C_{\mathrm{BSC}} =1+p\log_2p+(1-p)\log_2(1-p) }

BEC 容量

BEC 的輸出機率為

P(Y=0)=1−ρ2,P(Y=0)=\frac{1-\rho}{2}, P(Y=1)=1−ρ2,P(Y=1)=\frac{1-\rho}{2}, P(Y=Erasure)=ρ.P(Y=\text{Erasure})=\rho.

因此

H(Y)=−(1−ρ)log⁡21−ρ2−ρlog⁡2ρ.H(Y) =-(1-\rho)\log_2\frac{1-\rho}{2} -\rho\log_2\rho.

給定 XX 後,只有「正確接收」與「被抹除」兩種情況,所以

H(Y∣X)=h2(ρ).H(Y|X)=h_2(\rho).

故

I(X;Y)=H(Y)−H(Y∣X)=1−ρ.\begin{aligned} I(X;Y) &=H(Y)-H(Y|X)\\ &=1-\rho. \end{aligned}

因此 BEC 容量為

CBEC=1−ρ\boxed{C_{\mathrm{BEC}}=1-\rho}

其物理意義是:以機率 1−ρ1-\rho 成功傳送一個位元,以機率 ρ\rho 完全失去該位元。


(c) 容量比較

題目給定

0<p=ρ<0.5.0<p=\rho<0.5.

兩者容量分別為

CBSC=1−h2(p),C_{\mathrm{BSC}}=1-h_2(p), CBEC=1−p.C_{\mathrm{BEC}}=1-p.

對於 0<p<0.50<p<0.5,有

h2(p)>p.h_2(p)>p.

因此

1−p>1−h2(p),1-p>1-h_2(p),

得到

🔒

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

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

免費註冊

第 1 題10 分

Consider a signal xs(t)=x(t)s(t)x_s(t) = x(t)s(t), where x(t)x(t) is a bandlimited signal with bandwidth WBW_B Hz and s(t)s(t) is a periodic signal, expressed as s(t)=∑k=−∞∞Akδ(t−kT0)s(t) = \sum_{k=-\infty}^{\infty} A_k \delta(t - k T_0), where
Ak={1,if k is even0.8,if k is oddA_k = \begin{cases} 1, & \text{if } k \text{ is even} \\ 0.8, & \text{if } k \text{ is odd} \end{cases}
(a) Find the spectrum of s(t)s(t).
(b) What is the maximum value of T0T_0 such that x(t)x(t) can be fully recovered from xs(t)x_s(t) by passing xs(t)x_s(t) through a linear time-invariant filter?
(c) What is the impulse response of the linear filter in (b)?

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

這一題的完整詳解

核心觀念

本題考查:

  1. 週期脈衝串的傅立葉轉換。
  2. 時域相乘對應頻域摺積。
  3. 取樣頻譜的頻譜重複與避免混疊條件。
  4. 理想低通濾波器的設計與脈衝響應。

假設 x(t)x(t) 為低通限頻訊號,其頻譜滿足

X(f)=0,∣f∣>WBX(f)=0,\qquad |f|>W_B

因此原始訊號頻譜範圍為 [−WB,WB][-W_B,W_B]。


(a)求 s(t)s(t) 的頻譜

題目給定

s(t)=∑k=−∞∞Akδ(t−kT0)s(t)=\sum_{k=-\infty}^{\infty}A_k\delta(t-kT_0)

其中

Ak={1,k 為偶數0.8,k 為奇數A_k= \begin{cases} 1, & k\text{ 為偶數}\\ 0.8, & k\text{ 為奇數} \end{cases}

將 AkA_k 改寫成

Ak=0.9+0.1(−1)kA_k=0.9+0.1(-1)^k

因為:

  • kk 為偶數時,Ak=0.9+0.1=1A_k=0.9+0.1=1
  • kk 為奇數時,Ak=0.9−0.1=0.8A_k=0.9-0.1=0.8

所以

s(t)=0.9∑k=−∞∞δ(t−kT0)+0.1∑k=−∞∞(−1)kδ(t−kT0)s(t) = 0.9\sum_{k=-\infty}^{\infty}\delta(t-kT_0) + 0.1\sum_{k=-\infty}^{\infty}(-1)^k\delta(t-kT_0)

利用週期脈衝串的傅立葉轉換:

∑k=−∞∞δ(t−kT0)⟷1T0∑n=−∞∞δ(f−nT0)\sum_{k=-\infty}^{\infty}\delta(t-kT_0) \longleftrightarrow \frac{1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n}{T_0}\right)

而

(−1)k=ejπk(-1)^k=e^{j\pi k}

會使頻譜向右平移 1/(2T0)1/(2T_0),因此

∑k=−∞∞(−1)kδ(t−kT0)⟷1T0∑n=−∞∞δ(f−n+12T0)\sum_{k=-\infty}^{\infty}(-1)^k\delta(t-kT_0) \longleftrightarrow \frac{1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n+\frac12}{T_0}\right)

故

S(f)=0.9T0∑n=−∞∞δ(f−nT0)+0.1T0∑n=−∞∞δ(f−n+12T0)\boxed{ S(f)= \frac{0.9}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n}{T_0}\right) + \frac{0.1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n+\frac12}{T_0}\right) }

也可將頻率位置寫成

S(f)=0.9T0∑n=−∞∞δ(f−nT0)+0.1T0∑n=−∞∞δ(f−2n+12T0)\boxed{ S(f)= \frac{0.9}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{n}{T_0}\right) + \frac{0.1}{T_0} \sum_{n=-\infty}^{\infty} \delta\left(f-\frac{2n+1}{2T_0}\right) }

頻譜結構

s(t)s(t) 的週期為 2T02T_0,所以頻譜基本間隔為

12T0\frac{1}{2T_0}

其中:

  • 偶數次諧波位於 f=n/T0f=n/T_0,權重為 0.9/T00.9/T_0
  • 奇數次諧波位於 f=(2n+1)/(2T0)f=(2n+1)/(2T_0),權重為 0.1/T00.1/T_0

(b)T0T_0 的最大值

由時域相乘定理:

xs(t)=x(t)s(t)⟺Xs(f)=X(f)∗S(f)x_s(t)=x(t)s(t) \quad\Longleftrightarrow\quad X_s(f)=X(f)*S(f)

因此

Xs(f)=0.9T0∑n=−∞∞X(f−nT0)+0.1T0∑n=−∞∞X(f−n+12T0)X_s(f) = \frac{0.9}{T_0} \sum_{n=-\infty}^{\infty} X\left(f-\frac{n}{T_0}\right) + \frac{0.1}{T_0} \sum_{n=-\infty}^{\infty} X\left(f-\frac{n+\frac12}{T_0}\right)

原始頻譜 X(f)X(f) 位於 [−WB,WB][-W_B,W_B]。由於 S(f)S(f) 中最近的非零頻譜線位於

f=±12T0f=\pm\frac{1}{2T_0}
🔒

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

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

免費註冊

第 2 題10 分

Assume the DSB demodulation process is modeled as below, where the received signal is expressed as xr(t)=m(t)cos⁡(2πfct)+n(t)x_r(t) = m(t) \cos(2\pi f_c t) + n(t) with the noise n(t)n(t) being modeled as white, Gaussian noise with two-sided power spectral density of N0/2=1×10−10N_0/2 = 1 \times 10^{-10} W/Hz.

🖼️【此處有附圖,請對照原卷】
The figure shows a block diagram of a demodulator. The received signal xr(t)x_r(t) is first passed through a bandpass filter H1(f)H_1(f) centered at fcf_c. The output of the bandpass filter is then multiplied by cos⁡(2πfdt)\cos(2\pi f_d t). The result is then passed through a lowpass filter H2(f)H_2(f) to obtain the output y(t)y(t). The magnitude responses of H1(f)H_1(f) and H2(f)H_2(f) are given in the figure.

Assume the demodulation output y(t)y(t) can be decomposed into two parts, the signal part ys(t)y_s(t) and the noise part yn(t)y_n(t).

Case 1: fc=fd=200,000f_c = f_d = 200,000 Hz.
(a) Find ys(t)y_s(t).
(b) Find the total power of yn(t)y_n(t).

Case 2: fc=200,000f_c = 200,000 Hz but fd=200,100f_d = 200, 100 Hz.
(c) Plot the spectrum of ys(t)y_s(t).
(d) Find the auto-correlation of yn(t)y_n(t).

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

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

這一題的完整詳解

核心觀念

本題考查:

  1. DSB-SC 訊號的頻譜平移:
    xs(t)=m(t)cos⁡(2πfct)x_s(t)=m(t)\cos(2\pi f_ct)
    Xs(f)=12[M(f−fc)+M(f+fc)]X_s(f)=\frac{1}{2}\left[M(f-f_c)+M(f+f_c)\right]

  2. 乘上載波 cos⁡(2πfdt)\cos(2\pi f_dt) 所造成的頻譜平移:

F{v(t)cos⁡(2πfdt)}=12[V(f−fd)+V(f+fd)]\mathcal{F}\{v(t)\cos(2\pi f_dt)\} =\frac{1}{2}\left[V(f-f_d)+V(f+f_d)\right]
  1. 白色高斯雜訊通過濾波器後,其功率頻譜密度為
    Sn(f)=N02∣H(f)∣2S_n(f)=\frac{N_0}{2}|H(f)|^2

圖中資訊為:

  • M(f)M(f) 為三角形頻譜,支撐於 ∣f∣≤500|f|\leq 500 Hz,且 M(0)=1M(0)=1。
  • H1(f)H_1(f) 在 ±fc\pm f_c 附近通過,頻寬各為 10001000 Hz,增益為 22。
  • H2(f)H_2(f) 通過 ∣f∣≤500|f|\leq 500 Hz,增益為 11。
  • 雜訊雙邊 PSD:
    N02=10−10 W/Hz\frac{N_0}{2}=10^{-10}\ \text{W/Hz}

解題方法

先利用 H1(f)H_1(f) 保留 DSB 訊號的兩個旁帶,再乘上 cos⁡(2πfdt)\cos(2\pi f_dt) 將頻譜搬移至低頻,最後由 H2(f)H_2(f) 保留基頻部分。

對訊號而言,通過 H1(f)H_1(f) 後,由於 H1(f)H_1(f) 增益為 22,恰好抵消 DSB 調變頻譜中的係數 1/21/2:

X1(f)=M(f−fc)+M(f+fc)X_1(f)=M(f-f_c)+M(f+f_c)

再乘上 cos⁡(2πfdt)\cos(2\pi f_dt):

Ys(f)=12[X1(f−fd)+X1(f+fd)]Y_s(f)=\frac{1}{2}\left[X_1(f-f_d)+X_1(f+f_d)\right]


Case 1:fc=fd=200000f_c=f_d=200000 Hz

(a) 求 ys(t)y_s(t)

將 fc=fdf_c=f_d 代入:

Ys(f)=12[M(f)+M(f−2fc)+M(f+2fc)+M(f)]Y_s(f)=\frac{1}{2} \left[ M(f)+M(f-2f_c)+M(f+2f_c)+M(f) \right]

經過低通濾波器 H2(f)H_2(f) 後,只保留基頻部分:

Ys(f)=12[M(f)+M(f)]=M(f)Y_s(f)=\frac{1}{2}[M(f)+M(f)]=M(f)

因此

ys(t)=m(t)\boxed{y_s(t)=m(t)}

此時為理想同步解調,載波頻率完全匹配,所以原始訊號可完整恢復。


(b) 求 yn(t)y_n(t) 的總功率

H1(f)H_1(f) 的增益為 22,因此通過 H1(f)H_1(f) 後:

Sn1(f)=10−10×∣2∣2=4×10−10 W/HzS_{n1}(f)=10^{-10}\times |2|^2 =4\times10^{-10}\ \text{W/Hz}

在乘法器之後,每一個頻譜平移分量乘上 1/21/2,所以其 PSD 乘上 1/41/4。在低通濾波器通過的 ∣f∣≤500|f|\leq500 Hz 範圍內,正、負頻帶均搬移至基頻,因此共有兩個雜訊分量:

Syn(f)=2(14)(4×10−10)=2×10−10 W/HzS_{yn}(f) =2\left(\frac14\right) \left(4\times10^{-10}\right) =2\times10^{-10}\ \text{W/Hz}

低通頻寬為 10001000 Hz,因此雜訊總功率為

Pyn=∫−500500Syn(f) dfP_{yn} =\int_{-500}^{500}S_{yn}(f)\,df Pyn=(2×10−10)(1000)=2×10−7 WP_{yn} =(2\times10^{-10})(1000) =2\times10^{-7}\ \text{W}

因此

Pyn=2×10−7 W\boxed{P_{yn}=2\times10^{-7}\ \text{W}}


Case 2:fc=200000f_c=200000 Hz,fd=200100f_d=200100 Hz

此時解調頻率比載波頻率高 100100 Hz:

fd−fc=100 Hzf_d-f_c=100\ \text{Hz}

因此解調後的訊號頻譜會偏移 ±100\pm100 Hz。

(c) 繪出 ys(t)y_s(t) 的頻譜

低通範圍內的訊號頻譜為

Ys(f)=12[M(f−100)+M(f+100)],∣f∣≤500Y_s(f) =\frac12\left[M(f-100)+M(f+100)\right], \qquad |f|\leq500

其中

🔒

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

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

免費註冊

第 3 題13 分

Assume X(t) and Y(t) are two wide-sense stationary random processes with the autocorrelation functions RX(τ)=10e−2∣τ∣+2R_X(\tau) = 10e^{-2|\tau|} + 2 and RY(τ)=5e−10∣τ∣+5R_Y(\tau) = 5e^{-10|\tau|} + 5, respectively.

(a) Find the mean functions E[X(t)]E[X(t)] and E[Y(t)]E[Y(t)].
(b) Find the variance functions Var[X(t)]Var[X(t)] and Var[Y(t)]Var[Y(t)].
(c) Which random process has a wider bandwidth? Why?
(d) Which random process has a larger power? Why?
(e) Which random process is more like a Gaussian random process? Why?

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

這一題的完整詳解

核心觀念

本題考查寬平穩隨機程序(WSS)的自相關函數、平均功率、變異數、功率譜密度與頻寬,以及高斯隨機程序的判定。

對寬平穩隨機程序:

RX(τ)=E[X(t)X(t+τ)]R_X(\tau)=E[X(t)X(t+\tau)]

若自相關函數可分解為

RX(τ)=CX(τ)+mX2R_X(\tau)=C_X(\tau)+m_X^2

其中 mX=E[X(t)]m_X=E[X(t)],CX(τ)C_X(\tau) 為自共變異函數,則

CX(τ)=RX(τ)−mX2C_X(\tau)=R_X(\tau)-m_X^2

且

Var⁡[X(t)]=CX(0)=RX(0)−mX2\operatorname{Var}[X(t)]=C_X(0)=R_X(0)-m_X^2

平均功率為

PX=E[X2(t)]=RX(0)P_X=E[X^2(t)]=R_X(0)

解題方法

題目中的指數項會隨著 ∣τ∣|\tau| 增大而趨近於零,因此常數項代表平均值平方:

mX2=lim⁡∣τ∣→∞RX(τ)m_X^2=\lim_{|\tau|\to\infty}R_X(\tau) mY2=lim⁡∣τ∣→∞RY(τ)m_Y^2=\lim_{|\tau|\to\infty}R_Y(\tau)

採用隨機程序平均值為非負的慣例,取正平方根。

(a) 平均函數

對 X(t)X(t):

RX(τ)=10e−2∣τ∣+2R_X(\tau)=10e^{-2|\tau|}+2

因此

mX2=lim⁡∣τ∣→∞RX(τ)=2m_X^2=\lim_{|\tau|\to\infty}R_X(\tau)=2 E[X(t)]=mX=2E[X(t)]=m_X=\sqrt{2}

對 Y(t)Y(t):

RY(τ)=5e−10∣τ∣+5R_Y(\tau)=5e^{-10|\tau|}+5

因此

mY2=lim⁡∣τ∣→∞RY(τ)=5m_Y^2=\lim_{|\tau|\to\infty}R_Y(\tau)=5 E[Y(t)]=mY=5E[Y(t)]=m_Y=\sqrt{5}

嚴格而言,自相關函數只能決定平均值的平方,因此平均值也可能是負值;本題依一般教科書慣例取非負平均值。


(b) 變異數函數

由於兩個程序皆為寬平穩,平均值與自相關函數均與時間 tt 無關,所以變異數也是常數。

對 X(t)X(t):

RX(0)=10e0+2=12R_X(0)=10e^0+2=12 Var⁡[X(t)]=RX(0)−mX2=12−2=10\operatorname{Var}[X(t)] =R_X(0)-m_X^2 =12-2 =10

對 Y(t)Y(t):

RY(0)=5e0+5=10R_Y(0)=5e^0+5=10 Var⁡[Y(t)]=RY(0)−mY2=10−5=5\operatorname{Var}[Y(t)] =R_Y(0)-m_Y^2 =10-5 =5

因此

Var⁡[X(t)]=10\boxed{\operatorname{Var}[X(t)]=10} Var⁡[Y(t)]=5\boxed{\operatorname{Var}[Y(t)]=5}

(c) 哪一個隨機程序具有較寬的頻寬?

先扣除平均值所造成的常數項,得到自共變異函數:

CX(τ)=10e−2∣τ∣C_X(\tau)=10e^{-2|\tau|} CY(τ)=5e−10∣τ∣C_Y(\tau)=5e^{-10|\tau|}

利用傅立葉轉換公式:

F{e−a∣τ∣}=2aa2+ω2\mathcal{F}\{e^{-a|\tau|}\} =\frac{2a}{a^2+\omega^2}

可得連續部分的功率譜密度:

SX(ω)=10⋅44+ω2=404+ω2S_X(\omega)=10\cdot\frac{4}{4+\omega^2} =\frac{40}{4+\omega^2}
🔒

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

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

免費註冊

第 4 題10 分

Consider a bandlimited communication system with the overall response of the transmitter, channel, and receiver filters being represented by the pulse shaping waveform p(t)p(t). The Nyquist pulse shaping criterion
∑n=−∞∞P(f+nT)=T\sum_{n=-\infty}^{\infty} P\left(f + \frac{n}{T}\right) = T
provides a rule to design p(t)p(t) to avoid inter-symbol interference (ISI) in the bandlimited system, where TT is the sampling interval and P(f)P(f) is the Fourier transform of the overall pulse p(t)p(t).

(a) A commonly used pulse shaping waveform is the raised cosine pulse defined by
P(f)={a,∣f∣≤W−βb+ccos⁡(π(∣f∣−W+β)2β),W−β≤∣f∣≤W+β0,∣f∣≥W+βP(f) = \begin{cases} a, & |f| \le W - \beta \\ b + c \cos\left(\frac{\pi( |f|-W+\beta)}{2\beta}\right), & W - \beta \le |f| \le W + \beta \\ 0, & |f| \ge W + \beta \end{cases}
where W=12TW = \frac{1}{2T} and β\beta is the excessive bandwidth with 0≤β≤W0 \le \beta \le W. Find the constants a,b,ca, b, c and dd such that the above raised cosine spectrum satisfies the Nyquist pulse shaping criterion.

(b) Consider another pulse shaping spectrum P(f)P(f) shown below, where the vertical axis of the plot corresponds to the value of 1TP(f)\frac{1}{T}P(f) and the labeled frequencies correspond to the corner points of the spectrum. Find one signaling interval TT such that the bandlimited system with the pulse p(t)p(t) results in zero ISI.

🖼️【此處有附圖,請對照原卷】
The figure shows the spectrum of a pulse shaping filter. The horizontal axis represents frequency ff, and the vertical axis represents 1TP(f)\frac{1}{T}P(f). The spectrum is non-zero for ∣f∣≤7/8|f| \le 7/8. It has a flat top for ∣f∣≤1/8|f| \le 1/8, a cosine-shaped roll-off from ∣f∣=1/8|f|=1/8 to ∣f∣=7/8|f|=7/8, and is zero elsewhere. The plot shows constant values a′a', b′b', c′c' on the vertical axis corresponding to different frequency ranges.

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

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

這一題的完整詳解

此題考驗對 Nyquist 準則、 Raised Cosine 濾波器頻譜的理解,以及如何從給定的頻譜圖推導出抽樣間隔 TT。

核心觀念:

  1. Nyquist 準則:避免符號間干擾 (ISI) 的充分必要條件。
  2. Raised Cosine 濾波器:一種常見的 Raised Cosine 濾波器頻譜定義。
  3. 傅立葉轉換與頻譜的性質。
  4. 頻譜圖的解讀。

解題步驟:

Nyquist Pulse Shaping Criterion:
∑n=−∞∞P(f+nT)=T\sum_{n=-\infty}^{\infty} P\left(f + \frac{n}{T}\right) = T
其中 W=12TW = \frac{1}{2T},所以 1T=2W\frac{1}{T} = 2W。
∑n=−∞∞P(f+n(2W))=T\sum_{n=-\infty}^{\infty} P(f + n(2W)) = T

Raised Cosine Pulse Spectrum Definition:
P(f)={a,∣f∣≤W−βb+ccos⁡(π(∣f∣−W+β)2β),W−β≤∣f∣≤W+β0,∣f∣≥W+βP(f) = \begin{cases} a, & |f| \le W - \beta \\ b + c \cos\left(\frac{\pi( |f|-W+\beta)}{2\beta}\right), & W - \beta \le |f| \le W + \beta \\ 0, & |f| \ge W + \beta \end{cases}
其中 0≤β≤W0 \le \beta \le W。

(a) 尋找 a,b,ca, b, c

根據 Nyquist 準則,在 f=0f=0 處:
∑n=−∞∞P(n(2W))=T\sum_{n=-\infty}^{\infty} P(n(2W)) = T。
由於 P(f)P(f) 在 ∣f∣≥W+β|f| \ge W+\beta 為零,且 W+β≤2WW+\beta \le 2W (因為 β≤W\beta \le W),所以只有 n=0n=0 的項 P(0)P(0) 是非零的。
P(0)=TP(0) = T。
根據定義,P(0)=aP(0) = a (因為 ∣0∣≤W−β|0| \le W-\beta 成立)。
所以,a=T=12Wa = T = \frac{1}{2W}。

在 f=Wf=W 處:
∑n=−∞∞P(W+n(2W))=T\sum_{n=-\infty}^{\infty} P(W + n(2W)) = T。
考慮 n=0n=0 (P(W)P(W)) 和 n=−1n=-1 (P(W−2W)=P(−W)P(W - 2W) = P(-W))。
由於 P(f)P(f) 是偶函數,P(−W)=P(W)P(-W)=P(W)。
P(W)+P(−W)=2P(W)=TP(W) + P(-W) = 2 P(W) = T。
WW 落在 W−β≤W≤W+βW-\beta \le W \le W+\beta 的範圍內(除非 β=0\beta=0)。

🔒

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

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

免費註冊

第 5 題24 分

Consider a communication system with the vector representation (all with dimension N×1N \times 1):
r=s+n\mathbf{r} = \mathbf{s} + \mathbf{n}
where r\mathbf{r} is the received signal, s∈{S1,S2,…,SM}\mathbf{s} \in \{S_1, S_2, \dots, S_M\} is the transmitted symbol and the noise n\mathbf{n} has i.i.d. Gaussian components each with zero mean and variance N0/2N_0/2. The transmitted signal s\mathbf{s} and the noise n\mathbf{n} are assumed to be independent. All the possible transmitted symbols S1,S2,…,SMS_1, S_2, \dots, S_M have equal prior probability, but may not have equal energy. The signaling interval for each symbol s\mathbf{s} is TT seconds.

Maximum-likelihood (ML) decision rule is often employed at the receiver for symbol decoding and can be represented by
S^ML=arg⁡max⁡smp(r∣s=sm)\hat{S}_{ML} = \arg \max_{\mathbf{s}_m} p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m)
where p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) is the likelihood function associated with sm\mathbf{s}_m.

(a) What is the average signal energy of s\mathbf{s}? What is the transmission bit rate of this system?
(b) Find the likelihood function p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m). (Note: You can express the likelihood function using vector 2-norm ∣∣⋅∣∣2||\cdot||_2, or define appropriate notations by yourself. The probability density function of a N(0,σ2)N(0, \sigma^2) Gaussian random variable ZZ is pZ(z)=12πσ2e−z22σ2p_Z(z) = \frac{1}{\sqrt{2\pi\sigma^2}} e^{-\frac{z^2}{2\sigma^2}}.)
(c) A decoding rule based on correlation of the received signal r\mathbf{r} and each sm\mathbf{s}_m for m=1,2,…,Mm=1, 2, \dots, M can be devised as follows:
"If rTsm≥rTsn\mathbf{r}^T \mathbf{s}_m \ge \mathbf{r}^T \mathbf{s}_n for all n≠mn \neq m, then the decoder decides that sm\mathbf{s}_m is the transmitted symbol."
Is this decoder optimum in terms of symbol error probability? Please explain your answers.
(d) One vector representation of QPSK signaling is designed as S1=[2,0]TS_1 = [\sqrt{2}, 0]^T, S2=[0,2]TS_2 = [0, \sqrt{2}]^T, S3=[2,22]TS_3 = [\sqrt{2}, 2\sqrt{2}]^T and S4=[22,2]TS_4 = [2\sqrt{2}, \sqrt{2}]^T with allocation of two information bits to the four symbols being arranged as follows: (00)→S1(00) \rightarrow S_1, (01)→S2(01) \rightarrow S_2, (10)→S3(10) \rightarrow S_3 and (11)→S4(11) \rightarrow S_4.
i. Given that the transmitted symbol is S1S_1. Find the probability that the ML decoder erroneously gives the decision S^ML=S4\hat{S}_{ML} = S_4 when N0=2N_0 = 2 (Express your answer either in terms of the Gaussian Q-function or the erfc function).
ii. What would be wrong with this design? What would you do to improve any of the system performance in a reasonable way? Please explain.

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

這一題的完整詳解

此題考驗對數位通訊系統中最大概似接收器(ML receiver)、符號錯誤機率、最佳解碼準則、QPSK 調變、以及高斯雜訊下的錯誤機率計算。

核心觀念:

  1. 最大概似接收器 (ML Receiver) 的原理。
  2. 高斯雜訊下的似然函數 (Likelihood Function)。
  3. Correlation Receiver 與 ML Receiver 的關係。
  4. Symbol Error Probability (SEP) 的計算。
  5. QPSK 調變的符號表示與能量。
  6. 高斯 Q-function 和 erfc function。

解題步驟:

系統模型:r=s+n\mathbf{r} = \mathbf{s} + \mathbf{n}
s∈{S1,…,SM}\mathbf{s} \in \{S_1, \dots, S_M\} 是傳送符號,具有相等先驗機率 P(Si)=1/MP(S_i) = 1/M。
n\mathbf{n} 是零均值、對角協方差矩陣為 (N0/2)IN(N_0/2)I_N 的高斯雜訊。
N0/2N_0/2 是單邊功率譜密度。

(a) 平均信號能量與傳輸位元率

平均信號能量 EsE_s 定義為所有可能符號能量的平均值:
Es=E[∣∣s∣∣22]=1M∑m=1M∣∣Sm∣∣22E_s = E[||\mathbf{s}||_2^2] = \frac{1}{M} \sum_{m=1}^M ||S_m||_2^2
題目提到符號能量可能不相等。

傳輸位元率 RbR_b:
假設每個符號 SmS_m 代表 kk 個位元。則 M=2kM = 2^k。
傳輸速率(符號速率)是 Rs=1/TR_s = 1/T (符號/秒)。
傳輸位元率 Rb=k×Rs=k/TR_b = k \times R_s = k/T (位元/秒)。
題目沒有給出 MM 或 kk,也沒有給出 TT。
但是,題目有提到 TT 秒是信號間隔。
所以 Rs=1/TR_s = 1/T。
如果 MM 個符號代表 kk 個位元,則 M=2kM = 2^k。
位元率 Rb=k×(1/T)R_b = k \times (1/T)。
如果 MM 和 kk 未知,則無法計算。
重新閱讀題目: 題目沒有給出 MM 或 kk。所以這部分可能需要用符號表示。
平均信號能量: Es=1M∑m=1M∣∣Sm∣∣22E_s = \frac{1}{M} \sum_{m=1}^M ||S_m||_2^2。
傳輸位元率: Rb=log⁡2MTR_b = \frac{\log_2 M}{T} bits/sec。

答案 (a):
平均信號能量:Es=1M∑m=1M∣∣Sm∣∣22E_s = \frac{1}{M} \sum_{m=1}^M ||S_m||_2^2。
傳輸位元率:Rb=log⁡2MTR_b = \frac{\log_2 M}{T}。

(b) 尋找似然函數 p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m)

似然函數是給定傳送符號 sm\mathbf{s}_m 時,接收到信號 r\mathbf{r} 的條件機率密度函數。
我們有 r=sm+n\mathbf{r} = \mathbf{s}_m + \mathbf{n}。
雜訊 n\mathbf{n} 是 NN 維高斯向量,其各分量獨立,均值為 0,方差為 N0/2N_0/2。
所以 n∼N(0,N02IN)\mathbf{n} \sim N(\mathbf{0}, \frac{N_0}{2} I_N)。
則 r−sm=n\mathbf{r} - \mathbf{s}_m = \mathbf{n}。
r−sm\mathbf{r} - \mathbf{s}_m 是一個均值為 sm\mathbf{s}_m、協方差為 N02IN\frac{N_0}{2} I_N 的高斯向量。
一個 NN 維高斯隨機向量 x\mathbf{x},其均值為 μ\boldsymbol{\mu},協方差矩陣為 Σ\Sigma,其 PDF 為:
p(x)=1(2π)Ndet⁡(Σ)exp⁡(−12(x−μ)TΣ−1(x−μ))p(\mathbf{x}) = \frac{1}{\sqrt{(2\pi)^N \det(\Sigma)}} \exp\left(-\frac{1}{2}(\mathbf{x}-\boldsymbol{\mu})^T \Sigma^{-1} (\mathbf{x}-\boldsymbol{\mu})\right)
在這裡,x=r\mathbf{x} = \mathbf{r},μ=sm\boldsymbol{\mu} = \mathbf{s}_m,Σ=N02IN\Sigma = \frac{N_0}{2} I_N。
det⁡(Σ)=(N02)N\det(\Sigma) = (\frac{N_0}{2})^N。
Σ−1=2N0IN\Sigma^{-1} = \frac{2}{N_0} I_N。
所以,似然函數為:
p(r∣s=sm)=1(2π)N(N02)Nexp⁡(−12(r−sm)T(2N0IN)(r−sm))p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{\sqrt{(2\pi)^N (\frac{N_0}{2})^N}} \exp\left(-\frac{1}{2}(\mathbf{r}-\mathbf{s}_m)^T \left(\frac{2}{N_0} I_N\right) (\mathbf{r}-\mathbf{s}_m)\right)
p(r∣s=sm)=1(N02)N/2(2π)Nexp⁡(−1N0(r−sm)T(r−sm))p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\frac{N_0}{2})^{N/2} \sqrt{(2\pi)^N}} \exp\left(-\frac{1}{N_0}(\mathbf{r}-\mathbf{s}_m)^T (\mathbf{r}-\mathbf{s}_m)\right)
p(r∣s=sm)=1(N02)N/2(2π)N/2exp⁡(−1N0∣∣r−sm∣∣22)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\frac{N_0}{2})^{N/2} (2\pi)^{N/2}} \exp\left(-\frac{1}{N_0} ||\mathbf{r}-\mathbf{s}_m||_2^2\right)
p(r∣s=sm)=1(πN0)N/2exp⁡(−1N0∣∣r−sm∣∣22)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\pi N_0)^{N/2}} \exp\left(-\frac{1}{N_0} ||\mathbf{r}-\mathbf{s}_m||_2^2\right)
其中 ∣∣r−sm∣∣22=(r−sm)T(r−sm)=∑i=1N(ri−sm,i)2||\mathbf{r}-\mathbf{s}_m||_2^2 = (\mathbf{r}-\mathbf{s}_m)^T (\mathbf{r}-\mathbf{s}_m) = \sum_{i=1}^N (r_i - s_{m,i})^2。

ML 解碼器是最大化這個似然函數。由於指數函數是單調遞增的,最大化似然函數等價於最小化指數的指數項。
最大化 p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) 等價於最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
s^ML=arg⁡min⁡sm∣∣r−sm∣∣22\hat{\mathbf{s}}_{ML} = \arg \min_{\mathbf{s}_m} ||\mathbf{r}-\mathbf{s}_m||_2^2
這也稱為最小距離準則 (Minimum Distance Criterion)。

答案 (b):
p(r∣s=sm)=1(πN0)N/2exp⁡(−1N0∣∣r−sm∣∣22)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) = \frac{1}{(\pi N_0)^{N/2}} \exp\left(-\frac{1}{N_0} ||\mathbf{r}-\mathbf{s}_m||_2^2\right)。

(c) 相關性解碼器是否最佳(在符號錯誤機率方面)

相關性解碼器規則:
若 rTsm≥rTsn\mathbf{r}^T \mathbf{s}_m \ge \mathbf{r}^T \mathbf{s}_n 對所有 n≠mn \neq m,則選擇 sm\mathbf{s}_m。
我們知道 ML 解碼器是最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
∣∣r−sm∣∣22=(r−sm)T(r−sm)=rTr−2rTsm+smTsm||\mathbf{r}-\mathbf{s}_m||_2^2 = (\mathbf{r}-\mathbf{s}_m)^T (\mathbf{r}-\mathbf{s}_m) = \mathbf{r}^T \mathbf{r} - 2\mathbf{r}^T \mathbf{s}_m + \mathbf{s}_m^T \mathbf{s}_m
∣∣r−sm∣∣22=∣∣r∣∣22−2rTsm+∣∣sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 = ||\mathbf{r}||_2^2 - 2\mathbf{r}^T \mathbf{s}_m + ||\mathbf{s}_m||_2^2
最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 等價於最大化 2rTsm−∣∣sm∣∣222\mathbf{r}^T \mathbf{s}_m - ||\mathbf{s}_m||_2^2。
或者,如果所有符號的能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 都相等,則最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 等價於最大化 rTsm\mathbf{r}^T \mathbf{s}_m。
此時,ML 解碼器就變成了相關性解碼器。

然而,題目明確指出「符號能量可能不相等」。
如果符號能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 不相等,那麼最大化 rTsm\mathbf{r}^T \mathbf{s}_m 並不一定等同於最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
ML 解碼器是最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
相關性解碼器是最大化 rTsm\mathbf{r}^T \mathbf{s}_m。

考慮以下情況:
令 r=[1,0]T\mathbf{r} = [1, 0]^T。
令 S1=[1,0]TS_1 = [1, 0]^T,則 ∣∣s1∣∣22=1||\mathbf{s}_1||_2^2 = 1。rTs1=1\mathbf{r}^T \mathbf{s}_1 = 1。∣∣r−s1∣∣22=0||\mathbf{r}-\mathbf{s}_1||_2^2 = 0。
令 S2=[2,0]TS_2 = [2, 0]^T,則 ∣∣s2∣∣22=4||\mathbf{s}_2||_2^2 = 4。rTs2=2\mathbf{r}^T \mathbf{s}_2 = 2。∣∣r−s2∣∣22=∣∣[1,0]T−[2,0]T∣∣22=∣∣[−1,0]T∣∣22=1||\mathbf{r}-\mathbf{s}_2||_2^2 = ||[1,0]^T - [2,0]^T||_2^2 = ||[-1,0]^T||_2^2 = 1。
ML 解碼器會選擇 S1S_1 (因為 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 最小,為 0)。
相關性解碼器會選擇 S2S_2 (因為 rTsm\mathbf{r}^T \mathbf{s}_m 最大,2 > 1)。
在這種情況下,相關性解碼器不是最佳的。

解釋:
ML 解碼器選擇 sm\mathbf{s}_m 使得 p(r∣s=sm)p(\mathbf{r}|\mathbf{s} = \mathbf{s}_m) 最大。
這等價於最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2。
∣∣r−sm∣∣22=∣∣r∣∣22−2rTsm+∣∣sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 = ||\mathbf{r}||_2^2 - 2\mathbf{r}^T \mathbf{s}_m + ||\mathbf{s}_m||_2^2。
ML 解碼器選擇 sm\mathbf{s}_m 使得 rTsm−12∣∣sm∣∣22\mathbf{r}^T \mathbf{s}_m - \frac{1}{2}||\mathbf{s}_m||_2^2 最大。
相關性解碼器選擇 sm\mathbf{s}_m 使得 rTsm\mathbf{r}^T \mathbf{s}_m 最大。

只有當所有符號能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 都相等時,最大化 rTsm\mathbf{r}^T \mathbf{s}_m 才等價於最大化 rTsm−12∣∣sm∣∣22\mathbf{r}^T \mathbf{s}_m - \frac{1}{2}||\mathbf{s}_m||_2^2。
由於題目明確指出符號能量可能不相等,所以相關性解碼器不是最佳的。

答案 (c):
否。這個解碼器不是最佳的。
解釋:最大概似解碼器是最小化 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2,這等價於最大化 rTsm−12∣∣sm∣∣22\mathbf{r}^T \mathbf{s}_m - \frac{1}{2}||\mathbf{s}_m||_2^2。相關性解碼器則是最大化 rTsm\mathbf{r}^T \mathbf{s}_m。只有當所有符號的能量 ∣∣sm∣∣22||\mathbf{s}_m||_2^2 都相等時,這兩者才等價。由於題目說明符號能量可能不相等,因此相關性解碼器不保證是最佳的(即能獲得最小的符號錯誤機率)。

(d) QPSK 符號表示與錯誤機率計算

QPSK 符號:
S1=[2,0]TS_1 = [\sqrt{2}, 0]^T (00)
S2=[0,2]TS_2 = [0, \sqrt{2}]^T (01)
S3=[2,22]TS_3 = [\sqrt{2}, 2\sqrt{2}]^T (10)
S4=[22,2]TS_4 = [2\sqrt{2}, \sqrt{2}]^T (11)
N=2N=2 (向量維度)。
N0=2N_0 = 2。所以 N0/2=1N_0/2 = 1 (雜訊方差)。

(d) i. P(S^ML=S4∣S1 transmitted)P(\hat{S}_{ML} = S_4 | S_1 \text{ transmitted})

ML 解碼器選擇 sm\mathbf{s}_m 使得 ∣∣r−sm∣∣22||\mathbf{r}-\mathbf{s}_m||_2^2 最小。
我們知道 r=s1+n\mathbf{r} = \mathbf{s}_1 + \mathbf{n}。
ML 解碼器錯誤地選擇 S4S_4 而不是 S1S_1 當:
∣∣r−s4∣∣22<∣∣r−s1∣∣22||\mathbf{r}-\mathbf{s}_4||_2^2 < ||\mathbf{r}-\mathbf{s}_1||_2^2
代入 r=s1+n\mathbf{r} = \mathbf{s}_1 + \mathbf{n}:
∣∣s1+n−s4∣∣22<∣∣s1+n−s1∣∣22||\mathbf{s}_1 + \mathbf{n} - \mathbf{s}_4||_2^2 < ||\mathbf{s}_1 + \mathbf{n} - \mathbf{s}_1||_2^2
∣∣n−(s4−s1)∣∣22<∣∣n∣∣22||\mathbf{n} - (\mathbf{s}_4 - \mathbf{s}_1)||_2^2 < ||\mathbf{n}||_2^2

計算向量差:
s4−s1=[22,2]T−[2,0]T=[2,2]T\mathbf{s}_4 - \mathbf{s}_1 = [2\sqrt{2}, \sqrt{2}]^T - [\sqrt{2}, 0]^T = [\sqrt{2}, \sqrt{2}]^T。
令 Δ41=s4−s1=[2,2]T\Delta_{41} = \mathbf{s}_4 - \mathbf{s}_1 = [\sqrt{2}, \sqrt{2}]^T。
不等式變為:
∣∣n−Δ41∣∣22<∣∣n∣∣22||\mathbf{n} - \Delta_{41}||_2^2 < ||\mathbf{n}||_2^2
(n−Δ41)T(n−Δ41)<nTn(\mathbf{n} - \Delta_{41})^T (\mathbf{n} - \Delta_{41}) < \mathbf{n}^T \mathbf{n}
nTn−2nTΔ41+Δ41TΔ41<nTn\mathbf{n}^T \mathbf{n} - 2 \mathbf{n}^T \Delta_{41} + \Delta_{41}^T \Delta_{41} < \mathbf{n}^T \mathbf{n}
−2nTΔ41+∣∣Δ41∣∣22<0- 2 \mathbf{n}^T \Delta_{41} + ||\Delta_{41}||_2^2 < 0
∣∣Δ41∣∣22<2nTΔ41||\Delta_{41}||_2^2 < 2 \mathbf{n}^T \Delta_{41}

計算 ∣∣Δ41∣∣22||\Delta_{41}||_2^2:
∣∣Δ41∣∣22=(2)2+(2)2=2+2=4||\Delta_{41}||_2^2 = (\sqrt{2})^2 + (\sqrt{2})^2 = 2 + 2 = 4。
所以,4<2nTΔ414 < 2 \mathbf{n}^T \Delta_{41}。
2<nTΔ412 < \mathbf{n}^T \Delta_{41}。
nTΔ41=n1(2)+n2(2)=2(n1+n2)\mathbf{n}^T \Delta_{41} = n_1 (\sqrt{2}) + n_2 (\sqrt{2}) = \sqrt{2} (n_1 + n_2)。
所以,2<2(n1+n2)2 < \sqrt{2} (n_1 + n_2)。
2<n1+n2\sqrt{2} < n_1 + n_2
令 Z=n1+n2Z = n_1 + n_2。

🔒

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

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

免費註冊

其他考古題