关于LMS算法的一些问题，希望大神指导

您好，您能把您所用的LMS算法的代码提供一下嘛？

%LMS组合滤波器的仿真
clear all;
close all;
M=1000; %迭代次数
N=100; %运行次数
u1=0.05;  %两个滤波器的步长
u2=0.005;
MM=3;
p =zeros(M+3,1);
q=zeros(M+3,1);
pq=zeros(M+3,1);

w0=[0.8212,0.1620,-0.5987];
abc=zeros(M+3,1);
% x = randn(M+3,1);%输入信号
% d=w0'*x+n0(n);

for m=1:N
  
    w1 =0.00005*ones(MM,1); %加权矢量初始化
    w2=0.00005*ones(MM,1);

    n0=0.01 * randn(M+3,1); %噪声
    c=0.5;
    a=0;
   
    x =0.1* randn(M+3,1);%输入信号
   

    for n =MM:(M+3)
        
     d=w0*x(n:-1:n-MM+1)+n0(n-MM+1);
         
     %滤波器一
     y1 = x(n:-1:n-MM+1)' * w1;     %滤波器输出
     e1 = d-y1;%误差函数
     w1 = w1+u1*e1*x(n:-1:n-MM+1);%权矢量更新
     p(n) = p(n)+e1.^2;
     
     %滤波器二
     y2 = x(n:-1:n-MM+1)'*w2;     %滤波器输出
     e2 = d - y2;%误差
     w2 = w2+u2*e2*x(n:-1:n-MM+1); %权矢量更新
     q(n) = q(n)+e2.^2;
   
     %组合滤波器
       y=c*y1+(1-c)*y2;
       e=d-y;
       a=a+10*e*(y1-y2)*c*(1-c);
       c=1/(1+exp(-a));
       abc(n)=c;
       pq(n)=pq(n)+e.^2;

    end
   
end

p = p./N;
p = 10*log10(p);%取dB
HL = fir1(20,0.1);
p = filter(HL,1,p);
figure;
n = 1:M;
plot(n,p(n),'r');
grid on;hold on;

q=q./N;
q= 10*log10(q);
HL = fir1(20,0.1);
q = filter(HL,1,q);
n = 1:M;
plot(n,q(n),'g');
grid on;

pq = pq./N;
pq = 10*log10(pq);
HL = fir1(20,0.1);
pq = filter(HL,1,pq);
n = 1:M;
plot(n,pq(n),'b');
grid on;hold off;

legend('u1-LMS','u2-LMS','c-LMS');
xlabel('iterations');
ylabel('MSE(dB)');



都在这了，希望能得到你的指导

您好，就我的问题而言，您能给我点建议吗，多谢了

抱歉啊，我也是刚接触这一块，你的问题我尝试着试了一下，还是没有办法得到你想要的结果，最主要的是不知道您这个算法的原理

ddddddddddddddddddddddddddddddddddddd