小波包获得某个节点信号的几个细节问题-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 考研考博>>

考研

>>

小波包获得某个节点信号的几个细节问题

小波包获得某个节点信号的几个细节问题

发布:jingweimo | 分类:考研

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

CDA数据分析师证书

谁适合考CDA证书?年龄18-40周岁,专业是计算机、工商管理、统计学、管理科学类,学历本科及以上,行业是金融、信息技术、电信等,岗位是数据、产品、运营等。

谁适合考CDA证书?年龄18-40周岁,专业是计算机、工商管理、统计学、管理科学类,学历本科及以上,行业是金融、信息技术、电信等,岗位是数据、产品、运营等。

http://cn.qzs.qq.com/ac/b.gifhttp://cn.qzs.qq.com/ac/b.gif小波包分解主要用了这四个函数wpdec,wprec,wprcoef,wpcoef。1)Wpcoef:求解某个节点的小波包系数,数据长度是L/(2^n)(n表示分解的层数)比如:原始信号 ...
坛友互助群


扫码加入各岗位、行业、专业交流群


http://cn.qzs.qq.com/ac/b.gif
http://cn.qzs.qq.com/ac/b.gif

小波包分解主要用了这四个函数wpdec,wprec,wprcoef,wpcoef。

1)Wpcoef:
求解某个节点的小波包系数,数据长度是L/(2^n)(n表示分解的层数)
比如:原始信号长度为1024,经过3层分解,则使用cfs=wpcoef(wpt,[3 0])得到128个数据点,用以描述原始信号的近似低频信息。
具体代码如下:

clc;

% clear all;

close all;

load noisdopp; x = noisdopp;

% Decompose x at depth 3 with db1 wavelet packets

% using Shannon entropy.

wpt = wpdec(x,3,'db1');

plot(wpt); % Plot wavelet packet tree wpt.

cfs = wpcoef(wpt,[3 0]);

figure;

subplot(211);

plot(x); title('Original signal');

subplot(212);

plot(cfs); title('Packet (3,0) coefficients');

运行效果图如下所示:

http://ugc.qpic.cn/adapt/0/f0d2080e-35b3-ce6b-4ab1-c1d331428526/800
原始信号长度为1024,经过3层分解后,原始信号的近似信号长度为128。但是,为什么c(3,0)包数据的幅值会比原始信号大,后面会作具体研究。
2)Wprcoef:
把某个节点的小波包系数重构,得到的是和原信号一样长度的信号。
具体代码如下:

% clear all;

close all;

load noisdopp; x = noisdopp;

% Decompose x at depth 3 with db1 wavelet packets

% using Shannon entropy.

t = wpdec(x,3,'db1','shannon');

plot(t); % Plot wavelet packet tree.

rcfs = wprcoef(t,[3 0]);

figure(2);

subplot(311);

plot(x); title('Original signal');

subplot(312);

plot(rcfs); title('Reconstructed packet 3,0)');

运行效果图如下所示:

http://ugc.qpic.cn/adapt/0/84f54015-52a8-71de-3357-4ce99a21def9/800
原始信号长度为1024,使用rcfs = wprcoef(t,[3 0])得到的近似信号长度为1024,不同于cfs = wpcoef(wpt,[3 0])。
3)分析wprcoef与wpcoef的异同
在本例中,原始信号长度为1024。使用wprcoef得到rcfs,其长度为1024;而使用wpcoef得到cfs,其长度为128。
为了方便显示,将上述组数据导入LabVIEW软件,通过波形图直观地查看这几组数据之间的区别。
将上述原始信号,rcfs,cfs信号存至Excel中,则使用如下命令:
>> xlswrite('C:\compare.xls',noisdopp','sheet1','A1');

>> xlswrite('C:\compare.xls',rcfs','sheet1','B1');

>> xlswrite('C:\compare.xls',cfs','sheet1','C1');

在LabVIEW软件中使用Report Generation Toolkit导入这三组数据,分别命令:original data,rcfs data,cfs data。
http://ugc.qpic.cn/adapt/0/91e2fe44-b3c9-2668-7660-70d34c9f6c29/800
http://ugc.qpic.cn/adapt/0/deaffe0e-2b66-f5f9-8dcf-7bec6f3212e9/800
在上图中,将cfs data和rcfs data(div8)作商,得到的商的值均为2.82843(恒值)。
结论:
wprcoef和wpcoef使用不同的方式对原始信号做小波包分解,均可实现对原始信号的分解,某一个节点处的信号的提取。 不同的是使用wprcoef得到rcfs,其长度为1024;而使用wpcoef得到cfs,其长度为128。将rcfs信号作8分频,其大致曲线一致,只是幅值不同。
疑问:

为什么使用cfs = wpcoef(wpt,[3 0])得到的分解信号,相对于原始信号,其幅值均增加至原始幅值的2.82843倍,有待进一步研究。

附录:

几组信号的实际数据对比图:

http://ugc.qpic.cn/adapt/0/6ae0cd6e-2d4e-f8a3-f854-a816e925aa0f/800

扫码或添加微信号:坛友素质互助


「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-2668779-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
数据分析师 人大经济论坛 大学 专业 手机版