jiaquanokzishiying.m

加权波束形成形成

%***************** 波束形成 ******************
         %   圆形阵加权波束形成
clc;
clear;
f=5000;
fs=100*f;
N=8;
ns=256;  
c=1500;
lmda=c/f;
r=0.5*lmda;
b=2*pi/N;
SNR=10;
st0=pi/6; %入射角度
st=(-180:180);
%********各基元接收信号 *******
s=zeros(N,ns);
for i=1:N
    for n=1:ns
     tao(i)=r*cos(st0-(i-1)*b)/c;
    s(n)=exp(j*2*pi*f*n/fs);
    end
end
X=zeros(N,ns);
for i=1:N
 noise(i,:)=randn(1,ns);
end
for i=1:N
    for n=1:ns
    a(i)=exp(-j*2*pi*f*tao(i));
    A=(2*10^(SNR/10))^0.5;
    X(i,:)=A*a(i)*s(n)+noise(i,:);
    end
end
R=1/ns*X*X';  R1=inv(R);  %求逆
%  ***** 加权处理 ****
stm=(-180:180)*pi/180;
for m=1:length(stm)
   for i=1:N
    c(i)=exp(-j*2*pi*r*cos(st(m)*pi/180-(i-1)*b)/lmda);
    end
   w=c(:);
 % Pbf(st1)=w'*R*w/(w'*w);
  Pbf(m)=1/(w'*R1*w);
end
p=max(Pbf);
D=Pbf/p;
zz=10*log10(D);
% **** 画图 ****
figure(9);
subplot(1,2,1);
plot(st,D);
title('归一化指向性');
xlabel('jiaodu');ylabel('归一化值');
grid;
subplot(1,2,2);
plot(st,zz);
title('dB指向性');
xlabel('jiaodu');ylabel('dB')
grid;