u7.txt

很好的二维阵元域内的MUSIC算法源代码

clear
zhenyuuanshu=16; % Number of array elements
kuaipai=20; % Number of snapshots
Doa=[10 -20 45];%*pi/180; % Sources DOA's in radian
S_number=length(Doa); % Number of sources
Lambda=0.1; % Wavelength of Sources
d=0.5*Lambda; % Sensor spacings
S_x=[0:zhenyuanshu-1]*d; % Sensor position for ULA
A=exp(-j*2*pi/Lambda*S_x'*sin(Doa*pi/180));
%A=exp(j*2*pi/Lambda*S_x'*sin(Doa)); % Array manifold
Snr=5; % Signal-to-noise ratio in dB
Snr=sqrt(10.^(Snr/10));
%fs=3000000000;
%f=[300000000,200000000,400000000,500000000,100000000]';
%f1=f/fs;
%Sig=randn*exp(2*j*pi*f1*(1:Snap));
Sig=sqrt(0.5)*(randn(S_number,Snap)+j*randn(S_number,Snap));
% Signal waveform (unit power)
Noise=sqrt(0.5)*(randn(Sensor,Snap)+j*randn(Sensor,Snap));
% Additive noise (unit power)
X=Snr*A*Sig+Noise; % Complete data generation
R=X*X';%Snap; % Array covariance matrix estimation
[U S V]=svd(R);
En=U(:,S_number+1:Sensor); % Estimation of noise subspace
theta=[-90:0.1:90]*pi/180;
for num=1:length(theta)
A_search=exp(-j*2*pi/Lambda*S_x'*sin(theta(num)));
P_music(num)=abs(1/(A_search'*En*En'*A_search));
end
h=[-90:0.1:90];
%plot(h,10*log10(P_music));
plot(h,10*log10(P_music/max(P_music)));
axis([-100 100 -50 +10])
xlabel('入射角/（度）','fontsize',10)
ylabel('Pmusic','fontsize',10)
%text(80,10*log10(P_music/max(P_music)),'S3','fontsize',10)
grid
%plot(theta,10*log10(P_music/max(P_music)));