对小波函数wenergy(C,L) 的计算方法的分析

通过Matlab 自带的小波函数[C,L] = wavedec(X,N,'name');可以目标分析数据

进行小波分解。并通过函数X = wrcoef('type',C,L,'wname',N)对小波低频数据a 和

高频数据d 进行重构。以下对该命令的计算逻辑进行简要分析:

主要结论:

该公式可以理解为计算低频信号a 能量和各高频信号d ,d ,…d 信号的能

n 1 2 n

量与总能量的比值。

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

设向量a=[2 4 5 8 6 7 8 9 1 5 8 7];

1. 对该向量采用db1 小波1 层分解,得到

a1=[ 3 3 6.5 6.5 6.5 6.5 8.5 8.5 3 3 7.5 7.5]

d 1= [-1 1 -1.5 1.5 -0.5 0.5 -0.5 0.5 -2 2 0.5 -0.5]

(1)手动计算能量

∑12 2

a1 能量= (对a1 中所有数据求平方和)=462

=1

∑12 2

d 1 能量= (对d1 中所有数据求平方和)=16

=1

a1 能量占比= a1 能量/( a1 能量+ d 1 能量)= 97.7901%

d1 能量占比= a1 能量/( a1 能量+ d 1 能量)= 2.2099%

(2 )通过命令计算能量

而直接通过命令[Ea,Ed] = wenergy(c,l)计算得到 a1, d1 能量占比为

97.7901%,2.2099%

两种方法计算结果相等。

相关命令如下:

clear all;

a=[2 4 5 8 6 7 8 9 1 5 8 7];

[c,l] = wavedec(a,1,'db1');

a1=wrcoef('a',c,l,'db1',1);

d 1=wrcoef('d',c,l,'db1',1);

[Ea,Ed] = wenergy(c,l)

2. 对该向量采用db1 小波2 层分解,得到

a2=[ 4.75 4.75 4.75 4.75 7.5 7.5 7.5 7.5 5.25 5.25 5.25 5.25]

d1= [-1 1 -1.5 1.5 -0.5 0.5 -0.5 0.5 -2 2 0.5 -0.5]

d2=[ -1.75 -1.75 1.75 1.75 -1 -1 1 1 -2.25 -2.25 2.25 2.25]

(1)手动计算能量

∑12 2

a2 能量= (对a2 中所有数据求平方和)=425.5

=1

∑12 2

d 1 能量= (对d1 中所有数据求平方和)=16

=1

∑12 2

d2 能量= (对d2 中所有数据求平方和)=36.5

=1

a2 能量占比= a2 能量/( a2 能量+ d 1 能量+d2 能量)= 89.0167%

d1 能量占比= d 1 能量/( a2 能量+ d 1 能量+d2 能量)= 3.3473%

d2 能量占比= d2 能量/( a2 能量+ d 1 能量+d2 能量)= 7.6360%

(2 )通过命令计算能量

而直接通过命令[Ea,Ed] = wenergy(c,l)计算得到 a1, d1, d2 能量占比为

用matlab求小波系数能量,Matlab小波能量计算函数wenergy(C-L).pdf相关推荐

  1. 用matlab求小波系数能量,Matlab小波能量计算函数wenergy(C,L)

    对小波函数wenergy(C,L)的计算方法的分析 通过Matlab 自带的小波函数[C,L] = wavedec(X,N,'name');可以目标分析数据进行小波分解.并通过函数X = wrcoef ...

  2. Matlab 求纳什效率系数(NSE)

    Matlab求NSE系数 function nse = NSE(sim, obs) %输入参数分别为模拟值.实测值(序列) ave_obs = sum(obs(:))/numel(obs); %实测数 ...

  3. 小波系数等值线图和小波方差图绘制教学

    目录 1.数据格式的转化 2. 边界效应的消除或减小 3. 计算(复)小波系数 4. 计算小波系数的实部.模和模方 5. 绘制小波系数实部等值线图 6. 绘制小波系数模和模方等值线图 7. 绘制小波方 ...

  4. matlab求球的体积,【matlab计算不规则物体体积资讯】matlab计算不规则物体体积足球知识与常识 - 足球百科 - 599比分...

    matlab,求不规则三维图形体积? 我咋看也是个光滑曲面呀,哪来的体积?当然假如xy面到曲面之间全被填充为实体,应该能用循环或积分计算每个xy单位面积上的曲顶柱体的体积,然后求和,我用循环尝试了下: ...

  5. 如何用matlab求线性方程的通解,用matlab求线性方程的解

    用矩阵除法,求线性方程组的特解 利用矩阵除法求线性方程组的形式为:AX=b,其中A为系数矩阵,b为方程组右侧的一列常数. 由AX=b,得X=A\b,在这种条件下,首先要保证A为满秩矩阵,否则无法求解. ...

  6. matlab求被21整除,用matlab求[100,999]之间能被21整除的数的个数,还有建立一个字符串向量,删除大写字母._...

    怕散!|用户:的问题 用 matlab求[100,999]之间能被21整除的数的个数,还有建立一个字符串向量,删除大写字母. Promise丶陌影|用户:问问回答: clc;clear; m=0; f ...

  7. matlab求三维图形表面积,基于MATLAB的不规则平面立体表面积计算

    ·12· 设计与研究 机械 2012 年第 2 期 总第 39 卷 --------------- 收稿日期:2011-06-10 作者简介:和蕴锋(1988-),女,河南漯河人,硕士研究生,主要研究 ...

  8. matlab求系统全响应,利用matlab求LTI连续系统的响应

    <利用matlab求LTI连续系统的响应>由会员分享,可在线阅读,更多相关<利用matlab求LTI连续系统的响应(3页珍藏版)>请在人人文库网上搜索. 1.实验3 利用mat ...

  9. matlab求两方程交点坐标,matlab求交点坐标

    b xi ? 0(i ? 1, 2, , n) 三.线性规划问题的求解方法二元线性规划问题的图解法 线性规划问题的理论解法 线性规划问题的MATLAB软件解法 线性规划问题的图解...... MATL ...

最新文章

  1. 2018.11.12
  2. 聊一聊多源最短路径问题(只有5行代码哦)
  3. windows如何实现视屏自动定时、全屏、轮播 播放
  4. 2.3.1 为超参数选择合适的范围
  5. C++实现仿射加密法
  6. LeetCode之Move Zeroes
  7. 广州市城市智能交通大数据体系研究与实践
  8. 网页嵌入pdf、在线预览pdf工具及插件(转)
  9. 2018年内蒙古孕产妇、婴儿死亡率实现双下降
  10. [leetcode] Nim Game
  11. Veritas Backup Exec 21安装
  12. 【平面设计】ACDSee5.0软件安装教程
  13. java 使用itextpdf分割pdf文件,pdf分页
  14. 寻找春天nbsp;九宫格日记-2013.06.01
  15. 2.6.2 QDate类与QDateEdit控件介绍(日期处理)
  16. 10-不用加减乘除做加法
  17. 偏门赚钱项目:公众号打赏引流日赚500元
  18. 交换机与IP数据报、ICMP协议、ARP协议
  19. Spring入门示例
  20. HTML写一个登录框样式

热门文章

  1. 如何在微信公众号做广告不让人反感?
  2. Linux---MISC杂项驱动
  3. 全国各省市县 人口密度 数据 下载 空间数据 高精度 空间分布数据 多年 人口热力图 地理信息 GIS...
  4. Python代码写得丑怎么办?推荐几个神器拯救你
  5. cs go社区服务器维护,《CS:GO》即日起全面禁止社区服务器武器皮肤MOD
  6. Gerrit version 2.12.9 is now available
  7. WinSCP 的使用
  8. 沈阳师范大学 计算机专业分数低,辽宁省录取分数不高的3所师范类高校、就业率不低...
  9. java 回合制_java – 用于j2me“回合制”游戏的Gameloop
  10. iOS 讯飞语音的使用