pc_mud.asv

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function SINR=pc_mud()
%
%2005.12.6
load r;
load goldseq;
goldseq = 1/sqrt(31)*goldseq;
s1=goldseq(1,:)';
s2=goldseq(2,:)';
s3=goldseq(3,:)';
s4=goldseq(4,:)';
s5=goldseq(5,:)';
s6=goldseq(6,:)';

A2=10;
A3=10;
A4=10;
A5=10;
A6=10;
%[dim,N]=size(r);
N=31;

e=0.001;

x(k)=sqrt(e)*r(:,t);
for i=1;L
    eta=U(:,k-1)'*x(k);
    y=abs(eta);
    z=y/gama(k-1);
end
for i=1:L
    y(k)=(1-e)*gama(k-1)+y^2;
end
v=x(k);
w=zeros(L,1);
for i=1:L
    v=v-y*U(:,k-1);
    U(:,k)=U(:,k-1)+z*v-y*w;
    U(:,k)=U(:,k)/norm(U(:,k));
    w=w+z*U(:,k-1);
end