基于matlab 从接收脉冲中提取波形参数

news2024/9/25 20:02:22

一、前言

现代飞机通常随身携带雷达警告接收器 (RWR)。RWR检测到雷达发射,并在雷达信号照射到飞机上时警告飞行员。RWR不仅可以检测雷达发射,还可以分析截获的信号并编目信号来自哪种雷达。此示例显示了 RWR 如何估计截获脉冲的参数。该示例模拟了地面监视雷达(发射器)和配备 RWR 的飞行飞机(目标)的场景。RWR拦截雷达信号,从截获的脉冲中提取波形参数,并估计发射器的位置。提取的参数可以被飞机用来采取对策。

此示例需要图像处理工具箱

二、介绍

RWR是一种被动电子战支持系统[1],可及时向飞行员提供有关其RF信号环境的信息。RWR拦截撞击信号,并使用信号处理技术提取有关截获波形特性以及发射器位置的信息。此信息可用于调用对策,例如干扰以避免被雷达检测到。下图描述了雷达和飞机之间的相互作用。

此示例模拟具有地面监视雷达和具有 RWR 的飞机的场景。RWR检测雷达信号,并从截获的信号中提取以下波形参数:

  1. 脉冲重复间隔

  1. 中心频率

  1. 带宽

  1. 脉冲持续时间

  1. 到达方向

  1. 发射器的位置

RWR链由相控阵天线、信道化接收器、包络检测器和信号处理器组成。截获信号的频带由信道化接收器和包络检测器估计,然后将检测到的子带状信号馈送到信号处理器。将波束转向应用于该子带信号的到达方向,并使用伪Wigner-Ville变换和Hough变换估计波形参数。使用到达角和单基线方法,还可以估计发射器的位置。

三、场景设置

假设地面监视雷达在L波段工作,发射3μs脉冲重复间隔为 15 时的持续时间μs.发送的线性调频带宽为30 MHz,载波频率为1.8 GHz。监视雷达位于原点并且是静止的,飞机以200 m / s(~0.6 Mach)的恒定速度飞行。

雷达的发射天线是8×8均匀矩形相控阵,间距L/2在其元素之间。信号从雷达传播到飞机,并由RWR拦截和分析。为简单起见,选择波形为峰值功率为 100 W 的线性 FM 波形。

地面监视雷达不知道目标的方向,因此,它需要扫描整个空间来寻找飞机。通常,雷达会在移动到下一个方向之前在每个方向传输一系列脉冲。因此,在不失去通用性的情况下,此示例假设雷达向零度方位角和仰角传输。下图显示了四脉冲列车到达飞机的时间频率表示。请注意,尽管脉冲序列到达特定的延迟,但第一个脉冲到达的时间延迟与RWR无关,因为它没有知识传输时间,并且必须不断监控其环境。

RWR 配备一个 10 x 10 的均匀矩形阵列,间距为L/2在其元素之间。它在整个 L 波段工作,中心频率为 2 GHz。RWR 侦听环境,并不断将收集的数据馈送到处理链中。

四、RWR 包络检测器

RWR 中的包络检测器负责检测是否存在任何信号。当RWR持续接收数据时,接收器链缓冲并将接收到的数据截断为50μs段。

由于RWR不知道发射波形中使用的确切中心频率,因此它首先使用一组滤波器,每个滤波器调谐到略有不同的RF中心频率,将接收的数据划分为子频段。然后在每个频段应用包络检测器以检查信号是否存在。在本例中,信号被划分为100 MHz带宽的子带。这种操作的另一个好处是,每个子带中的信号可以下采样至100 MHz的采样频率,而不是对RWR覆盖的整个带宽进行采样。

此图显示了滤波器组创建的前四个波段。

接收到的数据 具有三个维度。第一维表示快速时间,第二维表示子带,第三维对应于接收数组的接收元素。对于此示例中使用的 RWR 的 10 x 10 天线配置,您有 100 个接收元件。由于发射功率低,接收噪声高,雷达信号与噪声无法区分。因此,在这些元件上对接收功率求和,以增强信嗅比(SNR),并更好地估计每个子带中的功率。具有最大功率的频段是雷达使用的频段。

五、RWR 信号处理器

尽管所选频段的功率高于相邻频段,但频段内的SNR仍然较低,如下图所示。

现在是一个二维矩阵。第一维表示快速时间样本,第二维是跨 100 个接收天线通道的数据。计算检测到的子带起始频率以找到检测到的信号的载波频率。subData

RWR的下一步是找到无线电波到达的方向。该到达角信息将用于将接收天线波束引导到发射器的方向,并使用单基线方法将发射器定位在地面上。RWR使用二维音乐估计器估计到达方向。波束控制是使用相移波束成形器完成的,以实现信号的最大SNR,从而有助于波形参数提取。

假设地平面是平坦的,并且平行于坐标系的 xy 平面。然后,RWR可以使用飞机高度计读数中的高度信息以及到达方向对发射器的位置进行三角测量。

该图清楚地显示了发射器的位置。

应用波束转向后,天线在信号到达的方位角和仰角方面具有最大增益。这进一步提高了截获信号的信噪比。接下来,使用一种称为伪Wigner-Ville变换与Hough变换耦合的时频分析技术在信号处理器中提取信号参数,如[2]中所述。

首先,使用Wigner-Ville变换推导出截获信号的时间频率表示。

使用人眼,即使由此产生的时间频率表示是嘈杂的,将信号与背景分开并不难。每个脉冲在时间频率平面上显示为一条线。因此,使用时频线的开始和结束,您可以推导出脉冲宽度和脉冲带宽。同样,来自不同脉冲的线路之间的时间为您提供脉冲重复间隔。

要在不依赖人眼的情况下自动执行此操作,请使用 Hough 变换从图像中识别这些线条。Hough变换在存在噪声的情况下表现良好,该变换是对时频信号分析方法的增强。

要使用霍夫变换,必须将时频图像转换为二进制图像。下面对图像执行一些数据平滑,然后使用该函数进行转换。转换阈值可以根据接收器的信号噪声特性和工作环境进行修改。

使用霍夫变换,首先将二进制伪维格纳-维尔图像变换为峰值。这样,您只需要检测图像中的峰值,而不是检测图像中的线。

使用函数提取峰值位置。

使用这些位置,该函数可以重建原始二进制图像中的线条。然后如前所述,这些行的开头和结尾可帮助您估计波形参数。

下面列出了提取的波形特性。它们与事实非常吻合。然后,这些估计可用于对雷达进行编目,并在必要时准备对策。

五、总结

本示例展示了RWR如何使用信号处理和图像处理技术估计截获的雷达脉冲的参数。

六、参考文献

[1] Electronic Warfare and Radar Systems Engineering Handbook 2013, Naval Air Warfare Center Weapons Division, Point Mugu, California.

[2] Stevens, Daniel L., and Stephanie A. Schuckers. “Detection and Parameter Extraction of Low Probability of Intercept Radar Signals Using the Hough Transform.” Global Journal of Researches in Engineering, Jan. 2016, pp. 9–25. DOI.orgDOI.org (Crossref), https://doi.org/10.34257/GJREJVOL15IS6PG9.

七、程序

使用Matlab R2022b版本,点击打开。(版本过低,运行该程序可能会报错)

打开下面的“XXX.mlx”文件,点击运行,就可以看到上述效果。

程序下载:https://download.csdn.net/download/weixin_45770896/87666438

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/495301.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

10倍速度开发贪吃蛇游戏之AI辅助

今天就来聊聊AI代码辅助神器,即便是零基础也能上手,因为实在是太强了,这以后叫程序员们怎么活啊!话不多说,直接上神器 我用的是cursor,其实目前AI辅助代码不止cursor,还有微软家的copilot X这个根植于gith…

JVM-01-JVM知识

1-JVM内存模型 Java开发人员一般情况下,使用过程中,不用关注内存的申请和释放,得益于JVM自动内存分配机制,但是其实是个双刃剑,这可以提升Java开发的效率,但是弱化了开发人员内存管理意识,系统容…

四元数快速入门【Quaternion】

四元数(Quaternion)是用于旋转和拉伸向量的数学运算符。 本文提供了一个概述,以帮助理解在空间导航等应用程序中对四元数的需求。 推荐:用 NSDT场景设计器 快速搭建3D场景。 可以通过多种方式在空间中准确定位、移动和旋转物体。 …

U盘在电脑上读不出来怎么办?详细解决方法在这!

案例:u盘在电脑上读不出来 【不知道为什么,我把u盘插入电脑后电脑完全读不出。我也不知道到底哪里出现了问题,有人可以帮我解答一下吗?】 U盘作为我们生活中经常使用的存储工具,其能帮我们存储大量的文件&#xff0c…

BM54-三数之和

题目 给出一个有n个元素的数组S,S中是否有元素a,b,c满足abc0?找出数组S中所有满足条件的三元组。 数据范围:0≤n≤1000,数组中各个元素值满足 ∣val∣≤100。 空间复杂度:O(n^2),时间复杂度 O(n^2)。 注…

React 第三方插件 —— Cron 表达式生成器(qnn-react-cron)

qnn-react-cron 可以看做 react-cron-antd 的升级版(具体“渊源”可见文档),现有功能如下: 🎉 全面支持 cron:秒、分、时、日、月、周、年 🎉 日及周条件互斥,自动改变响应值 &…

8年测试开发,写给1-3年功能测试的几点建议,满满硬货指导

从15年毕业到现在也从业八年了,普通本科毕业,现在一家互联网公司担任测试部门总监,摸爬打滚,坑坑洼洼也经历了不少。思绪很久决定还是写下这篇,希望对后进的小伙子少走一点弯路。 很多人把职场想得太美好,其…

学node写接口!!!

fs 可以读取文档 fs.readFild() 用于读取文件 第一个参数 路径 第二个参数 "utf8"(默认值) 第三个参数 函数 function(err , dataStr ){ 第一个参数是错误, 第二个参数是正确的可以拿到读取文件里面的值 } fs.writeFile() 用于创建文件添加内容 …

【LeetCode中等】1419.数青蛙

给你一个字符串 croakOfFrogs,它表示不同青蛙发出的蛙鸣声(字符串 “croak” )的组合。由于同一时间可以有多只青蛙呱呱作响,所以 croakOfFrogs 中会混合多个 “croak” 。 请你返回模拟字符串中所有蛙鸣所需不同青蛙的最少数目。…

如何根据参考文献查找原文及详细的文献信息

当我们已知参考文献想要查看下载原文以及相关信息时可以用下面的方法: 例如这篇参考文献:Alsamhi S H, Almalki F, Ma O, et al. Predictive estimation of optimal signal strength from drones over IoT frameworks in smart cities[J]. IEEE Transac…

FPGA实现SDI视频解码PCIE传输 提供工程源码和QT上位机源码加技术支持

目录 1、前言2、我已有的SDI编解码方案3、我已有的PCIE方案4、总体设计思路和方案SDI摄像头Gv8601a单端转差GTX解串SDI解码VGA时序恢复YUV转RGB图像缓存PCIE发送通路SDI同步输出通路 5、vivado工程详解6、驱动安装7、QT上位机软件8、上板调试验证SDI同步HDMI输出验证PCIE输出验…

neo4j结合gds实现最短路径算法

背景:Neo4j自带的cypher语句中的 shortestpath allShortestPaths 返回值内容非常有限,不易处理, 在实际生产环境中可用性极低, 且若带where条件查询时,查询效率极低 因此,使用Neo4j自带的插件如apoc来进行最短路径查询…

使用@Autowired、@Qualifier、@Primary注解自动装配组件

1.Autowired、Qualifier、Primary注解 1.1.Autowired注解 Autowired注解可以对类成员变量、方法和构造函数进行标注,完成自动装配的工作。 package org.springframework.beans.factory.annotation;import java.lang.annotation.Documented; import java.lang.ann…

hexo stellar主题添加运行时间-利用不蒜子

在网站的页脚想添加博客运行的时间以及,查了一下资料,根据下面的这个博主的代码进行了修改。 打造一个舒服的写作环境(Hexo) 下面我们开始修改: stellar/layout/_partial/main/footer.ejs 找到//footer的 这一段 // …

基于logback 实现springboot的日志配置

目录 一、前言 二、使用详解 2.1、打印到文件中 2.2、打印级别控制 2.3、logback 详细配置 2.4、logback 配置文件的组成 2.4.1、<root>标签 2.4.2、<contextName>标签 2.4.3、<property>标签 2.4.4、<appender>标签 2.4.5、<logger&g…

Mysql索引(2):索引结构

1 概述 MySQL的索引是在存储引擎层实现的&#xff0c;不同的存储引擎有不同的索引结构&#xff0c;主要包含以下几种&#xff1a; 索引结构描述BTree索最常见的索引类型&#xff0c;大部分引擎都支持 B 树索引 Hash索引 底层数据结构是用哈希表实现的, 只有精确匹配索引列的…

ShardingProxy分库分表实战

目录 ShardingProxy简述 快速使用 ShardingProxy部署 ShardingProxy使用 ShardingProxy的服务治理 Shardingproxy的其他功能 ShardingProxy的SPI扩展 ShardingSphere总结 ShardingProxy简述 ShardingProxy的功能同样是分库分表&#xff0c;但是他是一个独立部署的服务…

林曦讲座 | 关于“无用之美”的11个闪光

林曦老师的新书《无用之美》即将和大家见面啦~    究竟什么是“无用之美”呢&#xff1f;早在2013年&#xff0c;林曦老师就曾在“一席”做过一场关于“无用之美”的分享。2019年又在微软做过一次《中国人的闲情逸致——无用之美》的演讲。这次将两次演讲的内容精华整理出来…

【LeetCode】14. 最长公共前缀

1. 问题 编写一个函数来查找字符串数组中的最长公共前缀。 如果不存在公共前缀&#xff0c;返回空字符串 “”。 示例 1 输入&#xff1a;strs [“flower”,“flow”,“flight”] 输出&#xff1a;“fl” 示例 2 输入&#xff1a;strs [“dog”,“racecar”,“car”] 输出…

虹科动态 | 2023成都国际工业博览会精彩瞬间回顾

2023成都国际工业博览会&#xff08;CDIIF&#xff09;于中国国际西部博览城圆满落下帷幕。本届成都工博会紧紧围绕“工业引领&#xff0c;赋能产业新发展”主题&#xff0c;聚焦工业自动化、数控机床与金属加工、机器人、新一代信息技术与应用、工业互联网、新材料、节能与工业…