过载lms.m

信号过载情况下的SVM波束形成, 比其他常用的方法效果好很多.

close all;
clear all;
N=9;
M=8;
j=sqrt(-1);
len=100;
% rp.TNRdB=10;
% rp.INRdB=10;
% rp.mean_v  = 0; % mean of complex-valued AWGN
% rp.var_v   = 1; % variance of complex-valued AWGN
theta=[-80*pi/180,-60*pi/180,-40*pi/180,-20*pi/180,0*pi/180,10*pi/180,20*pi/180,45*pi/180,60*pi/180];
s0=2*(rand(1,len)>0.5)-1;
s1=2*(rand(1,len)>0.5)-1;
s2=2*(rand(1,len)>0.5)-1;
s3=2*(rand(1,len)>0.5)-1;
s4=2*(rand(1,len)>0.5)-1;
s5=2*(rand(1,len)>0.5)-1;
s6=2*(rand(1,len)>0.5)-1;
s7=2*(rand(1,len)>0.5)-1;
s8=2*(rand(1,len)>0.5)-1;
% t=0:0.01:0.2;
% x=rectpuls(t,1);
% for i=1:len
%     x_bpsk1((i-1)*20+1:i*20)=s0(1,i)*x(1:20);%基带信号
%     x_bpsk2((i-1)*20+1:i*20)=s1(1,i)*x(1:20);
% end
for i=1:len
   x_bpsk1(i)=sqrt(10.^(-5/10))*s0(i);
    x_bpsk2(i)=sqrt(10.^(6/10))*s1(i);
     x_bpsk3(i)=sqrt(10.^(5/10))*s2(i);
      x_bpsk4(i)=sqrt(10.^(-5/10))*s3(i); 
       x_bpsk5(i)=sqrt(10.^(10/10))*s4(i);
        x_bpsk6(i)=sqrt(10.^(-8/10))*s5(i);
         x_bpsk7(i)=sqrt(10.^(-5/10))*s6(i);
          x_bpsk8(i)=sqrt(10.^(2/10))*s7(i);
           x_bpsk9(i)=sqrt(10.^(5/10))*s8(i);
end
s=[x_bpsk1;x_bpsk2;x_bpsk3;x_bpsk4;x_bpsk5;x_bpsk6;x_bpsk7;x_bpsk8;x_bpsk9];
for i=1:M
    for k=1:N
        A1(i,k)=exp(-j*pi*(i-1)*sin(theta(k)));%array responce
    end
end

rp.mean_v  = 0; % mean of complex-valued AWGN
rp.var_v   = 1; % variance of complex-valued AWGN
vr = sqrt(rp.var_v/2) * randn(M,len) + rp.mean_v;
vi = sqrt(rp.var_v/2) * randn(M,len) + rp.mean_v;
v  = vr + j*vi;
X_nn=A1*s;
X=X_nn+v;
d=s4;
y=zeros(1,len);
e=zeros(len,1);
 w=zeros(M,1);
   u=0.01; 
   for i=1:len
    y(i)=w'*X(:,i);
    e(i)=d(i)-y(i);
    w=w+u*X(:,i)*conj(e(i));
   end
 

sita1=[-90:2:90];
sita=sita1*pi/180;
for i=1:length(sita)
a_sita_i=[1;exp(-j*pi*sin(sita(i)));exp(-j*pi*2*sin(sita(i)));exp(-j*pi*3*sin(sita(i)));exp(-j*pi*4*sin(sita(i)));exp(-j*pi*5*sin(sita(i)));exp(-j*pi*6*sin(sita(i)));exp(-j*pi*7*sin(sita(i)))];
G_sita(i)=a_sita_i'*w;
end;
y_out=e;
figure;
plot(sita1,10*log10(abs(G_sita)));
grid on;