[小波1]小波变换公式由来的深入理解

1.内积,点乘

b向量调小夹角之后和a向量更像,因此点乘可以用于衡量一个东西和另一个东西的相似度(Similarity)。点乘出来的数值越大,相似度越高

2c96ce7062474397a4d15bf0abb56d71.png

d26e7174ed4a48b0bdafa4341734f2cc.png 03a7cc021dd447f5beeacf5e48e882b5.png

4467193ed3e948d8861d02cab54a953a.png 33f82563833c44d088d27856638cbbd9.png

①归一化可以保证各个观察者能量一致的。

②方便后续数据处理

其中i和j的长度被归一化成1,是该方向的单位向量,是归一化过的观察者。

4和3是投影系数,是不同方向上观察者看到的目标的量是多少。

下面这张图展示了谁归一化都是等价的

7725058ed3624b2db3b7c61fa3aefabe.png


 2.连续函数点乘

直角坐标系中点乘计算很简洁,就是两个分量相乘后相加

3d24513dcb384f1faf976f431552d294.png65c29d434bf04362bf688d1c81c90129.pnge4cce737817e4044bd77d0f1ca674676.png

把n维向量这些数均匀撒到时间轴t上,范围是[t1, t2],即定义域是[t1, t2]的离散函数,n变大,在纵轴上面的数值越来越多,直至这个区间段都被占满,刚刚的离散函数就变成连续函数了,向量可以直接表示一条连续线段。

96bc97da6ffb4b5cac065e3c8e251b07.png1b3bd313e7c940bea10e69bc9e5a3daa.png

一个连续函数可以看做一个连续的无穷维向量。两个无穷维向量的点乘就是对应点的函数值相乘,当他们是连续函数时,就直接用两个函数表达式相乘即可。

c4ac4c5de3694e878fc41101589f7535.png

eq?X%28t%29Y%28t%29%3D%5Cint_%7Bt_%7B1%7D%7D%5E%7Bt_%7B2%7D%7DX%28t%29Y%28t%29dt 

上式为连续函数点乘的定义。这么积分就是求X(t)Y(t)和x轴围成的面积,上方的面积多一些,点乘的结果就是正的。

a1e950483e954b179616fc51b1611539.pngeq?%5Csin%20xeq?%5Ccos%202x围成的面积上下侧相等,点乘为0

从数学上看sinx和cos2x在[-eq?%5Cpi, eq?%5Cpi]范围内一点也不“像”


3.傅里叶级数的再解读

傅里叶级数以一系列简单的正余弦函数作为基,这些基(eq?%5Ccos%20nteq?%5Ccos%202teq?%5Ccos%203t)是一系列的观察者,这些观察者通过不同频率来区分的,有不同频率的观察者(1rad/s、2rad/s)。这些观察者可以看到目标函数在自己频率世界中的投影/影响是多少(目标函数和观察者自己本身的相似度是多少)

e551097577694a57aa231e04711d5573.png

下面是其中一个观察者和目标函数衡量有多“像”的计算步骤:eq?%5Ccos%20nt是周期函数,拿其中一个积分即可

2f8d8795c55d4695b0e6fb2a3749f13e.png

 要进一步看“相似度”——f(t)在eq?%5Ccos%20nt上面的投影是多少,需要求投影系数,则离不开归一化的操作

被观察者=eq?%5Csum归一化的投影系数eq?%5Ctimes观察者


如何看eq?%5Ccos%20nt这个观察者是否已经是归一化的?如果是向量有没有归一化,直接计算模长是不是等于1即可

                      78f1824eb2b54ebb9422d4c3fec54f19.png54c2717ce813407f853f28e508abe9e4.png

ba06283ab0c9444896a94f2f318632b8.png

连续函数=连续无穷维向量,连续函数的模长(2-范数)就是自己和自己点乘,然后再积分,后开方,发现cosnt模长不是1,那么需要用到前面的归一化系数公式进行计算:

3f6ac25ffd6e4f02a83b7cee4efa0e73.png

得到的这个公式和傅里叶级数求An的公式相同(对应到不同频率基的振幅),这样求得原函数在不同频率下的投影是多少:

 eq?%5Cfrac%7B1%7D%7B%5Cpi%20%7D%5Cint_%7B-%5Cpi%20%7D%5E%7B%5Cpi%20%7Df%28t%29cos%28nt%29dt


4.复指数形式的傅里叶变换

观察者由原先的余弦函数变为eq?%5Cexp%20%28iwt%29,把原来余弦函数观察者的观察范围由实数域拓展到带虚数的复数域。

在复平面和时间轴组成的空间里面是这样螺旋前进的:

72937df173b443788fa6f7631d032ddd.png

其中w为旋转的角频率(w越大旋转的越快),旋转的半径都是1,不同的观察者根据旋转的快慢进行分类。

对于一个负数计算模长的方法:

区别于实函数的模,复数是自己和自己的共轭相乘再开方。类似的,两个复函数的点乘是取后者的共轭函数进行计算的。

89e981dbff98481aa4e96471892707ac.pngdf437db7756e45c39fb66f0785388491.png

 观察者表示为半径为1的旋转,模长用1,用定义计算也是这样。

常用的连续函数的傅里叶变换公式:

eq?%5Cint_%7B-%5Cinfty%20%7D%5E%7B+%5Cinfty%20%7Df%28x%29e%5E%7B-i%5Clambda%20x%7Ddx

f16954bd58c04077b079340a4f39840e.png

 积分出来的结果是关于频率w的复函数eq?%5Cchi%20%28w%29eq?%5Cleft%20%7C%20%5Cchi%20%28w%29%5Cright%20%7C投影数,即该频率下震动的幅度。


以这个螺旋前进的函数在实数域的投影为例,由于eq?%5Cexp%20%28iwt%29 这个观察者在时间轴上是无穷无尽的,因此有缺点。派一个高频观察者对下面的信号进行观察,会发现如果问在什么时间段上和信号最近呢?这个观察者会说什么时间?因为观察者是在整个时间轴无穷无尽的!观察者一边,整个时间轴一起变。所以实际上在整条时间轴上一起看,就算只有部分像,我都说像。

f439ff956e464d4caeeba1fe6ac5da72.png766796ac556c41ca885b04308e74a149.png


 5.观察者局域化大改造

把我们想要的观察者局域化到一个小段时间内,把eq?%5Ccos%20nt乘以一个衰减系数eq?e%5E%7B-at%5E%7B2%7D%7D,这个衰减过的观察者在时间轴上就不是无穷无尽的了。

cb0913f322934102a78dbcbf620f2ea2.png

让观察者顺着时间轴去巡逻,去观察目标函数,然后把不同时间点观察到的结果反馈,那么高频的观察者知道时间轴前面的信号和它像,后面低频的和他不像。

b36126d145ed4d898eca40f346f90bc3.png

 观察者如何顺着时间轴进行巡逻呢?只需要给观察者的t减去一个数字b,当b变化时整个观察者就沿着时间轴跑起来了,b为平移因子,不同的b值对应不同的时间点。

9c2606177a5a4e15a5b84b240a0d52cc.png


 观察者的高低频是怎么控制的?因为想让不同频率的观察者去巡逻,胖胖的上图就是中心频率低,瘦瘦的上图就是中心频率高。里面的系数除以一个伸缩因子a,eq?%5CPhi%20%28%5Cfrac%7Bt-b%7D%7Ba%7D%29,可以控制观察者变胖变瘦,那么控制a的大小就能控制频率的高低,a又称尺度因子

若高频的观察者说目标函数和我长得更像,那么目标函数高频成分更多。

550423d291944de8ad5c7b0fd1d27108.png


 ①第一个观察者是高频的信号,一路走来没有发现目标函数和自己有一丁点相似的,右面的记录就是一条数值为0的直线。

9ad4f4af0d9141fc974079fc456aeb2b.png

 ②频率稍低一点的观察者上场,和这个频率的观察者相似度高的部分数值波动就高,如果不怎么像的数值就接近于0。

0eac0854ad414728b6e8f5f8f9059cf0.png

③得到了一个既有时间信息,又有频率信息的函数图像,这是一个关于平移因子和尺度因子的二元函数图像,这个图像是三维的,可以在这个图中得到任意频率在任意时间点的数值。

32c049f43a9a4e12abb7026fc8da61b5.png

0b5cc89f20f348c8a9d09ea16102f171.png


小波变换能够同时得到时间和频率的信息,是不是解决了时间不确定性的制约呢?——并没有

小波——一系列在时域上做了局域化的观察者

小波函数是一个非常庞大的家族,根据应用的场景和所需的特性不同,所选的小波不同:

191ee0a48e104b529e92d3e2f28b607e.png

 eq?f%28t%29为目标函数,变着胖瘦(a),顺着时间轴和目标函数进行点乘,再求积分eq?%5Cint_%7B-%5Cinfty%20%7D%5E%7B+%5Cinfty%20%7Df%28t%29%5CPhi%20%28%5Cfrac%7Bt-b%7D%7Ba%7D%29dx,得到的结果是关于时间和频率的二元函数,要想得到具体的数值,把a(频率)和b(时间)带进去即可。

问题:不同观察者eq?%5CPhi%20%28%5Cfrac%7Bt-b%7D%7Ba%7D%29和x轴围成的面积不一样,就要进行归一化,有一个eq?%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cleft%20%7C%20a%20%5Cright%20%7C%7D%7D的归一化因子

5b1cdad8d646435a941f02df78e0fa69.png3fcede995f6842c7b6eac408d97e8faa.png

归一化之后,变瘦的同时右面的图会变高,这样保证了围成的面积相同,能量一致。


6.小波变换的公式:

eq?C_%7Ba%2Cb%7D%3D%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cleft%20%7C%20a%20%5Cright%20%7C%7D%7D%5Cint_%7B-%5Cinfty%20%7D%5E%7B+%5Cinfty%20%7Df%28t%29%5CPhi%5E%7B*%7D%20%28%5Cfrac%7Bt-b%7D%7Ba%7D%29dx

7782fcd5be384739b58aa2f9ac0c0af0.png在时间轴上滚动着去和别人相乘求积分?实际上就是卷积!小波函数就是卷积核。

  • 21
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值