NTMYW.rar_SER_空间谱

clear all clc %----------------SER算法与LMS算法的比较--------------------------------% %----------------SER算法-q0=10-----------------------------------------% N=16; m=0.01; M=100; k=0:M; x=sin(2*pi*k/N)+sqrt(m)*randn(1,M+1); d=2*cos(2*pi*k/N); q0=10; alpha=0.93; I=eye(2); Q=q0*I; W=[6,-9]'; u=0.05; p=0.51; for k=2:M S=Q*x(k:-1:k-1)'; r=alpha+x(k:-1:k-1)*S; Q=(1/alpha)*(Q-(1/r)*S*S'); y=x(k:-1:k-1)*W; e=d(k)-y; W=W+2*u*p*(1-alpha^(k+1))*Q*e*x(k:-1:k-1)'/(1-alpha); t0(k-1)=W(1); t1(k-1)=W(2); end k=1:M; plot(t0,t1,'g'); hold on %----------------------SER算法-q0=1---------------------------------% M=100; k=0:M; x=sin(2*pi*k/N)+sqrt(m)*randn(1,M+1); d=2*cos(2*pi*k/N); q0=10; alpha=0.93; I=eye(2); Q=q0*I; W=[6,-9]'; u=0.05; p=0.51; for k=2:M S=Q*x(k:-1:k-1)'; r=alpha+x(k:-1:k-1)*S; Q=(1/alpha)*(Q-(1/r)*S*S'); y=x(k:-1:k-1)*W; e=d(k)-y; W=W+2*u*p*(1-alpha^(k+1))*Q*e*x(k:-1:k-1)'/(1-alpha); t0(k-1)=W(1); t1(k-1)=W(2); end k=1:M; plot(t0,t1,'r'); hold on % %------------------LMS 算法----------------------------------------% M=100; u=0.05; W=[6,-9]'; r0(1)=W(1);r1(1)=W(2); k=0:M; x=sin(2*pi*k/N)+sqrt(m)*randn(1,M+1); d=2*cos(2*pi*k/N); for k=2:1:M+1 %迭代次数循环 y=x(k:-1:k-1)*W; e=d(k)-y; W=W+2*u*e*x(k:-1:k-1)'; r0(k)=W(1); r1(k)=W(2); end plot(r0,r1); hold on legend('SER q0=10','SER q0=1','LMS'); title('SER算法与LMS算法的比较'); xlabel('w0');ylabel('w1'); %-----------------等高线----------------------------------------------% N=16 ; %定义x(k)的抽样频率 m=0.01; %定义随机信号的功率 [w0,w1]=meshgrid(-2:.1:8,0:-.1:-10) ; %定义权向量和坐标轴 q=(0.5+m)*(w0.^2+w1.^2)+w0.*w1.*cos(2*pi/N)+2*w1.*sin(2*pi/N)+2; contour(w0,w1,q,[37.0,23.8,13.6,6.3,1.9,0.78,0.49,0.40]); xlabel('w0');ylabel('w1'); hold on