CAPON参数估计

%************CAPON算法**************%
%%
clear all;
close all;
clc;
tic;
%% 算法
K=3;%信号源的数量
L=10000;%采样点数（脉冲积累数）   脉冲积累越大，则估计经度越高，分别率也高
M=10;%发射阵元数
N=10;%接收阵元数
thetaa=[-30 -10 20];%三个目标的真实波达角
snr=0;%dB
f0=10e6;%信号发射频率
C=3e8;
fs=200;%采样频率
Ts=1/fs;
lamda0=C/f0;%发射波长
dt=0.5*lamda0*M;
dr=0.5*lamda0;
%**************产生信号**************%
u0=10;u1=10;u2=10;
t=（0:L-1）*Ts;
s0=exp（j*2*pi*（f0*t+0.5*u0*t.^2））;      %此处默认每个发射阵对于同一目标发射相同的信号
s1=exp（j*2*pi*（f0*t+0.5*u1*t.^2））;
s2=exp（j*2*pi*（f0*t+0.5*u2*t.^2））;
S=[s0;s1;s2];
%*************生成联合导向矢量矩阵************%
at=[0:M-1]‘;
ar=[0:N-1]‘;
at_theta0=exp（-j*2*pi*dt*sin（thetaa（1）*pi/180）/lamda0*at）;
at_theta1=exp（-j*2*pi*dt*sin（thetaa（2）*pi/180）/lamda0*at）;
at_theta2=exp（-j*2*pi*dt*sin（thetaa（3）*pi/180）/lamda0*at）;
ar_theta0=exp（-j