%%三信號基于MUSIC的均勻矩陣平面空間功率譜估計（信號采用AM調制）
tic;
clear?all;
clc;
close?all;
img=sqrt（-1）;
L=15;????????????????????????????????%矩行陣面長寬上的陣元數
W=8;?????????????????????????????????????
element_num=2*（L+W）-4;????????????????%陣元個數;???
sn=4;????????????????????????????????%信號個數
d_lamda=0.6;??????????????????????????%間距波長比???
snr_db=[5101520];
snr=10.^（snr_db/20）;???????????????????%信噪比

f_carrier=[2e6?3e6?4e6??6e6?];????????%載波頻率MHz
f_signal=[2e3?3e3?4e3?6e3?];?????????%信號頻率KHz
fs=20e6;?????????????????????????????%采樣頻率MHz

N=1024;??????????????????????????????%采樣點數
n=1:N;?
theta=0:0.5:180;
fai=0:1:90;


theta0=[70?150?50?120]?;???%信號方位角
fai0=?[25?45?65?65]?;?????%信號仰角
direction_vector=zeros（element_numsn）;%%均勻矩形陣的方向矢量

%%矢量矩陣
L1=zeros（L2）;
L1（:2）=0:L-1;
W2=zeros（W-12）;
W2（:1）=1:W-1;
W2（:2）=L-1;
L3=zeros（L-12）;
L3（:1）=W-1;
for?i=1:L-1
L3（i2）=L-1-i;
end
W4=zeros（W-22）;
for?i=1:W-2
W4（i1）=W-1-i;
end
L_W=[L1;W2;L3;W4];

for?i?=?1:sn
?????direction_vector（:i）?=exp（img*2*pi*d_lamda*（L_W*[cos（theta0（i）*pi/180）*cos（fai0（i）*pi/180）;?sin（the