SAR成像的位置误差校正

%%
%说明：场景和模型如参考文献所示。参数设置参见表6.2。只修正雷达运动位置。相对于第二版，将场景元素看成一个整体，加快速度。
%设计思想：
%输入参数：
%输出参数：
%参考文献：《压缩感知成像雷达》（第六章，LFM）《Compressed Sensing Radar Imaging With Compensation of Observation Position Error》

%（C）中国科学技术大学 电磁场与微波技术 编写：李博 2015/10/06

clc;clear;close all;
%%
%表6.2 LFM波形仿真参数
fc=10e9; %中心频率，载频
B=512e6; %带宽
T=1e-6; %脉冲持续时间，信号持续时间
gamma=5.12e14; %调频率
N=1024; %距离采样数，距离向采样数
M=256; %观测数，方位向采样数
Ns=128; %压缩感知选取的距离采样数
Ms=32; %压缩感知选取的观测数
Y=40; %场景方位点数
X=40; %场景距离点数
K=Y*X; %场景的网格总数
R0=30; %场景中心到雷达轨迹的距离
c=3e8; %光速

%%
%场景构建
disp（‘场景构建‘）
X_all= linspace（28.531.5X）; %场景距离方向的网格点
Y_all= linspace（-1.5+1.5Y）; %场景方位方向的网格点
X_k=zeros（K1）;
Y_k=zeros（K1）;
for j_k=1:Y*X
    x=ceil（j_k/Y）;
    y=mod（j_kY）;
    if y==0
        y=Y;
    end
    X_k（j_k）=X_all（x）;
    Y_