一维时间序列信号的小波时间散射变换(MATLAB 2021)

2024-06-06 02:52

本文主要是介绍一维时间序列信号的小波时间散射变换(MATLAB 2021),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

小波散射变换的目的在于获取第一层次的特征信息,即免疫平移、轻微形变的信息。而低通的滤波器能够获取输入信号的概貌,获取反映其整体大尺度特征的信息,以图像为例,由低通滤波器选取的信号对于图像的平移、伸缩、旋转等局部变化有良好的不变性。因此,小波散射变换获取了这一部分的信息。

相应的,高频部分的信号则对于信号的平移、伸缩、轻微形变类的变化具有协变性。换而言之,带通滤波器选取的信号成分受平移、伸缩、旋转等变化的影响较大,越是高频处的滤波选取越易受这些变化的影响。高频信息中包含着对提取共性特征不利的成分。但是,仅获取低频信息就意味着把包含重要信息成分的高频部分丢弃掉,这对于特征提取来说是一种损失。可以认为,散射变换最重要的目的就是在获取高频信息的同时,丢弃掉其对平移、形变敏感的部分。散射变化对此采取了如下的形式:

图片

总的来说,小波散射变换的形式可以归结为,对信号做尺度不同的小波变换,对其进行取模,继而对变换结果做一个低通滤波,所得即是所提取的特征。从数学上的理解来看,由于同一小波的正交特性,不同特征间应是不具有相关性,每一个各代表了图像一部分的特征信息,所有层特征的组合即可认为是散射网络从中获取的特征之和。不同于卷积神经网络,散射网络每一层的特征提取器是复数个而非单一特征提取器,且每一层用的都是同一套特征提取器。这与其他的多层架构网络也不同。同时,不同于卷积神经网络只用最后一层的输出,即将经过所有特征提取器的信息来作为最后的特征输出,散射网络在所有层都有输出,最后将这些输出的特征组合来作为总的输出特征。最后也是最显著的一点,散射网络不是前馈式网络,它的所有卷积核都是预先选定的,换句话说,在对图像进行特征提取之前,它所要提取的目标特征就是确定的。这样做的好处在于,省略了训练的过程,可以根据理论提取到具有目标性质的特征,计算量大大减少的同时也具有了结果的理论解释,这是一种介于传统图像识别和深度学习之间的算法。

鉴于此,采用小波时间散射变换对一维时间序列信号进行分析,运行环境为MATLAB 2021B。

图片

图片

图片

图片

图片

for nl = 1:numel(fb)phift = self.filters{fb(nl)}.phift;psift = self.filters{fb(nl)}.psift;% Obtain wavelet filterspositiveMag2Psift = psift.*conj(psift);negativeMag2Psift = circshift(flip(positiveMag2Psift),1);positiveMag2Psift = sum(positiveMag2Psift,2);negativeMag2Psift = sum(negativeMag2Psift,2);Mag2Phift = phift.*phift;lpsum(:,nl) = Mag2Phift+1/2*(positiveMag2Psift+negativeMag2Psift);
完整代码https://mbd.pub/o/bread/mbd-ZJyYmZxq
end

工学博士,担任《Mechanical System and Signal Processing》《中国电机工程学报》《控制与决策》等期刊审稿专家,擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。

这篇关于一维时间序列信号的小波时间散射变换(MATLAB 2021)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1034892

相关文章

服务器集群同步时间手记

1.时间服务器配置(必须root用户) (1)检查ntp是否安装 [root@node1 桌面]# rpm -qa|grep ntpntp-4.2.6p5-10.el6.centos.x86_64fontpackages-filesystem-1.41-1.1.el6.noarchntpdate-4.2.6p5-10.el6.centos.x86_64 (2)修改ntp配置文件 [r

uva 10131 最长子序列

题意: 给大象的体重和智商,求体重按从大到小,智商从高到低的最长子序列,并输出路径。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vect

POJ1631最长单调递增子序列

最长单调递增子序列 import java.io.BufferedReader;import java.io.InputStream;import java.io.InputStreamReader;import java.io.PrintWriter;import java.math.BigInteger;import java.util.StringTokenizer;publ

MiniGPT-3D, 首个高效的3D点云大语言模型,仅需一张RTX3090显卡,训练一天时间,已开源

项目主页:https://tangyuan96.github.io/minigpt_3d_project_page/ 代码:https://github.com/TangYuan96/MiniGPT-3D 论文:https://arxiv.org/pdf/2405.01413 MiniGPT-3D在多个任务上取得了SoTA,被ACM MM2024接收,只拥有47.8M的可训练参数,在一张RTX

leetcode105 从前序与中序遍历序列构造二叉树

根据一棵树的前序遍历与中序遍历构造二叉树。 注意: 你可以假设树中没有重复的元素。 例如,给出 前序遍历 preorder = [3,9,20,15,7]中序遍历 inorder = [9,3,15,20,7] 返回如下的二叉树: 3/ \9 20/ \15 7   class Solution {public TreeNode buildTree(int[] pr

批处理以当前时间为文件名创建文件

批处理以当前时间为文件名创建文件 批处理创建空文件 有时候,需要创建以当前时间命名的文件,手动输入当然可以,但是有更省心的方法吗? 假设我是 windows 操作系统,打开命令行。 输入以下命令试试: echo %date:~0,4%_%date:~5,2%_%date:~8,2%_%time:~0,2%_%time:~3,2%_%time:~6,2% 输出类似: 2019_06

matlab读取NC文件(含group)

matlab读取NC文件(含group): NC文件数据结构: 代码: % 打开 NetCDF 文件filename = 'your_file.nc'; % 替换为你的文件名% 使用 netcdf.open 函数打开文件ncid = netcdf.open(filename, 'NC_NOWRITE');% 查看文件中的组% 假设我们想读取名为 "group1" 的组groupName

利用matlab bar函数绘制较为复杂的柱状图,并在图中进行适当标注

示例代码和结果如下:小疑问:如何自动选择合适的坐标位置对柱状图的数值大小进行标注?😂 clear; close all;x = 1:3;aa=[28.6321521955954 26.2453660695847 21.69102348512086.93747104431360 6.25442246899816 3.342835958564245.51365061796319 4.87

【MRI基础】TR 和 TE 时间概念

重复时间 (TR) 磁共振成像 (MRI) 中的 TR(重复时间,repetition time)是施加于同一切片的连续脉冲序列之间的时间间隔。具体而言,TR 是施加一个 RF(射频)脉冲与施加下一个 RF 脉冲之间的持续时间。TR 以毫秒 (ms) 为单位,主要控制后续脉冲之前的纵向弛豫程度(T1 弛豫),使其成为显著影响 MRI 中的图像对比度和信号特性的重要参数。 回声时间 (TE)

C# double[] 和Matlab数组MWArray[]转换

C# double[] 转换成MWArray[], 直接赋值就行             MWNumericArray[] ma = new MWNumericArray[4];             double[] dT = new double[] { 0 };             double[] dT1 = new double[] { 0,2 };