115 一般警察特考二等 刑事警察人員類別犯罪分析組 數位訊號處理(DSP) 考古題詳解

115 年・一般警察特考・二等・刑事警察人員類別犯罪分析組
試題代碼 數位訊號處理(DSP)

第 一 題25 分

假設一個週期為 TT 的連續時間訊號 x(t)x(t),其基本角頻率(radian frequency)為

ω0=2πT\omega_0=\frac{2\pi}{T}

該訊號可以表示為複數指數形式的傅立葉級數(Fourier Series):

x(t)=∑k=−∞∞akejω0ktx(t)=\sum_{k=-\infty}^{\infty}a_ke^{j\omega_0kt}

(一)請推導出上式之傅立葉係數(Fourier coefficients){ak}\{a_k\} 計算公式為:(15 分)

ak=1T∫0Tx(t)e−jω0kt dta_k=\frac{1}{T}\int_0^T x(t)e^{-j\omega_0kt}\,dt

(二)當訊號 x(t)=sin⁡3(4πt)x(t)=\sin^3(4\pi t),請計算出該訊號之傅立葉係數 {ak}\{a_k\}。(10 分)

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

這一題的完整詳解

核心觀念

複數指數傅立葉級數以 ejω0kte^{j\omega_0kt} 作為基底,其中 ω0=2π/T\omega_0=2\pi/T。利用不同基底在一個週期上的正交性,可以從級數中分離出特定索引 kk 的係數:

1T∫0Tejω0(n−k)t dt={1,n=k,0,n≠k.\frac{1}{T}\int_0^T e^{j\omega_0(n-k)t}\,dt = \begin{cases} 1, & n=k,\\ 0, & n\ne k. \end{cases}

解題方法

(一)從傅立葉級數定義出發:

x(t)=∑n=−∞∞anejω0nt.x(t)=\sum_{n=-\infty}^{\infty}a_ne^{j\omega_0nt}.

兩邊同乘 e−jω0kte^{-j\omega_0kt},並在一個週期 [0,T][0,T] 上積分:

∫0Tx(t)e−jω0kt dt=∫0T(∑n=−∞∞anejω0nt)e−jω0kt dt.\int_0^T x(t)e^{-j\omega_0kt}\,dt = \int_0^T\left(\sum_{n=-\infty}^{\infty}a_ne^{j\omega_0nt}\right)e^{-j\omega_0kt}\,dt.

交換加總與積分後:

∫0Tx(t)e−jω0kt dt=∑n=−∞∞an∫0Tejω0(n−k)t dt.\int_0^T x(t)e^{-j\omega_0kt}\,dt = \sum_{n=-\infty}^{\infty}a_n\int_0^T e^{j\omega_0(n-k)t}\,dt.

當 n=kn=k 時,積分為 TT。當 n≠kn\ne k 時:

∫0Tejω0(n−k)t dt=ejω0(n−k)tjω0(n−k)∣0T=ej2π(n−k)−1jω0(n−k)=0,\int_0^T e^{j\omega_0(n-k)t}\,dt = \left.\frac{e^{j\omega_0(n-k)t}}{j\omega_0(n-k)}\right|_0^T = \frac{e^{j2\pi(n-k)}-1}{j\omega_0(n-k)} =0,

因為 n−kn-k 是非零整數,ej2π(n−k)=1e^{j2\pi(n-k)}=1。因此加總中只留下 n=kn=k 的項:

∫0Tx(t)e−jω0kt dt=akT.\int_0^T x(t)e^{-j\omega_0kt}\,dt=a_kT.

除以 TT,得到傅立葉係數公式:

ak=1T∫0Tx(t)e−jω0kt dt.\boxed{a_k=\frac{1}{T}\int_0^T x(t)e^{-j\omega_0kt}\,dt}.

(二)先利用三倍角恆等式展開訊號:

sin⁡3θ=3sin⁡θ−sin⁡3θ4,\sin^3\theta=\frac{3\sin\theta-\sin3\theta}{4},

所以

🔒

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

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

免費註冊

第 二 題25 分

一個簡單的移動平均(moving average)濾波器系統定義如下:

y[n]=1M∑k=0M−1x[n−k]y[n]=\frac{1}{M}\sum_{k=0}^{M-1}x[n-k]

(一)求該系統之脈衝響應(impulse response)h[n]h[n]。(5 分)

(二)請推導該系統之頻率響應(frequency response)H(ejω)H(e^{j\omega}) 表示式。(10 分)

(三)假設濾波器 M=11M=11,若輸入訊號 x[n]=ejωnx[n]=e^{j\omega n},請找出 ω\omega 值,−π≤ω≤π-\pi\leq\omega\leq\pi,使得 y[n]=0y[n]=0,∀n\forall n 成立。(10 分)

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

這一題的完整詳解

核心觀念

移動平均濾波器是有限長度的線性時不變系統(LTI system),會將目前輸入樣本與前面 M−1M-1 個樣本等權平均。本題使用的基本關係為:

  • 脈衝響應:h[n]h[n] 是系統對單位脈衝 δ[n]\delta[n] 的輸出。
  • 頻率響應:H(ejω)=∑n=−∞∞h[n]e−jωnH(e^{j\omega})=\sum_{n=-\infty}^{\infty}h[n]e^{-j\omega n}。
  • 複指數 ejωne^{j\omega n} 是 LTI 系統的特徵函數,輸出為 H(ejω)ejωnH(e^{j\omega})e^{j\omega n}。

解題方法

(一)脈衝響應 h[n]h[n]

令輸入為單位脈衝 x[n]=δ[n]x[n]=\delta[n],代入系統定義:

h[n]=y[n]=1M∑k=0M−1δ[n−k]h[n]=y[n] =\frac{1}{M}\sum_{k=0}^{M-1}\delta[n-k]

當 n=0,1,…,M−1n=0,1,\ldots,M-1 時,和式中恰有一項為 11;其他 nn 則所有項皆為 00。因此:

h[n]={1M,0≤n≤M−10,其他h[n]= \begin{cases} \frac{1}{M}, & 0\leq n\leq M-1\\ 0, & \text{其他} \end{cases}

亦可寫成:

h[n]=1M∑k=0M−1δ[n−k]h[n]=\frac{1}{M}\sum_{k=0}^{M-1}\delta[n-k]

(二)頻率響應 H(ejω)H(e^{j\omega})

由脈衝響應的離散時間傅立葉轉換:

H(ejω)=∑n=−∞∞h[n]e−jωn=1M∑n=0M−1e−jωnH(e^{j\omega}) =\sum_{n=-\infty}^{\infty}h[n]e^{-j\omega n} =\frac{1}{M}\sum_{n=0}^{M-1}e^{-j\omega n}

這是等比級數。當 e−jω≠1e^{-j\omega}\neq 1 時,利用等比級數公式:

H(ejω)=1M1−e−jMω1−e−jωH(e^{j\omega}) =\frac{1}{M}\frac{1-e^{-jM\omega}}{1-e^{-j\omega}}

分子與分母各以半角形式整理:

1−e−jMω=e−jMω/2(ejMω/2−e−jMω/2)=2je−jMω/2sin⁡(Mω2)1-e^{-jM\omega} =e^{-jM\omega/2}\left(e^{jM\omega/2}-e^{-jM\omega/2}\right) =2j e^{-jM\omega/2}\sin\left(\frac{M\omega}{2}\right) 1−e−jω=2je−jω/2sin⁡(ω2)1-e^{-j\omega} =2j e^{-j\omega/2}\sin\left(\frac{\omega}{2}\right)

所以頻率響應也可表示為:

H(ejω)=1Me−jω(M−1)/2sin⁡(Mω/2)sin⁡(ω/2)H(e^{j\omega}) =\frac{1}{M}e^{-j\omega(M-1)/2} \frac{\sin(M\omega/2)}{\sin(\omega/2)}

此式在 ω=2πr\omega=2\pi r(rr 為整數)時呈現 0/00/0,須回到原始和式求值;此時 H(ejω)=1H(e^{j\omega})=1。

(三)M=11M=11 時使輸出為零的 ω\omega

輸入為 x[n]=ejωnx[n]=e^{j\omega n} 時,直接代入系統:

🔒

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

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

免費註冊

第 三 題25 分

考慮下面的差分方程式的系統:

y[n]=y[n−1]−y[n−2]+2x[n]+2x[n−1]y[n]=y[n-1]-y[n-2]+2x[n]+2x[n-1]

(一)找出系統函數(System function)

H(z)=Y(z)X(z)H(z)=\frac{Y(z)}{X(z)}

以 zz 的正冪次方多項式表示。(10 分)

(二)求該系統的脈衝響應(impulse response)h[n]h[n],以純實數表示。(15 分)

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

這一題的完整詳解

核心觀念

本題考查以差分方程式求系統函數,以及由系統函數的極點與留數求因果系統的脈衝響應。假設系統初始靜止,使用單邊或雙邊 zz 轉換時,延遲一個取樣的訊號對應為乘上 z−1z^{-1}:

Z{y[n−k]}=z−kY(z)\mathcal{Z}\{y[n-k]\}=z^{-k}Y(z)

系統函數定義為零初始條件下的輸出與輸入之比:

H(z)=Y(z)X(z)H(z)=\frac{Y(z)}{X(z)}

解題方法

(一)求系統函數

對差分方程式

y[n]=y[n−1]−y[n−2]+2x[n]+2x[n−1]y[n]=y[n-1]-y[n-2]+2x[n]+2x[n-1]

兩邊取 zz 轉換:

Y(z)=z−1Y(z)−z−2Y(z)+2X(z)+2z−1X(z)Y(z)=z^{-1}Y(z)-z^{-2}Y(z)+2X(z)+2z^{-1}X(z)

將 Y(z)Y(z) 項移到左側、X(z)X(z) 項留在右側:

Y(z)(1−z−1+z−2)=2X(z)(1+z−1)Y(z)(1-z^{-1}+z^{-2})=2X(z)(1+z^{-1})

因此

H(z)=Y(z)X(z)=2(1+z−1)1−z−1+z−2H(z)=\frac{Y(z)}{X(z)} =\frac{2(1+z^{-1})}{1-z^{-1}+z^{-2}}

題目要求以 zz 的正冪次方多項式表示,分子、分母同乘 z2z^2:

H(z)=2z2+2zz2−z+1H(z)=\frac{2z^2+2z}{z^2-z+1}

(二)求脈衝響應

先將系統函數改寫成便於部分分式展開的形式:

H(z)=2z(z+1)z2−z+1H(z)=\frac{2z(z+1)}{z^2-z+1}

分母因式分解為

z2−z+1=(z−ejπ/3)(z−e−jπ/3)z^2-z+1 =\left(z-e^{j\pi/3}\right)\left(z-e^{-j\pi/3}\right)

令 p=ejπ/3p=e^{j\pi/3},則另一根為 p∗=e−jπ/3p^*=e^{-j\pi/3}。將 H(z)H(z) 寫成部分分式:

H(z)=Azz−p+A∗zz−p∗H(z)=A\frac{z}{z-p}+A^*\frac{z}{z-p^*}

由留數計算

🔒

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

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

免費註冊

第 四 題25 分

對於任何實數序列 x[n]x[n],其傅立葉轉換為

X(ejω)=∑nx[n]e−jωnX(e^{j\omega})=\sum_n x[n]e^{-j\omega n}

(一)請證明其傅立葉轉換之大小(Magnitude)∣X(ejω)∣\left|X(e^{j\omega})\right| 為偶函數,而相位(Phase)∠X(ejω)\angle X(e^{j\omega}) 為奇函數,即:(10 分)

∣X(ejω)∣=∣X(e−jω)∣\left|X(e^{j\omega})\right|=\left|X(e^{-j\omega})\right|

∠X(ejω)=−∠X(e−jω)\angle X(e^{j\omega})=-\angle X(e^{-j\omega})

(二)當實數序列 x[n]=0.5nu[n]x[n]=0.5^n u[n],請驗證(一)之等式成立。(15 分)

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

這一題的完整詳解

核心觀念

對實數序列,x[n]∈Rx[n]\in\mathbb{R},其離散時間傅立葉轉換(DTFT)為

X(ejω)=∑n=−∞∞x[n]e−jωn.X(e^{j\omega})=\sum_{n=-\infty}^{\infty}x[n]e^{-j\omega n}.

實數序列的 DTFT 具有共軛對稱性:

X(e−jω)=X∗(ejω).X(e^{-j\omega})=X^*(e^{j\omega}).

因此,正負頻率的頻譜大小相同,相位則互為相反數;相位敘述須以相位採用一致分支,且 X(ejω)≠0X(e^{j\omega})\neq 0 為前提。

(一)證明大小為偶函數、相位為奇函數

由定義,

X(e−jω)=∑nx[n]ejωn.X(e^{-j\omega}) =\sum_n x[n]e^{j\omega n}.

因為 x[n]x[n] 為實數,x∗[n]=x[n]x^*[n]=x[n],故

X∗(ejω)=(∑nx[n]e−jωn)∗=∑nx[n]ejωn=X(e−jω).X^*(e^{j\omega}) =\left(\sum_n x[n]e^{-j\omega n}\right)^* =\sum_n x[n]e^{j\omega n} =X(e^{-j\omega}).

所以

X(e−jω)=X∗(ejω).X(e^{-j\omega})=X^*(e^{j\omega}).

取絕對值:

∣X(e−jω)∣=∣X∗(ejω)∣=∣X(ejω)∣,|X(e^{-j\omega})| =|X^*(e^{j\omega})| =|X(e^{j\omega})|,

因此 ∣X(ejω)∣|X(e^{j\omega})| 對 ω\omega 為偶函數。

令 X(ejω)=∣X(ejω)∣ejϕ(ω)X(e^{j\omega})=|X(e^{j\omega})|e^{j\phi(\omega)}。共軛會反轉相位符號,因此

X(e−jω)=∣X(ejω)∣e−jϕ(ω),X(e^{-j\omega}) =|X(e^{j\omega})|e^{-j\phi(\omega)},

從而

∠X(e−jω)=−∠X(ejω).\angle X(e^{-j\omega})=-\angle X(e^{j\omega}).

這就是相位的奇對稱性。相位是模 2π2\pi 定義;若採主值相位,等式在相位不連續處應理解為模 2π2\pi,並需避開頻譜為零、相位未定義的點。

(二)驗證 x[n]=0.5nu[n]x[n]=0.5^n u[n]

由於 u[n]=1u[n]=1(n≥0n\ge 0)、u[n]=0u[n]=0(n<0n<0),所以

X(ejω)=∑n=0∞(0.5)ne−jωn=∑n=0∞(0.5e−jω)n.X(e^{j\omega}) =\sum_{n=0}^{\infty}(0.5)^n e^{-j\omega n} =\sum_{n=0}^{\infty}(0.5e^{-j\omega})^n.

這是等比級數,公比大小為 0.5<10.5<1,故

X(ejω)=11−0.5e−jω=11−0.5cos⁡ω+j0.5sin⁡ω.X(e^{j\omega}) =\frac{1}{1-0.5e^{-j\omega}} =\frac{1}{1-0.5\cos\omega+j0.5\sin\omega}.

令分母為 a+jba+jb,其中 a=1−0.5cos⁡ωa=1-0.5\cos\omega、b=0.5sin⁡ωb=0.5\sin\omega。因 a≥0.5>0a\ge 0.5>0,分母不為零。故頻譜大小為

🔒

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

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

免費註冊

其他考古題