Burg算法估計(jì)模型參數(shù)_第1頁(yè)
Burg算法估計(jì)模型參數(shù)_第2頁(yè)
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、已知AR(4)過(guò)程:x(n)=2.760x(n-1)-3.809x(n-2)+2.654x(n-3)-0.924x(n-4)+w(n)試用Burg算法估計(jì)模型參數(shù)、階數(shù)和功率譜。實(shí)驗(yàn)仿真50次,給出估計(jì)結(jié)果的方差和偏差(用估計(jì)結(jié)果均值與真值的差來(lái)估計(jì)偏差)。用Burg算法對(duì)該AR過(guò)程估計(jì),50次實(shí)驗(yàn)得平均階數(shù)P=4.8,取整數(shù)即為P=5,階數(shù)方差var(p)=12.53,階數(shù)偏差為0.8。AR參數(shù)(平均):-2.06882.2674-1.16920.3671-0.00640.0082-0.01020.00290.0067-0.00550.00090.0007aa0.83604210AR系數(shù)估計(jì)

2、方差6ARinJ50.5521012ar索數(shù)估計(jì)如差6ARin:-D.5Burg最大熵譜估計(jì):3uiy大誼估計(jì)珀率諾0.611.622.633.5頻率matlab代碼:u=1;whileu<=50%50次仿真實(shí)驗(yàn)n=10000;w二randn(l,n);%生成零均值,方差為0.2的白噪聲w=w-mean(w);w=w/std(w);w=sqrt(0.2)*w;x(l)=l+w(l);x(2)=2+w(2);x(3)=3+w(3);x(4)=4+w(4);%由ar模型生成含噪信號(hào);fori=5:nx(i)=2.760*x(i-1)-3.809*x(i-2)+2.654*x(i-3)-0.9

3、24*x(i-4)+w(i);endP(l+l)=l/n*sum(x."2);%計(jì)算預(yù)測(cè)誤差功率的初始值g(1,:)=x;f(1,:)=g(1,:);m=1;k=0;l=0;a(1,1)=1;P(1)=1;whileabs(P(m+1)-P(m)/P(m)>=0.01forj=m+1:nk=k+f(m,j)*g(m,j-1);l=l+f(m,j廠2+g(m,jl廠2;endK(m)=-k/(1/2*l);%求反射系數(shù)a(1,1)=1;ifm>=2forj=1:m-1%計(jì)算前向預(yù)測(cè)濾波器系數(shù)a(m,j)=a(m-1,j)+K(m)*a(m-1,m-j);a(m,m)=K(m

4、);endendP(m+l+l)=(l+abs(K(m)廠2)*P(m+1);%計(jì)算預(yù)測(cè)誤差功率forj=l:n%計(jì)算濾波器輸出f(m+1,j)=f(m,j)+K(m)*g(m,n-1);g(m+1,j)=K(m)*f(m,j)+g(m,n-1);endm=m+1;endM(u)=m-1;%記錄每次實(shí)驗(yàn)的階數(shù)u=u+1;endfori=1:50ar二arburg(x,M(i);%估計(jì)AR系數(shù)ar_coef(i,1:length(ar)=ar;Px,fx=pburg(x,M(i);%估計(jì)功率譜Pxx(i,1:length(Px)=10*log10(Px');fxx(i,1:length(

5、fx)=fx'end%求階數(shù)和AR系數(shù)的方差和偏差var_p=var(M);fori=1:max(M)var_ar(i)=var(ar_coef(:,i+1);endbias_p=mean(M)-4;ar_r(max(M)=0;ar_r(1)=-2.76;ar_r(2)=3.809;ar_r(3)=-2.654;ar_r(4)=0.924;fori=1:max(M)mean_ar(i)=mean(ar_coef(:,i+1);bias_ar(i)=mean(ar_coef(:,i+1)-ar_r(i);end%50次實(shí)驗(yàn)功率譜平均fori=1:length(Pxx(1,:)mean_P

6、xx(i)=mean(Pxx(:,i);end%作圖plot(fxx(l,:),mean_Pxx);%作功率譜圖title('burg最大熵估計(jì)功率譜');xlabel('頻率');ylabel('功率譜密度');figure(2);plot(var_ar);%AR系數(shù)估計(jì)方差title('AR系數(shù)估計(jì)方差');xlabel('AR(n)');ylabel('方差');figure(3);plot(bias_ar);%AR系數(shù)估計(jì)偏差title('AR系數(shù)估計(jì)偏差');xlabel('AR(n)');ylabel('偏差');disp

溫馨提示

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

最新文檔

評(píng)論

0/150

提交評(píng)論