Keystone变换仿真程序

包含一阶keystone和二阶keystone的matlab仿真程序，帮助初学者快速上手。 若不能运行，尽快联系wangwu163@163.com

clear all;
close all;
clc
%mex BPimaging.cpp
%% 仿真SAR模式下动目标成像
% By WW
% 20160904
% 注意接收方式为deChirp接收，不是正交解调
%% LFM信号基本参数设置
paraSig = setSig（1）;

%% 成像参数设置
[paraAzm paraRag paraAirp] = setImg（1 paraSig）;

%% 目标参数设置
paraTar = setTar（1 paraRag.Rc）;

%% 回波仿真 deChirp接收
flag = 0; % 0不加误差，1加误差
flagV = 0;
[echoSig kr ksi hDelta yDelta] = setEcho（flag flagV paraAzm paraRag paraAirp paraSig paraTar）; 
figure imagesc（paraRag.tr paraAzm.ta abs（echoSig）） axis xy
title（‘原始回波‘） xlabel（‘快时间‘） ylabel（‘慢时间‘）

dksi = ksi（2） - ksi（1）;
x = linspace（-1/dksi/2 1/dksi/2 length（ksi））*2*pi;

%% 去除RVP相位
dtr = paraRag.tr（2）-paraRag.tr（1）;
fi  = linspace（-1/dtr/2 1/dtr/2 paraRag.Nr）;
r   = fi*3e8/2/paraSig.kai;
rvpH = ones（paraAzm.Na 1）*exp（-1i*pi*fi.^2/paraSig.kai）;

echoSig = ifftshift（ifft（echoSig [] 2） 2）.*rvpH;  clear rvpH

figure
imagesc（r x 20*log（abs（echoSig）./max（abs（echoSig（:）））） [-60 0]） axis xy colormap jet;
xlabel（‘距离‘）ylabel（‘慢时间‘）title（‘脉压结果‘）;


%% keystone变换 interp2
t = cputime
flagK = 1; % flagK= 0， 一阶keystone； flagK= 1， 二阶keystone。
echoSig = myKeystone（kr paraAzm paraSig echoSig flagK）;
t = cputime-t

figure
imagesc（r x 20*log（abs（echoSig）./max（abs（echoSig（:）））） [-60 0]） axis xy colormap jet;
xlabel（‘距离‘）ylabel（‘慢时间‘）title（‘一阶keystone变换之后脉压结果‘）;
if flagK
    title（‘二阶keystone变换之后脉压结果‘）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-04-08 11:03  keystone\
     文件        1521  2018-04-08 11:08  keystone\MTI_SAR.m
     文件         873  2018-04-08 10:37  keystone\myKeystone.m
     文件        1717  2018-04-08 11:09  keystone\setEcho.m
     文件        1776  2018-04-08 10:27  keystone\setImg.m
     文件         545  2018-04-08 10:26  keystone\setSig.m
     文件         475  2018-04-08 11:09  keystone\setTar.m