clc;
clear?all;
%=============基本參數設置=======================
c=3.0e+8;???????????????????%光速
r0=100;???????????????????????%雷達到目標旋轉中心的距離為100m?
w=0.125;?????????????????????%目標旋轉角速度，單位（rad）??????????????????????
Mr=128;?????????????????????%距離向采樣點數
Ma=128;?????????????????????%方位向采樣點數
%=============基本參數設置=======================

%=============調頻信號參數設置====================
f0=10e+9;????????????????????%載頻10GHz
lamda=c/f0;??????????????????%波長
pulsewide=25.6e-6;?????%發射信號的脈沖寬度為25.6微秒
B=400e+6;???????????????????%發射信號的帶寬為400MHz
PRF=400;??????????????????????%?脈沖重復頻率?400Hz
PRT=1/PRF;??????????????????%?脈沖重復周期為2.5毫秒
K=B/pulsewide?;???????????%?調頻斜率
%=============調頻信號參數設置====================

%=====================Load?Target?===============
load?Target_Plane_SeventyFour_Point.mat?posxy
x=posxy（:1）;
y=posxy（:2）;
[M_Target_Plane_SeventyFour_PointN_Target_Plane_SeventyFour_Point]=size（posxy）;
%=====================?Load?Target?===============
%=====================Load?Target?===============
%?%?load?Target_Plane_SeventyFour_Point.mat?posxy
%?x=[-9?-8?-7?-7?-6?-6?-5?-5?-4?-4?-4?-3?-3?-3?-2?-2?-2?-1?-1?-1?...?
%?????0?0?0?0?0?0?0?0?0?0?1?1?2?2?3?3?3?3?3];
%?y=[0?0?0.5?0.5?1?1?1.5?0?-1.5?2?0?-2?2.5?0?-2.5?3?0?-3?3.5?...
%?????0-3.5?4?3.6?2.8?2?1?-1?-2?-2.8?-3.6?-4?0.5?-0.5?1?-1?2?1?0?-1?-2];
%?[M_Target_Plane_SeventyFour_PointN_Target_Plane_SeventyFour_Point]=size（x）;
%?%=====================?Load?Target?===============

sigma=5*ones（1M_Target_Plane_SeventyFour_Point）;??????????????%?為什么全是5？為什么不是目標存在點的σ為1，其他為0？
j=sqrt（-1）;?
t=-pulsewide/2:pulsewide/（Ma-1）:pulsewide/2;
Hamming_Windows=hamming（Mr）;
F=4.0*pi*（f0+K*t）/c;

%=====================開辟矩陣便于保存數據?==========
echos=[MaMr];?
echos_real=[MaMr];
echos_final=[MaMr];
%=====================開辟矩陣便于保存數據?==========

%=====================仿真回波數據?=================成像區域化為Ma*Mr
for?m=1:Ma
????time=（m-Ma/2）*PRT;?????%慢時間
????for?n=1:Mr
????????s1=0.0;
????????for?k=1:M_Target_Plane_SeventyFour_Point
????????????s2=sigma（k）*exp（-j*（（F（n）*（x（k）*cos（w*time）+y（k）*sin（w*time）））））;??%此處為回波
????????????s1=s1+s2;
????????end
????????s1=exp（j*F（n））*s1;??%匹配濾波
????????s1=s1*Hamming_Windows（n）;??%?對行（距離向）加海明窗?將距離向的波形控制在n個采樣點內，
????????echos（mn）=s1;
????end
end
%=====================仿真回波數據?=================

%=====================輸出顯示回波數據?==================
figure（1）imagesc（abs（echos））;colormap（1-gray）;
xlabel（‘距離向采樣點（個）‘）ylabel（‘回波幀數（個）‘）;title（‘仿真回波數據‘）;
set（get（gca‘xlabel‘）‘FontSize‘16）;
set（get（gca‘ylabel‘）‘FontSize‘16）;
set（get（gca‘title‘）‘FontSize‘16）;
set（gca‘FontName‘‘Times?New?Roman‘‘FontSize‘16）;
save?echos???echos??
%=====================輸出顯示回波數據?==================

%=====================將回波數據處理為實測數據?=============
for?m=1:Ma
????for?n=1:Mr
????????echos_real（n）=echos（mn）;
????end
????echos_real=ifftshift（ifft（echos_real））;?????????????%?對行（距離向）傅氏逆變換
????for?n=1:Mr
????????echos_final（mn）=echos_real（n）;
????end
end
%=

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4809??2017-09-30?11:41??點目標的飛機模型ISAR仿真程序\ISAR_Imaging.m
?????文件?????????438??2017-09-27?17:47??點目標的飛機模型ISAR仿真程序\Target_Plane_SeventyFour_Point.mat
?????文件?????????673??2016-01-15?20:43??點目標的飛機模型ISAR仿真程序\Target_Plane_SeventyFour_Point1.m
?????文件??????253804??2017-09-30?11:41??點目標的飛機模型ISAR仿真程序\echos.mat
?????文件????????2192??2017-09-30?11:41??點目標的飛機模型ISAR仿真程序\echos_real.mat
?????目錄???????????0??2017-09-30?16:32??點目標的飛機模型ISAR仿真程序\