lfm1_mp_f_1.m

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

tic
clear;
format short;
c=3*10.^8;
L=1;
N=100;
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;
  for  nn=1%:2;   
  snr=-10%2*(nn-2)
  Amp=sqrt(2*10^(snr/10));
t=1:N;
sig1=Amp.*exp(j*2*pi*(Km1/2*(t.^2)/N+fm1*t+Phia));
f1=1*Km1*t/N+fm1;
s=[sig1];
num_rig_f=0;
num_rig_K=0;
  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.'];
  time=10;
 for times=1:time
     times   
 noise=randn(L,N)+j*randn(L,N);
 z=S+noise;
 [gamma,rr,atom]=mp_f_K(N,z(1,:)); 
 r_f(times)= gamma(1,1);
 r_K(times)= gamma(1,2);
    if  gamma(1,1)==fm1
        num_rig_f=num_rig_f+1;
    end
        if  round(gamma(1,2).*1000)/1000==Km1
        num_rig_K=num_rig_K+1;
    end
end
r_f;
r_K;
n_err_f(nn)=num_rig_f/time;
n_err_K(nn)=num_rig_K/time;
end
n_err_f
n_err_K

toc