數(shù)字信號(hào)處理課程設(shè)計(jì)報(bào)告(編輯修改稿)
2025-03-14 09:01 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 c os( ) c os( )x n n n???? 作如下譜分析： (1)截取 ()xn 使 ()xn 成為有限長(zhǎng)序列 N(0??n N1 )， (長(zhǎng)度 N 自己選 )寫(xiě)程序計(jì)算出()xn 的 N 點(diǎn) DFT ()Xk ,并畫(huà)出相應(yīng)的幅頻圖 ( ) ~X k k 。 （ 假設(shè) N 取 15，即 0≤ n≤ 14 時(shí) , 編寫(xiě)程序，計(jì)算出 X(n)的 15 點(diǎn) DFT Xk） n = 0:14。 xn=cos(*pi*n)+cos(*pi*n)。 Xk = fft (xn, 15)。 subplot(2,1,1)。 stem(n, xn)。 grid。 subplot(2,1,2)。 stem(n, abs(Xk))。 grid。 (2)將 (1)中 ()xn 補(bǔ)零 加長(zhǎng)至 M 點(diǎn)（長(zhǎng)度 M 自己選），編寫(xiě)程序計(jì)算 ()xn 的 M 點(diǎn) DFT 1()Xk ,并畫(huà)出相應(yīng)的圖 1( ) ~X k k 。 (2)將 (1)中 ()xn 補(bǔ)零加長(zhǎng)至 M 點(diǎn)（長(zhǎng)度 M 自己選），編寫(xiě)程序計(jì)算 ()xn 的 M 點(diǎn) DFT 6 1()Xk ,并畫(huà)出相應(yīng)的圖 1( ) ~X k k 。 取 M 為 32 時(shí)，因?yàn)榫帉?xiě)的 dft 函數(shù)在長(zhǎng)度 N 大于序列長(zhǎng)度時(shí)自動(dòng)補(bǔ)零，故直接調(diào)用函數(shù)即可。 程序如下： for n=0:9 xn(n+1)=cos(*pi*n)+cos(*pi*n) end dft32=dft(xn,32)。 stem(0:31,dft32)。 運(yùn)行結(jié)果如下 (3)利用補(bǔ)零 DFT 計(jì)算 (1)中 N 點(diǎn)有限長(zhǎng)序列 ()xn 頻譜 ()jXe? 并畫(huà)出相應(yīng)的幅頻圖( ) ~jXe? ? 。 當(dāng) dft 的 N 取值足夠大時(shí)，可以用 plot 函數(shù)來(lái)畫(huà)平滑的圖像，就是所要的幅頻圖 程序如下： for n=0:9 xn(n+1)=cos(*pi*n)+cos(*pi*n) end dft1024=dft(xn,1024)。 plot(0:1023,dft1024)。 7 3. 研究高密度譜與高分辨率頻譜。 對(duì)連續(xù)確定信號(hào) 3( ) c os( 2 6. 5 10 ) * tax t t e ?? ?? ? ? 令 ? =103 ， 104 ， 105 (1)對(duì)不同的 ? 考慮恰當(dāng)?shù)牟蓸宇l率進(jìn)行取樣分析，說(shuō)明理由 (2) 對(duì)某一個(gè) ? 采樣數(shù)據(jù)選取適當(dāng)長(zhǎng)度，計(jì)算 DFT，畫(huà)出相應(yīng)的 X(k)和 ()jXe? ~? (1)對(duì)不同的 ? 考慮恰當(dāng)?shù)牟蓸宇l率進(jìn)行取樣分析，說(shuō)明理由 當(dāng) ? =103 時(shí)，抽樣頻率可以選擇為 32khz，因?yàn)闈M(mǎn)足 32Khz 大于兩倍最高頻率 for n=0:15 xn(n+1)=cos(2*pi**1000*n/32021)*exp(1*1000*n/32021) end subplot(2,1,1) stem(0:15,xn) dft16=dft(xn,16)。 subplot(2,1,2)。 stem(0:15,dft16)。 運(yùn)行結(jié)果為 8
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號(hào)-1