小波变换与其分析.ppt

例1-1　在某工程实际应用中，有一信号的主要频率成分是由50 Hz和300 Hz的正弦信号组成，该信号被一白噪声污染，现对该信号进行采样，采样频率为1000 Hz。通过傅里叶变换对其频率成分进行分析。解　该问题实质上是利用傅里叶变换对信号进行频域分析，其MATLAB程序如下：　　t＝0：0.001：1.3；　　%时间间隔为0.001说明采样频 率为1000 Hz　　x＝sin(2*pi*50*t)＋sin(2*pi*300*t)；%产生主要频率为50 Hz和300 Hz的信号;　f＝x＋3.5*randn(1，length(t))；%在信号中加入白噪声　subplot(321)；plot(f)；　%画出原始信号的波形图　　　Ylabel(￠幅值￠)；　　　Xlabel(￠时间￠)；　　　title(￠原始信号￠)；　　　y＝fft(f，1024)；　%对原始信号进行离散傅里叶变换，参加DFT的采样点个数为1024　　　p＝y.*conj(y)/1024；　%计算功率谱密度;　　　ff＝1000*(0：511)/1024；%计算变换后不同点所对应的频率值　　　subplot(322)；plot(ff，p(1：512))；%画出信号的频谱图　　　Ylabel(￠功率谱密度￠)；　　　Xlabel(￠频率￠)；　　　title(￠信号功率谱图￠)；　　程序输出结果如图1.1所示。;　;　　从图1.1(a)中我们看不出任何频域的性质，但从信号的功率谱图(图1.1(b))中，我们可以明显地看出该信号是由频率为50 Hz和300 Hz的正弦信号和频率分布广泛的白噪声信号组成的，也可以明显地看出信号的频率特性。;7 小波变换简介 ; 1. 连续小波变换（CWT） 像傅立叶分析一样，小波分析就是把一个信号分解为将母小波经过缩放和平移之后的一系列小波，因此小波是小波变换的基函数。小波变换可以理解为用经过缩放和平移的一系列小波函数代替傅立叶变换的正弦波和余弦波进行傅立叶变换的结果。 图7-13表示了正弦波和小波的区别，由此可以看出，正弦波从负无穷一直延续到正无穷，正弦波是平滑而且是可预测的， 而小波是一类在有限区间内快速衰减到0的函数，其平均值为0， 小波趋于不规则、不对称。 ;图7-13 正弦波和小波 （a） 正弦波曲线； (b) 小波曲线 ; 从小波和正弦波的形状可以看出，变化剧烈的信号，用不规则的小波进行分析比用平滑的正弦波更好，即用小波更能描述信号的局部特征。 连续小波变换（Continuous Wavelet Transform， CWT）用下式表示： ; 基本小波函数ψ()的缩放和平移操作含义如下：  (1) 缩放。简单地讲， 缩放就是压缩或伸展基本小波， 缩放系数越小， 则小波越窄，如图7-14所示。 ; (2) 平移。简单地讲，平移就是小波的延迟或超前。在数学上， 函数f(t)延迟k的表达式为f(t-k)，如图7-15所示。 ; CWT计算主要有如下五个步骤： 第一步： 取一个小波， 将其与原始信号的开始一节进行比较。? 第二步： 计算数值C， C表示小波与所取一节信号的相似程度，计算结果取决于所选小波的形状， 如图7-16所示。 第三步：向右移动小波，重复第一步和第二步，直至覆盖整个信号，如图7-17所示。 第四步： 伸展小波， 重复第一步至第三步， 如图7-18所示。 ;图7-16 计算系数值C ;图7-17 计算平移后系数值C ;图7-18 计算尺度后系数值C ; 第五步：对于所有缩放，重复第一步至第四步。 小波的缩放因子与信号频率之间的关系是：缩放因子scale越小，表示小波越窄，度量的是信号的细节变化，表示信号频率越高；缩放因子scale越大， 表示小波越宽，度量的是信号的粗糙程度，表示信号频率越低。 ;常用的基本小波 ;常用的基本小波 ;;常用的基本小波 ;;4)Coiflet（coifN）小波系;5)SymletsA（symN）小波系;常用的基本小波 ;常用的基本小波 ;常用的基本小波 ;常用的基本小波 ;常用的基本小波 ;常用的基本小波 ; 2. 离散小波变换（DWT） 在每个可能的缩放因子和平移参数下计算小波系数，其计算量相当大， 将产生惊人的数据量，而且有许多数据是无用的。如果缩放因子和平移参数都选择为2j（j0且为整数）的倍数， 即只选择部分缩放因子和平移参数来进行计算， 就会使分析的数据量大大减少。使用这样的缩