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('匹配濾波分析--方位向');