radar.m

一种基于匹配滤波的对地面目标成像的算法 属于雷达成像部分的内容

T=9.2e-6;                %脉冲宽度
fs =73.33e6;             %距离向采样频率
B=60.e6;                 %发射信号带宽
c=3.0e8;                 %光速
lamd=0.0308;             %发射信号波长
fc =9.75e9;              %载频
va=105;                  %载机飞行速度
PRF=735;                 %方位向采样频率
Nr=675;                  %距离向采样点数
Na=2200;                 %方位向采样点数
r0=15000;                %斜距
% k=B/T;
% N=round(T*fs);
% n=ceil(-1*N/2):floor(1*N/2);
% f=exp(j*0.5*k*pi*(n/fs).^2);
% f=[zeros(1,167),f,zeros(1,167)];
% subplot(5,1,1),plot(real(f));
% title('原始信号');
% ff=fft(f);
% subplot(5,1,2),plot(abs(ff));
% title('原始信号频谱');
% h=exp(-1*j*0.5*k*pi*((ceil(N/2)-n)/fs).^2);
% h=[zeros(1,167),h,zeros(1,167)];
% y=conv(f,h);
% subplot(5,1,3),plot(abs(y));
% title('匹配过后的信号频谱');
% m=1:Na;
% hh=exp(-1*j*2*pi*va^2*((Na/2-m)/PRF).^2/(lamd*r0));
% yh=fft(hh);
% subplot(5,1,4),plot(abs(yh));
% title('方位向频谱')
% fh=exp(j*2*pi*va^2*(Na/4/PRF-(Na/2-m)/PRF).^2)/(lamd*r0);
% yy=conv(hh,fh);
% subplot(5,1,5),plot(abs(yy));
% title('方位向匹配后的信号');

k=B/T;
t_r=1:Nr;
t_r=(t_r-Nr/2)/fs;                  %距离向采样向量-T/2～T/2，采样频率fs
t_a=1:Na;
t_a=(t_a-Na/2)/PRF;                 %方位向采样向量-Ta/2～Ta/2，采样频率PRF
Rt=sqrt(va^2*t_a.^2+r0^2);          %点目标与相位中心的实时距离
input=zeros(Nr,Na);
for n=1:Na;                         %产生数据
    delta_R=Rt(n)-r0;
    delta=ceil(delta_R/(c/(2*fs)));
    %if delta==1;
%         for m=delta:Nr;
%             input(m,n)=exp(j*(k*pi*(t_r(m))^2-4*pi*Rt(n)/lamd));
%         end
%     else
        for m=(delta+1):Nr;
            input(m,n)=exp(j*(k*pi*(t_r(m-delta))^2-4*pi*Rt(n)/lamd));
        end
%     end
 end
%n=(round(Nr/2)-Nr):floor(Nr/2);
n=1:Nr;
h=exp(-j*k*pi*((-n)/fs).^2);
fh=fft([h,zeros(1,674)]);
for m=1:2200;
   fout(:,m)=(fft([input(:,m)',zeros(1,674)]))';
   fhfout(:,m)=fh'.*fout(:,m);
   out(:,m)=(ifft(fhfout(:,m)'))';
end
% for m=1:2200;
%     out(:,m)=conv(input(:,m)',h)';
% end
t=1:Na;
hh=exp(j*2*pi*va^2*(t/PRF).^2/lamd/r0);
fhh=fft([hh,zeros(1,2199)]);
for m=1:1349
    ffout(m,:)=fft([out(m,:),zeros(1,2199)]);
    fhhfout(m,:)=fhh.*ffout(m,:);
    output(m,:)=ifft(fhhfout(m,:));
end
% for m=1:1349
%     output(m,:)=conv(out(m,:),hh);
% end
%contour(abs(output));
subplot(3,1,1),contour(abs(output));  %找出要求成像点的位置。
title('成像质量图');
subplot(3,1,2),plot(real(output(:,1100)));
title('匹配滤波分析--距离向');
subplot(3,1,3),plot(real(output(675,:)));
title('匹配滤波分析--方位向');