【经验模态分解】4.信号由时域向频域的转换

news2024/9/20 22:29:51
/**
 * @poject          经验模态分解及其衍生算法的研究及其在语音信号处理中的应用
 * @file            傅里叶变换与小波变换
 * @author			jUicE_g2R(qq:3406291309)
 * 
 * @language        MATLAB
 * @EDA				Base on matlabR2022b
 * @editor			Obsidian(黑曜石笔记软件)
 * 
 * @copyright		2023
 * @COPYRIGHT	    原创学习笔记:转载需获得博主本人同意,且需标明转载源
 */
  • E M D EMD EMD 是 基于 傅里叶变换 与 小波变换 的改进
  • E M D EMD EMD离散小波变换 产生的背景是 基于解决 傅里叶变换 与 小波变换时间尺度 上存在的缺陷

文章目录

  • 1 傅里叶变换与短时傅里叶变换
    • 1-1 (传统的)傅里叶变换
      • 1-1-1 优点
      • 1-1-2 缺陷
    • 1-2 短时傅里叶变换
    • 1-3 信号的频谱图和时频图实验
  • 2 离散小波变换
    • 2-1 尺度因子
    • 2-2 平移因子

1 傅里叶变换与短时傅里叶变换

1-1 (传统的)傅里叶变换

  • 傅里叶变换(Fourier Transform)是一种将 时域上的信号转换到频域上 的数学工具

1-1-1 优点

优点详细说明
可时频转换傅里叶变换可以将时域上的信号转换为频域上的信号,便于对信号的频率特征进行分析和处理
可分解傅里叶变换可以将信号分解成多个不同频率的正弦和余弦函数,这些分量可以提供有关信号频率和振幅的信息
线性性傅里叶变换具有线性性,可以对多个信号进行叠加处理
可逆性可以将频率域上的信号转换回时域上的信号

1-1-2 缺陷

缺陷详细说明
实时性差傅里叶变换是一种理论分析工具,不适合对信号进行实时处理,因此需要使用快速傅里叶变换(FFT)等算法来加速计算
普适性差傅里叶变换需要将整个信号在时间上进行全局分析,计算复杂度较高,不适合对长时间信号进行处理
局部性差傅里叶变换对信号的局部特征无法进行分析,对于非周期性信号和突发性信号的频谱分析效果较差
处理特殊情况效果差傅里叶变换对于非线性信号、带噪声信号等情况处理效果较差
  • 傅里叶变换 是 对 一段信号 不分时间先后 来观察其中的频率成分 及 各个频率成分的贡献,即没有时间分辨率不可通过傅里叶变换知道一个特定时间点的频率成分

1-2 短时傅里叶变换

  • 短时傅里叶变换 是 基于 (传统的)傅里叶变换时间分辨率的缺陷 给出了相应的改进方法

短时:即将 长段信号源 拆分成 多段 短的信号源(但是 窗口大小是固定的【固定窗长】,这也是它的缺陷)

Ts = 0:0.001:2; 						% 采样频率1000次每秒,时长为2s
%% 两个源信号用的是扫频函数chirp()
Sig_Low2High = chirp(Ts, 50, 2, 300); 	% 低频到高频:从 50Hz  线性升至 300Hz
Sig_High2Low = chirp(Ts, 300, 2, 50);	% 高频到地频:从 300Hz 线性升至 50Hz

1-3 信号的频谱图和时频图实验

步骤:

1-3-1 打开 信号分析器APP

1-3-2 面板设置

在这里插入图片描述

1-3-3 拖入信号

在这里插入图片描述

1-3-4 设置 采样率 S a m p l e R a t e SampleRate SampleRate

在这里插入图片描述
注:两个信号要分别修改对应的 S a m p l e R a t e SampleRate SampleRate

1-3-5 对每个面板选中相应的信号

在这里插入图片描述

1-3-6 选中面板并点击上方工具栏的频谱获得相应的频谱图

在这里插入图片描述

  • 分析上述频谱图
    左面板是 高频 到 低频,右面板是 低频 到 高频:
    两边的频谱图是一模一样的,一个原因是两个信号都是扫频函数(频率线性变化),另一个是频谱工具箱得到的频谱图是通过 傅里叶变换 得到的(傅里叶变换是全局的,两个全局确实一样,但局部不一样)

  • 选中面板并点击上方工具栏的时谱获得相应的时谱图
    在这里插入图片描述
    时频图是根据 短时傅里叶变换 得到的,很明显的可以看到两个信号的区别,同时体现了 短时傅里叶变换 具有 时间与频率的分辨率

1-3-7 结论

傅里叶变换 相较于 短时傅里叶变换 存在明显的 时频分辨率 上的差异

2 离散小波变换

  • 频域 上来看,离散小波变换 相当于一个 带通滤波器,相当于通过 尺度因子 与 平移因子
因子作用
尺度因子控制 小波基 的 伸展程度(控制 频率
平移因子控制 时间
  • 短时傅里叶变换 而言,离散小波变换 具有平衡 频率分辨率 与 时间分辨率 的优势

2-1 尺度因子

  • 小波母函数
    C e n t F r e q ( 中心频率 ) = S c a l e F a c t o r ( 尺度因子 ) ∗ f CentFreq(中心频率)=ScaleFactor(尺度因子)*f CentFreq(中心频率)=ScaleFactor(尺度因子)f

2-2 平移因子

  • 能使得 时间和频率之间 达到 平衡
    离散小波变换 是 对 短时傅里叶变换固定窗长 这一特性出现 “低频精度不高,高频分辨率太低” 缺陷的一种解决办法。

  • C ( 常数 ) = W d n L e n ( 窗长 ) ∗ f C(常数)=WdnLen(窗长)*f C(常数)=WdnLen(窗长)f
    高频时:具有很高的 时间分辨率,窗长小
    低频时:具有很高的 频率精度
    ,窗长大

本节参考视频源:NO.11 两步搞定经验模态分解与离散小波变换-哔哩哔哩

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

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

相关文章

Git的进阶操作,在idea中部署gie

🏅我是默,一个在CSDN分享笔记的博主。📚📚 ​​ 🌟在这里,我要推荐给大家我的专栏《git》。🎯🎯 🚀无论你是编程小白,还是有一定基础的程序员,这…

CCIA数安委等组织发布PIA星级标识名单,合合信息再次通过数据安全领域权威评估

近期,“中国网络安全产业联盟(CCIA)数据安全工作委员会”、“数据安全共同体计划(DSC)”等组织共同发起“个人信息保护影响评估专题工作(简称“PIA专题工作”)”,并为入围企业颁发了…

交通信号标志识别系统 python 深度学习 YOLOv5

[毕业设计]2023-2024年最新最全计算机专业毕设选题推荐汇总 感兴趣的可以先收藏起来,还有大家在毕设选题,项目以及论文编写等相关问题都可以给我留言咨询,希望帮助更多的人 。 1、项目介绍 本系统基于YOLOv5,采用登录注册进行用…

查找或替换excel换行符ctrl+j和word中的换行符^p,^l

一、excel中 直接上图。使用ctrlh调出替换,查找内容里按ctrlj(会出现一个闪的小点),即为换行符。 二、word中 在word中,^p和^l分别代表换行符(enter)和手动换行符(使用shiftenter&…

Spring高手之路16——解析Spring XML配置的BeanDefinition源码

文章目录 1. BeanDefinition阶段的分析2. 加载xml配置文件2.1 XML配置文件中加载bean的代码示例2.2 setConfigLocations - 设置和保存配置文件路径2.3 refresh - 触发容器刷新,配置文件的加载与解析2.4 loadBeanDefinitions - 具体的BeanDefinition加载逻辑2.5 load…

java数据结构--阻塞队列

目录 一.概念 二.生产者消费者问题 三.阻塞队列接口BlockingQueue 四.基于数组实现单锁的阻塞队列 1.加锁方式 2.代码实现 3.解释说明 (1).offer添加元素 (2)poll取出元素 4.timeout超时时间 5.测试 五.基于数组实现双锁的阻塞队列 1.问题 …

数据的使用、表关系的创建、Django框架的请求生命周期流程图

目录 一、数据的增删改查 1. 用户列表的展示 2. 修改数据的逻辑分析 3. 删除功能的分析 二、如何创建表关系 三、Django的请求生命周期流程图 一、数据的增删改查 1. 用户列表的展示 把数据表中得用户数据都给查询出来展示在页面上 查询数据 def userlist(request):&qu…

luckysheet的使用——14.开启表格只读模式(所有单元格无法编辑)

开启只读模式后,所有的单元格都无法编辑,与非编辑模式做区分。 1.在src/global/api.js文件中,新增开启只读模式的方法: /*** 开启工作表只读模式(所有单元格无法编辑)*/ export function setWorkBookReadOnly() {Store.allowEdi…

Docker部署ubuntu1804镜像详细步骤

Docker部署ubuntu1804镜像详细步骤 ubuntu镜像库地址:https://hub.docker.com/_/ubuntu/tags?page1&ordering-name 拉取镜像(默认为最新版本): docker pull ubuntu或,拉取指定版本镜像: docker pull…

P1529 [USACO2.4] 回家 Bessie Come Home 题解

文章目录 题目描述输入格式输出格式样例样例输入样例输出 提示完整代码 题目描述 现在是晚餐时间,而母牛们在外面分散的牧场中。 Farmer John 按响了电铃,所以她们开始向谷仓走去。 你的工作是要指出哪只母牛会最先到达谷仓(在给出的测试数…

【文件IO】认识文件

文章目录 认识文件文件的结构和目录文件路径 认识文件 我们先来认识狭义上的文件(file),针对硬盘这种持久化存储的I/O设备,当我们想要进行数据保存时,往往不是保存一个整体,而是独立成一个个单位进行保存,这个独立的单…

人工智能技术的高速发展,普通人如何借助AI实现弯道超车?

人工智能技术的高速发展,普通人如何借助AI实现弯道超车? 随着互联网信息传播的爆炸,人类科技文明的快速发展“人工智能”成为新的话题,科技的进步也让普通人觉得自己与社会脱节,找工作越来越难,创业越来越难…

前端开发学习指南

前端是一个看似入门门槛不高,但要学好很难的领域。前端的知识体系庞杂又松散,技术演进快,如果摸不清脉络的话很容易陷入盲人摸象的困境甚至跑偏。 其实只要掌握了正确的方法,学习前端和学好前端就只是个时间问题,希望下…

阻塞队列和定时器的使用

阻塞队列 谈到队列,大家就能想到队列的先进先出原则,但有些特殊的队列,虽然也是先进先出的,但是带有阻塞功能,我们把这种队列叫做阻塞队列. ★如果队列为空,执行出队操作就会阻塞,阻塞到另外一个线程往队列里添加元素(队列不为空)为止. ★如果队列满了,执行入队操作时,也会阻…

网工内推 | 运维工程师,软考认证优先,全额社保

01 北京中科网威信息技术有限公司 招聘岗位:运维工程师 职责描述: 1 熟悉网络安全标准,等级保护管理制度 2 负责等级保护管理制度的的企业管理要求编写; 3 熟系网络组网和相关安全产品; 4 负责用户需求挖掘、分析和…

关于锁策略

常见的锁策略悲观锁乐观锁读写锁轻量级锁、重量级锁自旋锁公平锁和非公平锁可重入锁 vs 不可重入锁synchronized是什么锁呢? 常见的锁策略 锁策略不仅仅限制于Java;其它锁相关的也是会涉及这些策略;这些特性主要是在实现锁的时候运用的。虽然我们的工作可能就是把轮…

中断处理程序的延迟可能导致中断标志位仍然被置位

当中断处理程序的执行时间超过了中断事件的频率时,可能出现中断标志位仍然被置位的情况。让我们来详细解释一下这种情况。 在一个典型的系统中,中断处理程序会在中断事件发生时被触发执行。中断处理程序负责处理中断事件,并可能执行一系列操…

sjvisualizer,一个超强的Python数据可视化动画库

大家好,今天给大家介绍一个非常棒的数据可视化库,sjvisualizer。 根据时间序列数据制作动态图表,包含条形图、饼图、堆叠条形图、折线图、堆叠面积图。 可以先看一下官方的示例~ 只需几行代码,就可以制作电脑浏览器发展史的动态…

使用jmeter进行简单压力测试

前言 最近项目要上线,需要项目进行简单的压力测试,本次使用的是jmeter来进行的,由于本人不是专业测试,只是对本次使用过程进行简单的记录. 一.jemeter的下载与安装 我这个已经安装很久了,具体过程这个可以查询下其他博客(偷个懒). 二.使用过程 1.测试计划右击-添加(add)-线…

Kubernetes 创建pod的yaml文件-简单版-nginx

apiVersion: v1 #api文档版本 kind: Pod # 资源类型 Deployment,StatefulSet之类 metadata: #pod元数据 描述信息 name: nginx-demo labels: type: app #自定义标签 version: 1.0.0 # 自定义pod版本 namespace: default spec: #期望Pod按照这里的描述创建 cont…