lms_s.m

比较LMS

%LMS_S.m
%LMS algorithm
%static status
%synchronous CDMA
%channel: White Gaussis Noise

function [P1,P_i_n,y,correct,E_ex] = LMS_S(b,N,K,step,yN_step,S,SNR)

%multiuser detection
x=zeros(1,N);
%x(ceil(N/2))=1;
mu=1e-5;
lamda=1.01;
for i=1:step
    if mu>1e-8
        mu=mu/lamda;
    end  
    Z=(S(1,:)+x)*yN_step(:,i);
    Z_MF=S(1,:)*yN_step(:,i);
    x=x-mu*Z*(yN_step(:,i).'-Z_MF*S(1,:));
    xx(i,:)=x;
    c=S(1,:)+x;
    P1(i)=(c*S(1,:).')^2;
    P_i_n(i)=( c*(yN_step(:,i)-b(1,i)*S(1,:).') )^2;
    %SINR(i)=P1/P_i_n;
    y(i)=sign(c*yN_step(:,i));
    E=(c*yN_step(:,i)/N)^2/10;
    E_min=1;
    E_ex(i)=E-E_min;
end
correct=(y==b(1,:));

%plot(SINR);