psvm.m

近似支持响亮机用来分类, 在这个程序中用来做波束形成,效果很好.

close all;
clear all;
N=5;
M=8;
j=sqrt(-1);
len=100;

theta=[-45*pi/180,-30*pi/180,0*pi/180,10*pi/180,40*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;

for i=1:len
   x_bpsk1(i)=sqrt(10.^(-5/10))*s0(i);
    x_bpsk2(i)=sqrt(10.^(8/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.^(-8/10))*s4(i);
end
s=[x_bpsk1;x_bpsk2;x_bpsk3;x_bpsk4;x_bpsk5];
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;
L=100;

%PSVM法
a=zeros(L,1);
y=s2;
w=zeros(M,1);
e=ones(L,1);
D=zeros(L,L);
for i=1:L
    D(i,i)=y(i);
end
I=eye(L,L);
V=0.1;
a=inv(I/V+D*(X'*X+e*e')*D)*e;%目录函数中输入条件的参数
w=X*D*a;%得到的权值
r=-e'*D*a;%r决定了超平面相对与原点的位置
%分类判断
f=zeros(1,L);
K=0;
for i=1:L
    f(i)=sign(X(:,i)'*w-r);
    if f(i)=y(i)
        K=K+1;%BPSK码判断正确的个数
    end
end
    

sita1=[-90:2:90];
sita=sita1*pi/180;
G_sita=zeros(1,length(sita));
a_sita_i=zeros(M,length(sita));
for i=1:length(sita)
a_sita_i(:,i)=[1;exp(-j*pi*sin(sita(i)));exp(-j*2*pi*sin(sita(i)));exp(-j*3*pi*sin(sita(i)));exp(-j*4*pi*sin(sita(i)));exp(-j*5*pi*sin(sita(i)));exp(-j*6*pi*sin(sita(i)));exp(-j*7*pi*sin(sita(i)))];
G_sita(i)=a_sita_i(:,i)'*w;
end


figure;
plot(sita1,10*log10(abs(G_sita)));
grid on ;