【正文】 ,L]=wavedec(b,2,39。db439。)。q=。p = C(126:492)。V = length(p)。I = (1:V)39。cVID = 1。cVN = sum(1./(1:V))。pID = p(max(find(p=I/V*q/cVID)))。pN = p(max(find(p=I/V*q/cVN)))。figureplot(C(126:493))E1=wden(b,39。rigrsure39。,39。h39。,39。one39。,2,39。db439。)。E2=wden(b,39。heursure39。,39。h39。,39。one39。,2,39。db439。)。E3=wden(b,39。sqtwolog39。,39。h39。,39。one39。,2,39。db439。)。E4=wdencmp(39。gbl39。,b,39。db439。,2,39。s39。,1)figureplot(E1)。figureplot(E2)。figureplot(E3)。figureplot(E4)。y1=0。y2=0。y10=0。y20=0。y11=0。y21=0。y12=0。y22=0。y13=0。y23=0。for i=1:479 y10=y10+(E1(i)c(i))^2。endfor i=1:479 y20=y20+c(i)^2。endSEM1=(y10/y20)^。for i=1:479 y11=y11+(E2(i)c(i))^2。endfor i=1:479 y21=y21+c(i)^2。endSEM2=(y11/y21)^。for i=1:479 y12=y12+(E3(i)c(i))^2。end for i=1:479 y22=y22+c(i)^2。endSEM3=(y12/y22)^。for i=1:479 y13=y13+(E4(i)c(i))^2。endfor i=1:479 y23=y23+c(i)^2。endSEM4=(y13/y23)^。 disp(SEM1)。注：其中a、b、c是在通過測試二氧化硫在給定波長的光譜吸光度，a為給定的波長數(shù)組，b是含噪的光譜吸光度數(shù)組，c是純凈的光譜吸光度數(shù)組。是給定的數(shù)據(jù)。參 考 文 獻(xiàn)[1] 楊福生,小波分析的工程分析與應(yīng)用[M].北京:科學(xué)出版社,1999:1145[2] [J].兵工自動化, 2002,(06) [3] [M].北京:[4] 龍瑞麟,[M].北京:[5] [M].西安:[6] 楊族橋,方志雄,[J].黃岡師范學(xué)院學(xué)報,:6062.[7] 胡昌華，張軍波等，基于MATLAB 的系統(tǒng)分析與設(shè)計——小波分析，第二版，西安：西安電子科技大學(xué)出版社，2003：前言[8] 檀承志,胡貴軍,[J].光子學(xué)報, 2007, 36(4):577580.[9] 秦前清,[M].西安:[10] 劉貴忠, [M].西安:[11] 譚利娜, 高協(xié)平, [J]. 光子學(xué) 報,2007,36（8）:15521556[12] Shao X, Cai W, Sun P, Zhang M, Zhao G. Quantitative Determination of the Components in Overlapping Chromatographic Peaks Using Wavelet Transform[J].Analytical Chemistry,1997,69(9):1722[13] , E. Hoogendam, Wavelet transform for the evaluation of peak intensities inflowinjection analysis, Anal. Chim. Acta, 1992, 267: 7380.[14] Mittermayer, ., Nikolov, ., Hutter, H., Grasserbauer, M., Wavelet Denoising of Guassian Peak: A Comparative Study,Chemometrics and Intelligent Laboratory Systems, 1996, 34(2): 187202.[15] Alsberg Bk, Woodward AM, Winson MK, Wavelet denoising of infrared spectra,Analyst, 1997, 122: 645652.[16]JouanRimbaud, Noord, , Application of wavelet transform to extract the Relevant ponent from Spectral data for multivariate calibration, Anal Chem, 1997, 69(21): 43174323.[17] Ingrid Daubechies 著，李建平，楊萬年譯，小波十講，北京，國防工業(yè)出版社，2003：1313。[18] 楊黎,[J].色譜,1998， 16(5）:386389[19] Donoho , Denoising by soft thresholding. IEEE Trans. On IT, 1995, 41(3):613627.[20] Mallat S. Multiresolution approximations and wavelet orthonormal based of L2(R).Trans. Aamer. Maath. Soc., 1989, 315: 6987.[21]張麗，羅超，張志沛，柴油機(jī)強噪聲信號的小波去噪辨識，內(nèi)燃機(jī)學(xué)報，2001，[22]Donoho DL, Johnstone IM, Wavelet Shrinkage: Asymptoia, Jounal of the Royal Satistical Society Series(B), 1995, 57: （4）：365368。[23] Donoho DL, Johnstone IM, Threshold Selection For Wavelet Shrinkage of Noisy Data,Engineering In Medicine And Biology Society, 1994, Engineering Advances。 New Opportunities for Biomedical Engineers. Proceedings of The 16th Annual International Coference of the IEEE,36Nov 1994, 16(1), A2425.[24] 周忠來，施聚生，小波分析去噪方法在聲目標(biāo)識別系統(tǒng)中的應(yīng)用研究，現(xiàn)代引信，1998，（4）：1925。[25]趙海英， Mallab中的實現(xiàn)[J]．兵工 自動化，2006，25(2)：54～55 第 25 頁 共 25 頁