小波变换主篇(2)Wavelet Transforms in One Dimension

1.The Wavelet Series Expansions
结合上一篇的知识点,我们知道一个函数f(x)可做如下分解:

f(x)=kcj0(k)φj0,k(x)+j=j0dj(k)ψj,k(x)(1) f ( x ) = ∑ k c j 0 ( k ) φ j 0 , k ( x ) + ∑ j = j 0 ∞ d j ( k ) ψ j , k ( x ) ( 1 )

其中
cj0=f(x)φj0,k(x)dx(2) c j 0 = ∫ f ( x ) φ j 0 , k ( x ) d x ( 2 )

dj=f(x)ψj,k(x)dx(3) d j = ∫ f ( x ) ψ j , k ( x ) d x ( 3 )

注意:j代表尺度当前空间的尺度,k代表当前空间的平移。
以上便是一维连续小波变换,本来应该给出一个例子,但是码代码太麻烦了,这里上传照片:
这里写图片描述
这里的基函数使用的是harr函数,计算比较复杂的就是通过上一篇的公式从harr函数求得其对应的 ψ(x) ψ ( x ) 函数,不过上述例子直接给出了 ψ(x) ψ ( x ) 函数,只用计算系数即可。由于第j0阶函数比较宽,所以其只有计算k=0的系数其余系数为0,而j1阶以上就得计算k=0,1…等多阶数据。

2.The Discrete Wavelet Transform
离散的小波变换相比于连续的小波变换主要是系数的计算方式发生了变化:

Wφ(j0,k)=1Mnf(n)φj0,k(n)(4) W φ ( j 0 , k ) = 1 M ∑ n f ( n ) φ j 0 , k ( n ) ( 4 )

Wψ(j0,k)=1Mnf(n)ψj,k(n)(5) W ψ ( j 0 , k ) = 1 M ∑ n f ( n ) ψ j , k ( n ) ( 5 )

其中,M为离散点的个数, φj0,k(n)=φj0,k(xs+nΔxs) φ j 0 , k ( n ) = φ j 0 , k ( x s + n Δ x s ) ,n=0,1,2,…,M-1.即将 φ(x) φ ( x ) 平均分割为M块,然后分别与f(n)相乘。
这个地方是与连续小波变换求系数有差别的,因为连续小波变换的系数是f(x)和 φ(x) φ ( x ) 的内积他们的交叠部分才进行相乘处理,而离散小波变换,f(x)中的x可能取1,2,3,4,5.但是其要与 φ(x) φ ( x ) x(0,1) x ∈ ( 0 , 1 ) 将其进行M等分进行相乘即f(1)* φ(1/M) φ ( 1 / M ) …..。即他们在x轴上的范围并不相交。
f(n)1MkWφ(j0,k)φj0,k(n)+1Mj=j0kWψ(j,k)φj,k(n)(6) f ( n ) 1 M ∑ k W φ ( j 0 , k ) φ j 0 , k ( n ) + 1 M ∑ j = j 0 ∞ ∑ k W ψ ( j , k ) φ j , k ( n ) ( 6 )

同样给一个例子,同样拍照:
这里写图片描述
下一篇将正式开始熟悉的小波变换。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值