- 1、本文档共10页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
由于在的区间内,,所以上式又可写为:(9.7.2)由该式可以看出,小波变换可看作是和的卷积后的累加所得到的结果,卷积的中间变量是t,卷积后的变量为a及b。MATLAB中的cwt.m即是按此思路来实现的。从时间轴坐标的起点开始求积分,02先由指定的小波名称得到母小波及其时间轴上的刻度,假定刻度长为;01求和所选中的积分值序列的卷积,然后再作差分,即完成()式。04由尺度a确定对上述积分值选择的步长,a越大,上述积分值被选中的越多;03小波变换的大致过程:方法的不足:在a变化时,)式中括号内的积分、差分后的点数不同,也即和卷积后的点数不同。解决的方法:是在不同的尺度下对作插值,使其在不同的尺度下,在其有效支撑范围内的点数始终相同。有关CWT快速计算的方法还可借助于CZT及梅林变换等方法。例题:例9.7.1令为一正弦加噪声信号,它取自MATLAB中的noissin.mat。对该信号作CWT,a分别等于2和128,a=2时,小波变换的结果对应信号中的高频成份,a=128时,小波变换对应信号中的低频成份。其原始信号及变换结果见图9.7.1(a),(b)和(c)。图9.7.1信号“noissin”的小波变换原信号x(t),(b)a=2,(c)a=128例9.7.2仍然使用例的信号“noissin”,对其作CWT时a分别取10,30,60,90,120及150。所得到的图是在各个尺度下的小波系数的灰度图。颜色越深,说明在该尺度及该位移(水平轴)处的小波系数越大。此例旨在说明对小波变换的结果具有不同的表示方式。图9.7.2多尺度下小波变换的灰度表示9.8尺度离散化的小波变换及小波标架对同一个信号,在“时-频平面”a-b上,给出几种不同的表示形式:STFT:()Gabor变换:()WVD:()小波变换:()01020304目前通用的对a离散化的方法是按幂级数的形称为“半离散化二进小波”,而(9.8.6)式逐步加大a,即令。若取,则(9.8.5)称为二进小波变换。9.8.1尺度离散化的小波变换设:母小波的中心频率:,带宽:,当时,的中心频率变为,带宽。若时,的中心频率和带宽分别是:,。从对信号作频域分析的角度,希望当a由变成时,和在频域对应的分析窗和能够相连。。这样,当j由0变至无穷时,的傅里叶变换可以覆盖整个轴。设是的对偶小波,并令和取类似的形式,即(9.8.7)这样,通过对偶小波,我们希望能重建:(9.8.8)对上式作如下变换:123由恢复:由()和()式,有(9.8.9)显然,若(9.8.10)则()式的右边变成的傅里叶反变换,自然就是。01020304对于满足容许条件的小波,当时,其二进制小波对应的傅里叶变换应满足()式的稳定性条件。这样,结合()和()式,我们可由下式得到对偶小波:(9.8.11)由于()式的分母满足()式,因此有(9.8.12)这样,对偶小波也满足稳定性条件,也即,总可以找到一个“稳定的”对偶小波由()式重建出。定理9.4:如果存在常数,使得(9.8.13)则(9.8.14)如果满足(9.8.15)则(9.8.16)证明:由()式小波变换的定义,有将()式代入该式,有此即()和()式。余。的。这一结论告诉我们可以用
文档评论(0)