lms_bmud.m

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

function lms_bmud()
%
%2005.12.1
%加载数据
load goldseq;
load b;
%gold码归一化
goldseq = 1/sqrt(31)*goldseq;
s1=goldseq(1,:)';
P=eye(31)-s1*s1';
u=10^-4;


r1=goldseq(1,:)'*b(1,:);
r1=awgn(r1,10);
r2=goldseq(2,:)'*b(2,:);
%r2=awgn(r2,SNR);
r3=goldseq(3,:)'*b(3,:);
%r3=awgn(r3,SNR);
r4=goldseq(4,:)'*b(4,:);
%r4=awgn(r4,SNR);
r5=goldseq(5,:)'*b(5,:);
%r5=awgn(r5,SNR);
r6=goldseq(6,:)'*b(6,:);
%r6=awgn(r6,SNR);

r=r1+r2+r3+r4+r5+r6;
 x=zeros(31,1);
 z(1)=(s1+P*x)'*r(:,1);
for i=2:1000-1        
    z(i)=(s1+P*x)'*r(:,i);
    b1(i+1)=sign(z(i)*z(i-1));
    for j=1:100
    x=x-u*z(i)*(r(:,i)-(s1'*r(:,i))*s1);     
    m1(:,j)=s1+x;    
    end
    SINR(i+1)=mean((m1'*s1).^2)/mean((m1'*(r(i)-b(1,i)*s1)).^2);
    %SINR(i+1)=((m1'*s1)^2)/()
end
i=1:1000;
semilogy(i,SINR);
xlabel('number of iteration');
ylabel('SIR');