这两天在做用小波阈值法进行信号去噪的问题,需要将信号进行五层小波分解,每进行一层分解时都将该层得到的细节系数CDi进行一次阈值化处理。这可以通过多次调用matlab中的dwt函数来实现。但是,在信号进行重构的过程中遇到了问题。

我原来的想法是,多次调用idwt函数来实现逐层重构。然而,遇到了系数长度不匹配的问题。比如,对于长度为1000点的数据,进行五层小波分解时,各层系数的长度依次是502,253,129,67,36。在用长度为36的CA5和CD5,采用idwt函数进行重构时,得到的CA4长度是68!而CD4长度是67!如果将CA4进行补零,则会导致最后重构的信号比原来长;如果将CD4截断,则会导致引入误差。即使取信号长度为2的整数次幂512,得到分解的长度依次是258,131,68,36,20,中间也会得到CA3长度为132,而CD3长度为131的情况。

这该如何解决呢?

我想到了MATLAB中小波重构函数waverec。它能够实现对信号的重构,并且能够解决上述系数长度不匹配的问题(虽然我还没搞懂它是怎么解决这一问题的,呵呵)。想要使用这个函数,就必须先弄清waverec需要的矩阵C和L中存储的是什么东西。然后将得到的一系列阈值化后的CA和CD组合成C,并由分解过程得到矩阵L。

经过多次实验,我发现,C是列向量,其中存储的是各层小波分解后得到的细节系数CDi和最后一层的近似系数CA。以5层分解为例,它们的存储结构是,C=[CA5;CD5;CD4;CD3;CD2;CD1]。L也是列向量,其中存储的是C中各组系数的长度,以5层分解为例,L=[len5;len5;len4;len4;len2;len1]。了解了这一点后,便可将阈值化后的系数组合成C,然后直接调用waverec进行信号重构了!

哈哈~
(2011-05-26)

MATLAB小波分解/重构wavedec/waverec函数中L和C的意义相关推荐

  1. matlab 小波的分解与重构

    小波的分解和重构在很多领域都有运用. 第一: matlab中相关小波函数 1. wavedec:  Multilevel 1-D wavelet decomposition    小波分解 [C,L] ...

  2. matlab 小波阈值去噪函数简介

    文章目录 目录 文章目录 前言 一.基于小波的各函数简介 1.小波分解函数dwt和wavedec 1)dwt:单尺度(单级)离散一维小波变换 2)wavedec:多尺度(多级)一维小波分解 2.提取各 ...

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

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

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

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

  5. matlab 小波启发式阈值滤波,小波阈值去噪

    小波阈值去噪 小波阈值去噪 目录 1.概念 2.原理 3.影响降噪效果的因素 3.1小波基的选择 3.2分解层数的选择 3.3阈值的选择 3.4阈值函数的选择 4.MATLAB代码 参考文献 小波分析 ...

  6. python小波分解与重构_小波分解和重构

    小波变换能够很好地表征一大类以低频信息为主要成分的信号, 小波包变换可以对高频部分提供更精细的分解 详见(http://www.cnblogs.com/welen/articles/5667217.h ...

  7. matlab dwt 多层分解,多尺度小波分解Matlab/Python实现与原理分析

    小波分解的理解 Matlab中实现小波分解方式 在matlab中实现小波分解的方式比较简单,通过[C,L] = wavedec(X,N,'wname')得到C和L,如下图,C表示分解得到的各个层的小波 ...

  8. Matlab小波工具箱常用函数

    1.Cwt :一维连续小波变换 格式:coefs=cwt(s,scales,'wavename') coefs=cwt(s,scales,'wavename','plot') scales:尺度向量, ...

  9. 一维信号小波去噪算法C语言,[转载]一维小波分解与去噪重构

    对随机一维信号实现多尺度小波分解,长度M=256,层数N=3. 解:一维随机信号是用nelec函数作为信号源,产生0~256范围内的信号. 实现小波按层分解的函数是: [C,L]=wavedec(s, ...

最新文章

  1. [高数][高昆轮][高等数学上][第二章-导数与微分]01.导数的定义
  2. 机器学习中的异常检测手段
  3. OSS上传图片无法在线预览的解决方案
  4. word文档图标变成白纸_挽救你的文件 修复变成乱码的Word文档
  5. 作为一个程序员,进步完全取决于自己
  6. libjpeg: cannot open
  7. 【POJ - 3272】Cow Traffic(dp,建反向图,DAG拓扑图)
  8. 手动安装lzop压缩工具 - JerryMo06的专栏 - 博客频道 - CSDN.NET
  9. sqlsrv 过去一小时_包裹分拣神器!每小时8.6万件包裹,是整个车间人数的两倍...
  10. ghost linux pe,winPE下安装系统ghost的详细教程
  11. 磁力链接地址生成教程
  12. 【基因组学】系统发育分析-进化树的相关知识点
  13. 作者序:互联网的负能量之声
  14. oracle10g没有行列转换函数的替代方法(转)
  15. SQL Server evaluation period has expired
  16. 微信公众号二次开发实现自动回复文字,图片,图文功能
  17. 基于steam的游戏销量预测 — PART 2 — 文本分析
  18. mac环境下搭建hexo+github pages+next个人博客
  19. 程序员必备:拯救颈椎,只需三个步骤
  20. 4个步骤,搭建接口测试自动化框架

热门文章

  1. sniffer技术原理及应用,包括编程方法和工具使用(2)
  2. Springboot项目集成Consul后,怎样不启用Consul呢?
  3. 计算机英语z,英文字母_计算机词汇(R-Z)_沪江英语
  4. 计算机主机启动不了系统安装,电脑开机后进不了系统只是在运行怎么办 电脑开机进不了系统原因【图文】...
  5. iFunk,科技享乐新主张
  6. 测试面试问道MySQL,面试---mysql相关
  7. Sketch-Based Skeleton-Driven 2D Animation and Motion Capture 翻译
  8. 盘点2021,拥抱变化,不负韶华
  9. 华为又来黑马招聘了!薪资待遇令人眼馋!
  10. “哈啰出行” 进入公共交通,哈罗的情商摩拜比不了!