%************CAPON算法**************%
%%
clear?all;
close?all;
clc;
tic;
%%?算法
K=3;%信號(hào)源的數(shù)量
L=10000;%采樣點(diǎn)數(shù)（脈沖積累數(shù)）???脈沖積累越大，則估計(jì)經(jīng)度越高，分別率也高
M=10;%發(fā)射陣元數(shù)
N=10;%接收陣元數(shù)
thetaa=[-30?-10?20];%三個(gè)目標(biāo)的真實(shí)波達(dá)角
snr=0;%dB
f0=10e6;%信號(hào)發(fā)射頻率
C=3e8;
fs=200;%采樣頻率
Ts=1/fs;
lamda0=C/f0;%發(fā)射波長(zhǎng)
dt=0.5*lamda0*M;
dr=0.5*lamda0;
%**************產(chǎn)生信號(hào)**************%
u0=10;u1=10;u2=10;
t=（0:L-1）*Ts;
s0=exp（j*2*pi*（f0*t+0.5*u0*t.^2））;??????%此處默認(rèn)每個(gè)發(fā)射陣對(duì)于同一目標(biāo)發(fā)射相同的信號(hào)
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];
%*************生成聯(lián)合導(dǎo)向矢量矩陣************%
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