改进的MUSIC算法对信号的DOA估计

% 改进的程序
clc 
clear all 
format long                 %将数据显示为长整型科学计数
N = 200;                    %快拍数
doa = [20 60]/180*pi;       %信号到达角
w = [pi/4 pi/4]‘;           %信号频率
M = 10;                     %阵元数
P = length（w）;              %信号个数
lambda = 150;               %波长
d = lambda/2;               %阵元间距 
snr = 20;                   %信噪比
B = zeros（PM）;             %创建一个P行M列的0矩阵

for k = 1:P 
B（k:） = exp（-j*2*pi*d*sin（doa（k））/lambda*[0:M-1]）; %矩阵赋值
end 
B = B‘; 
xx = 2*exp（j*（w*[1:N]））;    %仿真信号
x = B*xx;  
x = x + awgn（xsnr）;        %加入高斯白噪声
R = x*x‘;