ADSP仿真报告-维纳滤波器.doc

Forpersonaluseonlyinstudyandresearch;mercialuseADSP仿真报告---维纳滤波器姓名:×××院系:×××学号:×××指导老师:杨灵日期:(1)信号模型AR(1)信号模型:其中,是均值为0,。(2)测量模型测量模型:其中,是均值为0,方差为1的高斯白噪声信号。图(1)(1)因果IIR维纳滤波器已知:,可知解得故待求IIR维纳滤波器的传输函数为:(2)FIR维纳滤波器设FIR维纳滤波器冲激响应序列的长度为M,冲激矢量为滤波器数据矢量滤波器输出令,;closeall;clc;%******************建立信号模型s(n)=a*s(n-1)+w(n)***************N=2000;%数据长度w=sqrt()*randn(1,N);%噪声w(n)均值为0,(1)=w(1);fori=2:Ns(i)=*s(i-1)+w(i);end%******************建立测量模型x(n)=s(n)+w(n)*******************v=randn(1,N);%噪声v(n)均值为0,方差为1x=s+v;%接收信号序列%******************IIR因果维纳滤波器**********************%计算得到滤波器参数G=,f==[1,-];y_IIR=filter(,A,x);%IIR因果维纳滤波器输出MSE_IIR=mean((s-y_IIR).^2)%IIR因果维纳滤波器均方误差%******************FIR维纳滤波器*************************%设FIR滤波器阶数为M,即冲激响应序列长度为M+1%M=30;rx=xcorr(x);rsx=xcorr(s,x);forM=1:500%求不同阶数下的滤波输出Rxx=toeplitz(rx(N:N+M))/N;%求x(n)自相关,(M+1)*(M+1)矩阵P=(rsx(N:N+M)/N)';%求s(n)与x(n)互相关,(M+1)*1向量h_opt=inv(Rxx)*P;y_FIR=filter(h_opt,1,x);%FIR滤波器输出MSE_FIR(M)=mean((s-y_FIR).^2);%FIR滤波器均方误差end%********************滤波结果************************%IIR因果维纳滤波器输出结果figuren=[1:100];plot(s(1:100),'b-',v(1:100),'g:',y_IIR(1:20:2000),'r--');legend('原始数据s(n)','噪声数据v(n)','IIR滤波输出');ylabel('振幅');xlabel('时间(n)');title('IIR因果维纳滤波器输出结果y_IIR(1:100)');%FIR维纳滤波器输出结果figuren=[1:100];plot(n,s(1:100),'b-',n,v(1:100),'g:',n,