SAR signal processing with MATLAB 书籍 对应代码 matlab

SAR signal processing with MATLAB 书籍和对应代码 matlab

        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        %        CROSS-RANGE IMAGING         %
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


colormap（gray（256））
cj=sqrt（-1）;
pi2=2*pi;
%
c=3e8;              % propagation speed
fc=200e6;           % frequency
lambda=c/fc;        % Wavelength
k=pi2/lambda;       % Wavenumber
Xc=1.e3;            % Range distance to center of target area
%
% Case 1:
L=400;              % synthetic aperture is 2*L
Y0=100;             % target area in cross-range is within [Yc-Y0Yc+Y0]
Yc=0;               % Cross-range distance to center of target area
%
% Case 2:
% L=150;
% Y0=200;
% Yc=.5e3;
%
theta_c=atan（Yc/Xc）;  % squint angle to center of target area
Rc=sqrt（Xc^2+Yc^2）;   % squint range to center of target area
kus=2*k*sin（theta_c）; % Doppler frequency shift in ku domain due to squint
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Program performs slow-time compression to save PRF   %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
Xcc=Xc/（cos（theta_c）^2）;       % redefine Xc by Xcc for squint processing
%
du=（Xcc*lambda）/（4*（Y0+L））;    % sample spacing in aperture domain
duc=（Xcc*lambda）/（4*Y0）;       % sample spacing in aperture domain
                               % for compressed signal
DY=（Xcc*lambda）/（4*L）;         % Cross-range resolution
%
L_min=max（Y0L）;               % Zero-padded aperture
%
% u domain parameters and arrays for compressed signal
%
mc=2*ceil（L_min/duc）;             % number of samples on aperture
uc=duc*（-mc/2:mc/2-1）;            % synthetic aperture array
dkuc=pi2/（mc*duc）;                % sample spacing in ku domain
kuc=dkuc*（-mc/2:mc/2-1）;          % kuc array
%
dku=dkuc;                         % sample spacing in ku domain
%
% u domain parameters and arrays for Synthetic aperture signal
%
m=2*ceil（pi/（du*dku））;            % number of samples on aperture
du=pi2/（m*dku）;
u=du*（-m/2:m/2-1）;                % synthetic aperture array
ku=dku*（-m/2:m/2-1）;              % ku array
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%          SIMULATION      %%%%%%%          
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
ntarget=4;             % Number of targets
%
% Targets‘ coordinates and reflectivity

yn（1）=0;             fn（1）=1;
yn（2）=.7*Y0;         fn（2）=1;
yn（3）=yn（2）-4*DY;    fn（3）=1.;
yn（4）=-.8*Y0;        fn（4）=1;
%
s=zeros（1mc）;         % Measured Echoed Signal （baseband）
for i=1:ntarget;
 dis=sqrt（Xc^2+（Yc+yn（i）-uc）.^2）;
 s=s+fn（i）*exp（-cj*2*k*dis）.*（abs（uc） <= L）;
end;
%
s=s.*exp（-cj*kus*uc）;     % Slow-time baseband conversion for squint
fs=fty（s）;                % aliased aperture signal spectrum
%
plot（ucreal（s））
xlabel（‘Synthetic Aperture u meters‘）
ylabel（‘Real Part‘）
title（‘Aliased Aperture Signal‘）
axis（[uc（1） uc（mc） 1.1*min（real（s）） 1.1*max（real（s））]）;
axis（‘square‘）
print P2.1a.ps
pause（1）
%
plot（kuc+kusabs（fs））
xlabel（‘Synthetic Aperture Frequency ku rad/m‘）
ylabel（‘Magnitude‘）
title（‘Aliased Aperture Signal Spectrum‘）
axis（‘square‘）
prin

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6377  1999-02-19 09:10  SAR signal processing with MATLAB\soumekh\crange.m

     文件        182  1999-02-19 09:10  SAR signal processing with MATLAB\soumekh\ftx.m

     文件        188  1999-02-19 09:10  SAR signal processing with MATLAB\soumekh\fty.m

     文件        183  1999-02-19 09:10  SAR signal processing with MATLAB\soumekh\iftx.m

     文件        190  1999-02-19 09:10  SAR signal processing with MATLAB\soumekh\ifty.m

     文件      10037  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.10a.ps

     文件      35201  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.1a.ps

     文件      36632  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.2a.ps

     文件      29708  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.3a.ps

     文件      29193  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.4a.ps

     文件      29566  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.5a.ps

     文件      10045  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.6a.ps

     文件      30692  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.7a.ps

     文件       8679  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.8a.ps

     文件      35943  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P1.9a.ps

     文件      13981  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.10a.ps

     文件       8568  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.1a.ps

     文件       8957  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.2a.ps

     文件       8633  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.3a.ps

     文件       9687  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.4a.ps

     文件      16384  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.5a.ps

     文件      13895  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.6a.ps

     文件      14380  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.7a.ps

     文件      13736  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.8a.ps

     文件      10078  1999-02-19 10:44  SAR signal processing with MATLAB\soumekh\P2.9a.ps

     文件      27918  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P3.1a.ps

     文件      11300  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P3.2a.ps

     文件      27231  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P3.3a.ps

     文件      13699  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P3.4a.ps

     文件      26504  1999-02-19 10:45  SAR signal processing with MATLAB\soumekh\P3.5a.ps

............此处省略34个文件信息