请高手帮俺看下LMS的程序！！！

bravecapta

谢谢你的回复，我昨天对程序进行了修改，改成100次计算的平均，结果波形不平坦，和理论计算波形差别很大，单次计算的波形又不正确，还望高手出手相助啊！！

1. a1=-0.195;a2=0.95;p=-10.22;%功率是-10.22dbw
   
2. x=zeros(1,10000);
   
3. w1=zeros(1,10000);
   
4. w2=zeros(1,10000);
   
5. a=zeros(1,100);
   
6. b=zeros(1,100);
   
7. v=wgn(1,1000,p);
   
8. x(1)=v(1);
   
9. x(2)=v(2)-a1*x(1);
   
10. e(1)=x(1);
    
11. e(2)=x(2);
    
12. w1(1)=0;w2(1)=0;
    
13. w1(2)=0;w2(2)=0;
    
14. w1(3)=0.02*x(1)*e(2);
    
15. w2(3)=0;
    
16. for n=3:1000;
    
17. for m=1:100;
    
18. e(n)=x(n)-w1(n)*x(n-1)-w2(n)*x(n-2);
    
19. w1(n+1)=w1(n)+ 0.02*x(n-1)*e(n);
    
20. w2(n+1)=w2(n)+ 0.02*x(n-2)*e(n);
    
21. x(n)=v(n)-a1*x(n-1)-a2*x(n-2);
    
22. a(m)=w1(n);
    
23. b(m)=w2(n);
    
24. end;
    
25. for i=1:100;
    
26. w1(n)=w1(n)+a(i);
    
27. w2(n)=w2(n)+b(i);
    
28. end
    
29. w1(n)=w1(n)/100;
    
30. w2(n)=w2(n)/100;
    
31. end
    
32. n=1:1000;
    
33. plot(n,w1(n),'r',n,w2(n),'b');

复制代码

simon21

你的程序一点注解都没有，读起来很来
给你个例子，你自己参考改改吧

matlab lms 算法实例

1. function [wn,en]=lms(xn,dn,M,delt,varargin)
   
2. % LMS Algorithm ,返回滤波器加权系数矩阵和误差向量
   
3. %
   
4. % 调用格式
   
5. % [wn,en]=lms(xn,dn,M,delt,itr)
   
6. % en=滤波器输出和d(n)的误差序列,为列向量
   
7. % wn=滤波器的加权参量序列,为一矩阵,其每行代表一个加权参量,每列代表一次迭代;初始化值设为0
   
8. % xn=输入列向量信号
   
9. % dn=期望列向量信号
   
10. % M=滤波器阶数
    
11. % delt=步长
    
12. % itr=迭代次数
    
13. %
    
14. % [wn,en]=lms(xn,dn,M,delt)
    
15. % 迭代次数为默认值,即等于x(n)的点数
    
16. 
    
17. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
18. % 输入输出参数的检查
    
19. %--------------------------------------------------------------------------
    
20. vin=length(varargin); Item=vin+4;
    
21. error(nargchk(4,Item,nargin)); % 检查输入变量数目是否合适,其中前四个参数必须输入
    
22. if nargout>2 % 检查输出变量数目是否合适
    
23. error('Too many output arguments');
    
24. end
    
25. %------------------------------------------------------------------------
    
26. Nx=length(xn); % x(n)的长度
    
27. if Nx~=length(dn) % 检查x(n)和d(n)长度是否相等
    
28. error('The length of x(n) is not equal to that of d(n)');
    
29. end
    
30. %------------------------------------------------------------------------
    
31. sizex=size(xn); % 检查输入向量是否为列向量
    
32. if sizex(1)<sizex(2)
    
33. xn=xn.';
    
34. end
    
35. sizedn=size(dn); % 检查期望信号向量是否为列向量
    
36. if sizedn(1)<sizedn(2)
    
37. dn=dn.';
    
38. end
    
39. %-------------------------------------------------------------------------
    
40. itr=Nx; % 迭代次数取默认值
    
41. %-------------------------------------------------------------------------
    
42. % 当输入变量为5个时
    
43. if Item==5 % 确定迭代次数
    
44. itr=varargin{1};
    
45. if itr>Nx | itr<M
    
46. error('Too many or too few iterations');
    
47. end
    
48. end
    
49. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
50. % 参数的初始化
    
51. en=zeros(itr,1);
    
52. wn=zeros(M,itr); % 每行代表一个加权参量,每列代表一次迭代
    
53. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
54. % 迭代计算
    
55. for k=M:itr % 第k次迭代
    
56. x_tap=xn(k:-1:k-M+1);
    
57. en(k)=dn(k)-wn(:,k-1)'*x_tap;
    
58. wn(:,k)=wn(:,k-1)+2*delt*en(k,1)*x_tap;
    
59. end

复制代码

yang_miya

另外你的信号x(n)产生有点问题，用二阶自回归模型产生的，使用循环产生即可。

bravecapta

谢谢大家的回复！！不胜感激！