ISAR成像仿真程序和算法说明

ISAR成像的简单仿真，包括回波建模和仿真，理想情况下的成像，越距离单元校正，运动补偿

clc
clear all
close all

j = sqrt（-1）;
pi2 = 2*pi;
c = 3e8;

%% radar parameters
fc = 10e9;
lamda = c / fc;
Tp = 2e-6;
B = 800e6;
K = B / Tp;
PRF = 1000;

%% target parameters
L = 13;
x = [0 5 10 15... 
	-5 -10 -15...
	5 10 15...
	-5 -10 -15];
y = [0 5 10 15...
	5 10 15...
	-5 -10 -15...
	-5 -10 -15];
% L = 1;
% x = [15];
% y = [15];

scatters = x + j*y;
y0 = 100e3;
fai0 = 10/180*pi;
v = 7.5e3;
x0 = y0*tan（fai0） + （-500:v/PRF:500）;
% bulk location
bulk_pos = x0 + j*y0;
M = length（bulk_pos）;

%% receiving design
% ideal reference distances
Rref = abs（bulk_pos）;
taoref = 2*Rref/c;

% receiving gate
ts_min = -2*15*sqrt（2）/c - Tp/2;
ts_max = 2*15*sqrt（2）/c + Tp/2;
osv = 1.2;				
fs = （ts_max - ts_min - Tp）*K*osv;
ts = ts_min:1/fs:ts_max;
N = length（ts）;

%% backwave simulation
raw = zeros（N M）;
for i = 1:M
	for l = 1:L
		R = abs（bulk_pos（i）+scatters（l））;
		t = ts + taoref（i）;
		tao = 2*R/c;
		win = （t - tao >= -Tp/2 & t - tao < Tp/2）;
		raw（: i） = raw（: i） + ...
			win.‘ .* exp（j*pi2*（fc*（t-tao）+1/2*K*（t-tao）.^2））.‘;
	end
	% dechirp
	sref = exp（j*2*pi*fc*（t-taoref（i））+j*pi*K*（t-taoref（i））.^2）;
	raw（: i） = raw（: i） .* sref‘;		%‘
end

save data1


%% 

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

     文件     902796  2018-03-04 09:44  ISAR成像报告.pdf

     文件       1277  2018-01-24 04:36  code\backwave.m

     文件     647150  2018-01-27 12:15  code\data1.mat

     文件     649394  2018-01-27 12:27  code\data2.mat

     文件       2783  2018-01-24 03:59  code\ideal_case.m

     文件        831  2018-01-24 03:35  code\imagesc_data.m

     文件       1339  2018-01-27 12:23  code\movebackwave.m

     文件       3520  2018-01-27 12:22  code\move_case.m

     文件        333  2018-01-27 12:38  code\readme.txt

     文件        355  2018-01-21 13:30  code\sinc_interpl.m

     目录          0  2018-01-27 12:38  code

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

              2209778                    11