High Resolution Tactical Synthetic Aperture Radar (TSAR) - MATLAB Simulations for Radar Systems Design - page 618

Q_Temp = sin(PHI) + Q_Temp;
end;
IQ_Raw(J,I) = I_Temp + i*Q_Temp;
I_Temp = 0.0;
Q_Temp = 0.0;
end;
end;
%%%%%%%%%% End target return section %%%%%
%%%%%%%%%% Range Compression %%%%%%%%%%%%%
Num_RB = 2*Num_Pulse;
WR = hamming(Num_Pulse);
for I = 1:Num_U
Range_Compressed(:,I) = fftshift(ifft(IQ_Raw(:,I).*WR,Num_RB));
end;
%%%%%%%%%% Focus Range Compressed Data %%%%
dn = (1:Num_U)*du - L/2;
PHI_Focus = -2*pi*(dn.^2)/(Lambda*xn(1));
for I = 1:Num_RB
Temp = angle(Range_Compressed(I,:)) - PHI_Focus;
Focused(I,:) = abs(Range_Compressed(I,:)).*exp(i*Temp);
end;
%Focused = Range_Compressed;
%%%%%%%%%% Azimuth Compression %%%%%%%%%%%%
WA = hamming(Num_U);
for I = 1:Num_RB
AZ_Compressed(I,:) = fftshift(ifft(Focused(I,:).*WA'));
end;
SAR_Map = 10*log10(abs(AZ_Compressed));
Y_Temp = (1:Num_RB)*(c/(2*Num_RB*df));
Y = Y_Temp - max(Y_Temp)/2;
X_Temp = (1:length(IQ_Raw))*(Lambda*xn(1)/(2*L));
X = X_Temp - max(X_Temp)/2;
image(X,Y,20-SAR_Map); %
%image(X,Y,5-SAR_Map); %
axis([-25 25 -25 25]); axis equal; colormap(gray(64));
xlabel('Cross Range (m)'); ylabel('Down Range (m)');
grid
%print -djpeg .jpg