rls_s.m

比较LMS

%RLS_S.m
%RLS algorithm
%static status
%synchronous CDMA
%channel: White Gaussis Noise

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

%multiuser detection
yyN_step=(yN_step);
delta=1e-2;
P=1/delta*eye(N);%P=1/R
lamda=0.997;
for i=1:step
    KK=P*yyN_step(:,i)/(lamda+yyN_step(:,i)'*P*yyN_step(:,i));
    P=1/lamda*(P-KK*yyN_step(:,i)'*P);
    h=P*S(1,:).';
    c=h/(S(1,:)*h);
    
    P1(i)=abs( (c.'*S(1,:).')^2 );
    P_i_n(i)=abs( ( c.'*(yyN_step(:,i)-b(1,i)*S(1,:).') )^2 );
    %SINR(i)=P1/P_i_n;
    y(i)=sign(real(c.'*yyN_step(:,i)));
    
    E=abs( (c.'*yyN_step(:,i))^2 );
    E_min=1;
    E_ex(i)=E-E_min;
end
correct=(y==b(1,:));

%plot(SINR);