数字信号处理（ii）

理想滤波器的h(n)Р h1=zeros(N:1);Р for i=1:NР h1(i)=0.238*0.724^(i-1)*u(i);Р endР %利用公式,计算Rxx和rxsР Rxx=zeros(N,N);Р rxs=zeros(N,1);Р for i=1:NР for j=1:NР m=abs(i-j);Р tmp=0;Р for k=1:(L-m)Р tmp=tmp+x(k)*x(k+m);Р endР Rxx(i,j)=tmp/(L-m);Р endР endР Р for m=0:N-1Р tmp=0;Р for i=1: L-mР tmp=tmp+x(i)*s(m+i);Р endР rxs(m+1)=tmp/(L-m);Р endР %产生FIR维纳滤波器的h(n)Р h2=zeros(N,1);Р h2=Rxx^(-1)*rxs;Р %绘出理想和维纳滤波器h(n)的曲线图Р Р i=1:N;Р subplot(2,2,2);Р plot(i,h1(i),i,h2(i),'r:');Р title('h(n) & h~(n)');Р legend('h(n) ','h~(n)');Р %计算SiР Si=zeros(1,L);Р Si(1)=x(1);Р for i=2:LР Si(i)=0.724*Si(i-1)+0.238*x(i);Р endР %绘出Si(n)和s(n)曲线图Р Р i=L-100:L;Р subplot(2,2,3);Р plot(i,s(i),i,Si(i),'r:');Р title('Si(n) & s(n)');Р legend('Si(n) ','s(n)');Р Р %计算SrР Sr=zeros(1,L);Р for i=1:LР tmp=0;Р for j=1:N-1Р if(i-j<=0)Р tmp=tmp;