現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7

現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7

ID:37141523

大?。?.14 MB

頁(yè)數(shù):31頁(yè)

時(shí)間:2019-05-18

現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7_第1頁(yè)
現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7_第2頁(yè)
現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7_第3頁(yè)
現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7_第4頁(yè)
現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7_第5頁(yè)
資源描述:

《現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)7》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)

1、專業(yè)資料現(xiàn)代數(shù)字信號(hào)處理仿真作業(yè)1.仿真題3.17仿真結(jié)果及圖形:圖1基于FFT的自相關(guān)函數(shù)計(jì)算word完美格式專業(yè)資料圖2基于式3.1.2的自相關(guān)函數(shù)的計(jì)算圖3周期圖法和BT法估計(jì)信號(hào)的功率譜word完美格式專業(yè)資料圖4利用LD迭代對(duì)16階AR模型的功率譜估計(jì)16階AR模型的系數(shù)為:a1=-0.402637623107952-0.919787323662670i;a2=-0.013530139693503+0.024214641171318i;a3=-0.074241889634714-0.088834852915013i;a4=0.0278810223539

2、97-0.040734794506749i;a5=0.042128517350786+0.068932699075038i;a6=-0.0042799971761507+0.028686095385146i;a7=-0.048427890183189-0.019713457742372i;a8=0.0028768633718672-0.047990801912420ia9=0.023971346213842+0.046436389191530i;a10=0.026025963987732+0.046882756497113i;a11=-0.03392939778

3、4767-0.0053437929619510i;a12=0.0082735406293574-0.016133618316269i;a13=0.031893903622978-0.013709547028453i;a14=0.0099274520678052+0.022233240051564i;a15=-0.0064643069578642+0.014130696335881i;a16=-0.061704614407581-0.077423818476583i.仿真程序(3_17):clearallclc%%產(chǎn)生噪聲序列N=32;%基于FFT的樣本長(zhǎng)度wor

4、d完美格式專業(yè)資料%N=256;%周期圖法,BT法,AR模型功率譜估計(jì)的長(zhǎng)度vn=(randn(1,N)+1i*randn(1,N))/sqrt(2);%%產(chǎn)生復(fù)正弦信號(hào)f=[0.150.170.26];%歸一化頻率SNR=[303027];%信噪比A=10.^(SNR./20);%幅度signal=[A(1)*exp(1i*2*pi*f(1)*(0:N-1));%復(fù)正弦信號(hào)A(2)*exp(1i*2*pi*f(2)*(0:N-1));A(3)*exp(1i*2*pi*f(3)*(0:N-1))];%%產(chǎn)生觀察樣本un=sum(signal)+vn;%%利用3.

5、1.1的FFT估計(jì)Uk=fft(un,2*N);Sk=(1/N)*abs(Uk).^2;r0=ifft(Sk);r1=[r0(N+2:2*N),r0(1:N)];%%利用3.1.2估計(jì)Rr2=xcorr(un,N-1,'biased');%畫(huà)圖k=-N+1:N-1;figure(1)subplot(1,2,1)stem(k,real(r1))xlabel('m');ylabel('實(shí)部');subplot(1,2,2)stem(k,imag(r1))xlabel('m');ylabel('虛部');figure(2)subplot(1,2,1)stem(k,re

6、al(r2))xlabel('m');ylabel('實(shí)部');subplot(1,2,2)stem(k,imag(r2))xlabel('m');ylabel('虛部');%%周期圖法NF=1024;Spr=fftshift((1/NF)*abs(fft(un,NF)).^2);kk=-0.5+(0:NF-1)*(1/(NF-1));Spr_norm=10*log10(abs(Spr)/max(abs(Spr)));%%BT法M=64;r3=xcorr(un,M,'biased');BT=fftshift(fft(r3,NF));word完美格式專業(yè)資料BT_

7、norm=10*log10(abs(BT)/max(abs(BT)));figure(3)subplot(1,2,1)plot(kk,Spr_norm)xlabel('w/2pi');ylabel('歸一化功率譜/DB');title('周期圖法')subplot(1,2,2)plot(kk,BT_norm)xlabel('w/2pi');ylabel('歸一化功率譜/DB');title('BT法')%%LD迭代算法p=16;r0=xcorr(un,p,'biased');r4=r0(p+1:2*p+1);%計(jì)算自相關(guān)函數(shù)a(1,1)=-r4(2)/r4(1)

8、;sigma(1)=r4

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫(huà)的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無(wú)此問(wèn)題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫(kù)負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無(wú)法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。