adapbf_snr.m

在无缘雷达探测里雷达天线接收的自适应波束形成算法

%adaptive beamforing using max SNR law.
clear all;
%参数赋值
lmda=0.24;             %L波段
N=8;                     %阵元数
d=0.12;                   %阵元间距
DOA_sig=45;                   %目标源来波方向与阵法线方向夹角（单位：度）
DOA_inte=-20;                 %干扰方向（单位：度）         
u=1;                     %任意常数
Pinte=100;                  %干扰源功率
deltn=1;                  %噪声均方根
%产生干扰方向矢量(steering vector)
sv_inte=[];
for i=1:N
    sv_inte=[sv_inte;exp(j*2*pi*d*(i-1)*sin(DOA_inte*pi/180)/lmda)/sqrt(N)];
end
%产生目标方向矢量(steering vector)
sv_sig=[];
for i=1:N
    sv_sig=[sv_sig;exp(j*2*pi*d*(i-1)*sin(DOA_sig*pi/180)/lmda)/sqrt(N)];
end
Ri=N*Pinte*sv_inte*sv_inte';       %干扰协方差矩阵 
Rn=deltn^2*eye(N)+Ri;   %（干扰＋噪声）协方差矩阵
Wopt=u*inv(Rn)*sv_sig;     %最优权值
%自适应阵列的方向图
arf=[];
for sita=-pi/2:1/100:pi/2; 
    sv_sita=[];
    for i=1:N
        sv_sita=[sv_sita;exp(j*2*pi*d*(i-1)*sin(sita)/lmda)/sqrt(N)];
    end
    arf=[arf;Wopt'*sv_sita];
end  
%绘图
sita=-pi/2:1/100:pi/2;
sita_deg=sita*180/pi;
arf_log=20*log10(abs(arf));
plot(sita_deg,arf_log)
axis([-90,90,-60,0])
title('自适应阵列方向图')
xlabel('方位角（度）')
ylabel('增益（dB）')
hold on
stem(DOA_inte,-60,'*')
stem(DOA_sig,-60,':')