利用充分光滑的S形函數(shù)構(gòu)造Meyer小波

邵云虹 鄧彩霞 賀鵬

摘 要：為了在信號或圖像的重構(gòu)中獲得較好的平滑效果，必須盡量增大小波的正則性或者連續(xù)可微性。在Meyer小波構(gòu)造中S形函數(shù)的選取影響著Meyer小波的可微性、光滑性和衰減速度等性質(zhì)，所以S形函數(shù)的選取至關(guān)重要。給出一種構(gòu)造充分光滑的S形函數(shù)的方法，并以一個充分光滑的非多項(xiàng)式S形函數(shù)為例，將其作為BP神經(jīng)網(wǎng)絡(luò)中的激勵函數(shù)進(jìn)行函數(shù)逼近得到好的逼近效果且訓(xùn)練次數(shù)少。然后通過充分光滑的S形函數(shù)得到Meyer小波的尺度函數(shù)，給出相應(yīng)的具有充分光滑、高階消失矩且無窮次可微性的頻譜有限的Meyer小波。最后把充分光滑的Meyer小波與剪切波變換結(jié)合進(jìn)行圖像去噪，與傳統(tǒng)的Meyer小波剪切波變換去噪相比較，峰值信噪比高于傳統(tǒng)的Meyer剪切波變換且去噪后的圖像紋理和邊緣信息保留更加完整。

關(guān)鍵詞：小波分析;Meyer小波;S形函數(shù);尺度函數(shù);高階消失矩

DOI：10.15938/j.jhust.2019.02.019

中圖分類號： O174.22

文獻(xiàn)標(biāo)志碼： A

文章編號： 1007-2683（2019）02-0127-08

Abstract：In order to obtain better smooth effect in signal or image reconstruction， the regularity or continuous differentiability of wavelet must be increased as much as possible. The selection of sigmoid function in Meyer wavelet construction affects the differentiability， smoothness and attenuation speed of Meyer wavelet， so it is exceedingly essential to select sigmoid function. Firstly， a method of constructing fully smooth sigmoid function is given， and a sufficiently smooth nonpolynomial sigmoid function is taken as an example， which is used as the excitation function in BP neural network for function approximation to obtain good approximation results. Then， the scale functions of Meyer wavelets are obtained by fully smooth sigmoid functions. Meyer wavelets with limited spectrum are given with sufficient smooth， highorder vanishing moments and infinitesimal differentiability. Finally， the full smooth Meyer wavelet and the shearlet transform are combined to denoise the image. Compared with the traditional Meyer wavelet shearlet transform， the peak signaltonoise ratio is higher and the denoised image texture and edge information are more complete.

Keywords：wavelet analysis; Meyer wavelet; sigmoid function; scaling function; high order vanishing moment

收稿日期： 2018-09-04

基金項(xiàng)目： 國家自然科學(xué)基金（11871181）.

作者簡介：

鄧彩霞（1965—），女，博士，教授，碩士研究生導(dǎo)師.

通信作者：

邵云虹（1994—），女，碩士研究生，Email：18846141926@163.com.

0 引 言

小波分析是80年代中期新興的一門學(xué)科，它是在Fourier分析的基礎(chǔ)上創(chuàng)建起來的[1]。由于它具有良好的時頻特性，因此可以更好地實(shí)現(xiàn)信號的處理[2]、小波降噪[3-4]、模式識別[5-6]等。1988年，Mallat提出多分辨率分析概念，為正交小波的構(gòu)造提供可行性方法[7]。正交小波按照其特點(diǎn)可以分成兩類：一類是頻譜有限函數(shù)，通常稱為Meyer型小波，例如Meyer小波;另一類是正交小波在時間域上具有緊支集，稱為Daubechies小波[8]。Daubechies小波是工程中常用的小波，然而許多能量有限信號是頻譜有限的，Daubechies小波在應(yīng)用中受到一定限制，因此研究在頻域內(nèi)具有緊支撐的Meyer型小波是十分重要的。

經(jīng)典的Shannon小波，它有理想的頻域性能，但是它在時域中的局部性卻很差。因此，法國數(shù)學(xué)家Meyer于1986年建議修正頻域上的Shannon尺度函數(shù)，既平滑了尖銳的邊緣值，又保留了其正交性條件，得到Meyer小波的尺度函數(shù)，進(jìn)而得到Meyer小波函數(shù)[9]。Meyer小波具有一定的可微性、衰減速度快和高階消失矩等性質(zhì)，因此Meyer小波有著廣泛的應(yīng)用。如Meyer小波可以用于解決方程的不適定問題[10]，實(shí)現(xiàn)對脈沖噪聲的信號處理[11]和利用Meyer小波的尺度函數(shù)進(jìn)行采樣[12]等。由于Meyer小波這些良好的性質(zhì)取決于它的尺度函數(shù)的性質(zhì)，而尺度函數(shù)的性質(zhì)與它的S形函數(shù)的性質(zhì)密切相關(guān)，所以研究不同的S形函數(shù)可以得到性質(zhì)不同的Meyer小波的尺度函數(shù)，再根據(jù)多分辨率分析方法給出相應(yīng)的Meyer小波函數(shù)。Meyer通過概率的形式引入了多項(xiàng)式型S形函數(shù)[13]，理論上可以獲得在（-∞，+∞）上任意有限階可導(dǎo)連續(xù)的多項(xiàng)式型S形函數(shù)，但隨著其可導(dǎo)階數(shù)的增加，它在[0，1]上的多項(xiàng)式次數(shù)也增加。例如常用的多項(xiàng)式型S形函數(shù)在[0，1]上為一個七次多項(xiàng)式x4（35-84x+70x2-20x3），但它在（-∞，+∞）上是三階可導(dǎo)連續(xù)的函數(shù)?？梢奙eyer給出的多項(xiàng)式型S形函數(shù)在[0，1]上的多項(xiàng)式次數(shù)越高越不便于數(shù)值計(jì)算。雖然存在三角函數(shù)型S形函數(shù)[14]，如在[0，1]上為sin2π2x，但它在（-∞，+∞）上僅為一階可導(dǎo)連續(xù)函數(shù)，光滑性較差。而S形函數(shù)的光滑性直接影響Meyer小波函數(shù)的可微性，可微性越高，Meyer小波的消失矩也越高。小波的消失矩階數(shù)越高，光滑函數(shù)在小波展開式中的零元就越多，也就是說小波的消失矩的大小決定了小波逼近光滑信號的能力，越大的消失矩將使高頻系數(shù)越小、小波分解后的圖像能量也就越集中、壓縮比例就越高，這一點(diǎn)也可以用來進(jìn)行圖像壓縮[15]和圖像分解[16]，所以考慮具有高階可微的S形函數(shù)是必要的。本文給出充分光滑的S形函數(shù)的一種構(gòu)造方法，得到相應(yīng)的Meyer尺度函數(shù)和小波函數(shù)，進(jìn)一步提升Meyer小波函數(shù)的可微性、衰減速度和消失矩等性質(zhì)。Meyer小波可以與剪切波變化結(jié)合進(jìn)行圖像去噪，剪切波是一種新型多尺度幾何分析工具，以普通小波作為基函數(shù)，對基函數(shù)進(jìn)行剪切、平移和伸縮變換而生成具有多方向的多分辨分析函數(shù)[17]，剪切波被廣泛應(yīng)用于圖像去噪[18-19]。剪切波為了去噪處理后獲取更佳的視覺效果，在選取基函數(shù)時，需要增加小波的光滑性或連續(xù)可微性，所以利用本文構(gòu)造的充分光滑的Meyer小波作為剪切波的基函數(shù)進(jìn)行圖像去噪，提高圖像的去噪效果。

1 多分辨分析及正交小波函數(shù)

通過仿真實(shí)驗(yàn)結(jié)果可以看出，本文選取的充分光滑的Meyer小波作為剪切波的基函數(shù)得到的剪切波在圖像去噪中峰值信噪比相對較高，說明本文的算法實(shí)現(xiàn)是有實(shí)際意義的，從圖6-8可以看出本文充分光滑的Meyer剪切波去噪效果更好，在邊緣和紋理的處理上也更加完整。

4 結(jié) 論

由于不同的S形函數(shù)構(gòu)成的尺度函數(shù)不同，通過改進(jìn)S形函數(shù)進(jìn)一步改進(jìn)由尺度函數(shù)構(gòu)成的Meyer小波函數(shù)的性質(zhì)。本文給出一種充分光滑的S形函數(shù)構(gòu)造方法，以一個充分光滑的非多項(xiàng)式型S形函數(shù)為例，將其應(yīng)用在BP神經(jīng)網(wǎng)絡(luò)函數(shù)逼近中作為一個新的激勵函數(shù)進(jìn)行仿真應(yīng)用，得到較好地逼近效果且用時相對較少。同時基于充分光滑的S形函數(shù)得到相應(yīng)的Meyer小波的尺度函數(shù)，根據(jù)Fourier變換及小波分析理論得到了Meyer小波函數(shù)時域和頻域的函數(shù)表達(dá)式，給出Meyer小波函數(shù)具有無窮次可微性、高階消失矩和快速衰減性等良好性質(zhì)。同時把構(gòu)造的充分光滑的Meyer小波作為剪切波的基函數(shù)，提高剪切波圖像去噪的效果，更好地保留圖像的紋理和邊緣等信息。

參 考 文 獻(xiàn)：

[1] GROSSMANN A， MORLET J. Decomposition of Hardy Functions into Square Integrable Wavelets of Constant Shape [J]. SIAM Journal on Mathematical Analysis， 1984， 15（4）： 723.

[2] MALLATS. A Wavelet Tour of Signal Processing [J]. Wavelet Analysis & Its Applications， 1999：74.

[3] 丁鋒， 秦峰偉. 小波降噪及Hilbert變換在電機(jī)軸承故障診斷中的應(yīng)用[J]. 電機(jī)與控制學(xué)報(bào)， 2017， 21（6）： 89.

[4] 孫鵬躍， 唐小妹， 黃仰博， 等. 利用實(shí)時小波降噪的抗電離層閃爍載波跟蹤[J]. 國防科技大學(xué)學(xué)報(bào)， 2017， 39（5）： 38.

[5] 尤波， 李忠杰， 黃玲， 等. 手部抓取動作特征提取算法研究[J]. 電機(jī)與控制學(xué)報(bào)， 2017， 21（12）： 75.

[6] 劉玉琪. 基于隨機(jī)森林算法的人體運(yùn)動模式識別研究[D]. 北京：北京郵電大學(xué)， 2018.

[7] MALLAT S. Multiresolution Approximations and Wavelet Orthonormal Bases of L2（R） [J]. ?Transactions of the American Mathematical Society， 1989， 315（1）： 69.

[8] 樊啟斌. 小波分析[M]. 武漢大學(xué)出版社， 2008： 92.

[9] MEYER Y. Uncertainty Principle， Hilbertian Bases and Operator Algebras [J]. Paris Bourbaki Seminar， 1986， 145-146（4）： 209.

[10]張占利. 一類橢圓型方程不適定問題的Meyer小波正則化方法[J]. 數(shù)學(xué)進(jìn)展， 2016， 45（3）： 390.

[11]JUN QIN， PENGFEI SUN. Applications and Comparison of Continuous Wavelet Transforms on Analysis of Awave Impulse Noise [J]. Archives of Acoustics， 2016， 40（4）： 84.

[12]SELVAN A A， RADHA R. An Optimal Result for Sampling Density in Shiftinvariant Spaces Generated by Meyer Scaling Function [J]. Journal of Mathematical Analysis & Applications， 2017， 451（1）： 197.

[13]BRANI VIDAKOVIC. Statistical Modeling by Wavelets [M]. Beijing：Higher Education Press， 2007： 52.

[14]彭玉華. 小波變換與工程應(yīng)用[M]. 北京：科學(xué)出版社， 1999： 21.

[15]QIAN L I， WANG Z. Image Compression Based on Wavelet Analysis [J]. Communications Technology， 2010， 25（2）： 1778.

[16]汪可， 張書琦， 李金忠， 等. 基于灰度圖像分解的局部放電特征提取與優(yōu)化[J]. 電機(jī)與控制學(xué)報(bào)， 2018， 22（5）： 25.

[17]EASLEY G， LABATE D， LIM W Q. Sparse Directional Image Representations Using the Discrete Shearlet Transform[J]. Applied & Computational Harmonic Analysis， 2008， 25（1）： 25.

[18]程義平. 基于平移不變性的剪切波變換醫(yī)學(xué)超聲圖像去噪算法[D]. 杭州：浙江工業(yè)大學(xué)， 2016.

[19]ABAZARI R， LAKESTANI M. A Hybrid Denoising Algorithm Based on Shearlet Transform Method and Yaroslavskys Filter [J]. Multimedia Tools & Applications， 2018， 77（14）： 17829.

[20]潘慶先， 董紅斌， 韓啟龍， 等. 一種基于BP神經(jīng)網(wǎng)絡(luò)的屬性重要性計(jì)算方法[J]. 中國科學(xué)技術(shù)大學(xué)學(xué)報(bào)， 2017（1）： 18.

（編輯：關(guān) 毅）