我采用经典MUSIC法进行频率估计,得不到正确结果。谐波模型为:x(n)=0.1*sin(2*pi*1*n)+0.06*sin(2*pi*1.03*n)+0.05*sin(2*pi*2*n)+0.03*sin(2*pi*2.06*n)+0.01×w(n+1*loop),请高手帮忙!程序见附件。
M=100; %等距阵线的接收点个数
N=100; %快拍次数
end_f=2.5;%频率搜索的范围
count=0;
for delta_f =0:0.001:end_f
count=count+1;
end
loops=5;
Pw=zeros(count,loops);
for loop=1:1:loops
p=4;
% while p~=4
%生成观测数据
%w = wgn(1,2000,0);
w = randn(2000,1);
for n=1:1:300
x(n)=0.1*sin(2*pi*1*n)+0.06*sin(2*pi*1.03*n)+0.05*sin(2*pi*2*n)+0.03*sin(2*pi*2.06*n)+0.01*w(n+1*loop);
end
X=zeros(M,N);
for Nn=1:1:N
X(:,Nn)=x(Nn:Nn+M-1)';
end
%求相关矩阵Rxx
Rxx=zeros(M,M);
for Nn&