✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
⛄ 内容介绍
随着科学技术的发展进步,武器系统受到越来越多的电子干扰威胁,而计算机仿真具有安全,高效,经济等特点,已成为目前及未来雷达与电子对抗领域研究中的重要手段.
⛄ 完整代码
clear all;clc;close all;
fc=3e9; %载波频率
PRF=2000;
Br=5e6; %带宽
fs=10*Br; %采样频率
Tp=5e-6; %脉宽
Kr=Br/Tp; %频率变化率
c=3e8; %光速
lamda=c/fc; %波长
Tr=1/PRF; %脉冲重复周期
N_mc=1.5/60*PRF; %脉冲个数
t=0:1/fs:15*Tp+Tp; %采样时间
N_r=length(t); %采样点数
N_target=5; %目标个数
Rmax=c/2*15*Tp; %目标最大距离
R_t=Rmax*abs(rand(1,N_target)); %目标的距离
RCS_t=10*(exp(i*2*pi*rand(1,N_target))); %目标RCS,幅度为10,相位在(0,2pi)之间随机分布
Vmax=lamda*PRF/2; %目标最大速度
v=Vmax*((1+rand(1,N_target))/2); %目标速度
%% 生成目标矩阵
sr=zeros(N_mc,N_r);
for i=1:N_mc
ta=(i-1)*Tr;
sri=0;
for k=1:N_target
tao=2*(R_t(k)-v(k).*(ta+t))/c;
srj=RCS_t(k).*rectpuls(t-tao-Tp/2,Tp).*exp(-1j*2*pi*fc*tao+1j*pi*Kr.*(t-tao-Tp/2).^2);
sri=sri+srj;
end
sr(i,:)=sri;
end
%% 距离压缩前的回波
tm=(1:N_mc)/PRF;
R=c*t/2;
figure(1);
image(R,tm,255*abs(sr)/max(max(abs(sr))))
figure(2);
plot(t*c/2,abs(sr(1,:)))
%% 距离压缩
st=rectpuls(t-Tp/2,Tp).*exp(1i*pi*Kr*(t-Tp/2).^2);
stf=conj(fft(st));
for i=1:N_mc
sr(i,:)=ifft(fft(sr(i,:)).*stf);
end
figure(3);
image(R,tm,255*abs(sr)/max(max(abs(sr))))
figure;
plot(t*c/2,abs(sr(1,:)))
sr=fft(sr,[],1);
V=linspace(0,PRF,50)*lamda/2;
figure;image(R,V,255*abs(sr)/max(max(abs(sr))))
⛄ 运行结果
⛄ 参考文献
[1]李斌, 王磊, 张少华. 脉冲多普勒雷达Matlab建模与仿真[J]. 光电技术应用, 2020, 35(4):9.