版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、現(xiàn)代信號(hào)處理大作業(yè)姓名:潘曉丹學(xué)號(hào):0140349045班級(jí):A1403492作業(yè)1LD算法實(shí)現(xiàn)AR過(guò)程估計(jì)1.1 AR模型 p階AR模型的差分方程為:,其中是均值為0的白噪聲。AR過(guò)程的線性預(yù)測(cè)方法為:先求得觀測(cè)數(shù)據(jù)的自相關(guān)函數(shù),然后利用Yule-Walker方程遞推求得模型參數(shù),再根據(jù)公式求得功率譜的估計(jì)。Yule-Walker方程可寫(xiě)成矩陣形式:1.2 LD算法介紹Levinson-Durbin算法可求解上述問(wèn)題,其一般步驟為:1) 計(jì)算觀測(cè)值各自相關(guān)系數(shù);i=1;2) 利用以下遞推公式運(yùn)算:3) i=i+1,若ip,則算法結(jié)束;否則,返回(2)。1.3 matlab編程實(shí)現(xiàn)以AR模型:
2、xn=12xn-1-12xn-2+w(n)為例,Matlab 程序代碼如下:clear; clc;var = 1;noise = var*randn(1,10000);p = 2;coefficient = 1 -0.5 0.5;x = filter(1,coefficient,noise);divide = linspace(-pi,pi,200);for ii = 1:200 w = divide(ii); S1(ii) = var/(abs(1+coefficient(2:3)*exp(-j*w*(1:2)2;enda_p var_p=Levinson_Durbin(x,p);for i
3、i = 1:200 w = divide(ii); Sxx(ii) = var_p/(abs(1+a_p(2:p+1)*exp(-j*w*(1:p)2;endfigure;subplot(2,2,1); plot(divide,S1,b);grid onxlabel(w); ylabel(功率); title(AR 功率譜);subplot(2,2,2); plot(divide,Sxx,r-);grid onxlabel(w); ylabel(功率); title(L-D算法估計(jì));subplot(2,2,3); plot(divide,S1,b);hold onplot(divide,Sx
4、x,r-);hold offgrid onxlabel(w); ylabel(功率); title(AR功率譜和算法比較);子函數(shù):Levinson_Durbin.mfunction a_p var_p = Levinson_Durbin(x,p)N = length(x);for ii=1:N Rxx(ii) = x(1:N-ii+1)*(x(ii:N)/N; enda(1)=1; a(2)=-Rxx(2)/Rxx(1);for k=1:p-1 % Levinson-Durbin algorithm var(k+1) = Rxx(0+1)+a(1+1:k+1)*Rxx(1+1:k+1); r
5、eflect_coefficient(k+1+1) = -a(0+1:k+1)*(fliplr(Rxx(2:k+1+1)/var(k+1); var(k+1+1) = (1-(reflect_coefficient(k+1+1)2)*var(k+1); a_temp(1) = 1; for kk=1:k a_temp(kk+1) = a(kk+1)+reflect_coefficient(k+1+1)*a(k+1-kk+1); end a_temp(k+1+1) = reflect_coefficient(k+1+1); a = a_temp;enda_p = a; % prediction
6、coeffecientsvar_p = var(p+1); % prediction error power1.4 仿真結(jié)果1)p=2時(shí),仿真結(jié)果圖如下預(yù)測(cè)系數(shù):a20,a21,a22=1,-0.5068,0.5031誤差功率:var_p=1.01942)p=20時(shí),仿真結(jié)果圖如下預(yù)測(cè)系數(shù):a20,a21,a22,a23,a24,=1,-0.5098,0.4999,-0.0066,0.0060,-0.0179,0.0193,誤差功率:var_p=0.99983)p=50時(shí),仿真結(jié)果圖如下預(yù)測(cè)系數(shù):a20,a21,a22,a23,a24,=1,-0.4951,0.5178,-0.0145,0.0
7、117,-0.0169,0.0141,誤差功率:var_p=0.99551.5 結(jié)果分析由不同階數(shù)(P值)得到的仿真結(jié)果可得:當(dāng)P的階數(shù)較低時(shí),L-D算法估計(jì)AR模型對(duì)功率譜估計(jì)的分辨率較低,有平滑的效果,從P=2的仿真結(jié)果可以看出估計(jì)得到的功率譜與原始功率譜基本吻合,且曲線平滑沒(méi)有毛刺;隨著階數(shù)增大,采用L-D算法進(jìn)行估計(jì)后,得到的功率譜會(huì)產(chǎn)生振蕩,從仿真可以看到,當(dāng)階數(shù)P較高為50時(shí),估計(jì)得到的功率譜與原始功率譜基本吻合,但估計(jì)得到的功率譜曲線不平滑,有急劇的振蕩。從LD算法得到的預(yù)測(cè)系數(shù)可得,不論階數(shù)P(2)取何值,通過(guò)該算法得到的預(yù)測(cè)參數(shù)與原始目標(biāo)函數(shù)中的一致,其余各個(gè)參數(shù)均接近0。因
8、此,L-D算法得到可計(jì)算得到較為精確的預(yù)測(cè)值或估計(jì)值。作業(yè)二非平穩(wěn)信號(hào)由兩個(gè)高斯信號(hào)疊加而成,為其中,分別求出的WV分布及其模糊函數(shù),畫(huà)出二者的波形圖,指出并分析其信號(hào)項(xiàng)和交叉項(xiàng)。2.1 WV分布由WV分布的定義知:若記,則 ,其中,由此,我們計(jì)算得到:信號(hào)項(xiàng):交叉項(xiàng):其中,2.2 模糊函數(shù)由模糊函數(shù)的定義知:若記,則 ,其中,由此,我們計(jì)算得到:信號(hào)項(xiàng):交叉項(xiàng):其中,2.3 Matlab 編程實(shí)現(xiàn)取,進(jìn)行matlab編程如下clear all;format long;alpha1=20; t1=10;t2=4;w1=8;w2=4;a=alpha1/2; td=t1-t2; omegad=w1
9、-w2;tm=0.5*(t1+t2); omegam=0.5*(w1+w2); m=1;n=1;for t=0:0.1:8for omega=-6:0.1:12W_auto(m,n)=2*(exp(-a*(t-t1)2-1/a*(omega-w1)2)+exp(-a*(t-t2)2-1/a*(omega-w2)2); W_cross(m,n)=4*exp(-a*(t-tm)2-1/a*(omega-omegam)2)*cos(omega-omegam)*td+omegad*t)W(m,n)=W_auto(m,n)+W_cross(m,n); n=n+1;end m=m+1; n=1;endfi
10、gure; mesh(-6:0.1:12,0:0.1:8,W);xlabel(time);ylabel(frequency);title(WV分布);figure; mesh(-6:0.1:12,0:0.1:8,W_auto);xlabel(time);ylabel(frequency);title(WV分布信號(hào)項(xiàng));figure;mesh(-6:0.1:12,0:0.1:8,W_cross);xlabel(time);ylabel(frequency);title(WV分布交叉項(xiàng));format long; %模糊函數(shù)a=10; t1=6;t2=2;w1=6;w2=2;td=t1-t2; w
11、d=w1-w2;tm=0.5*(t1+t2); wm=0.5*(w1+w2);m=1;n=1;for t=-10:0.1:10for w=-10:0.1:10 A_auto(m,n)=abs(exp(-a/4*t2-1/(4*a)*w2)*(exp(i*w1*t-i*t1*w)+exp(i*w2*t-i*t2*w); A_cross(m,n)=abs(exp(i*wm*t+i*w*tm+i*wd*tm)*(exp(-1/(4*a)*(w+wd)2-a/4*(t-td)2)+exp(-1/(4*a)*(w-wd)2-a/4*(t+td)2);A(m,n)=A_auto(m,n)+A_cross(
12、m,n); n=n+1;end m=m+1; n=1;endfigure; mesh(-10:0.1:10,-10:0.1:10,A);xlabel(time);ylabel(frequency);title(模糊函數(shù));figure; mesh(-10:0.1:10,-10:0.1:10,A_auto);xlabel(time);ylabel(frequency);title(模糊函數(shù)信號(hào)項(xiàng));figure;mesh(-10:0.1:10,-10:0.1:10,A_cross);xlabel(time);ylabel(frequency);title(模糊函數(shù)交叉項(xiàng));2.4 仿真結(jié)果及分析
13、1)Z(t)函數(shù)的WV分布波形圖如下2)WV分布信號(hào)項(xiàng)波形圖如下3)WV分布交叉項(xiàng)波形圖如下根據(jù)結(jié)果可以看到,WV分布的兩個(gè)信號(hào)項(xiàng)是分開(kāi)的,分別以為中心;WV分布的交叉項(xiàng)則耦合在一起,以為中心。4)模糊函數(shù)波形圖5) 模糊函數(shù)信號(hào)項(xiàng)6) 模糊函數(shù)交叉項(xiàng)由結(jié)果可以看出,模糊函數(shù)的兩個(gè)信號(hào)項(xiàng)耦合在一起,以原點(diǎn)(0,0)為中心;而交叉項(xiàng)是分開(kāi)的,分別以為中心。作業(yè)三信號(hào)zt=0.25e-t2/2ejt2/2ejt是由xt=0.25e-t2/2與yt=ejt2/2ejt相乘得到。分別求出三者的WV分布,并畫(huà)出三維分布圖。該計(jì)算過(guò)程和作業(yè)二類(lèi)似,在此不再贅述。3.1 Matlab編程實(shí)現(xiàn)syms x t
14、a=10; b=2; o=2;zt=(a/pi)0.25*exp(-a*t2/2+1i*b*t2/2+1i*o*t);xt=(a/pi)0.25*exp(-a*t2/2);yt=exp(1i*b*t2/2+1i*o*t);x1=(a/pi)0.25*exp(-a*(t+x/2)2/2);x2=(a/pi)0.25*exp(-a*(t-x/2)2/2);x12=(a/pi)0.5*exp(-a*(t-x/2)2/2-a*(t-x/2)2/2);y1=exp(1i*b*(t+x/2)2/2+1i*o*(t+x/2);y2=exp(-1i*b*(t-x/2)2/2-1i*o*(t-x/2);z1=(
15、a/pi)0.25*exp(-a*(t+x/2)2/2+1i*b*(t+x/2)2/2+1i*o*(t+x/2);z2=(a/pi)0.25*exp(-a*(t-x/2)2/2+-1i*b*(t-x/2)2/2-1i*o*(t-x/2);z12=(a/pi)0.5*exp(-a*(t+x/2)2/2+1i*b*(t+x/2)2/2+1i*o*(t+x/2)-a*(t-x/2)2/2+-1i*b*(t-x/2)2/2-1i*o*(t-x/2);Wx=simple(fourier(x12)*(a/pi)(1/2);Wz=simple(fourier(z12)*(a/pi)(1/2);Wx=inline(Wx)Wz=inline(Wz) t=-6:0.01:6; f=-6:0.01:6;T,F=meshgrid(t,f);Wxwv= abs(Wx(T,2*pi*F); Wzwv= abs(Wz(T,2*pi*F);figure; mesh(T
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 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ì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 安徽省淮南市第二中學(xué)2022年高一物理第二學(xué)期期末調(diào)研模擬試題含解析
- 四年級(jí)語(yǔ)文上冊(cè)課件
- 2022年浙江省公立寄宿學(xué)校物理高一第二學(xué)期期末質(zhì)量檢測(cè)試題含解析
- 2022年云南民族大學(xué)附屬中學(xué)物理高一下期末達(dá)標(biāo)測(cè)試試題含解析
- 2022年物理高一第二學(xué)期期末經(jīng)典試題含解析
- 中醫(yī)護(hù)理課件教學(xué)
- 2024年傳熱液項(xiàng)目申請(qǐng)報(bào)告模稿
- 2024年高效包裝機(jī)項(xiàng)目申請(qǐng)報(bào)告模稿
- 2024年耐高溫超輕硅酸鈣隔熱保濕材料項(xiàng)目立項(xiàng)申請(qǐng)報(bào)告
- 2024年飼料級(jí)磷酸二氫鈣項(xiàng)目申請(qǐng)報(bào)告模稿
- JT-T-1088-2016公路工程噴射混凝土用無(wú)堿速凝劑
- 房地產(chǎn)開(kāi)發(fā)商買(mǎi)房合同模板
- YY-T1676-2020超聲內(nèi)窺鏡行業(yè)標(biāo)準(zhǔn)
- 康復(fù)醫(yī)學(xué)導(dǎo)論案例分析
- 中國(guó)的地理環(huán)境與國(guó)家安全
- 安全生產(chǎn)責(zé)任制(各工種)
- 第十二章 全等三角形 作業(yè)設(shè)計(jì)-2023-2024學(xué)年人教版八年級(jí)數(shù)學(xué)上冊(cè)
- 2024-2030年國(guó)內(nèi)油井水泥行業(yè)市場(chǎng)發(fā)展分析及發(fā)展前景與投資機(jī)會(huì)研究報(bào)告
- 工業(yè)級(jí)海綿鋯市場(chǎng)現(xiàn)狀研究分析與發(fā)展前景預(yù)測(cè)報(bào)告
- (正式版)JBT 7122-2024 交流真空接觸器 基本要求
- 發(fā)酵工程為人類(lèi)提供多樣的生物產(chǎn)品(教學(xué)課件)高二生物(浙科版2019選擇性必修3)
評(píng)論
0/150
提交評(píng)論