隨機過程的模擬與特征估計_第1頁
隨機過程的模擬與特征估計_第2頁
隨機過程的模擬與特征估計_第3頁
隨機過程的模擬與特征估計_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進行舉報或認領(lǐng)

文檔簡介

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)實驗2.1 隨機過程的模擬與特征估計實驗結(jié)果及分析：實驗2.1 （1）估計x(n)=0.8*x(n-1)+1+4.*randn(N,1)隨機序列的自相關(guān)函數(shù)和功率譜MATLAB仿真程序%估計x(n)=0.8*x(n-1)+1+4.*randn(N,1)隨機序列的自相關(guān)函數(shù)和功率譜%x(n)=0.8*x(n-1)+1+4.*randn(N,1)隨機序列的產(chǎn)生a=0.8;N=500;w=1+2.*randn(N,1);x(1)=w(1); for n=2:Nx(n)=a*x

2、(n-1)+w(n); endsubplot(3,2,1);plot(x);title(隨機序列x(n)=0.8*x(n-1)+1+4.*randn(N,1);grid on%估計自相關(guān)函數(shù)R=xcorr(x,coeff);subplot(3,2,2);axis(0 500 0 1);plot(R);title(自相關(guān)函數(shù));grid on%估計功率譜%周期圖功率譜估計subplot(3,2,3);periodogram(x,512,1000);axis(0 500 -50 0);title(周期圖功率譜估計)%加漢寧窗window=hann(500);subplot(3,2,4);perio

3、dogram(x,window,512,1000);axis(0 500 -50 10);title(漢寧周期功率譜估計)%相關(guān)函數(shù)法R=xcorr(x)/15000;Pw=fft(R);subplot(3,2,5);f=(0:length(Pw)-1)*1000/length(Pw);plot(f,10*log10(abs(Pw);axis(0 500 -50 10);title(BT功率譜估計)grid onsubplot(3,2,6);pwelch(x,128,64,1000);axis(0 500 -50 10);title(韋爾奇功率譜估計);grid on；實驗2.1 （2）x=s

4、in(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)隨機序列的自相關(guān)函數(shù)和功率譜N=256時的結(jié)果：N=1024時的結(jié)果：MATLAB仿真程序N=256:%估計x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)隨機序列的自相關(guān)函數(shù)和功率譜%x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)隨機序列的產(chǎn)生N=256; %N=256或1024w=randn(N,1); for n=1:N x(n)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+

5、w(n); endsubplot(3,2,1);plot(x);axis(0 260 -8 8);title(隨機序列x(N)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)/N=256);grid on%估計自相關(guān)函數(shù)R=xcorr(x,coeff);subplot(3,2,2);plot(R);axis(0 500 -1 1);title(自相關(guān)函數(shù)/N=256);grid on%估計功率譜%周期圖功率譜估計subplot(3,2,3);periodogram(x,512,1000);axis(0 500 -50 0);title(周期圖功率譜

6、估計/N=256)%加漢寧窗window=hann(256);subplot(3,2,4);periodogram(x,window,256,1000);axis(0 500 -50 10);title(漢寧周期功率譜估計)%相關(guān)函數(shù)法R=xcorr(x)/15000;Pw=fft(R);subplot(3,2,5);f=(0:length(Pw)-1)*1000/length(Pw);plot(f,10*log10(abs(Pw);axis(0 500 -50 10);title(BT功率譜估計/N=256)grid onsubplot(3,2,6);pwelch(x,128,64,1000

7、);axis(0 500 -50 10);title(韋爾奇功率譜估計/N=256);grid on;N=1024:%估計x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)隨機序列的自相關(guān)函數(shù)和功率譜%x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)隨機序列的產(chǎn)生N=1024; %N=256或1024w=randn(N,1); for n=1:N x(n)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+w(n); endsubplot(3,2,1);plot(x);axis(

8、0 1030 -8 8);title(隨機序列x(N)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)/N=1024);grid on%估計自相關(guān)函數(shù)R=xcorr(x,coeff);subplot(3,2,2);plot(R);axis(0 2000 -1 1);title(自相關(guān)函數(shù)/N=1024);grid on%估計功率譜%周期圖功率譜估計subplot(3,2,3);periodogram(x,1024,1000);axis(0 500 -50 0);title(周期圖功率譜估計/N=1024)%加漢寧窗window=hann(1024);subplot(3,2,4);periodogram(x,window,1024,1000);axis(0 500 -50 10);title(漢寧周期功率譜估計)%相關(guān)函數(shù)法R=xcorr(x)/15000;Pw=fft(R);subplot(3,2,5);f=(0:length(Pw)-1)*1000/length(Pw);plot(f,10*log10(abs(Pw

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論