当前位置: 开发笔记 > 编程语言 > 正文

matlab中morl小波,小波变换后的各频率分量的功率谱，

作者：mobiledu2502855463 | 来源：互联网 | 2023-09-24 11:55

该楼层疑似违规已被系统折叠隐藏此楼查看此楼下面程序是用来提取脑电波信号的，利用小波提取四种频率分量的波α、β、δ、θ，现在不明白如何求频率分量的功率谱&

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

下面程序是用来提取脑电波信号的&＃xff0c;利用小波提取四种频率分量的波α、β、δ、θ&＃xff0c;现在不明白如何求频率分量的功率谱&＃xff1f;

其中一段代码是&＃xff1a;

%dalt

fs1&＃61;fft(dalt,1024);%快速付氏变换

pp1&＃61;fs1.*conj(fs1)/1024;%计算功率谱

ff1&＃61;((0:511)/1024)*173.61;%计算各点对应的频率值

其中173.61是什么意思

完整代码如下&＃xff1a;

load singal.mat

s&＃61;sig;

fs&＃61;1024;%采样频率

dalt&＃61;1/fs;%采样间隔

t&＃61;0:length(s)-1*dalt;

figure(1);

subplot(311);plot(s);title(&＃39;原始信号&＃39;);

axis([-10 4000 -200 200]);

fs&＃61;fft(s,1024);%快速付氏变换

pp&＃61;fs.*conj(fs)/1024;%计算功率谱

ff&＃61;(0:511)/1024/dalt;%计算各点对应的频率值

subplot(312);plot(ff,pp(1:512));ylabel(&＃39;功率谱密度&＃39;);xlabel(&＃39;频率&＃39;);title(&＃39;信号功率谱&＃39;);axis([0 55 -10 110000]);

l&＃61;wmaxlev(s,&＃39;db4&＃39;);%取db4分解变量s的最大尺度l 实际分解尺度小于l

sd&＃61;wden(s,&＃39;minimaxi&＃39;,&＃39;s&＃39;,&＃39;mln&＃39;,real(l),&＃39;db4&＃39;);%除噪 dbN 对称 minimaxi用极大极小原理进行阈值选择 &＃39;s&＃39;为软阈值 &＃39;mln&＃39;在不同层估计噪声层并调整阈值

subplot(313);plot(sd);xlabel(&＃39;滤除噪声信号&＃39;);axis([-10 4000 -200 200]);

[LD,HD,LR,HR]&＃61;wfilters(&＃39;db4&＃39;);

figure(2);

subplot(221);stem(LD);title(&＃39;LD 低通分解滤波器&＃39;);grid;

subplot(222);stem(HD);title(&＃39;HD 低通重构滤波器&＃39;);grid;

subplot(223);stem(LR);title(&＃39;LR 高通分解滤波器&＃39;);grid;

subplot(224);stem(HR);title(&＃39;HR 高通重构滤波器&＃39;);grid;

[C,L]&＃61;wavedec(s,7,&＃39;db4&＃39;); %用小波db4对信号进行多尺度分解(7层)

C7&＃61;appcoef(C,L,&＃39;db4&＃39;,7);%尺度128

D7&＃61;detcoef(C,L,7);%细节系数

D6&＃61;detcoef(C,L,6);%尺度64

D5&＃61;detcoef(C,L,5);%尺度32

D4&＃61;detcoef(C,L,4);%尺度16

D3&＃61;detcoef(C,L,3);%尺度8

D2&＃61;detcoef(C,L,2);%尺度4

D1&＃61;detcoef(C,L,1);%尺度2

figure(3);

subplot(711);plot(D7);ylabel(&＃39;D7&＃39;);title(&＃39;细节系数&＃39;);

subplot(712);plot(D6);ylabel(&＃39;D6&＃39;);

subplot(713);plot(D5);ylabel(&＃39;D5&＃39;);

subplot(714);plot(D4);ylabel(&＃39;D4&＃39;);

subplot(715);plot(D3);ylabel(&＃39;D3&＃39;);

subplot(716);plot(D2);ylabel(&＃39;D2&＃39;);

subplot(717);plot(D1);ylabel(&＃39;D1&＃39;);

[C,L]&＃61;wavedec(s,6,&＃39;db4&＃39;);%近似系数尺度64

C6&＃61;appcoef(C,L,&＃39;db4&＃39;,6);

[C,L]&＃61;wavedec(s,5,&＃39;db4&＃39;);%尺度32

C5&＃61;appcoef(C,L,&＃39;db4&＃39;,5);

[C,L]&＃61;wavedec(s,4,&＃39;db4&＃39;);%尺度16

C4&＃61;appcoef(C,L,&＃39;db4&＃39;,4);

[C,L]&＃61;wavedec(s,3,&＃39;db4&＃39;);%尺度8

C3&＃61;appcoef(C,L,&＃39;db4&＃39;,3);

[C,L]&＃61;wavedec(s,2,&＃39;db4&＃39;);%尺度4

C2&＃61;appcoef(C,L,&＃39;db4&＃39;,2);

[C,L]&＃61;wavedec(s,1,&＃39;db4&＃39;);%尺度2

C1&＃61;appcoef(C,L,&＃39;db4&＃39;,1);

figure(4);

subplot(711);plot(C7);ylabel(&＃39;C7&＃39;);title(&＃39;近似系数&＃39;);

subplot(712);plot(C6);ylabel(&＃39;C6&＃39;);

subplot(713);plot(C5);ylabel(&＃39;C5&＃39;);

subplot(714);plot(C4);ylabel(&＃39;C4&＃39;);

subplot(715);plot(C3);ylabel(&＃39;C3&＃39;);

subplot(716);plot(C2);ylabel(&＃39;C2&＃39;);

subplot(717);plot(C1);ylabel(&＃39;C1&＃39;);

%带宽 0.5 Hz to 85 Hz;采样频率fs&＃61;173.61Hz

%δ-wave(1~4Hz);θ-wave(4~8Hz);α-wave(8~13Hz);β-wave(14~30Hz);

%******************************

[C,L]&＃61;wavedec(s,7,&＃39;db4&＃39;);%尺度128~64

C7&＃61;appcoef(C,L,&＃39;db4&＃39;,7);

D7&＃61;detcoef(C,L,7);

D6&＃61;detcoef(C,L,6);%尺度64

D5&＃61;detcoef(C,L,5);%尺度32

D4&＃61;detcoef(C,L,4);%尺度16

D3&＃61;detcoef(C,L,3);%尺度8

D2&＃61;detcoef(C,L,2);%尺度4

D1&＃61;detcoef(C,L,1);%尺度2

[C,L]&＃61;wavedec(s,7,&＃39;db4&＃39;);

SRD7&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,7);

SRD6&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,6);

SRD5&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,5);

SRD4&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,4);

SRD3&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,3);

SRD2&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,2);

SRD1&＃61;wrcoef(&＃39;d&＃39;,C,L,&＃39;db4&＃39;,1);

dalt&＃61;[SRD5];%δ-wave(1~4Hz)

thalt&＃61;[SRD4];%θ-wave(4~8Hz)

alpha&＃61;[SRD3];%α-wave(8~13Hz)

belta&＃61;[SRD2];%β-wave(14~30Hz)

%dalt

fs1&＃61;fft(dalt,1024);%快速付氏变换

pp1&＃61;fs1.*conj(fs1)/1024;%计算功率谱

ff1&＃61;((0:511)/1024)*173.61;%计算各点对应的频率值

%thalt

fs2&＃61;fft(thalt,1024);%快速付氏变换

pp2&＃61;fs2.*conj(fs2)/1024;%计算功率谱

ff2&＃61;((0:511)/1024)*173.61;%计算各点对应的频率值

%alpha

fs3&＃61;fft(alpha,1024);%快速付氏变换

pp3&＃61;fs3.*conj(fs3)/1024;%计算功率谱

ff3&＃61;((0:511)/1024)*173.61;%计算各点对应的频率值

%belta

fs4&＃61;fft(belta,1024);%快速付氏变换

pp4&＃61;fs4.*conj(fs4)/1024;%计算功率谱

ff4&＃61;((0:511)/1024)*173.61;%计算各点对应的频率值

figure(5);%特征波形提取

subplot(411);plot(ff1,pp1(1:512));ylabel(&＃39;功率谱密度&＃39;);title(&＃39;dalt信号功率谱&＃39;);%axis([0 30 -10 110000]);

subplot(412);plot(ff2,pp2(1:512));ylabel(&＃39;功率谱密度&＃39;);title(&＃39;thalt信号功率谱&＃39;);%axis([0 55 -10 110000]);

subplot(413);plot(ff3,pp3(1:512));ylabel(&＃39;功率谱密度&＃39;);title(&＃39;alpha信号功率谱&＃39;);%axis([0 55 -10 110000]);

subplot(414);plot(ff4,pp4(1:512));ylabel(&＃39;功率谱密度&＃39;);xlabel(&＃39;频率&＃39;);title(&＃39;belta信号功率谱&＃39;);%axis([0 55 -10 110000]);

figure(6);

subplot(311);xlabel(&＃39;时间 t/s&＃39;);ylabel(&＃39;频率 f/Hz&＃39;);zlabel(&＃39;幅值 A/uV&＃39;);

c &＃61; cwt(s,[64 32 16:-2:2],&＃39;morl&＃39;,&＃39;3Dlvl&＃39;); %morl只能进行连续小波变换 morl具有对称性&＃xff0c;不生成正交系&＃xff0c;不具有有限支集等性质

subplot(312);xlabel(&＃39;时间 t/s&＃39;);ylabel(&＃39;频率 f/Hz&＃39;);zlabel(&＃39;幅值 A/uV&＃39;);

c &＃61; cwt(s,[16 8 4 2 1],&＃39;morl&＃39;,&＃39;3Dlvl&＃39;); %morl进行连续小波变换 &＃39;lvl&＃39; By scale不同尺度对比

subplot(313);xlabel(&＃39;时间 t/s&＃39;);ylabel(&＃39;频率 f/Hz&＃39;);zlabel(&＃39;幅值 A/uV&＃39;);

c &＃61; cwt(s,[16 8 4 2 1],&＃39;mexh&＃39;,&＃39;3Dlvl&＃39;);%用mexihat进行连续小波变换不同小波对比

推荐阅读

7层
Vue 实现表格分页功能详解

本文将详细介绍如何在 Vue 中实现表格的分页功能，包括代码示例和具体实现步骤，帮助开发者更好地理解和应用这一技术。 ... [详细]

蜡笔小新 2024-11-14 16:00:58
case
Matlab学以致用模拟os任务装载情况

2012-06-0821:26:42 用matlab来建模，仿真不同时刻ostask在队列中的装载情况。输入参数如下作为初学者，M文件写的有点长。能实现功能就算学以致用了。cle ... [详细]

蜡笔小新 2024-11-12 19:40:34
format
使用多项式拟合分析淘宝双11销售趋势

根据天猫官方数据，2019年双11成交额达到2684亿元，再次刷新历史记录。本文通过多项式拟合方法，分析并预测未来几年的销售趋势。 ... [详细]

蜡笔小新 2024-11-12 16:57:26
main
JavaAwtSwing笔记之 JComponent JLabel 的 setBorder(Border 边框) 用例

JComponentJLabel的setBorder前言用例2205262241前言setBorder(Border边框)实现自JComponentjava.awt.Insets ... [详细]

蜡笔小新 2024-11-12 13:48:19
select
在AX2012中使用自定义查询在数据网格视图中显示数据

本文介绍了如何在AX2012中通过自定义查询在数据网格视图中显示所有记录的方法。 ... [详细]

蜡笔小新 2024-11-12 12:02:50
go
如何在Webpack项目中集成ECharts

本文将详细介绍如何在Webpack项目中安装和使用ECharts，包括全量引入和按需引入的方法，并提供一个柱状图的示例。 ... [详细]

蜡笔小新 2024-11-12 09:49:07
format
Python 序列图分割与可视化编程入门教程

本文介绍了如何使用 Python 进行序列图的快速分割与可视化。通过一个实际案例，详细展示了从需求分析到代码实现的全过程。具体包括如何读取序列图数据、应用分割算法以及利用可视化库生成直观的图表，帮助非编程背景的用户也能轻松上手。 ... [详细]

蜡笔小新 2024-11-11 07:14:26
main
Java项目中两个文件互相调用时出现函数错误，请求专业解决方案

在Java项目中，当两个文件进行互相调用时出现了函数错误。具体问题出现在 `MainFrame.java` 文件中，该文件位于 `cn.javass.bookmgr` 包下，并且导入了 `java.awt.BorderLayout` 和 `java.awt.Event` 等相关类。为了确保项目的正常运行，请求提供专业的解决方案，以解决函数调用中的错误。建议从类路径、依赖关系和方法签名等方面入手，进行全面排查和调试。 ... [详细]

蜡笔小新 2024-11-08 17:04:11
js
木偶评估函数拒绝将JSHandles和ElementHandles用作参数

在处理木偶评估函数时，我发现可以顺利传递本机对象（如字符串、列表和数字），但每当尝试将JSHandle或ElementHandle作为参数传递时，函数会拒绝接受这些对象。这可能是由于这些句柄对象的特殊性质导致的，建议在使用时进行适当的转换或封装，以确保函数能够正确处理。 ... [详细]

蜡笔小新 2024-11-07 18:34:55
js
利用 Python 中的 Altair 库实现数据抖动的水平剥离分析

利用 Python 中的 Altair 库实现数据抖动的水平剥离分析 ... [详细]

蜡笔小新 2024-11-02 12:39:54
js
基于POI的Java Maven项目中实现Excel数据的高效导入与导出功能

本项目在Java Maven框架下，利用POI库实现了Excel数据的高效导入与导出功能。通过优化数据处理流程，提升了数据操作的性能和稳定性。项目已发布至GitHub，当前最新版本为0.0.5。该项目不仅适用于小型应用，也可扩展用于大型企业级系统，提供了灵活的数据管理解决方案。GitHub地址：https://github.com/83945105/holygrail，Maven坐标：`com.github.83945105:holygrail:0.0.5`。 ... [详细]

蜡笔小新 2024-10-31 18:52:27
main
Java理论基础题库精编_必备|Java核心知识点练习题库

本题库精选了Java核心知识点的练习题，旨在帮助学习者巩固和检验对Java理论基础的掌握。其中，选择题部分涵盖了访问控制权限等关键概念，例如，Java语言中仅允许子类或同一包内的类访问的访问权限为protected。此外，题库还包括其他重要知识点，如异常处理、多线程、集合框架等，全面覆盖Java编程的核心内容。 ... [详细]

蜡笔小新 2024-10-31 16:12:02
format
C#编程指南：实现列表与WPF数据网格的高效绑定方法

C#编程指南：实现列表与WPF数据网格的高效绑定方法 ... [详细]

蜡笔小新 2024-10-31 10:46:47
testing
软件测试工程师必备：全面解析Web应用的浏览器兼容性测试方法与工具

随着市场上浏览器种类的不断增多，特别是在平板和移动设备上的普及，确保网站在各种浏览器中都能正常运行变得尤为重要。本文将深入探讨Web应用的浏览器兼容性测试方法与工具，帮助软件测试工程师全面掌握这一关键技能，以提高应用的用户体验和可靠性。 ... [详细]

蜡笔小新 2024-10-31 07:56:54
main
尽管存在唯一列，仍显示“当前选择不包含唯一列。网格编辑、复选框、编辑、复制和删除功能不可用”的消息。

尽管存在唯一列，仍显示“当前选择不包含唯一列。网格编辑、复选框、编辑、复制和删除功能不可用”的消息。 ... [详细]

蜡笔小新 2024-10-30 18:52:18

mobiledu2502855463

这个家伙很懒，什么也没留下！

Tags | 热门标签

RankList | 热门文章