fangzhen.m

基于matlab的目标回波信号的相干积累过程分析与仿真。其中包括目标的积累完成过程。

c=3e8;         % 光速
Te=2e-4;       % 发射脉宽
B=1e6;         %带宽
mu=B/Te;       %调频斜率
Ts=1/(2*B);
R0=70e3;        % 起始距离
f0=0;           % 中心频率
Vr=612;          %径向速度
t=0:Ts:Te-Ts;
W=exp(j*pi*mu*t.^2);
Wf=fft(W,1024);
nnn=fix((R0-30e3)/75);%采样的起始位置,从30km开始采样
R=0:75:15e3-75;    %在30km和45km之间采样，采样间隔75m
for i=1:200
    for k=1 :64
        Ri(k,i)=R(i)-Vr*Ts*(k-1);
    end
end
taoi=2*Ri/c;
echo=10*(randnr(64,1024)+j*randn(64,1024));
j=sqrt(-1); 
for i=1:64    %回波信号由噪声与信号组成，信号加入重频和多普勒引起的相位%  差，其中0.096只是换算系数，周期由24变为10(最大不模糊
%  速度)
  echo(i,nnn:nnn+199)=echo(i,nnn:nnn+199)+....
 exp(-j*2*pi*f0*taoi(i,:)+j*pi*mu*taoi(i,:).^2).*exp(j*2*pi*(2*Vr/0.03)*0.096*i/64);
end            
for i=1:64     %脉冲压缩
     sp2(i,:)=ifft(fft(echo(i,:),1024).*conj(Wf),1024);
 end
 for k=1:1024      % 相干积累，用FFT做相干积累，既加快运算速度，又使%  得重频引起的相位差得到补偿
     sct(:,k)=abs((fft(sp2(:,k),64)));
 end
 sct=sct./max(max(sct));%归一化 
 sp=sp2./max(max(sp2));%归一化
 figure(1)
 plot(20*log10(abs(sp')),'b')
 ylabel('dB')
 title('相干积累前')
 axis([1  1024 -30 0])
 figure(2)
 plot(20*log10(sct'),'b')
 ylabel('dB')
 title('相干积累输出')
 axis([1  1024 -30 0])
r=((1:1024)*75+30e3)./1e3;
dp=0:63;
figure(3)
mesh(r,dp,sct)
xlabel('距离 km')
ylabel('Doppler - KHz')
title('相干积累输出结果')
figure(4)
sct(find(sct<10^(-0.4)))=0;
contour(r,dp,sct);grid on
xlabel('距离 km')
ylabel('Doppler -  KHz')
title('R-fd 等高线')