中国地质大学现代数字信号处理选读.doc

现代数字信号处理上机报告 任课老师 姓名 学 号 专 业 第1题 （1）代码： clear all; close all; sample=1000; %取1000个样本点 %误差判定阈值 e1=0.0008; e2=0.002; %N=300;%滤波器阶数 n=1:sample; %产生x,y方向的期望信号 dnx=cos(0.002*pi*n); dny=sin(0.002*pi*n); %产生x,y方向的噪声 xn=sqrt(0.03)*randn(1,sample); yn=sqrt(0.02)*randn(1,sample); %产生x,y方向的信号 x=dnx+xn; y=dny+yn; %产生x方向上观测信号的自相关函数 rxx=xcorr(x); %产生x方向上观测信号与期望信号的互相关函数 rxd=xcorr(x,dnx); %产生y方向上观测信号的自相关函数 ryy=xcorr(y); %产生y方向上观测信号与期望信号的互相关函数 ryd=xcorr(y,dny); %x方向滤波 for Nx=1:10:999 for i=1:Nx for j=1:Nx Rxx(i,j)=rxx(sample-i+j); %x方向信号的自相关矩阵 end end for i=1:Nx Rxd(i)=rxd(sample-1+i); %x方向信号和期望信号的互相关矩阵 end hoptx=inv(Rxx)*Rxd;%由维纳-霍夫方程得到的x方向上的滤波器最优解 ax=[1 zeros(1,Nx-1)]; fx=filter(hoptx,ax,x);%滤波后x方向上的输出 ex=dnx-fx;%误差信号 Ex=sum(abs(ex).^2)/sample; %最小均方误差 if(Ex=e1) break; end end %y方向滤波 for Ny=1:10:999 for i=1:Ny for j=1:Ny Ryy(i,j)=ryy(sample-i+j); %y方向信号的自相关矩阵 end end for i=1:Ny Ryd(i)=ryd(sample-1+i); %y方向信号和期望信号的互相关矩阵 end hopty=inv(Ryy)*Ryd;%由维纳-霍夫方程得到的y方向上的滤波器最优解 ay=[1 zeros(1,Ny-1)]; fy=filter(hopty,ay,y);%滤波后y方向上的输出 ey=dny-fy;%误差信号 Ey=sum(abs(ey).^2)/sample; %最小均方误差 if(Ey=e2) break; end end %%%%%%%%%%%%绘制图形%%%%%%%%%%%%%% subplot(2,5,1) plot(dnx); axis([0 1000 -1.5 1.5]); title(x方向期望信号); subplot(2,5,2) plot(xn); axis([0 1000 -1.5 1.5]); title(x方向噪声信号); subplot(2,5,3) plot(x); axis([0 1000 -1.5 1.5]); title(x方向观测信号); subplot(2,5,4) plot(fx