【雷达LFM】矩形脉冲+高斯脉冲和LFM信号仿真【含Matlab源码 2893期】

💥💥💥💥💥💥💥💥💞💞💞💞💞💞💞💞💞Matlab武动乾坤博客之家💞💞💞💞💞💞💞💞💞💥💥💥💥💥💥💥💥
🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀
在这里插入图片描述
🔊博主简介:985研究生,Matlab领域科研开发者;

🚅座右铭:行百里者,半于九十。

🏆代码获取方式:
CSDN Matlab武动乾坤—代码获取方式

更多Matlab信号处理仿真内容点击👇
Matlab信号处理(进阶版)

⛳️关注CSDN Matlab武动乾坤,更多资源等你来!!

⛄一、简介

1 雷达工作原理
雷达是Radar(RAdio Detection And Ranging)的音译词,意为“无线电检测和测距”,即利用无线电波来检测目标并测定目标的位置,这也是雷达设备在最初阶段的功能。典型的雷达系统如图1.1,它主要由发射机,天线,接收机,数据处理,定时控制,显示等设备组成。利用雷达可以获知目标的有无,目标斜距,目标角位置,目标相对速度等。现代高分辨雷达扩展了原始雷达概念,使它具有对运动目标(飞机,导弹等)和区域目标(地面等)成像和识别的能力。雷达的应用越来越广泛。
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
2 线性调频(LFM)信号
脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。这种体制采用宽脉冲发射以提高发射的平均功率,保证足够大的作用距离;而接受时采用相应的脉冲压缩算法获得窄脉冲,以提高距离分辨率,较好的解决雷达作用距离与距离分辨率之间的矛盾。
脉冲压缩雷达最常见的调制信号是线性调频(Linear Frequency Modulation)信号,接收时采用匹配滤波器(Matched Filter)压缩脉冲。
LFM信号(也称Chirp 信号)的数学表达式为:
在这里插入图片描述
3 LFM脉冲的匹配滤波
信号的匹配滤波器的时域脉冲响应为:
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
线性调频脉冲雷达仿真是一种通过使用线性调频脉冲压缩算法来模拟雷达系统的工作原理和性能的方法。在仿真中,可以使用MATLAB中的phased.LinearFMWaveform函数来创建线性调频脉冲波形,并自定义波形特征,如采样率、脉宽、脉冲重复频率、扫频带宽和扫频方向等[1]。通过仿真,可以模拟雷达发射信号经过线性时不变系统后的回波信号,从而实现对雷达的性能评估和优化。

⛄二、部分源代码

%{
模糊函数
paraInput:
signal :1*N matrix
fs :double
Doppler cut: hz
指示多普勒维度上显示位置
默认为中间
dalay cut : ms
指示时延维度上的显示位置
默认为中间
Max_doppler : 多普勒维最大点数
Max_delay : 时延维最大点数
default_cut: 当信号长度大于时延维最大点数时对信号的处理方式
“cut” :直接阶段信号
“resample”:对信号重采样
returntype:

        ambiguity(signal, fs, varargin);
        ambi = ambiguity(signal, fs, varargin); 此时不绘图
        [ambi, timeax] = ambiguity(signal, fs, varargin); 此时不绘图
        [ambi, timeax,fd] ambiguity(signal, fs, varargin); 此时不绘图

        ambi:模糊函数值
        timeax:时延轴
        fd:多普勒频移轴

%}

function varargout = ambiguity(signal, fs, varargin)
p = inputParser;

addOptional(p, “dalay_cut”, 0);
addOptional(p, “Doppler_cut”,0);
addOptional(p, “Max_doppler”,1000,@(x)validateattributes(x,{‘numeric’},{‘nonnegative’}));
addOptional(p, “Max_delay”,1000,@(x)validateattributes(x,{‘numeric’},{‘nonnegative’}));
addOptional(p, “default_cut”,“cut”,@(s)isstring(s));
parse(p,varargin{:});

Max_doppler = floor(p.Results.Max_doppler);
Max_delay = floor(p.Results.Max_delay);

if Max_doppler <= 1 || Max_delay<= 1
ME = MException(‘myComponent:inputError’,‘最大点数太小’);
throw(ME);
end

[m, len] = size(signal);
if m > 1 && len > 1
ME = MException(‘myComponent:inputError’,‘signal must be 1 dimension’);
throw(ME);
elseif m > 1 && len == 1
signal = signal’;
len=m;
elseif m == 0
ME = MException(‘myComponent:inputError’,‘signal is empty’);
throw(ME);
end

if len > Max_doppler
Ndoppler = Max_doppler;
else
Ndoppler = len;
end
if len > Max_delay
Ndelay = Max_delay;
if strcmpi(p.Results.default_cut, “resample”)
signal = resample(signal, Max_delay, len);
elseif strcmpi(p.Results.default_cut, “cut”)
signal = signal(1:Max_delay);
end
else
Ndelay = len;
end
fd=linspace(-fs/2,fs/2,Ndoppler);
t2=linspace(0,len/fs,Ndelay);

ambi = abs(xcorr2(bsxfun(@times,signal,exp(1j2pi*fd’*t2) ),signal));
timeax=[fliplr(-t2),t2(2:end)] * 1e3;

isplot = false;
if nargout == 0
isplot = true;

if isplot
[~,DopplerPos] = min(abs(fd - p.Results.Doppler_cut));
[~,DalayPos] = min(abs(timeax - p.Results.dalay_cut));
freq = fd(DopplerPos);
delay = timeax(DalayPos);

figure(1);
mesh(timeax,fd,ambi),grid on;
xlabel('delay [ms]');
ylabel('frequency[Hz]');
zlabel('magnitude');
title("ambiguity function");

figure(2)
plot(timeax, ambi(DopplerPos,:) )
xlabel('delay [ms]');
ylabel('magnitude');
title(sprintf("dalay cut, Doppler shift freq=%.2fHz",freq) )

figure(3)
plot(fd, ambi(:,DalayPos) )
xlabel('frequency[Hz]');
ylabel('magnitude');  
title(sprintf("Doppler cut,delay=%fms",delay) )

figure(4)
contour(timeax,fd,ambi)
xlabel('frequency[Hz]');
ylabel('delay [ms]');
title("等高线图");

end
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 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值