二维MUSIC算法仿真.m

clear all;close all;clc;
sensor_number=8;
source_number=2;
snapshort_number=250;
source_theta=[ -3050];
source_fei=[20-60];
snr=10;
A=zeros（sensor_numbersource_number）;
k=0:sensor_number-1;
 for i=1:source_number
   A（:i）=exp（-j*pi*cos（2*pi*k‘/sensor_number-source_theta（i）*pi/180）*sin（source_fei（i）*pi/180））;
 end

length=500;
M1=2;
M2=4;
S1=randi（[0 M1-1]length1）;
S2=randi（[0 M2-1]length1）;
ps1=pskmod（S1M1）;
ps2=pskmod（S2M2）;
t=151:150+snapshort_number;
S=[ps1（t） ps2（t）]‘;
s=A*S;
noise=randn（sensor_numbersnapshort_number）+j*randn（sensor_numbersnapshort_number）;
ss=s+noise;
%ss=awgn（ssnr）;

Rs=ss*ss‘/snaps