MATLAB环境下使用二维高分辨时频分析方法提取波状分量

news2024/12/23 22:13:38

MATLAB环境下使用二维高分辨时频分析方法提取波状分量(分离混合地震数据)。

为了得到更高的时频分辨率,近年来涌现出了大量的新的时频分析方法。有些以线性和非线性时频分析为基础,有些则另辟蹊径,比如Hilbert-Huang变换HHT。HHT包括经验模态分解EMD和希尔伯特变换两部分,其中经验模态分解是其核心部分。由于经验模态分解存在模态混叠和端点效应,集合经验模态分解和互补集合经验模态分解被提出,集合经验模态分解在一定程度上克服了模态混叠问题,互补集合经验模态分解提高了集合经验模态分解的计算效率,并减少了重构误差。

希尔伯特黄变换虽然能得到瞬时频率信息,但是其数学理论基础不够完备,模态混叠问题难以消除。匹配追踪算法与希尔伯特黄变换都属于贪婪算法,主要用于储层预测,噪声压制,匹配追踪时频分辨率较高,但是计算效率较低。稀疏约束的谱反演方法通过加入正则化约束条件控制稀疏程度,得到高分辨率的时频分析结果,近年来受到了越来越多的关注,稀疏参数的选择非常重要,如果选取不合适,可能会把一些有效的弱信号稀疏掉。

还有一些高分辨率高聚焦的时频分析方法是建立在线性非线性时频分析方法基础上的,比如分数阶变换,谱重排,同步压缩变换。1980年Namias在量子力学中引入了分数傅里叶变换FRFT这一概念,1987年McBride和Kerr对其进行了完善,之后一些学者做了进一步的分析研究,并把这种分数阶变换推广到了Radon-Wigner、分数阶短时傅里叶变换STFRFT、分数阶Gabor变换FRGT,分数阶小波变换FRWT、分数阶S变换等。

傅里叶变换是将时间域信号转换到频率域,傅里叶反变换是把频率域转到时间域,而分数阶傅里叶变换可以将信号(时间域或者频率域)转换到时间和频率之间的域。由于分数阶傅里叶变换采用的是全局核,因此只能得到频谱信息,不能得到时间局部化信息,与短时傅里叶变换相似,分数阶短时傅里叶变换通过加窗分段处理,能在时间频率域联合表示信号,分数阶小波变换具有线性特性,既有小波变换多分辨率特性,又有类似于分数阶傅里叶变换分数域表征信号的能力,适用于频率域能量非最佳聚焦信号的处理。

时频谱重排是一种重新分配时频谱坐标的方法,通过瞬时频率和群延迟的局部估计,将数据映射到更接近信号真实支撑区间的时间-频率坐标。谱重排可以得到非常高的时频分辨率,在地震信号处理中也得到了应用,但是其最大的问题是建立在谱图的基础上的,把时间和频率两个方向的坐标都做了改变,不能重构信号。

同步压缩小波变换SST是一种新的时频能量排谱算法,与之前的谱重排方法不同,同步压缩小波变换是只对频率进行重排,可以重构原始信号,因此受到了广泛的欢迎。近年来,以同步压缩变换为核心发展了多种时频变换方法,包括同步压缩短时傅里叶变换和同步压缩S变换,同步压缩小波包变换等。随着对地震勘探精度要求的越来越高,这些高分辨率时频分析方法也在不同的地震处理问题上展现了自身的优势。同步压缩变换作为一种新发展起来的时频分析方法,将会在地球物理领域有更进一步的发展和应用。

曲波变换具有强大的多尺度分析和多方向分析的能力,在地震勘探领域得到了广泛的应用。可以利用曲波变换进行随机噪声和相干线性噪声衰减;可以利用自适应调整曲波阈值来压制随时间空间改变的非相干噪声;可以在曲波域进行稀疏反褶积去除随机噪声;可以在贝叶斯框架下利用曲波稀疏性压制面波;可以将曲波和奇异值分解相结合,对面波做滤波处理。众多研究结果表明,曲波变换面对地震数据的分析处理时具有非常大的优势。

本算法将同步压缩变换和曲波变换相结合,提出一种二维高分辨时频分析方法,并在MATLAB环境下使用二维高分辨时频分析方法提取波状分量(分离混合地震数据)。

部分代码如下:

is_real = 1;
sz = [N N]/8;
R_high = N/sqrt(2);
R_low = 0;
is_cos = 1;
is_unif = 1;
typeNUFFT = 1;
[xo yo] = ndgrid(0:1/N:(N-1)/N);
%ccc is the wave atom coefficients
ccc = gdct2_fwd(fff, is_real,is_unif,xo,yo,typeNUFFT,sz,R_high,R_low,rad,is_cos,t_sc,s_sc);
%aaa and bbb are the wave atom coefficients with the derivatives in b_1 and b_2
aaa = gdct2_ext_1(fff, is_real,is_unif,xo,yo,typeNUFFT,sz,R_high,R_low,rad,is_cos,t_sc,s_sc);
bbb = gdct2_ext_2(fff, is_real,is_unif,xo,yo,typeNUFFT,sz,R_high,R_low,rad,is_cos,t_sc,s_sc);

部分出图如下:

工学博士,担任《Mechanical System and Signal Processing》审稿专家,担任
《中国电机工程学报》优秀审稿专家,《控制与决策》,《系统工程与电子技术》,《电力系统保护与控制》,《宇航学报》等EI期刊审稿专家,担任《计算机科学》,《电子器件》 , 《现代制造过程》 ,《电源学报》,《船舶工程》 ,《轴承》 ,《工矿自动化》 ,《重庆理工大学学报》 ,《噪声与振动控制》 ,《机械传动》 ,《机械强度》 ,《机械科学与技术》 ,《机床与液压》,《声学技术》,《应用声学》,《石油机械》,《西安工业大学学报》等中文核心审稿专家。
擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。

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

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

相关文章

2.12日学习打卡----初学RocketMQ(三)

2.12日学习打卡 目录: 2.12日学习打卡一. RocketMQ高级特性(续)消息重试延迟消息消息查询 二.RocketMQ应用实战生产端发送同步消息发送异步消息单向发送消息顺序发送消息消费顺序消息全局顺序消息延迟消息事务消息消息查询 一. RocketMQ高级特…

Packet Tracer - Configuring ASA Basic Settings and Firewall Using CLI

Packet Tracer - 使用CLI配置ASA基本设置和防火墙 IP地址表 目标 验证连接并探索ASA设备使用CLI配置ASA的基本设置和接口安全级别使用CLI配置路由、地址转换和检查策略配置DHCP、AAA和SSH服务配置DMZ区域、静态NAT和访问控制列表(ACL) 场景 您的公司…

grafana配置钉钉告警模版(一)

1、配置钉钉告警模版 创建钉钉告警模版,然后在创建钉钉告警时调用模版。 定义发送内容具体代码 my_text_alert_list 是模版名称后面再配置钉钉告警时需要调用。 {{/* 定义消息体片段 */}} {{ define "my_text_alert_list" }}{{ range . }}告警名称&…

redis为什么使用跳跃表而不是树

Redis中支持五种数据类型中有序集合Sorted Set的底层数据结构使用的跳跃表,为何不使用其他的如平衡二叉树、b树等数据结构呢? 1,redis的设计目标、性能需求: redis是高性能的非关系型(NoSQL)内存键值数据…

在Postgresql 下安装QGIS

安装QGIS的前提是需要 在windows下安装Postgres,具体可以参考文章: Windows 安装和连接使用 PgSql数据库 安装GIS的具体步骤如下: 一.打开 Application Stack Builder 二.选择默认端口和安装目标 三.选择【Spatial Extensions】 四.选择安装…

链式结构实现队列

链式结构实现队列 1.队列1.1队列的概念及结构1.2队列的实现 2. 队列的各种函数实现3. 队列的全部代码实现 1.队列 1.1队列的概念及结构 队列:只允许在一端进行插入数据操作,在另一端进行删除数据操作的特殊线性表,队列具有先进先出 FIFO(Fi…

【大厂AI课学习笔记】【2.1 人工智能项目开发规划与目标】(5)数据管理

今天学习了数据管理,以及数据管理和数据治理的区别和联系。 数据管理:利用计算机硬件和软件技术对数据进行有效的收集、存储、处理和应用的过程其目的在于充分有效地发挥数据的作用。 实现数据有效管理的关键是数据组织。 数据管理和数据治理的区别&am…

无人驾驶控制算法LQR和MPC的仿真实现

1. LQR控制器 1.1 问题陈述 考虑一个质量为 m m m 的滑块在光滑的一维地面上运动。初始时,滑块的位置和速度均为 0 0 0。我们的目标是设计一个控制器,基于传感器测得的滑块位置 x x x,为滑块提供外力 u u u,使其能够跟随参考…

每日一题——LeetCode1455.检查单词是否为句中其他单词的前缀

方法一 js函数slice() 将字符串按空格符分割为单词数组&#xff0c;记searchWord的长度为n&#xff0c;分割每个单词的前n位看是否和searchWord匹配 var isPrefixOfWord function(sentence, searchWord) {let res sentence.split(" ")for(i 0 ; i < res.lengt…

七天入门大模型 :大模型LLM 训练理论和实战最强总结!

本文对于想入门大模型、面试大模型岗位、大模型实具有很强的指导意义。喜欢记得收藏、关注、点赞 文章目录 技术交流群用通俗易懂方式讲解系列总览介绍预训练范式如何确定自己的模型需要做什么训练&#xff1f;模型推理的一般过程PyTorch 框架设备PyTorch基本训练代码范例Trans…

【复现】cellinx摄像设备 未授权漏洞_50

目录 一.概述 二 .漏洞影响 三.漏洞复现 1. 漏洞一&#xff1a; 四.修复建议&#xff1a; 五. 搜索语法&#xff1a; 六.免责声明 一.概述 cellinx是一家韩国的摄像设备 二 .漏洞影响 通过未授权访问可以创建用户进入后台&#xff0c;可能造成系统功能破坏。 三.漏洞复…

CCF编程能力等级认证GESP—C++8级—20231209

CCF编程能力等级认证GESP—C8级—20231209 单选题&#xff08;每题 2 分&#xff0c;共 30 分&#xff09;判断题&#xff08;每题 2 分&#xff0c;共 20 分&#xff09;编程题 (每题 25 分&#xff0c;共 50 分)奖品分配大量的工作沟通 答案及解析单选题判断题编程题1编程题2…

GIS利用不舒适指数绘制地区的生物气候舒适度图

生物气候舒适度定义了最适宜的气候条件,在这种条件下,人们感到健康和充满活力。生物气候舒适度地图对城市规划研究特别有用。温度、相对湿度和风速等要素对评估生物气候舒适度非常重要。[1] 人们已经得出了许多不同的指数来确定生物气候舒适度。在本博文中,我们将使用广泛使…

基于SringBoot+Vue的大学生社团管理系统

末尾获取源码作者介绍&#xff1a;大家好&#xff0c;我是墨韵&#xff0c;本人4年开发经验&#xff0c;专注定制项目开发 更多项目&#xff1a;CSDN主页YAML墨韵 学如逆水行舟&#xff0c;不进则退。学习如赶路&#xff0c;不能慢一步。 目录 一、项目简介 1.1 研究背景 1.…

英文论文(sci)解读复现【NO.21】一种基于空间坐标的轻量级目标检测器无人机航空图像的自注意

此前出了目标检测算法改进专栏&#xff0c;但是对于应用于什么场景&#xff0c;需要什么改进方法对应与自己的应用场景有效果&#xff0c;并且多少改进点能发什么水平的文章&#xff0c;为解决大家的困惑&#xff0c;此系列文章旨在给大家解读发表高水平学术期刊中的 SCI论文&a…

leetcode hot100不同路径Ⅱ

本题和之前做的不同路径类似&#xff0c;区别是本题中加入了障碍&#xff0c;遇到障碍之后需要避开&#xff08;注意&#xff0c;这里依旧是只能向下向右移动&#xff09;&#xff0c;那么也就是说&#xff0c;有障碍的点是到达不了的&#xff0c;并且 &#xff0c;我在初始化的…

Java基于微信小程序的医院挂号小程序,附源码

博主介绍&#xff1a;✌程序员徐师兄、7年大厂程序员经历。全网粉丝12w、csdn博客专家、掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java技术领域和毕业项目实战✌ &#x1f345;文末获取源码联系&#x1f345; &#x1f447;&#x1f3fb; 精彩专栏推荐订阅&#x1f447;…

Atmel ATSHA204应用总结

1 ACES软件安装 Atmel Crypto Evaluation Studio (ACES) https://www.microchip.com/DevelopmentTools/ProductDetails/PartNO/Atmel%20Crypto%20%20Studio%20(ACES) 2 基本概念 ACES CE&#xff1a;Atmel Crypto Evalution Studio Configuration Environment&#xff08;基于加…

Intelij Terminal中文乱码解决

第一&#xff1a; &#xff08;重启Intelij生效&#xff09; -Dfile.encodingUTF-8 第二&#xff1a; &#xff08;重启Intelij生效&#xff09; 如果还不行&#xff0c;第三&#xff1a; 测试结果很ok&#xff1a;

红队打靶练习:IMF: 1

目录 信息收集 1、arp 2、nmap 3、nikto 目录探测 gobuster dirsearch WEB 信息收集 get flag1 get flag2 get flag3 SQL注入 漏洞探测 脱库 get flag4 文件上传 反弹shell 提权 get flag5 get flag6 信息收集 1、arp ┌──(root㉿ru)-[~/kali] └─# a…