高速目标距离走动校正 keystone 插值算法

高速目标距离走动校正，matlab工程代码，保证实用，毕设必备。假如有人需要，话可以上传keystone-DFT-IFFT算法实现的与keystone-CZT-IFFT算法资源。

clear;clc;close all
c=3e8;
f0=35e9;
Tr=100e-6;
fr=1/Tr;
Tp=20e-6;
B=20e6;
fs=2*B;
Ts=1/fs;
R0=4000;
v0=8000;
u=B/Tp;
lambda=c/f0;
SNR=-20;
Rmin=c*Tp/2;
Rmax=c*Tr/2;
a0=0;
delta_R=c/2/fs;
v_max=lambda/2*fr;
M1=fix（Tr*fs）;           
number=fix（Tp*fs）;  
M=2^（ceil（log2（M1+number-1）））;
N=64;

tt=-Tp/2:1/fs:（Tp/2-1/fs）;
sref=exp（1j*2*pi*（0.5*u*tt.^2））; 
ha=hamming（length（sref））;
s_fft=fft（conj（fliplr（sref））.*ha‘M）;
s_fft=repmat（s_fft.‘1N）;
t=-Tp/2+1/fs:1/fs:（Tr-Tp/2）;

tm=（0:N-1）‘*Tr;
vt=v0+a0*tm;
R=R0+vt.*tm+0.5*a0*tm.^2;

fd=2*v0*f0/c;
tau=2*R/c;
sr1=zeros（NM1）;

num_mohu=fix（fd/fr）;

for kk=1:N
    td=t-tau（kk）;
    huibo=10^（SNR/20）*（abs（td）<=Tp/2）.*exp（1j*pi*u*td.^2）.*exp（-1j*2*pi*f0*tau（kk））;
    sr1（kk:）=huibo+0*0.707*（randn（1length（t））+1j*randn（1length（t）））;
end

yt_fft=fft（sr1.‘M）;
yy=yt_fft.*s_fft;
y_out=（ifft（yy））.‘;
pc_out=y_out（:number:M1+number-1）;

figure;plot（abs（pc_out（1:）））;hold on;plot（abs（pc_out（N:））‘r--‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2423  2018-12-11 17:14  keystone_sinc.m