MATLAB music分解信号,MUSIC算法信号频率问题求解

clear all

close all

clc

%参数设定

M=10;                                %阵列数量

doa=[-20 40 50]/180*pi;              %波达方向

P=length(doa);

w=[pi/4 pi/3 pi/6].';                %设置角频率

%w=pi/4;

f=(w(1)+w(2)+w(3))/3/2/pi;           %信号频率

%f=1000;

c=3.0*10^8;                          %光速

lambda=c/f;d=lambda/2;               %lambda:波长    d:间距为半波长的等距均匀线阵

snr=[-10 5 20];                      %信噪比

N=128;                               %快拍数(全部阵元在时域的一次采样)

for ii=1:length(snr)

%阵列流型A

for i=1:P

A(:,i)=exp(-j*2*pi*d*[0:M-1]/lambda*sin(doa(i)));

end

%信源模型建立

S=sqrt(10.^(snr(ii)/10)*exp(j*w*[0:N-1]));

%接收信号模型建立

X=A*S+1/sqrt(2)*(randn(M,N)+j*randn(M,N));   %sqrt(2)是指白噪声功率开方

%协方