%%
%說明：場景和模型如參考文獻(xiàn)所示。參數(shù)設(shè)置參見表6.2。只修正雷達(dá)運(yùn)動(dòng)位置。相對(duì)于第二版，將場景元素看成一個(gè)整體，加快速度。
%設(shè)計(jì)思想：
%輸入?yún)?shù)：
%輸出參數(shù)：
%參考文獻(xiàn)：《壓縮感知成像雷達(dá)》（第六章，LFM）《Compressed?Sensing?Radar?Imaging?With?Compensation?of?Observation?Position?Error》

%（C）中國科學(xué)技術(shù)大學(xué)?電磁場與微波技術(shù)?編寫：李博?2015/10/06

clc;clear;close?all;
%%
%表6.2?LFM波形仿真參數(shù)
fc=10e9;?%中心頻率，載頻
B=512e6;?%帶寬
T=1e-6;?%脈沖持續(xù)時(shí)間，信號(hào)持續(xù)時(shí)間
gamma=5.12e14;?%調(diào)頻率
N=1024;?%距離采樣數(shù)，距離向采樣數(shù)
M=256;?%觀測數(shù)，方位向采樣數(shù)
Ns=128;?%壓縮感知選取的距離采樣數(shù)
Ms=32;?%壓縮感知選取的觀測數(shù)
Y=40;?%場景方位點(diǎn)數(shù)
X=40;?%場景距離點(diǎn)數(shù)
K=Y*X;?%場景的網(wǎng)格總數(shù)
R0=30;?%場景中心到雷達(dá)軌跡的距離
c=3e8;?%光速

%%
%場景構(gòu)建
disp（‘場景構(gòu)建‘）
X_all=?linspace（28.531.5X）;?%場景距離方向的網(wǎng)格點(diǎn)
Y_all=?linspace（-1.5+1.5Y）;?%場景方位方向的網(wǎng)格點(diǎn)
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_