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

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.mzph.cn/news/691906.shtml

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

相关文章

高并发系统中常见的问题

在当今的高并发系统中,常见的问题是多种多样的,这些问题往往会对系统的稳定性和性能产生重大影响。本文将详细介绍高并发系统中常见的问题,并探讨其产生原因和解决方案。 一、高并发系统概述 高并发系统是指在同一时间内有大量用户同时访问…

对视频进行分块,断点续传

分块测试 //分块测试Testpublic void testChunk() throws IOException {//源路径File sourceFile new File("D:\\BaiduNetdiskDownload\\Day1-00.项目导学.mp4");//分块文件存储路径String chunkFilePath "D:\\develop\\chunk\\";//分块文件大小int chun…

解锁跨境电商新纪元:大数据驱动下的精准营销攻略

随着全球互联网的快速发展,跨境电商已经成为企业拓展国际市场的重要途径。在这个竞争激烈的环境中,如何有效地进行精准营销成为企业成功的关键之一。大数据技术的兴起为跨境电商提供了丰富的信息和洞察,为精准营销提供了坚实的基础。本文Nox聚…

【八股文】面向对象基础

【八股文】面向对象基础 面向对象和面向过程的区别 面向过程把解决问题的过程拆成一个个方法,通过一个个方法的执行解决问题。面向对象会先抽象出对象,然后用对象执行方法的方式解决问题。 创建一个对象用什么运算符?对象实体与对象引用有何不同? …

安达发|APS排产软件的机台产线任务甘特图功能详解

在现代制造业中,高级计划与排产是制造业运营的关键环节。为了提高生产效率、降低成本并确保产品质量,企业需要对生产过程进行精细化管理。APS(高级计划与排产)系统作为一种先进的生产计划和调度工具,可以帮助企业实现这…

过了30岁了,一定要专注一件事情?视频号值得尝试!

经常说视频号下载助手, 但发现大多数的大佬都只是先专注一件事情。 小编初6就回来了,和一个大佬吃饭,虽然人家规模并不大,但日引客户上千也是基本的。 这里给大家揭秘一下,他的做法!!&#x…

3ds Max视频怎么渲染 3ds Max云渲染操作

在3ds Max软件中制作视频动画本质上是逐帧生成画面,并将这些连续帧串联起来创造出动态连贯的视觉效果。常见的视频帧率包括25 FPS(每秒帧数)、60 FPS、以及120 FPS等,帧率的提升可以使视频动画更加流畅。在实质上,视频渲染就是动画渲染&#…

MakerPlace by Michaels:与Etsy一争高下的手工艺品新平台

Michaels是北美规模很大的专业艺术与手工艺品零售商,在加拿大和美国有1290家店。作为一个新的手工艺品在线交易平台——MakerPlace by Michaels 充满无限机遇!所以各位跨境人们也不要错过这个平台,跟着我一起来了解一下吧! MakerP…

探究二维码技术:连接现实与数字世界的桥梁

title: 探究二维码技术:连接现实与数字世界的桥梁 date: 2024/2/19 13:15:36 updated: 2024/2/19 13:15:36 tags: 二维码技术数据编码纠错算法图像处理商业应用安全验证实时交互 引言: 二维码已经成为现代社会中广泛应用的一种技术工具。它不仅在商业领…

牛客网 OR141 密码检查

答案&#xff1a; #include <stdio.h> #include <string.h> #include <ctype.h> int main() {int n 0;int count1 0, count2 0, count3 0;scanf("%d", &n);while (n--){char ch[100];scanf("%s", ch);int len strlen(ch);if (…

基于SSM的在线教学质量评价系统(有报告)。Javaee项目。ssm项目。

演示视频&#xff1a; 基于SSM的在线教学质量评价系统&#xff08;有报告&#xff09;。Javaee项目。ssm项目。 项目介绍&#xff1a; 采用M&#xff08;model&#xff09;V&#xff08;view&#xff09;C&#xff08;controller&#xff09;三层体系结构&#xff0c;通过Spri…

远程办公利器,有哪些远程办公利器好用

远程办公利器&#xff0c;有哪些远程办公利器好用 有的时候需要紧急办公位&#xff0c;但是你的电脑、公司、办公室、酒店、网吧、家里的电脑等极为不方便&#xff0c;这时有远程办公利器就好了。在网上查了下&#xff0c;发下个好用的工具&#xff0c;这里分享下。 1、…

API关键词接口的应用范围说明

一、应用场景 适用商家批发选品&#xff0c;价格监控&#xff0c;品牌监控&#xff0c;实地商家进货 二、公共参数 请求地址: https://api-gw.onebound.cn/1688/item_search 三、请求参数 请求参数&#xff1a;q女装&start_price0&end_price0&page1&cat0&…

UE5 C++ UObject实例化

一.创建UObject C类 在MyObject中声明结构体FMyDataTableStruct 在MyPawn里面&#xff0c;先将头文件里包含 MyObject.h 在MyPawn中声明一个UMyObject类型的指针 TSubclassOf 是提供 UClass 类型安全性的模板类。例如您在创建一个投射物类&#xff0c;允许设计者指定伤害类型…

openai chatGPT 原理通俗介绍

引言 近年来&#xff0c;随着深度学习技术的不断发展&#xff0c;自然语言处理&#xff08;NLP&#xff09;领域取得了长足的进步。ChatGPT&#xff08;Generative Pre-trained Transformer&#xff09;作为一种先进的语言生成模型&#xff0c;在各类对话系统和智能助手中得到…

Stable Diffusion WebUI 常用命令行参数

本文收录于《AI绘画从入门到精通》专栏&#xff0c;专栏总目录&#xff1a;点这里。 大家好&#xff0c;我是水滴~~ 本文主要讲解 Stable Diffusion WebUI 的一些常用命令行参数&#xff0c;内容详细的介绍了每一种参数的使用&#xff0c;并配有截图&#xff0c;非常适合初学者…

[规章规划]国家数据局-《关于开展全国数据资源调查的通知》要点解析

去年10月份国家数据局成立以后&#xff0c;各地方省级政府数据局纷纷快速挂牌成立。 仅以2024年1月举例&#xff0c;1月5日&#xff0c;江苏省数据局&#xff08;江苏省政务服务管理办公室&#xff09;挂牌成立&#xff0c;成为国家数据局正式揭牌后&#xff0c;第一个成立的省…

Web3区块链游戏:创造虚拟世界的全新体验

随着区块链技术的不断发展&#xff0c;Web3区块链游戏正逐渐崭露头角&#xff0c;为玩家带来了全新的虚拟世界体验。传统游戏中的中心化结构和封闭经济体系已经被打破&#xff0c;取而代之的是去中心化的游戏环境和真实所有权的数字资产。本文将深入探讨Web3区块链游戏的特点、…

回答一下 Spring Bean 生命周期的执行流程?

一、问题解析 Spring 生命周期全过程大致分为五个阶段&#xff1a;创建前准备阶段、创建实例阶段、依赖注入阶段、容器缓存阶段和销毁实例阶段。 这张是 Spring Bean 生命周期完整流程图&#xff0c;其中对每个阶段的具体操作做了详细介绍&#xff1a; &#xff08;把这张图按…

用冒泡排序实现快速排序(qsort函数),指针进阶实例

目录 1、qsort函数是什么 2、冒泡排序实现指针进阶 2.1 主函数 2.2 功能函数声明​编辑 2.3 my_qsort函数介绍 2.4 Swap函数 总结 1、qsort函数是什么 qsort函数是c语言自带的函数&#xff0c;其功能是实现快速排序。我们来看一下他的参数和返回值&#xff1a; 以上就是q…