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

matlabmag函数,频谱分析函数

作者：嘉苹民智79 | 来源：互联网 | 2023-05-28 12:14

大神你懂这个代码吗，我从网上找的，我想请你解释一下。。。下面程序是用来提取脑电波信号的，利用小波提取四种频率分量的波α、β、δ、θ

大神你懂这个代码吗&＃xff0c;我从网上找的&＃xff0c;我想请你解释一下。。。

下面程序是用来提取脑电波信号的&＃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进行连续小波变换不同小波对比

推荐阅读

java
如何在datetimebox中进行赋值与取值操作

在 datetimebox 中进行赋值和取值操作时，可以通过以下方法实现：使用 `$('#j_dateStart').datebox('setValue', '指定日期')` 进行赋值，而通过 `$('#j_dateStart').datebox('getValue')` 获取当前选中的日期值。若需要清空日期值，可以使用 `$('#j_dateStart').datebox('clear')` 方法。这些操作能够确保日期控件的准确性和灵活性，适用于各种前端应用场景。 ... [详细]

蜡笔小新 2024-10-27 10:48:27
case
Matlab学以致用模拟os任务装载情况

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

蜡笔小新 2024-11-12 19:40:34
php
使用多项式拟合分析淘宝双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
require
如何在Webpack项目中集成ECharts

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

蜡笔小新 2024-11-12 09:49:07
main
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
java
木偶评估函数拒绝将JSHandles和ElementHandles用作参数

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

蜡笔小新 2024-11-07 18:34:55
java
利用Python实现高效语音识别技术

本文探讨了利用Python实现高效语音识别技术的方法。通过使用先进的语音处理库和算法，本文详细介绍了如何构建一个准确且高效的语音识别系统。提供的代码示例和实验结果展示了该方法在实际应用中的优越性能。相关文件可从以下链接下载：链接：https://pan.baidu.com/s/1RWNVHuXMQleOrEi5vig_bQ，提取码：p57s。 ... [详细]

蜡笔小新 2024-11-07 13:05:53
config
利用 Python 中的 Altair 库实现数据抖动的水平剥离分析

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

蜡笔小新 2024-11-02 12:39:54
php
宏基因组学经典文献重现（一）：利用ggplot2进行散点图可视化分析

宏基因组学经典文献重现（一）：利用ggplot2进行散点图可视化分析 ... [详细]

蜡笔小新 2024-10-30 16:50:53
replace
在单个图表中实现饼图与条形图的精准对齐

在单个图表中实现饼图与条形图的精准对齐 ... [详细]

蜡笔小新 2024-10-29 13:42:17
select
深入解析 C# 中 SqlCommand 与 SqlDataAdapter 的核心差异与应用场景

本文深入探讨了 C# 中 `SqlCommand` 和 `SqlDataAdapter` 的核心差异及其应用场景。`SqlCommand` 主要用于执行单一的 SQL 命令，并通过 `DataReader` 获取结果，具有较高的执行效率，但灵活性较低。相比之下，`SqlDataAdapter` 则适用于复杂的数据操作，通过 `DataSet` 提供了更多的数据处理功能，如数据填充、更新和批量操作，更适合需要频繁数据交互的场景。 ... [详细]

蜡笔小新 2024-10-28 12:25:41
main
在Java应用中实现只读模式的切换方法与技巧

在Java应用中实现只读模式的切换方法与技巧 ... [详细]

蜡笔小新 2024-10-27 19:40:35

嘉苹民智79

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

Tags | 热门标签

RankList | 热门文章