CBF算法matlab程序

clc;
clear all;
% ---------------------------参数初始化------------------------------------
M=32;                         % 阵元数目
L=1024;                      % 信号长度/快拍数 
lam=0.15;                    % 信号波长
d=1/2*lam;                   % 阵元间距
angle=[-pi/3 pi/6 35*pi/180];      % 信号入射角度
snr=-10;                             %信噪比
N=3;                         %信源数
% ---------------------------信号模型--------------------------------------
A=zeros（NM）;
for k=1:N
A（k:）=exp（-j*2*pi/lam*（[0:M-1]*d）*sin（angle（k）））;   
end
A=A‘;                                             % 导向向量
SS=2*（randn（1L）+j*randn（1L））;
Wn=[0.36 0.38];
[ba] = fir1（512Wn‘bandpass‘）;      % 设计带通滤波器
SS1=filter（baSS）; 
Wn=[0.38 0.4];
[ba] = fir1（512Wn‘bandpass‘）;      % 设计带通滤波器
SS2=filter（baSS）; 
Wn=[0.4 0.42];
[ba] = fir1（512Wn‘bandpass‘）;      % 设计带通滤波器
SS3=filter（baSS）; 
SS=zeros（NL）;   SS=[SS1;SS2;SS3];
RS=SS*SS‘/L; 
gls=trace（RS）/N;
R=A*RS*A‘+（gls*10^（-snr/10））*eye（M）;             %白噪声协方差矩阵为对角阵
% ---------------------------角度估计--------------------------------------
theta=-pi/2:pi/3600:pi/2;       % ULA估计角度变化的范围和频率选择 
for k2=1:length（theta）      % 角度估计
    AA=exp（-j*2*pi/lam*（[0:M-1]*d）*si