資源簡介
強力推薦的ISAR程序代碼,自己調試過,可以運行,對學習雷達成像的朋友很有幫助!

代碼片段和文件信息
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仿真程序\
評論
共有 條評論