用music方法进行信源数估计的MATLAB程序

clearallcloseallclcN=16;%阵元个数En_n=0.5;snap_num=64;%快拍数c=3e+8;%c为光速fc=500e6;%雷达工作频率wavelength=c/fc;%波长d=wavelength/2;%阵元间距theta_0=25*pi/180;%来波方向theta_1=pi/6;theta_2=pi/3;%theta_0=0;%theta_1=pi/9;%theta_2=pi/6;theta=-pi/2:pi/180:pi/2;steer_sig_0=exp(j*2*pi*(0:N-1)'*sin(theta_0)*d/wavelength);%信号导向矢量steer_sig_1=exp(j*2*pi*(0:N-1)'*sin(theta_1)*d/wavelength);steer_sig_2=exp(j*2*pi*(0:N-1)'*sin(theta_2)*d/wavelength);snap_sig_0=steer_sig_0*(randn(1,snap_num)+j*randn(1,snap_num))/sqrt(2);%snap_num个采样数据snap_sig_1=steer_sig_1*(randn(1,snap_num)+j*randn(1,snap_num))/sqrt(2);snap_sig_2=steer_sig_2*(randn(1,snap_num)+j*randn(1,snap_num))/sqrt(2);%input_x=snap_sig_0+snap_sig_1+snap_sig_2;%Rx=(input_x*input_x')/snap_num+En_n*eye(N);input_x=snap_sig_0+1.5*snap_sig_1+1.5*snap_sig_2+En_n*(randn(N,snap_num)+j*randn(N,snap_num))/sqrt(2);Rx=(input_x*input_x')/snap_num;r=rank(Rx);x=1:r;[V,D]=eig(Rx);D=diag(D);M=0;form=1:rifD(m)<1M=M+1;v(:,M)=V(:,m);endendM_0=r-M;Pn=v*v';steer=exp(j*2*pi*(0:N-1)'*sin(theta)*d/wavelength);%搜索矢量form=1:length(theta)S(m)=1/(steer(:,m)'*Pn*steer(:,m));endfiguresubplot(1,2,1)plot(x,D,'O-')title('EVD')xlabel('numberofeigenvalue')ylabel('eigenvalue/dB')gridonsubplot(1,2,2)%plot(theta*180/pi,abs(S))semilogy(theta*180/pi,abs(S))title('Music谱估计')xlabel(['信号源个数为',num2str(M_0),'来波方向分别为',num2str(theta_0*180/pi),'\circ',num2str(theta_1*180/pi),'\circ',num2str(theta_2*180/pi),'\circ'])ylabel('powerofMUSICmethod')gridon