bmbijiaofxt.m

常规波束形成与最小方差法的比较

clc
clear
%--------------------基本参数----------------------------------
M=8;                                                %阵元数
N=256;                                              %快拍数
SNR=50;                                             %信噪比
f=1500;                                             %信号频率
fs=5*f;                                             %采样频率                             
rsj=[0 20];                                       %信号源入射角
Q1=-3;                                             %扫描扇区
Q2=3;
z=0:M-1;
n=1:N;
%------------------产生两个独立信号----------------------------
s1=sin(2*pi*f/fs*n);
s2=cos(2*pi*f/fs*n);                 %s1 s2是两个独立源，取样256
s=[s1;s2];                            %s为信号数据（2x256)
A=exp(-j*pi*z'*sin(rsj*pi/180));      %导向矢量(8x2)
S=(10^(SNR/20))*A*s;                  %信号空间(8x256)
noise=hilbert(randn(M,N))/sqrt(2);    %噪声空间(8x256)
%-----------------构造协方差矩阵------------------------------
X=S+noise;                            %接受数据空间(8x256)
R=1/N*X*X';                           %数据协方差矩阵R(8x8)
%--------------------谱峰搜索（Bartlett、capon)-------------------
R1=inv(R);                            %求R的逆矩阵
k=length(Q1:0.1:Q2);                  %搜索步长为0.1，点数1800
for n=1:k;
    a=exp(-j*pi*z'*sin((Q1+(n-1)*0.1)*pi/180));
    Pcbf(n)=a'*R*a/(a'*a);                   %波束形成算法a（1x8)
    Pmvm(n)=1./(a'*R1*a);                    %最小方差算法
end
pu1=10*log10(abs(Pcbf));
pu2=10*log10(abs(Pmvm));
%－－－－－－－－－－－－－描绘图形－－－－－－－－－－－－－－
figure
subplot(211)
plot(Q1:0.1:Q2,abs(pu1));xlabel('入射角度');ylabel('空间譜')
title('波束形成法')
grid on
subplot(212)
plot(Q1:0.1:Q2,(pu2));xlabel('入射角度');ylabel('空间譜')
title('最小方差法')
grid on