合成孔径雷达成像matlab程序

clear;clc;close all; %%================================================================ %%常数定义 C=3e8; %光速 %%雷达参数 Fc=1e9; %载频1GHz lambda=C/Fc; %波长 %%目标区域参数 Xmin=0; %目标区域方位向范围[Xmin,Xmax] Xmax=50; Yc=10000; %成像区域中线 Y0=500; %目标区域距离向范围[Yc-Y0,Yc+Y0] %成像宽度为2Y0 %%轨道参数 V=100; %SAR的运动速度100 m/s H=5000; %高度 5000 m R0=sqrt(Yc2+H2); %最短距离 %%天线参数 D=4; %方位向天线长度 Lsar=lambdaR0/D; %SAR合成孔径长度，《合成孔径雷达成像——算法与实现》P.100 Tsar=Lsar/V; %SAR照射时间 %%慢时间域参数 Ka=-2V^2/lambda/R0; %多普勒频域调频率P.93 Ba=abs(KaTsar); %多普勒频率调制带宽 PRF=Ba; %脉冲重复频率，PRF其实为多普勒频率的采样率，又为复频率，所以等于Ba.P.93 PRT=1/PRF; %脉冲重复时间 ds=PRT; %慢时域的时间步长 Nslow=ceil((Xmax-Xmin+Lsar)/V/ds); %慢时域的采样数，ceil为取整函数，结合P.76的图理解 Nslow=2^nextpow2(Nslow); %nextpow2为最靠近2的幂次函数，这里为fft变换做准备 sn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);%慢时间域的时间矩阵 PRT=(Xmax-Xmin+Lsar)/V/Nslow; %由于Nslow改变了，所以相应的一些参数也需要更新，周期减小了 PRF=1/PRT; ds=PRT; %%快时间域参数设置 Tr=5e-6; %脉冲持续时间5us Br=30e6; %chirp频率调制带宽为30MHz Kr=Br/Tr; %chirp调频率 Fsr=2Br; %快时域采样频率，为3倍的带宽 dt=1/Fsr; %快时域采样间隔 Rmin=sqrt((Yc-Y0)2+H2); Rmax=sqrt((Yc+Y0)2+H2+(Lsar/2)^2); Nfast=ceil(2(Rmax-Rmin)/C/dt+Tr/dt);%快时域的采样数量 Nfast=2^nextpow2(Nfast); %更新为2的幂次，方便进行fft变换 tm=linspace(2Rmin/C,2Rmax/C+Tr,Nfast); %快时域的离散时间矩阵 dt=(2Rmax/C+Tr-2Rmin/C)/Nfast; %更新间隔 Fsr=1/dt; %%分辨率参数设置 DY=C/2/Br; %距离向分辨率 DX=D/2; %方位向分辨率 %%点目标参数设置 Ntarget=5; %点目标的数量 %点目标格式[x,y,反射系数sigma] Ptarget=[Xmin,Yc-50DY,1; %点目标位置，这里设置了5个点目标，构成一个矩形以及矩形的中心 Xmin+50DX,Yc-50DY,1; Xmin+25DX,Yc,1 Xmin,Yc+50DY,1 Xmin+50DX,Yc+50*DY,1]; disp(‘Parameters:’) %参数显示 disp(‘Sampling Rate in fast-time domain’);disp(Fsr/Br) disp(‘Sampling Number in fast-time domain’);disp(Nfast) disp(‘Sampling Rate in slow-time domain’);disp(PRF/Ba) disp(‘Sampling Number in slow-time domain’);disp(Nslow) disp(‘Range Resolution’);disp(DY) disp(‘Cross-range Resolution’);disp(DX) disp(‘SAR integration length’);disp(Lsar) disp(‘Position of targets’);disp(Ptarget)