wan.m

谱估计源程序：是谱估计中的经典算法esprit,MUSIC等算法的源程序

%music
clear all;
close all;
clc;

sensor_number=8;
source_number=2;
snapshort_number=250;
source_theta=[ -30,50];
source_fei=[20,-60];
snr=10;

A=zeros(sensor_number,source_number);
k=0:sensor_number-1;
 for i=1:source_number
   A(:,i)=exp(-j*pi*cos(2*pi*k/sensor_number-source_theta(i)*pi/180)*sin(source_fei(i)*pi/180));
 end

length=500;
M1=2;
M2=4;
S1=randint(length,1,M1);
S2=randint(length,1,M2);
ps1=pskmod(S1,M1);
ps2=pskmod(S2,M2);
t=151:150+snapshort_number;
S=[ps1(t) ps2(t)]';
s=A*S;
noise=randn(sensor_number,snapshort_number)+j*randn(sensor_number,snapshort_number);
ss=s+noise;
%ss=awgn(s,snr);

Rs=ss*ss'/snapshort_number;
[U,S,V]=svd(Rs);
Vs=U(:,1:source_number);
Vn=U(:,source_number+1:sensor_number);
 
 

for theta=-90:90
for fei=-90:90
k=[0:sensor_number-1]'; 
AA=exp(-j*pi*cos(theta*pi/180-2*pi*k/8)*sin(fei*pi/180));
%AA=exp(-j*pi*(sin(st1*degrad)*cos(dt1*degrad)*cos(2*pi*k/8)+sin(st1*degrad)*sin(dt1*degrad)*sin(2*pi*k/8)));
WW=AA'*Vn*Vn'*AA;
Pmusic(theta+91 ,fei+91)=abs((AA'*AA)./WW);%角谱
end 
end 
%t=0;
%h=0;
%for theta=0:1:90
    %for fei=0:1:90
      %  a=exp(-j*pi*cos(2*pi*k/sensor_number-theta*pi/180)*sin(fei*pi/180))';
       % t=t+1;
        %h=h+1;
       %music(t,h)=abs(1./(a'*Vn*Vn'*a));
   % end
%end
theta=-90:90;
fei=-90:90;
mesh(theta,fei,Pmusic);
title('MUSIC空间谱');
xlabel('方位角'); 
ylabel('俯仰角');
zlabel('角谱');
grid on;