基于MUSIC的均匀矩形阵DOT估计matlab仿真程序

%%三信号基于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