【雷达通信】基于matlab线性调频脉冲雷达仿真【含Matlab源码 2104期】

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

⛄一、线性调频脉冲压缩雷达仿真研究

1 LFM脉冲压缩雷达工作原理
LFM脉冲压缩雷达的工作原理图如图1所示。
如果将雷达天线和目标看作一个系统,可得到如图2的等效LTI(线性时不变)系统。
等效LTI系统的冲击响应可写成:

式中,M为目标的个数,σi为目标散射特性,τi是光速在雷达与目标之间往返一次的时间:

式中,Ri为第i个目标与雷达的相对距离。
雷达发射信号s(t)经过该LTI系统,得到输出信号(即雷达的回波信号)sr(t):

2 LFM信号
脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。采用宽脉冲发射以提高发射的平均功率,保证足够大的作用距离;而接收时,采用相应的脉冲压缩算法获得窄脉冲,以提高距离分辨率,较好地解决雷达作用距离与距离分辨率之间的矛盾。

脉冲压缩雷达最常见的调制信号是线性调频(Linear Frequency Modulation,LFM)信号,接收时采用匹配滤波器(Matched Filter,MF)压缩脉冲。

LFM(Chirp)信号的数学表达式为:

⛄二、部分源代码

%% ==============================================================================%
close all; %关闭所有图形
clear all; %清除所有变量
clc;
%% 雷达参数设置
% ===================================================================================%
C=3.0e8; %光速 (m/s)
RF=3.140e9/2; %雷达射频
Lambda=C/RF; %雷达工作波长
PulseNumber=16; %回波脉冲数
BandWidth=2.0e6; %发射信号带宽2G
TimeWidth=42.0e-6;%发射信号时宽42us
PRT=240e-6; %雷达发射脉冲重复周期(s)=240us，对应1/2240300=36000米
PRF=1/PRT;
Fs=2.0e6; %采样频率2M
NoisePower=-12;%(dB);%噪声功率（目标为0dB）
% ---------------------------------------------------------------%
SampleNumber=fix(FsPRT); %计算一个脉冲周期的采样点数480；
TotalNumber=SampleNumberPulseNumber;%总的采样点数48016=TotalNumber；
BlindNumber=fix(FsTimeWidth); %计算一个脉冲周期的盲区-遮挡样点数=84；

%=%
%% 目标参数设置
%=%
X=9;Y=2;Z=0; %这三个参数可控制目标的距离和速度

TargetNumber=4; %目标个数

SigPower(1:TargetNumber)=[1 1 0.25 1]; %目标功率,无量纲

TargetDistance(1:TargetNumber)=[2800 8025 8025 9000+(Y*10+Z)200]; %目标距离,单位m 9200需要改9000+(Y10+Z)*200

DelayNumber(1:TargetNumber)=fix(Fs2TargetDistance(1:TargetNumber)/C);% 把目标距离换算成采样点（距离门）

TargetVelocity(1:TargetNumber)=[50 -100 0 (200+X10+Y10+Z)]; %目标径向速度 单位m/s 230需要改为(200+X10+Y10+Z)

TargetFd(1:TargetNumber)=2*TargetVelocity(1:TargetNumber)/Lambda; %计算目标多卜勒

%%
%% 产生线性调频信号
%%
number=fix(Fs*TimeWidth);%回波的采样点数=脉压系数长度=暂态点数目+1=83+1
if rem(number,2)~=0
number=number+1;
end

for i=-fix(number/2):fix(number/2)-1
Chirp(i+fix(number/2)+1)= exp(j*(pi*(BandWidth/TimeWidth)(i/Fs)^{2));%exp(j*pi*ut}2)，u为调制斜率
freq(i+fix(number/2)+1) = (BandWidth/TimeWidth)(i/Fs);%f=ut
end
coeff=conj(fliplr(Chirp));%产生脉压系数，对称反转，复共轭，h(t)=s*(-t)

%% 线性调频信号图示
figure(1);
subplot(2,1,1);
plot(real(Chirp),‘r-’);title(‘线性调频信号的实部(红色)和虚部（蓝色）’);
hold on;
plot(imag(Chirp),‘b-’);
hold off;
subplot(2,1,2);plot(freq);title(‘频率变化曲线’);
%========================================================

%% 产生目标回波串
SignalAll=zeros(1,TotalNumber);%所有脉冲的信号,先填0
for k=1:TargetNumber% 依次产生各个目标1 2 3 4
SignalTemp=zeros(1,SampleNumber);% 一个脉冲
SignalTemp(DelayNumber(k):DelayNumber(k)+number-1)=sqrt(SigPower(k))Chirp;%一个脉冲的1个目标（未加多普勒速度）
Signal=zeros(1,TotalNumber);
for i=1:PulseNumber
Signal((i-1)SampleNumber+1:iSampleNumber)=SignalTemp;
end
FreqMove=exp(j2piTargetFd(k)(0:TotalNumber-1)/Fs);%目标的多普勒速度时间=目标的多普勒相移
Signal=Signal.*FreqMove;
SignalAll=SignalAll+Signal;
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 廖建国,李永,李继杰.线性调频脉冲压缩雷达仿真研究[J].空间电子技术. 2010,7(02)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合