The Applications of Wavelet Transforms in Analytical Chemistry Wu Renan #, Qian Dongyun (The School of Chemical Engineering of Wenzhou University #The Modern Education Technology Center of Wenzhou University Wenzhou 325027) Abstract The basic concepts of wavelet transform were introduced. The applications of wavelet transforms in analytical chemistry were classified into de-noising, data compressing and multiresolution analysis. Key words Wavelet transform, Analytical chemistry, Signal processing 摘要 介绍了小波变换的基本概念，将小波变换在分析化学信号处理中的应用概括为小波滤波、小波数据压缩和小波多尺度分析三个方面。 关键词 小波变换 分析化学 信号处理 小波变换在分析化学信号处理中的应用 吴仁安# 钱冬云 (温州大学化工学院 #温州大学现代教育技术中心 温州 325027) 小波变换及其应用是八十年代后期发展起来的应用数学分支，被称为“Fourier 分析方法的突破性进展[1]”。 1986 年 Meyer Y 构造了一个真正的小波基，十多年间小波分析及其应用得到了迅速发展，原则上传统的傅里叶分析可用小波分析方法取代[2],它能对几乎所有的常见函数空间给出通过小波展开系数的简单刻划，也能用小波展开系数描述函数的局部光滑性质，特别是在信号分析中，由于它的局部分析性能优越，因而在数据压缩与边缘检测等方面它比现有的手段更为有效[3-8]。 小波变换在图像压缩中的应用因它的高压缩比和好的恢复图像质量而引起了广泛的注意，且出现了各种基于小波变换的图像压缩方案。 小波变换自 1992 年 Bos M 等[9]首先应用于流动注射信号的处理，至今虽才 8 年时间，但由于小波变换其优良的分析特性而迅速渗透至分析化学信号处理的各个领域。本文介绍了小波变换的基本原理及其在分析化学中的应用情况。 1 基本原理 设 f(t)为色谱信号，其小波变换在 L2(R)中可表示为： 其中 a, b∈R,a≠0,参数 a 称为尺度因子 b 为时移因子，而(Wf)(b, a)称为小波变换系数，y(t)为基本小波。在实际分析化学信号检测中其时间是有限长度，f(t)通常以离散数据来表达，所以要采用 Mallat 离散算法进行数值计算，可用下式表示： fj+1=θ j + f j 其中:N为分解起始尺度;M为分解次数；fj 和 q j 可由下式求得： 此处:Φ j, m 为尺度函数;Ψ j, m 为小波函数;系数 Cmj ,dmj 可由下式表达： hk-2m , gk-2m 取决于小波母函数的选取。 用图表示小波分解过程如下： 图中 fN 、fN-1....fN-m 和 θ N-1、θ N-2....θ N-m 分别称为在尺度 N 上的低频分量和高频分量。上述分解过程的逆过程即是信号的重构过程。 2 分析化学中的应用 根据小波变换基本原理及其优良的多分辩分析特性，本文将小波变换在分析化学信号处理中的应用划归为以下三个方面： 2.1 信号的滤波 小波滤波方法目前在分析化学中应用主要是小波平滑和小波去噪两种方法。小波平滑是将某一信号先经小波分解，将在时间域上的单一信号分解为一系列不同尺度上的小波系数(也称不同频率上的信号), 然后选定某一截断尺度，使高于此尺度的小波系数全部为零，再重构信号，这样就完成了一个低通小波滤波器的设计；而小波去噪，则是在小波分解基础上选定一阈值，对所有尺度空间的小波系数进行比较，使小于此阈值的小波系数为零，然后重构信号[10]。 邵利民[11]等首次将小波变换应用于高效液相色谱信号的滤波，他们应用了 Haar 小波母函数，由三次小波分解后所得的低频部分重构色谱信号，结果成功地去除了噪声，明显地提高了色谱信号的信噪比，而色谱峰位保持一致，此法提高了色谱的最低检测量和色谱峰的计算精度。董雁适[12]等提出了基于色谱信号的小波自适应滤波算法，使滤波与噪声的频带分布，强度及信噪在频带上的交迭程度基本无关，具有较强的鲁棒性。 在光谱信号滤噪中的应用，主要为红外光谱和紫外光谱信号滤噪方面的应用，如 Bjorn K A[13]等将小波变换用于红外光谱信号的去噪，运用 6 种不同的小波滤噪方法(SURE，VISU，HYBRID，MINMAX，MAD 和 WP)对加噪后红外光谱图进行了去噪,针对加噪与不加噪的谱图，对 Fourier 变换、移动平均滤波与小波滤波方法作了性能比较研究，结果认为 Fourier 变换、移动平均滤波等标准滤波方法在信噪比很低时滤噪性能与小波滤波方法差不多，但对于高信噪比的信号用小波滤噪方法(特别是 HYBRID 和 VISU)则更有效 。闵顺耕[14]等对近红外漫反射光谱进行了小波变换滤波。顾文良[15]等对示波计时电信号进行了滤噪处理。王立世[16]等对电泳信号也做了小波平滑和去噪,都取得了满意的效果。邹小勇[17]等利用小波的时频特性去除了阶跃伏安信号中的噪音，并提出了样条小波多重滤波分析方法，即将过滤后的高频噪音信号当成原始信号进行滤波处理，使之对有用信号进行补偿。鲍伦军等[18]将样条小波和傅里叶变换联用技术应用于高噪音信号的处理。另外，程翼宇[19]等将紫外光谱信号的滤噪和主成分回归法进行了有机的结合，提出了小波基主成分回归(PCRW)方法，改善了主成分回归算法。 2.1 信号小波压缩 信号经小波分解之后，噪音信号会在高频部分出现，而对于有用的信号分量大部分在低频部分出现，据此可以将高频部分小波系数中低于某一阈值的系数去除，而对其余系数重新编码，只保留编码后的小波系数，这样可大大减少数据贮存量，达到信号压缩的目的。 在近代分析化学中分析仪器的自动化水平在不断提高，分析仪器所提供的数据量越来越大。寻找一种不丢失有效信息的数据压缩方法，节省数据的贮存量，或降低与分析化学信息处理有关的一些算法的处理量，已成为人们关心的问题。Chau F T 等[20]用快速小波变换对模拟和实验所得的紫外可见光谱数据进行了压缩，讨论了不同阶数的 Daubechies 小波基、不同的分解次数及不同的阈值对压缩结果的影响。Barclay V J 和 Bonner R F[10]对实验光谱数据作了压缩，压缩率可达 1/2~1/10,并指出在数据平滑和滤噪的同时，也能进行数据的压缩是小波有别与其他滤波方法的一大特点。王洪等[21]用 Daubechies 二阶正交小波基对聚乙烯红外光谱进行了成功的压缩，数据可压缩至原来的 1/5 以下。邵学广等[22]对一维核磁共振谱数据作了小波变换压缩，分别对常用的 Haar、Daubechies 以及 Symmlet 小波基作了比较，其结果表明准对称的 Symmlet 小波基对数据的复原效果最佳，而且在压缩到 64 倍时，均方差仍然较小。章文军等[23]提出了常用小波变换数据压缩的三种方法，将紧支集小波和正交三次 B-样条小波压缩 4-苯乙基邻苯二甲酸酐的红外光谱数据进行了对比，计算表明正交三次 B-样条小波变换方法效果较好，而在全部保留模糊信号及只保留锐化信号中数值较大的系数时，压缩比大而重建光谱数据与原始光谱数据间的均方差较小。邵学广等[24]将小波数据压缩与窗口因子分析相结合，在很大程度上克服了用窗口因子分析直接处理原始信号时人工寻找最佳窗口的困难，在压缩比高达 8:1 的情况下，原始信号中的有用信息几乎没有丢失，窗口因子分析的解析时间大为缩短。Bos M 等[25]用 Daubechies 小波对红外光谱数据进行压缩，压缩后的数据作为人工神经网络算法的输入接点，从而提高了人工神经网络的训练速度，预测的效果也比直接用光谱数据训练的要好。 2.3 小波多尺度分析 在多尺度分析方面的应用主要是对化学电信号进行小波分解，使原来单一的时域信号分解为系列不同频率尺度下的信号，然后对这些信号进行分析研究。 小波在色谱信号处理方面的应用，主要是对重叠色谱峰的解析。邵学广[26-27]等对苯、甲苯、乙苯三元体系色谱重叠峰信号小波变换后的某些频率段进行放大，然后重构色谱信号，使重叠色谱峰得到了分离，定量分析结果得到了良好的线性关系。此后邵学广[28]等利用了谱峰提取法对植物激素重叠色谱峰作了定量计算，此法表明，利用小波变换从重叠色谱信号中提取的各组分的峰高与浓度之间仍然具有良好的线性关系。 重叠伏安峰的分辨是电分析化学中一个长期存在的难题。当溶液中存在两种或更多的电活性物质，而这些物质的氧化(或还原)电位又很靠近时，就会不可避免地出现重叠峰的现象，而给进一步的定性、定量分析带来了很大困难。因此，人们做了较多的工作去解决这一难题。数学方法是目前处理重叠峰的重要手段，如 Fourier 变换去卷积以及曲线拟合。曲线拟合通常用来获得“定量”的信息，但这种方法有较多的人为因素，重叠峰包含的峰的个数，相对强度都是靠假设得来，因而可能引入严重的误差；去卷积方法则是一种频域分析手段，但该方法需先找出一个函数来描述伏安峰，然后再根据这个函数来确定去卷积函数，因此，去卷积函数的确定是比较麻烦的，尤其是对不可逆电极过程，无法找到一个合适的函数表达式，而且该方法还需经正、反 Fourier 变换，比较繁琐费时, 而小波分析的出现成了电分析化学家关注的热点。 陈洁等[29]用 DOG 小波函数处理差分脉冲实验数据，通过选择合适的伸缩因子，成功地延长了用 DPV 法测定 Cu2+的线性范围。郑建斌等[30-31]将小波变换用于示波计时电位信号的处理，在有用信息提取、重叠峰分辨等方面进行了系统的研究。王洪等[32]将小波边缘检测的思想用于电位滴定终点的确定，找到了一种判断终点准确的终点判断方法。郑小萍等[33]将样条小波变换技术用于分辨重叠的伏安峰，以选定的分辨因子作用于样条小波滤波器，构造了一个小波峰分辨器，用它来直接处理重叠的伏安峰，取得了较好的分离效果，被处理重叠峰可达到完全基线分离，且峰位置和峰面积的相对误差均较小。 对于红外光谱图，目前也是通过对红外谱图进行小波分解，以提高红外谱图的分辩率。陈洁[34]等对辐射合成的丙烯酰胺、丙烯酸钠共聚物水凝胶的红外光谱信号经小波处理后，使其特征吸收带较好地得到分离，成功地提高了红外光谱图的分辨率。谢启桃[35]等对不同晶型聚丙烯红外光谱图作了小波变换，也得到了可用以区分聚丙烯 a、b 两晶型的红外光谱图。 3 展望 小波变换由于其优良的局部分析能力，使其在分析化学信号的滤噪、数据压缩和谱峰的分离方面得到了很好的应用。本人通过对小波变换在化学中应用的探索，认为对于分析化学中各种电信号的平滑、滤波还有待作更深入的研究，以设计出更为合理有效的小波滤波器，以消除由于平滑而导至的尖锐信号的峰高及峰面积的变化或由于去噪而带来的尖锐信号附近的不应有的小峰的出现；对于重叠峰的分离及其定量计算，还应该探讨如色谱峰基线的确定方法以及待分离频率段的倍乘系数的确定方法；另外对于色谱峰的保留指数定性问题，由于不同化合物在某一确定的分析条件下有可能会出现保留值相同的情况，这将使在未知样中加标准的峰高叠加法定性或外部标准物对照定性变得困难，我们是否可能对色谱峰进行小波分解，然后在不同的尺度上对其进行考察，以寻求色谱峰的小波定性方法，这可能是个可以进一步研究的问题。 小波变换将在分析化学领域得到更加广泛的应用，特别对于分析化学中的多元定量分析法，如多元线性回归法(MLR)，主成分回归法(PCR)，偏最小二乘法(PLS)等方法及人工神经网络(ANN)将会同小波变换进行有机的结合，以消除各种噪声干扰对定量分析的影响；或对相关数据进行压缩以减少待分析数据的冗余，提高分析精度和大大减少计算量提高分析速度。小波变换将会成为分析化学中定量和定性分析的一种非常重要的工具。 4 参考文献 [1] 秦前清,杨宗凯.实用小波分析,西安:西安电子科技大学出版社,1993,1. [2] Meyer Y. Ondelettes. Hermann(Paris),1990:126. [3] Mallat S, Huang W L. Singularity detection and processing with wavelets. IEEE Trans, IT-38, 1992, 2:617-643. [4] Mallat S, Zhong S. Characterization of signals from multiscale edges. IEEE Trans, PAMI-14,1992,7:710-732. [5] Marr D C, Hildreth E. Theory of edge detection. Proc R Soc,1980,B207:187-217. [6] Canny J. A computational approach to edge detection. IEEE Trans, PAMI-8, 1986, 6:679-698. [7] 姜杭毅,蔡元龙.用于边缘检测的 Laplace 样条算子.中国科学A辑,1989,(10):1113-1120. [8] 王玉平,蔡元龙.多尺度 B 样条边缘检测算子.中国科学 A 辑,1995,25(4):426-437. [9] Bos M, Hoogendam E. Wavelet Transform for the Evaluation of Peak Intensities in Flow-Injection Analysis. Anal Chim Acta, 1992, 267:72-80. [10] Barclay V J, Bonner R F, Hamilton I P. Application of wavelet transforms to experimental spectra-smoothing, denoising and data set compression. Anal Chem, 1997,69(1),78-90. [11] 邵利民,唐兵,邵学广等.小波变换用于高效液相色谱的噪声滤除.分析化学,1997,25(1):15-18. [12] 董雁适,程翼宇.基于小波理论的化学谱图数据自适应滤波方法研究.高等学校化学学报,1999,20(3):378-382. [13] Bjorn K A, Andrew M W, Michael K W et al. Wavelet denoising of infrared spectra. Analyst, 1997,122:645-652. [14] 闵顺耕,谢秀娟,周学秋等.近红外漫反射光谱的小波变换滤波.分析化学,1998,26(1):34-37. [15] 顾文良,张书玲,张红权等.小波滤噪用于示波计时电位信号处理的研究.分析化学,1999,27 (2): 220-223. [16] 王立世,杨晓云,席绍峰等.毛细管电泳信号的小波平滑与去噪.高等学校化学学报,1999,20(3):383-386. [17] 邹小勇,莫金垣.阶跃伏安法信号的样条小波分析.高等学校化学学报,1996,17(10):1522-1527. [18] 鲍伦军,莫金垣,唐祖英.样条小波和 Fourier 变换联用技术在分析化学中的应用.高等学校化学学报,1998,19(12):193-197. [19] 程翼宇,陈闽军,方慧生.一类用于多元光度分析的小波基主成分回归法.分析化学,1999,27(2):170-173. [20] Chau F T, Shih T M, Gao J B et al. Application of the fast wavelet transform method to compress ultraviolet-visible spectra. Applied Spectroscopy, 1996, 50(3):339-348. [21] 王洪,司圣柱,熊建辉.红外光谱数据的小波压缩和重建.高等学校化学学报,1996,17(10)：1537-1540. [22] 邵学广,顾华,蔡文生等.一维核磁共振谱数据的小波变换压缩. 光谱学与光谱分析,1999,19(2)：139-141. [23] 章文军,许禄,刘胜雄.小波变换方法的比较-红外光谱数据的压缩.高等学校化学学报,1999,20(4)：534-538. [24] 邵学广,邵利民,侯树泉.基于小波变换的数据压缩在窗口因子分析中的应用.中国科学技术大学学报,1997,27(3):461-465. [25] Bos M, Vrielink J A M. The Wavelet Transform for Preprocessing IR-Spectra in the Identification of Mono-substituted and Disubstituted Benzenes. Chemometrics and Intelligent Laboratory System, 1994, 23(1):115-112. [26] 邵学广,孙培艳,蔡文生等.小波变换及其在色谱重叠峰解析中的应用.化学通报,1997,(8):59-62. [27] 邵学广,孙培艳,蔡文生等.小波变换用于色谱重叠峰的解析.分析化学,1997, 25(6):671-674. [28] 邵学广,侯树泉,方能虎等.小波变换用于植物激素重叠色谱峰中组分的定量分析.分析化学,1 998,26(1):107-110. [29] 陈洁,仲洪波,潘忠孝.小波变换在差分脉冲伏安法数据处理中的应用.分析化学,1996,24(9)：1002-1006. [30] 郑建斌,仲洪波等.小波变换用于二次微分交流示波计时电位信号中有用信息的提取.分析化学,1998，26(1)：25-28. [31] 仲洪波,郑建斌,潘忠孝等.小波变换用于示波信号中有用信息提取的研究.高等学校化学学报,1998,19(4)：547-549. [32] 王洪,潘忠孝,刘伟等.用小波变换确定电位滴定终点. 高等学校化学学报,1997,18(8):128691290. [33] ZHENG Xiaoping, MO Jinyuan, CAI Peixiang. Spline Wavelet in the Resolution of Overlapping Voltammetric Peaks, Science in China(Series B),1999,42(2),145-152. [34] 陈洁,单军,张懋森.共聚物红外光谱的小波变换.分析化学,1997,25(2):172-174. [35] 谢启桃,胡克良,潘忠孝等.小波变换用于不同晶型聚丙烯红外特征峰的鉴别.化学通报,1998,(2):61-64. 吴仁安 男， 31 岁，浙江瑞安，讲师，硕士，从事色谱及化学计量学研究。