2维ISAR成像基于FEKO仿真模型

FEKO仿真的模型回波用matlab程序进行二维FFT成像的基础程序，雷达成像专业可用

%% 2d FFT成像
clear all%#ok
close all
clc

%% 参数设置
c = 3e8;
%%f的离散设置=======================
f_Min = 8e9;
f_Max = 12e9;
f_Interval = 80e6;
fn = （f_Max-f_Min）/f_Interval+1;
fc = （f_Max-f_Min）/2+f_Min;
B = f_Max-f_Min;
pr=c/B/2;
r_width = c/2/f_Interval;
%%=====================================

%%Phi的离散设置================================方位角观测角度
Phi_Min = -4 / 180 * pi;
Phi_Max = 4 / 180 * pi;
Phi_Interval = 0.2 / 180 * pi;
pn=ceil（（Phi_Max-Phi_Min）/Phi_Interval+1）;
pa=c/fc/2/（Phi_Max-Phi_Min）;
p_width = c/2/Phi_Interval/fc;
%%================================

%% 读取数据
E = zeros（fnpn）;
data_real=load（‘real.txt‘）;
data_imag=load（‘imag.txt‘）;
data=data_real + 1i * data_imag;
for i=1:fn
    E（i:）=data（pn*（i-1）+1:i*pn）;
end

%% FFT成像
ISAR=（fft2（E.‘4*2^nextpow2（pn）4*2^nextpow2（fn）））;
cmx=abs（fftshift（ISAR））;
S=cmx/max（max（abs（cmx）））;
ff=linspace（-（fn-1）/2*pr（fn-1）/2*pr4*2^nextpow2（fn））;
pp=linspace（-（pn-1）/2*pa（pn-1）/2*pa4*2^nextpow2（pn））;
figure（1）;
imagesc（ffpp20*log10（S））
xlabel（‘距离向/m‘ ‘fontsize‘12 ）;                        
ylabel（‘方位向/m‘‘fontsize‘12）;
title（‘回波成像‘）;
colormap（jet）;
axis（[-r_width/2 r_width/2 -r_width/2 r_width/2]）;
axis square
clim = get（gca‘CLim‘）;
set（gca‘CLim‘clim（2） + [-20 0]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      32862  2018-11-26 22:41  2D_FEKO_IFFT\imag.txt

     文件        890  2018-05-18 15:28  2D_FEKO_IFFT\old\cx1.m

     文件        566  2018-10-29 15:54  2D_FEKO_IFFT\old\fft_isar.m

     文件      30142  2018-10-29 11:55  2D_FEKO_IFFT\old\imag.txt

     文件      30118  2018-10-29 11:55  2D_FEKO_IFFT\old\real.txt

     文件      32031  2018-11-26 22:41  2D_FEKO_IFFT\real.txt

     文件       1316  2018-11-26 18:40  2D_FEKO_IFFT\TWO_D_FFT.m

     目录          0  2018-10-29 15:46  2D_FEKO_IFFT\old

     目录          0  2018-10-30 18:35  2D_FEKO_IFFT

----------- ---------  ---------- -----  ----

               127925                    9