基于MATLAB 7.x的系统分析与设计——小波分析 第三版 教学课件 作者 胡昌华 第2章.pptVIP

　　　　0　　0　　0　　　　7　　0　　8参见：dyaddown。 　　　　title(￠小波函数sym4的近似值(iter从1到10)￠)；　　　　holdoff　　　输出结果如图2.5所示。参见：intwave，waveinfo，wfilters。 　　freq＝［1：length(Lo_D)］/length(Lo_D)；　　subplot(427)；plot(freq，abs(fftld))；　　title(￠转换模量：低通￠)；grid； 　　subplot(428)；plot(freq，abs(ffthd))；　　title(￠转换模量：高通￠)；grid；输出结果如图2.9所示。　　检查双正交性的结果，df*rf＝：　　id＝　　　　1.0000　－0.0000　　　－0.0000　　1.0000参见：biorfilt，qmf，wfilters。 　　例2-11：　　%装载与某一正交小波相关的尺度滤波器　　loaddb10；　　subplot(621)；stem(db10)；title(￠db10低通滤波器￠)；grid；　　%计算积分镜像滤波器　　qmfdb10＝qmf(db10)；　　subplot(622)；stem(qmfdb10)；title(￠QMF低通滤波器￠)；grid；　　%检查频率条件 　　disp(￠频率尺度￠)　　disp(［a￠f￠per￠］)输出结果：　　频率尺度　　　2.0000　　 3.5294　　 0.2833 4.0000 1.7647 0.5667 8.0000 0.8824 1.1333 16.0000 0.4412 2.2667 32.0000 0.2206 4.5333 64.0000 0.1103 9.0667 128.0000 0.05511 8.1333 　　　set(axTITL，￠String￠，str1)　　　set(axXLAB，￠String￠，［str4，￠－￠，str5］)　　　clc；disp(strvcat(￠￠，str1，￠￠，str2，str3，￠￠，str4，str5))输出结果如图2.12所示。 　　　　　　　0.905 7　　　　　　　0.792 8　　　　　　　0.704 9　　　　　　　0.633 10　　　　　　　0.576 11　　　　　　　0.528 12　　　　　　　0.487 13　　　　　　　0.452 14　　　　　　　0.422 15　　　　　　　0.396 16 　　　　Pseudofrequencies＝0.792　　0.317　　0.158　　　　相应的尺度＝82040参见：centfrq。 　　　⑦OUT1＝wavemngr(￠read￠，IN2)；　　　⑧OUT1＝wavemngr(￠read_asc￠)。　　说明：该函数是一个用于小波管理的函数，它可以对小波进行新建、添加、删除、存储和读取操作。对于格式①、②可以把一个新的小波函数添加到工具箱中。 FN是FamilyName的缩写，FSN是FamilyShortName的缩写，WT是WaveletType的缩写。当WT＝1时，为正交小波；当WT＝2时，为双正交小波；当WT＝3时，为带尺度函数的小波；当WT＝4时，为不带尺度函数的小波；当WT＝5时，为不带尺度函数的复小波。 　　Meyer meyr　　DMeyer dmey　　Gaussian gaus　　Mexican_hat mexh　　Morlet morl　　ComplexGaussian cgau　　Shannon shan　　FrequencyBSpline fbsp　　ComplexMorlet cmor　　＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝ 　　程序：　　　　wavemngr(￠add￠，￠Lemarie￠，￠lem￠，1，￠12345￠，￠lemwavf￠)；　　　　wavemngr(￠read￠)　输出结果：　　ans＝　　＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝　　Haar　　　　　　　　　　　　haar　　Daubechies db　　Symlets sym　　Coiflets coif 　　Symlets sym--------------------------------　　　　sym2　　 sym3　　 sym4　　 sym5 　　sym6　　 sym7　　 sym8　　 sym**　　＝＝＝＝＝＝＝＝＝＝＝＝