%%?ISAR回波模擬+基于壓縮感知的一位距離像
%?Author：潘小義
clear?all;
clc;
close?all;
format?long;
echo?off;
%?%?目標幾何中心是原點，且是剛體運動的轉(zhuǎn)動中心
Omiga?=?0.05*2*pi;%剛體旋轉(zhuǎn)速度
%?d?=?[3?3;-3?3;-3?-3;3?-3];%剛體點相對于旋轉(zhuǎn)中心的坐標?此處目標是定義在目標坐標系中
%?d?=?[-10?2;-9?2;-8?2;-7?2;-6?2;-5?2;-4?2;-3?2;-2?2;-2?3;-2?4;-2?5;-2?6;-1?6;0?3;0?4;0?5;0?6;0?2;1?2;2?2;2?3;2?4;2?5;2?6;2?7;2?8;2?9;2?10;3?10;4?3;4?4;4?5;4?6;4?7;4?8;4?9;4?10;4?2;5?2;6?2;7?2;8?2;9?2;10?2;...
%?????-9?1;-8?0;-7?-1;10?-1;10?0;10?1;...
%?????-6?-2;-5?-2;-4?-2;-3?-2;-2?-2;-1?-2;0?-2;1?-2;2?-2;3?-2;4?-2;5?-2;6?-2;7?-2;8?-2;9?-2;10?-2];
d?=?[-10?2;-8?2;-6?2;-4?2;-2?2;-2?4;-2?6;-1?6;0?4;0?6;0?2;2?2;2?4;2?6;2?8;2?10;3?10;4?4;4?6;4?8;4?10;4?2;6?2;8?2;10?2;...
????-9?1;-8?0;-7?-1;10?0;...
????-6?-2;-4?-2;-2?-2;0?-2;2?-2;4?-2;6?-2;8?-2;10?-2];
figure
plot（d（:1）d（:2）‘*‘）;xlabel（‘cross-range/m‘‘Fontsize‘20）;ylabel（‘down-range/m‘‘Fontsize‘20）;axis（[-15?15?-15?15?]）
set（gca‘Fontsize‘20）;
[RowgColg]?=?size（d）;
theta_o?=?zeros（1Rowg）;%剛體點和目標幾何中心相對初始角度

for?jj=1:Rowg
????if?d（jj1）>0
????????if?d（jj2）>0
????????????theta_o（jj）?=?atan（d（jj2）/d（jj1））;
????????elseif?d（jj2）<0
????????????theta_o（jj）?=?2*pi+atan（d（jj2）/d（jj1））;
????????elseif?d（jj2）==0
????????????theta_o（jj）?=?0;
????????end
????elseif?d（jj1）<0
????????if?d（jj2）>0
????????????theta_o（jj）?=?pi+atan（d（jj2）/d（jj1））;
????????elseif?d（jj2）<0
????????????theta_o（jj）?=?pi+atan（d（jj2）/d（jj1））;
????????elseif?d（jj2）==0
????????????theta_o（jj）?=?pi;
????????end
????elseif?d（jj1）==0
????????if?d（jj2）>0
????????????theta_o（jj）?=?pi/2;
????????elseif?d（jj2）<0
????????????theta_o（jj）?=?3*pi/2;
????????elseif?d（jj2）==0
????????????theta_o（jj）?=?0;
????????end
????end
end
%?以ISAR視線構(gòu)建坐標系
f0?=?10e9;%載頻
c?=?3e8;%
lamda?=?c/f0;
PW?=?100e-6;%單個脈沖持續(xù)時間
PRF?=?1000;
PRI?=?1/PRF;%重復周期
M?=?128;%方位向采樣數(shù)
Tp?=?PRI*M;%總時長
BW?=?500e6;%脈內(nèi)調(diào)制帶寬
k?=?BW/PW;%調(diào)制斜率
N?=?32;%每一個脈沖采樣點數(shù)（距離單元數(shù)）
deltat?=?PW/（N-1）;%采樣間隔
t_all?=?0:deltat:Tp;%全時間
t?=?0:deltat:PW;%快時間

R0R?=?40000;%單站：雷達距離目標幾何中心位置
R0J?=?40000;%?雷達發(fā)射機和接收機的距離
alpha?=?（0）/180*pi;%雙基地角
RRJ?=?sqrt（R0R^2+R0J^2-2*R0R*R0J*cos（alpha））;%接收站距離雷達發(fā)射機位置

det_r?=?c/2/BW;%距離分辨率
det_a?=?lamda/2/Omiga/Tp/cos（alpha/2）;%單站ISAR方位向分辨率
for?ii=1:M%剛體
?????tm?=?（ii-1）*PRI;%慢時間
????for?jj?=?1:Rowg
????????R（jjii）?=?R0R+sqrt（d（jj1）^2+d（jj2）^2）*sin（theta_o（jj）+Omiga*tm）+R0J-sqrt（d（jj1）^2+d（jj2）^2）*cos（3*pi/2-theta_o（jj）-Omiga*tm-alpha）;
????end??
end
figure
plot（R‘）;
%%?回波
sr?=?zeros（MN）;%全回波
A?=?1;%剛體散射點反射強度
for?jj=1:M
????%%?剛體點回波
????for?ii=1:Rowg
????????tdelay（jj:）?=?t-R（iijj）/c;
????????sr（jj:）?=?sr（jj:）+A*exp（sqrt（-1）*2*pi*（f0*tdelay（jj:）+0.5*k.*tdelay（jj:）.*tdelay（jj:）））;
????end????
end
figure;imagesc（abs（sr））
sr?=?sr/max（max（abs（sr）））;%歸一化
%%?Dechirping處理獲取一維距離像
%參考信號模擬
HRRP?=?zeros（MN）;
tref?=?t-（R0R+R0J）/c;%參考距離對應延時
sref??=?exp（sqrt（-1）*2*pi*（f0*tref+0.5*k.*tref.*tref））;%參考信號
for?ii=1:M
????sif（ii:）?=?sr（ii:）.*conj（sref）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4783??2013-01-28?17:10??ISARCSSim_dr.m

-----------?---------??----------?-----??----

?????????????????4783????????????????????1