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