SAR成像之TDC算法

clc;
clear all;
close all;
u=zeros(121,121);
z=5;
c=3e8;
o=(2-0)^2+(3-0)^2+(5-3)^2;
m=1;
 for x=-6:0.1:6
     n=1;
  for y=-6:0.1:6
      I=0; 
       for f=40e9:50e6:50e9           
           k=2*pi*f/c;          
           s=exp(-1i*2*k*sqrt(o));%收发对（0，0，3）物体（2，3，5）
           I=s*exp(1i*k*2*sqrt((x-0)^2+(y-0)^2+(z-3)^2))+I;

       end
      u(m,n)=(abs(I))^2; 
      n=n+1;
     
  end
   m=m+1;
 
 end
 
           X=-4:0.5:4; Y =4:-0.5:-4;
           figure;imagesc(X,Y,flipud(u));
           colormap gray
           colorbar('Fontsize',14);
           xlabel('X(m)','FontSize',14);
           ylabel('Y(m)','FontSize',14);