mud.asv

此程序为DS-CDMA链路中加入基于子空间的PAST多用户检测方法。系统中采用m序列阔频

function [ich10,qch10]=mud(ich9,qch9,MUDmode,attn,A,code)
switch MUDmode
    case 1
        ich10 = ich9;
        qch10 = qch9;
    case 2
        R=(code*code')/7;
        ich10 = inv(R)*ich9;
        qch10 = inv(R)*qch9;
    case 3
        R=(code*code')/7;
        ich10 = inv(R+attn(1)^2*inv(A)^2)*ich9
        qch10 = inv(R+attn(1)^2*inv(A)^2)*qch9
    case 4
        [n,v]=size(ich9');
        X=zeros(n,v+1);   %initialazing parameters 
        y=zeros(v);
        namuda=zeros(v);
        namudapre=zeros(v);
        U=zeros(n,v);
        Upre=zeros(n,v);
        X(:,1)=ich9(1,:);
        beta = 1;
        for k = 1:v             %begin to estimation
            fid=fopen('PASTdU.mat','r');
            [Upre,count]=fread(fid,[n,v],int8);
            fclose(fid);
            fid=fopen('PASTdna.mat','r')
            [namudapre,count]=fread(fid,v,int8);
            fclose(fid);
            y(k)=Upre(:,k)'*X(:,k);
            namuda(k)=beta*namudapre(k)+abs(y(k))^2;
            U(:,k)=Upre(:,k)+(X(:,k)-Upre(:,k)*y(k))*y(k)/namuda(k);
            X(:,k+1)=X(:,k)-U(:,k)y(k);
        end
        
        
        
        
end