lfm1_mp_2.m

信号与信息处理——阵列信号处理DOA估计的matlab算法

tic
  %  the  data of signal
clear;
format short;
c=3*10.^8;
L=6;
N=128*2;
fm1=0.25;%fc1=0.275;
Km1=0.025;
Phia=0;
f0=0.275;
lamta=c/f0;
len=lamta/2;
%len=1/2;
P1=10*pi/180;
snr1=-10;
Amp1=sqrt(2*10^(snr1/10));
t=1:N;
sig1=Amp1.*exp(j*2*pi*(Km1/2*(t.^2)/N+fm1*t+Phia));
f1=1*Km1*t/N+fm1;
s=[sig1];
  for i=1:N;
      l=1:L;  
      x1=exp(j*2*pi*(len*f1(1,i)*(l-1)*sin(P1))/c);
      a1=x1.';
      a=[a1];
      S(:,i)=a*s(:,i);
  end
  %a=[a1.';a2.'];
   
 noise=randn(L,N)+j*randn(L,N);
 z=S+noise;

 [gamma,rr,atom]=mp_f_K_Pi(N,z(1,:));   
 gamma

 doa_range1=[0 90 1];
 doa_doa1=mp_doa(z,len,gamma,doa_range1)
 
 doa_range2=[doa_doa1-1 doa_doa1+1 0.01];
 doa_doa2=mp_doa(z,len,gamma,doa_range2)
 toc