雷达成像SINC插值程序

function xb=rdsinc（xRsPRFFsvlambda）
%输入 x距离压缩后矩阵  Rs斜距   PRF重频  Fs距离采样率 lambda波长
%输出 弯曲校正后矩阵
%因为正侧视，没有距离走动项，只校正弯曲
[nrnnan]=size（x）;
C=3e8;
for n=1:nrn
    xa（n:）=fftshift（fft（fftshift（x（n:））））;  %对每个距离时间内回波fft 变换到频域 （tr-ta域到tr-fa域）
end
P=4;%4点sinc插值
temp=zeros（nrnnan）;
for n=1:nan
    for m=P:nrn
       r=Rs+（m-nrn/2）*C/2/Fs;
       f=（n-nan/2）*PRF/nan;
       delta_R = （1/8）*（lambda/v）^2*r*f^2;  %根据公式算出弯曲量ΔR
       num_n