哪位帮我看看RLS自适应滤波的程序

1990mixiaobai

1. clear
   
2. N=1000;Fs=500;
   
3. n=0:N-1;t=n/Fs;
   
4. xs=( sin(2*pi*5*t))';
   
5. figure(1);
   
6. subplot(211);
   
7. plot(t,xs);grid;
   
8. ylabel('fuzhi');
   
9. title('\it{输入周期性信号}');
   
10. xn=( 0.6*randn(1,length(t)))';
    
11. subplot(212);
    
12. plot(t,xn);grid;
    
13. ylabel('幅度');
    
14. xlabel('时间');
    
15. title('\it{随机噪声信号}');
    
16. d=xs;
    
17. x=xs+xn;
    
18. N=10;
    
19. w=(zeros(1,N))';
    
20. M=length(x);
    
21. p=0.1*eye(N,N);
    
22. a=0.1;
    
23. y=(zeros(1,M))';
    
24. e=(zeros(1,M))';
    
25. sum1=zeros(N,N);
    
26. sum2=zeros(N,1);
    
27. for n=N:M
    
28. x1=x(n:-1:n-N+1);
    
29. juzhen=x1*x1';
    
30. 
    
31. k=((1/a)*p*x1)/(1+(1/a)*x1'*p*x1);
    
32. e(n)=d(n)-w'*x1;
    
33. w=w+k*conj(e(n));
    
34. p=(1/a)*p-(1/a)*k*x1'*p;
    
35. y(n)=w'*x1;
    
36. end
    
37. figure(2);
    
38. plot(y);grid;
    
39. title('\it{自适应滤波器输出信号}')
    
40. figure(3);
    
41. plot(e);grid;
    
42. title('\it{误差}')
    

复制代码