用sinc插值和最近领域插值完成距离弯曲校正的完整程序和几篇经典弯曲校正文章

用sinc插值和最近领域插值完成距离弯曲校正的完整程序和几篇弯曲校正文章

%%=========================================================================
%% Filename:  
%% Help file: .doc
%% Project: 
%% Author:    
%% E-mail:      
%%=========================================================================
%%=========================================================================
clear; clc; close all;
%%=========================================================================
%% Parameter--constant  
C = 3e8;               %speed of light
Fc = 1e9;              %carrier frequency
lambda = C/Fc;          %wavelength
Rfactor = 2;         % 距离维采样倍数因子
Afactor = 1;         %
%%Parameter--resolution  距离维和方位向的分辨率要求是多少
DY = 5;                   %range resolution  DY = C/2/Br
DX = 1;                   %cross-range resolution  DX = D/2  分辨率都是两名
%%Parameter--target area 观察场景相关参数
Xmin = 0;              %Targets area in azimuth is within [Xmin Xmax]
Xmax = 200;
Yc = 10000;             %Center of imaged area 场景中心线 
Yw = 300;              %Target area in range is within[Yc-YwYc+Yw]
                       %imaged width is 2*Yw
%%Parameter--orbital information 载机平台相关参数
V = 100;               %SAR velocity 100m/s 载机速度
H = 6000;              %height of airborne 6000m  载机运行高度
R0 = sqrt（H^2+Yc^2）;   %载机的中心垂直距离 是常数   Yc是场景中心线距离
theta=90*pi/180;	   % squint angle  斜视角
%%Parameter--antenna   天线参数
D = 2*DX;              %antenna length in azimuth direction  由方向向分辨率决定
Lsar = lambda*R0/D;    %SAR integration length  合成孔径长度
Tsar = Lsar/V;         %SAR integration time    合成孔径时间

%%Parameter--fast-time domain  快时间域参数设置
Tr = 10e-6;                 %pulse duration 10us  发射脉冲宽度 
Br = C/（2*DY）;              %chirp frequency modulation bandwidth 30MHz LFM信号的带宽 由距离维分辨率决定
Kr = -Br/Tr;                 %chirp slope         LFM信号的调频率 这里为了看回波采用负调频
Fsr = Rfactor*Br;           %sample frequency in fast-time domain  距离维采样频率 由奈奎斯特原理知要大于等于带宽
dt = 1/Fsr;                 %sample spacing in fast-time domain  距离维的采样间隔


Rmin = sqrt（（Yc-Yw）^2+H^2）;               %最小斜距在距离维   
% Rmax = sqrt（（Yc+Yw）^2+H^2+（Lsar/2）^2）;  %最大斜距在距离维
Rmax = sqrt（（Yc+Yw）^2+H^2）;               %最大斜距在距离维
Rwid = Rmax-Rmin;                         %也就是波门的范围
Twid = 2*Rwid/C;
Nfast = ceil（Twid/dt+Tr/dt）;              %sample number in fast-time domain 雷达成像技术书P23页说明
Nfast = 2^nextpow2（Nfast）;                %Ready for fft
tm=linspace（2*Rmin/C2*Rmax/C+TrNfast）;  %discrete time array in fast-time domain    离散化
dt=（2*Rmax/C+Tr-2*Rmin/C）/Nfast;          %修正和更新采样间隔
Fsr=1/dt;

%%Parameter--slow-time domain 慢时间域参数设置
Ka = -2*V^2/lambda/R0;          %doppler frequency modulation rate 多普勒频率调制率
Ba = abs（Tsar*Ka）;             %doppler frequency modulation bandwidth 多普勒频率调制带宽
PRF = Afactor*Ba;              %pulse repetition frequency 以Afactor倍的带宽去采样  也就是采样频率
PRT = 1/PRF;                   %pulse repetition interval 也就是采样间隔 
ds =  PRT;                     %sample spacing in slow-time domain  采样间隔
                               %因为我们知道慢时间域采样频率为PRF  
   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-01-04 18:38  校正RCM完整程序和资料\
     文件        8402  2011-01-03 22:11  校正RCM完整程序和资料\range_migration.m
     文件      860955  2010-12-05 15:35  校正RCM完整程序和资料\SAR成像处理中插值算法的比较研究.pdf
     文件      330556  2010-12-05 15:36  校正RCM完整程序和资料\SAR距离_多普勒成像算法中的距离徙动及校正.pdf
     文件       94852  2010-12-05 15:35  校正RCM完整程序和资料\几种SAR距离弯曲校正算法比较分析.pdf
     文件      375465  2010-12-05 15:35  校正RCM完整程序和资料\距离_多普勒方法中的几种插值算法比较.pdf