连续子波变换的MATLAB实现分析.pdfVIP

维普资讯 Feb．2005，Volume2，No．2(SerialNo．3) 通飙和 算楗 JournalofCommunicationandComputer,ISSN1548．7709，IJSA 逵缱子波燮换的MATLAB童琨分析 趟元英 ，袁 晚，滕旭柬，魏永豪 (四川大学鼋于信息学院，成都 610064) 摘 要：本文目的在于考察序列信我的连绩子波燮换 (CWT)在MATLAB中的程序贲现。首先介绍 子波熊换的基本概念和MATLAB中的算法原理，封MATLAB中的cwt．m程序逛行深入解请，阈述了各佃 参数的功能；而後餮现差分连算的存在使子波滤波器系数含有误差，继而封此误差的座生原因逛行分析， 解释如何稠整cwt．m算法中的参敷，以克服差分遵算的影喾，改善子波滤波器系数的取值，并且提出一檀 新算法，改善了CWT燮换的递算效果；最後以频移键控FSK信魏属例，分别用雨獯算法封之连行连绩子 波燮换，并封结果追行分析和比较，镫明了本文的改造算法侄于MATLAB中的原算法。 明键匍i子波庶波器；尺度；morlet子波；FSK信虢 算法中由于引入差分迎算而弹致的 差遣行了解 1．引言 析， 而提出一棰新算法，改善了CWT的逢算效果。 MATLAB是凿今世界上虑用最魔泛的欺罩软件 2．MATLAB中的CWT算法 之一，其最引人人腾的功能是具有各领域尊用的工 具箱，逭些工具箱的出现火大增强了MATLAB的寅 2．1算法原理及程序寅现 甩陛。子波分析是在现代羽和分析的基磋上嚣展起 子波燮换是一檀多分辨分析方法引，封于平方 来的一 翠科，具有深刻的理谕和魔泛的鹰用…。本 可稹函敷：s(t)e (尺)，其连镇子波燮换定羲焉 文研究了子波丁具箱中的连缱子波燮换 (CWT)的寅 现，群 解滔了MATLAB中的cwt_II1程序，并封歆 去 ()水 ∽)， ㈣ 式|}1a和b分别稻禺尺度参数和平移参敷。 ( )植分的差分造行卷稍，其中的差分序列就 在 MATLAB ·1『’基于子波燮换的定蓑式，届 便于数值 算，将信虢和子波函数都遣行丁雄散 是子波燮换畴采用的子波滤波器系数，算法的僵劣 化。令s(／)= (足)，t∈l，k+11，则有式 (2o由此 取决于此系敷的邃取是否精碓。 可兑，算法的本管是信虢s(k)舆子波滤波器函数 ， 士￡)v／()士 ∽()=击 ( (]，](2) 在 MATLAB中，基于上述算法思想，利用 CWt 是系数 coefs的给圃模式，xlim用柬放置横坐 的 命令封信虢造行连绩子波燮换。CWt命令的格式届 袍圈，coefs即禺封信虢s(t)连缒子波燮换後的系敷 functioncoefs：cwt(signal，scales，wname，plotmode， 值。通遇查找 toolbox／wavelet／wavelet／cwt．111文什， xlim)，signal是输入信虢 (足)，scales焉子波燮换 得到其主斡程序如下： 的尺度箍斟，wrlame是子波函数的名稻，plotmode 趟元英。女，四川大学鼋子信息学院硕士生，研究方向：信虢舆信，虑理。E—mail：zyyflyworld@163．com 7l 维普资讯 遑绩子波燮换的MATLAB寅琨分析