1.刚刚接触ISAR,找的matlab代码,想问代码中的R和RT代表什么以及运行以后生成的图怎么理解?
2.代码
clear;
close all;
clc;
f0=10e9; %载频10GHZ
ft=1.6e6;
T=5.12e-3;
t=0:1/ft:(T-1/ft);
M=8192;
B=1.57e9/4;%带宽
c=3e8;
R0=1e3; %目标中心点距离
k=B/T; %调频率k,调频带宽1.57GHZ;
K=menu('选择:','默认演示','自定义演示');
if K==2
disp('输入每一点位置,第一列为距离向y,第二列为方位向x:[y1,x1;y2,x2;......]');
disp('分辨力0.1911m');
R=input('R=');
else
%R=[0,0;0.5,0.5;-0.5,-0.5;-0.5,0.5;0.5,-0.5;1,1;-1,-1;-1,1;1,-1;1.5,-1.5;-1.5,1.5];
R=[4,4;4,3;4,5;3,4;5,4];
R2=[0,0;0,1;1,0;0,-1;-1,0];
R3=[-4,-4;-4,-3;-4,-5;-5,-4;-3,-4];
end
%回波信号形式
N=64;
u=0;
s=zeros(128,64);%生成一个128乘64的0矩阵
for Q=0:0.035*pi/180:(0.035*(N-1)*pi/180)
%for Q=1:0.14*pi/180:(1+0.14*(N-1)*pi/180)
u=u+1;
Rt1=R0+R*[cos(Q),-sin(Q)]';
Rt2=R0+R2*[cos(Q),-sin(Q)]';
Rt3=R0+R3*[cos(Q),-sin(Q)]';
Rt=[Rt1;Rt2;Rt3];
st=sum(exp(j*4*pi/c*Rt*(B*t/T+f0)));
%每一固定转角的信号采样8192点
%hanning
hn=hamming(M);
st1=st.*hn'; %加窗,压低旁瓣;
% figure(1);
%plot(st1);
%s1=fft(st1);
%figure(2);
%f=0:ft/M:(M-1)*ft/M;
%plot(f,abs(s1));
st2=zeros(1,8192);
st=[st,st2];
L=2*M;
s2=fft(st,L); %每一列进行傅立叶变换------实现距离向高分辨
%figure(1);
%plot(l,s2);
%identify the frequency index corresponding to the center of the
%object range field;
s3=zeros(1,128);
for i=1:128
s3(i)=s2(i+5168);
end %取128点,对应距离范围6m;
%64个不同角度所对应的数据阵
s(:,u)=s3';
end
figure(1);
contour(abs(s));
xlabel('方位向分辨单元');
ylabel('距离向分辨单元');
view([0,90]);