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,一经查实,立即删除!

相关文章

高并发系统中常见的问题

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

外文文献查找的6个途径

寻找外文文献的渠道有很多,以下六个是一些常用的途径: 1、学术搜索引擎:像Google Scholar、PubMed、Web of Science、Scopus等学术搜索引擎是查找外文文献的常用工具,它们提供了广泛的学术资源和文献数据库。 2、学术数据库&…

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

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

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

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

【八股文】面向对象基础

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

rust gui fltk

FLTK 图形用户界面库的 Rust 绑定。 fltk crate 是一个跨平台的轻量级 GUI 库,可以静态链接以生成小型、独立且快速的 GUI 应用程序。 doc https://www.rust-lang.org/zh-CN/learn/get-started https://docs.rs/fltk/latest/fltk/ install $ curl --proto http…

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

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

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

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

Nginx返回502错误提示问题原因常用解决方法

Nginx返回502错误提示问题原因和常用解决方法 1. 因网站的访问量大,而php-cgi的进程数偏少,导致nginx 返回502错误 解决方法:针对这种情况的502错误,需增加php-cgi的进程数。 一般一个php-cgi进程占20M内存,你可以自…

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

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

力扣:45. 跳跃游戏 II

动态规划: 1.声明一个dp数组来表示到下标i时要跳转的最小次数,同时要个每个dp数组的值赋值为最大。之后初始化dp【0】0,用for循环来进行遍历下标i和遍历下标i之前的值,判断条件为如果i之前的某个下标跳转的距离>i时要根新dp【…

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

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

平台+低代码:中小企业数字化转型普惠之路

随着数字化转型的深入推进,中小企业面临着数字化转型的压力和挑战。如何在有限的资源和条件下,实现高效、便捷的数字化转型,成为中小企业亟待解决的问题。本文将以“平台低代码”为主题,探讨中小企业数字化转型的新模式&#xff0…

武汉灰京文化展望未来游戏产业,科技创新引领全面升级的游戏体验

随着科技的迅速发展,未来游戏产业的发展将迎来一个全新的纪元。科技创新将引领游戏体验的全面升级,让玩家不再仅仅是通过屏幕与游戏互动,而是能够亲身感受到游戏世界的存在。这种全新的游戏体验将推动游戏产业不断突破创新,吸引更…

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

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 (…

图形化编程:启蒙儿童的创意之窗

传统的编程语言对于初学者&#xff0c;尤其是儿童来说&#xff0c;可能会显得复杂和难以理解。为了解决这个问题&#xff0c;图形化编程应运而生。图形化编程是一种通过拖拽图形块来编写程序的方法&#xff0c;它让编程变得简单、直观和有趣&#xff0c;为儿童提供了一个启蒙创…

【Python笔记-设计模式】惰性评价模式

一、说明 将某些对象的创建或计算延迟到真正需要它们的时候&#xff0c;以减少不必要的资源消耗和提高性能。 惰性评价在Python中实现也成为生成器&#xff0c;一般通过yield关键字实现。 (一) 解决问题 在处理大量数据时&#xff0c;使用惰性加载可以避免一次性加载所有数…

spark为什么比mapreduce快?

spark为什么比mapreduce快&#xff1f; 首先澄清几个误区&#xff1a; 1&#xff1a;两者都是基于内存计算的&#xff0c;任何计算框架都肯定是基于内存的&#xff0c;所以网上说的spark是基于内存计算所以快&#xff0c;显然是错误的 2;DAG计算模型减少的是磁盘I/O次数&…

linux部署File Browser文件管理系统

File Browser是go语言编写的软件&#xff0c;可以通过浏览器对服务器上的文件进行管理 创建部署目录 /data/filebrowser mkdir -p /data/filebrowser && cd /data/filebrowser下载File Browser软件包 wget https://github.com/filebrowser/filebrowser/releases/do…