clc;
close all;
clear all;
c=3e8; %光速
f0=10e9; %载频10GHZ
lamd=c/f0; %波长
Tp=1e-6; %脉宽1us
B=100e6; %带宽100M
PRF=3600; %脉冲频率3600hz
fs=120e6; %采样120M
K=B/Tp; %调频斜率
v=150; %雷达速度150m/s
x1=0;y1=0; %目标1坐标（0,0）
x2=30;y2=0; %目标2坐标（30,0）
x3=0;y3=30; %目标3坐标（0,30）
Na=8192; %脉冲个数8192
bet=18/180*pi; %波束夹角18°
R0=1000; %雷达距离场景中心1000m
Nr=256; %距离向采样点数
Tr=Nr/fs; %距离向采样的总时间
Ta=Na/PRF; %方位向总的时间
tr=2*R0/c+linspace(-Tr/2,Tr/2,Nr); %回波的快时间
ta=linspace(-Ta/2,Ta/2,Na); %回波的慢时间
xr=v*ta; %方位向值
yr=R0; %距离向值
fr=fftshift(linspace(-fs/2,fs/2,Nr)); %距离向频率
fa=linspace(-PRF/2,PRF/2,Na); %方位向频率
r=tr*c/2; %径向距离
%产生回波
rec_signal=zeros(Nr,Na);
for i=1:Na
puls_signal=zeros(1,Nr); %回波脉冲
r1=sqrt((y1+yr).^2+(xr(i)-x1).^2); %距离目标1的距离
cita1=atan(abs(x1-xr(i))/(yr+y1)); %雷达到目标1的角度，与垂线的角度
if (abs(cita1)<bet/2)
puls_signal=exp(j*(pi*K*(tr-2*r1/c).^2-2*pi*f0*2*r1/c)).*(abs(tr-2*r1/c)<Tp/2); %目标1回波
end
r2=sqrt((y2+yr).^2+(xr(i)-x2).^2); %距离目标2的距离
cita2=atan(abs(x2-xr(i))/(yr+y2)); %雷达到目标2的角度，与垂线的角度
if (abs(cita2)<bet/2)
puls_signal=puls_signal+exp(j*(pi*K*(tr-2*r2/c).^2-2*pi*f0*2*r2/c)).*(abs(tr-2*r2/c)<Tp/2);%目标2回波
end
r3=sqrt((y3+yr).^2+(xr(i)-x3).^2); %距离目标3的距离
cita3=atan(abs(x3-xr(i))/(yr+y3)); %雷达到目标3的角度，与垂线的角度
if (abs(cita3)<bet/2)
puls_signal=puls_signal+exp(j*(pi*K*(tr-2*r3/c).^2-2*pi*f0*2*r3/c)).*(abs(tr-2*r3/c)<Tp/2);%目标3回波
end
rec_signal(:,i)=puls_signal.';
end
figure(1);imagesc(abs(rec_signal));
%距离向脉压
tran_signal=exp(-j*pi*fr.^2/K).*(abs(fr)<B/2);
for i=1:Na
rec_signal(:,i)=ifft(fft(rec_signal(:,i)).*tran_signal');
end
figure(2);imagesc(abs(rec_signal));
%插值法消除距离徙动
rec_signal=fftshift(fft(rec_signal,[],2),2); %方位向变为频域
% figure(3);imagesc(abs(rec_signal));
for i=1:Na
delt_r=lamd.^2*R0*fa(i).^2/8/v.^2;
rec_signal(:,i)=interp1(r,rec_signal(:,i),r+delt_r); %插值
end
figure(4);imagesc(abs(rec_signal));
%方位向脉压
for ii=1:Nr
Ka=2*v.^2/lamd/r(ii);
puls_a=exp(-j*pi*fa.^2/Ka).*(abs(fa)<Ka*Ta/2);
rec_signal(ii,:)=ifft(rec_signal(ii,:).*puls_a);
end
figure(5);imagesc(abs(rec_signal));
评论0