Daubechies小波介绍
Daubechies小波是正交、连续且紧支撑的。
正交条件下,$H(\omega)$必须满足下式:
$|H(\omega)|^2+|H(\omega + \pi)|^2 =1$
连续紧支撑条件下,$H(\omega)$必须满足下式:
$H(\omega)=\left( \frac{1+e^{-i\omega}}{2} \right) ^N S(\omega)$
$S(\omega)=\sum^A_{k=0}a_ke^{-ik\omega}$
从上式可以看出,$h[n]$可看做N个[1,1]的卷积再和$s[n]$卷积(不考虑归一化),而$s[n]$就是$S(\omega)$的系数:$s[n]=[a_0, a_1, a_2,...]$。
将紧支撑条件代入正交条件,简化后可得:
$(cos^2(\frac{\omega}{2}))^N|s(\omega)|^2+(sin^2(\frac{\omega}{2}))^N|S(\omega+\pi)|^2=1$
若将$s[n]$限定为实数,那么$|S(\omega)|^2$是偶函数,所以$|S(\omega)|^2$可写作下式:
$|S(\omega)|^2=\sum^A_{k=0}d_k cos^k (\omega) = \sum^A_{k=0} d_k ( 1 - 2sin^2( \frac{\omega}{2} ) )$
至此,我们将三角变量全部化成$\frac{\omega}{2}$,且上式只包含$cos$和$sin$的平方项。
我们暂时令$ y=sin^2 ( \frac{\omega}{2} ) $,$P(y)=\sum^A_{k=0}d_k(1-2y)^k$,就得到:
$(1-y)^NP(y)+y^NP(1-y)=1$
这种形式的方程的解为:
$P(y)=\sum^{N-1}_{k=0}\dbinom{2N-1}{k}y^k(1-y)^{N-1-k}$
然后将$y$和$P(y)$恢复成三角形式,以$\frac{1-cos \omega}{2}$代替$sin^2(\frac{\omega}{2})$,再将$cos \omega$以欧拉公式替换为复指数形式,最后以$z=e^{-i\omega}$替换复指数,最终得到:
$T(z)=|S(z)|^2 = \frac{1}{ 4^{N-1} } \sum^{N-1}_{k=0} \dbinom{2N-1}{k} (-z+2-z^{-1})^k (z+2+z^{-1})^{N-1-k}$
上式满足$T(z)=T(\frac{1}{z})$。于是我们就能得到一系列关于$T(z)$的零点的性质:
(a) 若$z_0 \neq 0$是$T(z)$的零点,那么$1/z_0$也是$T(z)$的零点
(b) 若$z_0$是$T(z)$的零点,那么$\overline{z_0}$也是$T(z)$的零点
考虑到将$z$的取值限定在单位圆上,我们有(考虑$|z|=|\overline{z}|=z\overline{z}=1$):
$|(z-w)(z-1/\overline{w})|=|(z-w)(z-z\overline{z}/\overline{w})|=|(z-w)z(1-\overline{z/w})|=|z-w||w|^{-1}|\overline{w-z}|=|w|^{-1}|z-w|^2$
这样,$T(z)=|S(z)|^2$就可以化为所有项均为平方的$z$的一次多项式连乘,且$z$的所有零点全部“翻转”到单位圆内成为最小相位系统。将$T(z)$开方(并归一化)后得到的$F(z)$代回到文章开头的$H(\omega)=\left( \frac{1+e^{-i\omega}}{2} \right) ^N S(\omega)$,即可求得Daubechies尺度滤波器系数。
下面是一个求六抽头Daubechies尺度滤波器系数的例子:
$N=3$
$P(y)=1+3y+6y^2$
$T(z)=1+3(\frac{1-z}{2})(\frac{1-1/z}{2})+6((\frac{1-z}{2})(\frac{1-1/z}{2}))^2$\
$=\frac{3}{8} z^2+\frac{9}{4} z+\frac{19}{4} z^{-1}+\frac{3}{8} z^{-2}$
$\alpha=\frac{3}{8}$
$T(z)$用计算机求解,取单位圆内零点:$z_1=0.28725+0.15289i, \overline{z_1}=0.28725+0.15289i$
$F(z)=\sqrt{\frac{3}{8}} \frac{1}{z_1}(z-z_1)(z-\overline{z_1}) $
于是$H(\omega)=(\frac{1+e^{-i\omega} }{2})^3 F(e^{-i\omega}) $展开后即得滤波器系数。(习惯上将求得的滤波器系数逆序排列)
Daubechies小波介绍相关推荐
- 用迭代法求Daubechies小波的尺度函数
p = [0.683013,1.183013,0.316987,-0.183013]; N = 20; phi_before = 1; for i = 1:N; phi_after = con ...
- 有关小波的几个术语及常见的小波基介绍
题目:有关小波的几个术语及常见的小波基介绍 本篇是这段时间学习小波变换的一个收尾,了解一下常见的小波函数,混个脸熟,知道一下常见的几个术语,有个印象即可,这里就当是先作一个备忘录,以后若有需要再深入研 ...
- 小波的几个术语及常见的小波基介绍
小波的几个术语及常见的小波基介绍 本篇是这段时间学习小波变换的一个收尾,了解一下常见的小波函数,混个脸熟,知道一下常见的几个术语,有个印象即可,这里就当是先作一个备忘录,以后若有需要再深入研究. 一. ...
- 【小波分析】小波的正交性,消失矩等常见概念及常见的小波基介绍
小波的几个术语及常见的小波基介绍 本篇是这段时间学习小波变换的一个收尾,了解一下常见的小波函数,混个脸熟,知道一下常见的几个术语,有个印象即可,这里就当是先作一个备忘录,以后若有需要再深入研究. 一. ...
- 基于小波的图像边缘检测,小波变换边缘检测原理
1.什么是"小波神经网络"?能干什么用呀 小波神经网络(Wavelet Neural Network, WNN)是在小波分析研究获得突破的基础上提出的一种人工神经网络.它是基于小波 ...
- 方波信号发生器电路仿真,小波神经网络算法原理
什么是"小波神经网络"?能干什么用呀 小波神经网络(Wavelet Neural Network, WNN)是在小波分析研究获得突破的基础上提出的一种人工神经网络. 它是基于小波分 ...
- 有关小波的几个术语及常见的小波基介
版权声明:本文为博主原创文章,转载请注明出处,谢谢! https://blog.csdn.net/jbb0523/article/details/42586749 目录(?)[+] 题目:有关小波的几 ...
- java 小波去噪原理_小波去噪的基本知识
本篇是这段时间学习小波变换的一个收尾,了解一下常见的小波函数,混个脸熟,知道一下常见的几个术语,有个印象即可,这里就当是先作一个备忘录,以后若有需要再深入研究. 一.小波基选择标准 小波变换不同于傅里 ...
- morlet小波matlab代码,morlet小波matlab程序
论-Matlab 7.0复Morlet小波分析操作实例.? 92? 2015年5 月 科技创新 中文科技期刊数据库 ( 文摘版 )自然科学 Matlab 7.0复MorletP]~ 波分析 操作 .. ...
最新文章
- php路由器怎么登录认证,PHP用户身份验证,如路由器登录
- java等号_java等号
- 打印 指定目录下和子目录下的的所有.java文件的路径. (使用FileFilter过滤器)
- Android应用程序进程启动过程
- python右对齐输出如果结果超过30个字符_Python3 输出字符串左对齐、右对齐、居中对齐...
- mysql不用left join_MySQL在右表数据不唯一的情况下使用left join的方法_MySQL - join
- 小程序 mathjs渲染公式_Mac 3D渲染和动画制作----KeyShot 9 Pro
- Docker的常用操作
- TopFreeTheme精选免费模板【20130703】
- 一个人生活可以多痛快
- 智能语音识别究竟是如何实现的?
- 《构建高可用Linux服务器 第3版》—— 2.7 系统维护时应注意的非技术因素
- 给扔物线 HenCoder Plus 学员的一次分享文字版
- 金万维异速联服务器重装,金万维异速联服务器配置说明
- 软件工程——软件计划
- 宏病毒实验讲解(含共享文件夹建立)
- 【基本功】Java魔法类:Unsafe应用解析
- 如何在Windows 10中打开设置?
- Excel数据的导入
- Mac node.js express-generator脚手架安装
热门文章
- 按键精灵--关于函数返回值实例
- 白鹭引擎之Scroller容器
- sinx/x的极限为什么是1_在sinx/x中当x→0时极限为什么为1
- 营销管理手册_营销管理学院新生入学讲座系列第四讲——学习红绿灯之学习学生手册...
- ug585_ch19 UART 导读
- 父母最应该看的10本教育经典书
- JAVA 根据汉字首字母排序
- small_arts:球球大作战刷(秒点/代点)棒棒糖网站地址
- 推荐贼好用的csdn浏览器助手插件!!
- 科技书籍与文章排版:Ctex 与 texlive 的中文编译模板