⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源: 【雷达通信】基于matlab雷达前视成像仿真【含Matlab源码 643期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab信号处理(初级版)
备注:
点击上面蓝色字体付费专栏Matlab信号处理(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab信号处理(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、雷达前视成像简介
⛄三、部分源代码
clear all;clc; close all;
% 前视线阵
%% 参数
C=3e8;
fc=10e9;
lambda=C/fc;
v=100;
Tp=1.5e-6;
N=60; %阵元个数
B=100e6;
Fs=120e6;
gama=B/Tp;
Da=0.2;
H=600;
PRF=1000; %阵元切换速度
phi=45/180pi; %俯视角
theta=24/180pi; %波束宽度
d=0.3;% 阵元间距
va=dPRF; %方位速度
nrn=200;
nan=N;
Rs=Hsqrt(2);
%% point targets
% point=[600 0 0; 650 0 0; 600 50 0 ; 600 -50 0; 550 0 0 ];
point=[600 0 0];
point_num=size(point,1);
%%
La=N*(Da+d);
ta=[-N/2:N/2-1]/PRF;
Tstart=2Rs/C-nrn/2/Fs; %起始时刻
Tend=2Rs/C+(nrn/2-1)/Fs; %终止时刻
tr=[Tstart:1/Fs:Tend].‘; %快时间
%%
x_pos=ta.v;
y_pos=ta.va;
z_pos=H;
R0=sqrt(z_pos2+(600-v.*ta).2);
Rr=sqrt(R0.2+(va.*ta).2);
%% echo
x=zeros(nrn,nan,‘single’);
for m=1:N
temp1=zeros(nrn,1);
for num=1:point_num
x_temp=abs(x_pos(m)-point(num,1));
y_temp=y_pos(m)-point(num,2);
z_temp=z_pos-point(num,3);
Rt=sqrt(x_temp2+y_temp2+z_temp^2);
alpha=atan(y_temp./sqrt(x_temp.2+z_temp.2));
if (alpha>=-theta/2 && alpha<=theta/2)
Rt=sqrt(x_temp2+y_temp2+z_temp^2);
Wr=abs(tr-2Rt/C)<=Tp/2;
temp2=Wr.exp(jpigama*(tr-2Rt/C).^2).exp(-j4piRt/lambda);
temp1=temp1+temp2;
end
end
x(:,m)=temp1;
end
figure;imagesc(abs(x));
ka=(-2/lambda).(va^2./R0); %方位调频率
fr=[-nrn/2:nrn/2-1]’/nrn*Fs;
del_ran=Rr-Rs;
%% range compress
fr=[-nrn/2:nrn/2-1]'/nrnFs;
% length=(fr/abs(gama)/Tp);
Hp=exp(jpi*(2v/lambda+(fr.^2)/gama));
%%
x=fftshift(fft(fftshift(x)));
for m=1:N
temp2=exp(j4pidel_ran(m)/C*fr);
temp1=x(:,m).Hp.temp2;
x(:,m)=fftshift(fft(fftshift(temp1)));
end
figure;imagesc(abs(x));
%%
for n=1:nrn
xa(n,:)=fftshift(fft(fftshift(x(n,:))));
end
delta_R=((va.ta).^2)./(2R0);
Curve=delta_R2Fs/C;
delta_RMC=Curve-floor(Curve);%小数部分
% sinc
P=4;
temp=zeros(nrn,nan);
for n=1:nan
for m=P:nrn
for i=-P/2:P/2
if m+Curve(n)+i>nrn
temp(m,n)=temp(m,n)+xa(nrn,n)*sinc((-i+Curve(n)));
else
temp(m,n)=temp(m,n)+xa(m+floor(Curve(n))+i,n)*sinc((-i+delta_RMC(n)));
end
end
end
end
⛄四、运行结果
⛄五、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1]李树锋.基于完全互补序列的MIMO雷达与5G MIMO通信[M].清华大学出版社.2021
[2]何友,关键.雷达目标检测与恒虚警处理(第二版)[M].清华大学出版社.2011
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除