Non-Linear Weighting Function for Non-stationary Signal Denoising

Similar documents
Adaptive Harmonic IIR Notch Filter with Varying Notch Bandwidth and Convergence Factor

A Robust Noise Spectral Estimation Algorithm for Speech Enhancement in Voice Devices

A soft decision decoding of product BCH and Reed-Müller codes for error control and peak-factor reduction in OFDM

Speech Enhancement using Temporal Masking and Fractional Bark Gammatone Filters

EQUALIZED ALGORITHM FOR A TRUCK CABIN ACTIVE NOISE CONTROL SYSTEM

FORWARD MASKING THRESHOLD ESTIMATION USING NEURAL NETWORKS AND ITS APPLICATION TO PARALLEL SPEECH ENHANCEMENT

APPLICATION OF THE FAN-CHIRP TRANSFORM TO HYBRID SINUSOIDAL+NOISE MODELING OF POLYPHONIC AUDIO

TESTING OF ADCS BY FREQUENCY-DOMAIN ANALYSIS IN MULTI-TONE MODE

Kalman Filtering for NLOS Mitigation and Target Tracking in Indoor Wireless Environment

Selective Harmonic Elimination for Multilevel Inverters with Unbalanced DC Inputs

NONLINEAR WAVELET PACKET DENOISING OF IMPULSIVE VIBRATION SIGNALS NIKOLAOS G. NIKOLAOU, IOANNIS A. ANTONIADIS

Cross-correlation tracking for Maximum Length Sequence based acoustic localisation

OTC Statistics of High- and Low-Frequency Motions of a Moored Tanker. sensitive to lateral loading such as the SAL5 and

Power Improvement in 64-Bit Full Adder Using Embedded Technologies Er. Arun Gandhi 1, Dr. Rahul Malhotra 2, Er. Kulbhushan Singla 3

An orthogonal multi-beam based MIMO scheme. for multi-user wireless systems

RAKE Receiver. Tommi Heikkilä S Postgraduate Course in Radio Communications, Autumn II.

Spectral analysis of biosignals. Biosignal processing I, S Autumn 2017

Estimating Travel Time Distribution under different Traffic conditions

DSI3 Sensor to Master Current Threshold Adaptation for Pattern Recognition

Reliable A posteriori Signal-to-Noise Ratio features selection

Phase Noise Modelling and Mitigation Techniques in OFDM Communications Systems

Comparison of Fourier Bessel (FB) and EMD-FB Based Noise Removal Techniques for Underwater Acoustic Signals

Using Adaptive Modulation in a LEO Satellite Communication System

Relation between C/N Ratio and S/N Ratio

Overlapping Signal Separation in DPX Spectrum Based on EM Algorithm. Chuandang Liu 1, a, Luxi Lu 1, b

SECURITY AND BER PERFORMANCE TRADE-OFF IN WIRELESS COMMUNICATION SYSTEMS APPLICATIONS

Improved Codebook-based Speech Enhancement based on MBE Model

ELEC2202 Communications Engineering Laboratory Frequency Modulation (FM)

Notes on Orthogonal Frequency Division Multiplexing (OFDM)

Efficient Non-linear Changed Mel-filter Bank VAD Algorithm

Windowing High-Resolution ADC Data Part 2

Novel half-bridge inductive DC-DC isolated converters for fuel cell applications

Compound quantitative ultrasonic tomography of long bones using wavelets analysis

COMBINED FREQUENCY AND SPATIAL DOMAINS POWER DISTRIBUTION FOR MIMO-OFDM TRANSMISSION

ROBUST UNDERWATER LOCALISATION OF ULTRA LOW FREQUENCY SOURCES IN OPERATIONAL CONTEXT

Track-Before-Detect for an Active Towed Array Sonar

Design and Implementation of Block Based Transpose Form FIR Filter

Additive Synthesis, Amplitude Modulation and Frequency Modulation

Evaluation of Steady-State and Dynamic Performance of a Synchronized Phasor Measurement Unit

This is an author-deposited version published in: Eprints ID: 5737

PREDICTING SOUND LEVELS BEHIND BUILDINGS - HOW MANY REFLECTIONS SHOULD I USE? Apex Acoustics Ltd, Gateshead, UK

Evolutionary Computing Based Antenna Array Beamforming with Low Probabality of Intercept Property

Fundamental study for measuring microflow with Michelson interferometer enhanced by external random signal

Mitigation of GPS L 2 signal in the H I observation based on NLMS algorithm Zhong Danmei 1, a, Wang zhan 1, a, Cheng zhu 1, a, Huang Da 1, a

Alternative Encoding Techniques for Digital Loudspeaker Arrays

Evolutionary Computing Based Antenna Array Beamforming with Low Probabality of Intercept Property

Quality-enhanced Voice Morphing using Maximum Likelihood Transformations

Iterative Receiver Signal Processing for Joint Mitigation of Transmitter and Receiver Phase Noise in OFDM-Based Cognitive Radio Link

Ruohua Zhou, Josh D Reiss ABSTRACT KEYWORDS INTRODUCTION

A Preprocessing Method to Increase High Frequency Response of A Parametric Loudspeaker

Analysis of Time-Frequency Energy for Environmental Vibration Induced by Metro

Precise Indoor Localization System For a Mobile Robot Using Auto Calibration Algorithm

Keywords Frequency-domain equalization, antenna diversity, multicode DS-CDMA, frequency-selective fading

Speech Enhancement for Nonstationary Noise Environments

Linear MMSE detection technique for MC-CDMA

PARAMETER OPTIMIZATION OF THE ADAPTIVE MVDR QR-BASED BEAMFORMER FOR JAMMING AND MULTIPATH SUPRESSION IN GPS/GLONASS RECEIVERS

Adaptive noise level estimation

A NEW APPROACH TO UNGROUNDED FAULT LOCATION IN A THREE-PHASE UNDERGROUND DISTRIBUTION SYSTEM USING COMBINED NEURAL NETWORKS & WAVELET ANALYSIS

MMSE STSA Based Techniques for Single channel Speech Enhancement Application Simit Shah 1, Roma Patel 2

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

Robust Acceleration Control of Electrodynamic Shaker Using µ Synthesis

ARCING HIGH IMPEDANCE FAULT DETECTION USING REAL CODED GENETIC ALGORITHM

SUBJECTIVE QUALITY OF SVC-CODED VIDEOS WITH DIFFERENT ERROR-PATTERNS CONCEALED USING SPATIAL SCALABILITY

Implementation of Adaptive Viterbi Decoder

Statistical Singing Voice Conversion with Direct Waveform Modification based on the Spectrum Differential

POWER QUALITY ASSESSMENT USING TWO STAGE NONLINEAR ESTIMATION NUMERICAL ALGORITHM

A comparison of LSF and ISP representations for wideband LPC parameter coding using the switched split vector quantiser

New Adaptive Linear Combination Structure for Tracking/Estimating Phasor and Frequency of Power System

Dictionary Learning with Large Step Gradient Descent for Sparse Representations

Signal Characterization in terms of Sinusoidal and Non-Sinusoidal Components

LUENBERGER ALGORITHM BASED HARMONICS ESTIMATOR FOR FRONT END RECTIFIER AND PWM-VSI

SIG: Signal-Processing

Performance of Multiuser MIMO System Employing Block Diagonalization with Antenna Selection at Mobile Stations

CONFIDENCE FEATURES EXTRACTION FOR WYNER-ZIV VIDEO DECODING

Radar Imaging of Non-Uniformly Rotating Targets via a Novel Approach for Multi-Component AM-FM Signal Parameter Estimation

Improving Poor GPS Area Localization for Intelligent Vehicles

A Novel NLOS Mitigation Approach for Wireless Positioning System

Characterization and Modeling of Underwater Acoustic Communications Channels for Frequency-Shift-Keying Signals

A Novel TDS-FDMA Scheme for Multi-User Uplink Scenarios

High Impedance Fault Detection in Electrical Power Feeder by Wavelet and GNN

Statistical Singing Voice Conversion based on Direct Waveform Modification with Global Variance

International Journal of Advanced Research in Electronics and Communication Engineering (IJARECE) Volume 3, Issue 9, September 2014

A Novel High Resolution Spectrum Sensing Algorithm for Cognitive Radio Applications

Mismatch error correction for time interleaved analog-to-digital converter over a wide frequency range

EXPERIMENTATION FOR ACTIVE VIBRATION CONTROL

Speech Signal Enhancement Techniques

Direct F 0 Control of an Electrolarynx based on Statistical Excitation Feature Prediction and its Evaluation through Simulation

2. Continuous-wave modulation

WIPL-D Pro: What is New in v12.0?

Keywords: Equivalent Instantaneous Inductance, Finite Element, Inrush Current.

Multicarrier Interleave-Division Multiple Access Communication in Multipath Channels

Two Dimensional Linear Phase Multiband Chebyshev FIR Filter

Optimal Modulation Index of the Mach-Zehnder Modulator in a Coherent Optical OFDM System Employing Digital Predistortion

COMPARISON OF TOKEN HOLDING TIME STRATEGIES FOR A STATIC TOKEN PASSING BUS. M.E. Ulug

Ultrasonic Beamforming with Delta-Sigma Modulators

Transmit Beamforming and Iterative Water-Filling Based on SLNR for OFDMA Systems

NINTH INTERNATIONAL CONGRESS ON SOUND AND VIBRATION, ICSV9 PASSIVE CONTROL OF LAUNCH NOISE IN ROCKET PAYLOAD BAYS

On-line Discrete Wavelet Transform in EMTP Environment and Applications in Protection Relaying

DIRECT MAPPING OVSF-BASED TRANSMISSION SCHEME FOR UNDERWATER ACOUSTIC MULTIMEDIA COMMUNICATION

Performance Evaluation of UWB Sensor Network with Aloha Multiple Access Scheme

Transcription:

Non-Linear Weighting Function for Non-stationary Signal Denoising Farès Abda, David Brie, Radu Ranta To cite this version: Farès Abda, David Brie, Radu Ranta. Non-Linear Weighting Function for Non-stationary Signal Denoising. IEEE. May 26, IEEE, pp.cdrom, 26. <hal-121584> HAL Id: hal-121584 https://hal.archives-ouvertes.fr/hal-121584 Subitted on 21 Dec 26 HAL is a ulti-disciplinary open access archive for the deposit and disseination of scientific research docuents, whether they are published or not. The docuents ay coe fro teaching and research institutions in France or abroad, or fro public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de docuents scientifiques de niveau recherche, publiés ou non, éanant des établisseents d enseigneent et de recherche français ou étrangers, des laboratoires publics ou privés.

NON-LINEAR WEIGHTING FUNCTION FOR NON-STATIONARY SIGNAL DENOISING Farès Abda, David Brie, Radu Ranta CRAN UMR 739 CNRS-UHP-INPL, B.P. 239, 5456 Vandœuvre-lès-Nancy, France ABSTRACT We propose in this paper a new strategy for non-stationary signals denoising based on designing a tie-varying filter adapted to the signal short ter spectral characteristics. The basic idea leading us to use a new paraetric nonlinear weighting of the easured signal short ter spectral aplitude (STSA) is exposed. The overall syste consists in cobining the estiated STSA and the coplex exponential of the noisy phase. The proposed technique results in a significant reduction of the noise for a variety of non-stationary signals including speech signals. 1. INTRODUCTION The denoising proble consisting in enhancing a signal degraded by uncorrelated additive noise when only the noisy easure is available, is a fundaental tas in signal processing. Because of the non-stationarity of the original signal, the short-tie Fourier analysis and synthesis concepts have been widely used [1 4]. The noisy signal is decoposed into spectral coponents by eans of the short-tie Fourier transfor. The advantages of the spectral decoposition results fro a good separation of the original signal and noise as well as the decorrelation of spectral coponents. Therefore, it is possible to treat the frequency bins independently and the estiation proble is siplified since the tie varying filtering procedure can be ipleented in the frequency doain where ultiplicative odifications are applied on the short-tie spectru and the enhanced signal is synthesized fro the odified short-tie spectru using the OverLap Add ethod (OLA). In general, it is significantly easier to estiate the STSA of the original signal than to estiate both aplitude and phase. There are a variety of denoising techniques that capitalize on this aspect by focusing on enhancing only the STSA of the noisy signal. We propose in this article a new strategy for the restoration of the STSA. We introduce a paraetric spectral gain forulation based on a continuous non-linear function, naely, the arctangent. The paper is organized as follows. In section 2 we introduce the necessary concepts of analysis and synthesis using the short-tie Fourier transfor, with a special attention to the OLA ethod and the effects of spectral odifications on the reconstructed signal. In section 3, we expose the overall schee of the denoising syste using the nonlinear spectral weighting ethod. In section 4, we discuss the epirical paraeters deterination. In section 5, we suarize the paper and draw conclusions. 2. SHORT-TIME FOURIER ANALYSIS AND SYNTHESIS Let x(n) and d(n) denote the original signal and the statistically independent additive noise respectively. The set of observations expressed using the odel y x(n) + d(n) is analysed using the short-tie discrete Fourier transfor defined as Y w(n )y()e j2π/l (1) where w(n ) is an appropriate sliding window (e.g. Haing window) of size N, and L N is the nuber of analysis frequencies ω = 2π/L for =,...,L 1. The denoising procedure ay be viewed as the application of a weighting rule, or nonnegative real-valued spectral gain G (n), to each bin of the observed short-tie spectru Y (n), in order to for an estiate ˆX (n) of the original signal shorttie spectru. We recall in the following the OLA ethod [1] for reconstructing the estiated signal ˆx(n) fro the odified short-tie spectru ˆX (n). Expressing a tie-varying ultiplicative odification to the short-tie spectru of y(n) as ˆX Y (n)g (n), (2) the synthesized signal is obtained by inverse Fourier transfor of ˆX (n), and suing over all the windowed fraes ˆx [ ] 1 Y ()G ()e j2πn/l. (3) L By defining the tie-varying ipulse response corresponding to G (n) as g n () = 1 G ()e j2πn/l, (4) L and using the definition of Y (n), we obtain fro (3) ˆx [ ] y(l) w( l)g (n l). (5) l

If we let r = n l or l = n r, the latter relation becoes ˆx [ ] y(n r) w( n + r)g (r). (6) r If we define ĝ(r n,r) = ĝ(q,r) = g (r)w( q) (7) w(n) y(n) F F T Y (λ). I F F T G (λ) Noise PSD estiation Y (λ) OverLap Add ˆx(n) then (5) becoes ˆx r y(n r)ĝ(r n,r). (8) We see fro equation (7) that the resulting tie-varying filter ĝ(r n,r) is a filtered version of g (n) by eans of the low-pass filter w(n). However, in our case the spectral tievarying filter G (n) is a nonlinearly obtained spectral weighting fro Y (n). Thus, tie-varying odifications applied to the easured short-tie spectru Y (n) are only dependent on the latter. Therefore, since Y (n) is a slowly tie-varying function (by taing into account the band liiting effect of the analysis window), the analysis window will not have any significant effect on the designed tie-varying filter. We do not discuss here the selection of the rate at which Y (n) should be sapled in both tie and frequency to provide an unaliased representation of Y (n). This point is detailed in [1] fro which it turns out that a properly sapled short-tie transfor using a Haing window requires on the order of four ties ore inforation as would be required relative to the original signal. In turn to this redundancy one obtains a very flexible signal representation for which extensive odifications can be ade. 3. NONLINEAR WEIGHTING OF THE SHORT-TIME SPECTRAL AMPLITUDE Most of the denoising techniques operate in the frequency doain by applying a frequency dependent gain function to the spectral coponents of the noisy signal, in an attept to heavily attenuate the noise only spectral coponents, while preserving those corresponding to the original signal. These denoising algoriths coonly consists of three ajor coponents: a spectral analysis/synthesis syste which was described in section 2, a noise estiation algorith, and a spectral gain coputation. The gain odifies only the Fourier agnitudes of an input frae. 3.1. Miniu statistics noise PSD estiation Noise estiation usually involves soe ind of SNR based signal activity detection which restrict the update of the noise estiate to periods of signal absence. However, these traditional noise estiation approaches are difficult to tune and their application to low SNR signals results often in distorted Fig. 1: Bloc diagra of a single channel signal denoising algorith. signals. Due to these reasons, we used a useful noise estiation ethod, nown as the iniu statistics [5 7], which consists in tracing the inia values of a soothed power estiate of the noisy signal. 3.2. The spectral gain calculation The tie-varying filter G (n) is specified by a pass region (G (n) 1) which covers the effective support region of the original signal x(n), a stop region (G (n) ) which has to suppress the undesired spectral coponents representative of the noise and an interediate region ( < G (n) < 1) where both signal and noise are present with values depending on the relative spectral values of the original signal and noise. Based on this intuitively appealing interpretation of the spectral weighting function G (n), and given an estiate of the noise PSD ˆσ 2 (n), we propose to design the tie-varying filter as G π arctan [β ( Y (n) αˆσ (n))] + 1 }, (9) 2 where > are design paraeters that have to be chosen in order to adjust at best the shape of G (n). The values of α and β control the suppression level of the noise as well as the resulting signal distorsion. This nonlinear gain function is depicted in figure (2). We see that this function taes very low values corresponding to sall spectral aplitudes of the easured signal copared to the estiated noise PSD, and taes values close to the unity for large values of the spectral aplitudes. Consequently, this function perfors a selective weighting of spectral coponents based on the noise PSD estiate. However, since the gain function (9) is sensitive to the dynaical range of the spectral aplitudes, the input signal has to be noralised with its axial value in order to obtain the sae variation range for all signals. An other solution to this proble consists in using the odified gain function expressed as G where [ ( )] π arctan β γ (n) α + 1 }, (1) 2 γ Y (n) 2 ˆσ 2 (n) (11)

is the a posteriori SNR as defined in [2]. The proposed short-ter spectral aplitude estiator is therefore given by ˆX (n) = = π arctan [ β ( γ (n) α )] + 1 2 } Y (n) G (n) Y (n). (12) Finally, by cobining this spectral aplitude estiate with the noisy phase of the easured signal, we obtain the shortter spectral estiate as ˆX G (n)y (n), (13) and the corresponding synthesized signal ˆx (n) is obtained using (3). Input SNR db 5dB 1dB 15dB 2dB Signal 1 Signal 2 Signal 3 Signal 4 α 3.9 2.54 2.25 1.68.56 β 3.18 2.58 1.98 1.45 1.48 α 3. 2.78 2.55 2.19 1.75 β 3.26 2.77 2.34 2.36 2.42 α 3.6 2.84 2.62 2.41 2.19 β 2.93 2.66 2.58 2.76 2.48 α 3.1 2.75 2.43 1.92 1.25 β 2.77 2.46 2.11 1.97 1.87 Table 1: Optial values of paraeters α and β obtained over 1 siulations for four different speech signals 1 and negative values indicates that there is a signal degradation. G (n).5 ασ (n) Fig. 2: Weighting function (9) for β = 3. 4. EMPIRICAL PARAMETER DETERMINATION We consider in this section the experiental deterination of the paraeters α and β using the Nelder-Mead algorith for the iniisation of the MSE between the noise free and the estiated signal calculated as 1 M [ M n=1 x(n) ˆx (n) ] 2, where M the length of the input signal. We used in a first study four speech signals corrupted by coputer-generated white Gaussian noises in order to have different SNRs. The obtained optial values of α and β with the corresponding output SNRs are given in table (1). Several siulations have been conducted and we found that reasonable values in the case of speech signals ust be chosen such that 2 < α 5 and 2 < β 5. However, if we choose α = 2.8 and β = 3, the algorith perfors high noise reduction in the ost cases. Furtherore, we found that this algorith, using the proposed values outperfors (in ter of SNR iproveents) the Wiener [8] and the MMSE-LSA [3, 9] algoriths both using the decision-directed approach for estiating the a priori SNR. Table (2) shows the obtained SNR iproveents for four different speech signals using the proposed algorith, the Wiener and the MMSE-LSA estiators for several input SNRs. The SNR iproveent is calculated as [ M ] n=1 (x(n) y(n))2 SNR = 1log 1 M n=1 (x(n) ˆx(n))2 (14) Input SNR db 5dB 1dB 15dB 2dB Signal 1 Signal 2 Signal 3 Signal 4 A 6.72 4.94 3.1.9-1.87 W 6.1 3.94 1.78 -.42-3.3 M 6.23 4.24 2.26.15-2.64 A 7.49 5.8 4.4 1.99 -.55 W 7.2 4.93 2.96.86-1.74 M 6.95 5.13 3.29 1.28-1.28 A 7.69 6.13 4.82 3.56 2.39 W 7.23 5.4 3.85 2.53 1.2 M 7.19 5.59 4.16 2.9 1.66 A 7.35 5.53 3.61 1.56-1.15 W 6.89 4.8 2.55.26-2.64 M 6.88 5.3 2.96.84-1.91 Table 2: SNR iproveents (in db) for the proposed algorith (A), the Wiener (W) and the MMSE-LSA (M) filters for four speech signals and different noise levels. As we can see fro table (2), the proposed ethod gives the best results. Moreover, the residual noise sounds very siilarly to that obtained using the MMSE-LSA algorith which is nown to provide less usical tones phenoenon [1]. We also tested the proposed algorith for denoising siulated sparse signals, i.e., a few large coefficients doinate the signal tie-frequency representation. In this case, we found it necessary to use a larger value for the paraeter α while aintaining the sae value for β, and the best results were obtained using the fixed value α = 5. Table (3) gives the SNR iproveents calculated using (14) for different input SNRs obtained for three test signals, naely, the linear chirp, the quadratic chirp and the two chirps signals corrupted with coputer-generated white Gaussian noise. The corresponding noise-free spectrogras are represented in Figure (3). The obtained results confir the efficiency of the pro-

posed ethod in denoising such signals..5.25 1 2 3 4 5 Tie (a) Linear Chirp 5 1 15.5.25 1 2 3 4 5 Tie 5 1 15 (b) Quadratic Chirp Fig. 3: Sparse test signals..5.25 1 2 3 4 5 Tie (c) Two Chirps Input SNR db 5dB 1dB 15dB 2dB Linear Chirp Quadratic Chirp Two Chirps A 16.9 14.92 13.7 1. 5.88 W 13.64 13. 11.82 9.3 5.6 M 11.62 11.23 1.44 8.46 5.22 A 15.8 13.96 11.83 8.13 3.74 W 13.35 12.17 1.78 7.57 3.57 M 11.5 1.66 9.65 7.1 3.34 A 12.24 11.97 11.22 9.51 6.47 W 11.64 11.21 1.39 8.76 6.14 M 1.3 1. 9.38 8.8 5.75 Table 3: SNR iproveents (in db) for the proposed algorith (A), the Wiener (W) and the MMSE-LSA (M) filters for three sparse test signals and different noise levels. 5 1 15 [2] R. J. McAulay and M. L. Malpass, Speech enhanceent using a soft-decision noise suppression filter, IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 28, no. 2, pp. 137 145, April 198. [3] Y. Ephrai and D. Malah, Speech enhanceent using a iniu ean-square error short-tie spectral aplitude estiator, IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 32, no. 6, pp. 119 1121, Deceber 1984. [4] P. J. Wolfe and S. J. Godsill, Siple alternatives to the Ephrai and Malah suppression rule for speech enhanceent, in Proc. 11th IEEE Worshop on Statistical Signal Processing, vol. II, pp. 496 499, 21. [5] Douglas B. Paul, The spectral envelope estiation vocoder, IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 29, no. 4, pp. 786 794, August 1981. [6] R. Martin, Spectral subtraction based on iniu statistics, in Proc. Eur. Signal Processing Conf., pp. 1182 1185, 1994. [7] R. Martin, Noise power spectral density estiation based on optial soothing and iniu statistics, IEEE Transactions on Speech and Audio Processing, vol. 9, no. 5, pp. 54 512, July 21. 5. CONCLUSION We proposed a siple and intuitively appealing ethod for non-stationary signal denoising consisting in a non-linear weighting of the STSA coefficients. We believe that the full potential of the proposed approach is not yet fully exploited, although very encouraging results were obtained. A future wor will concern the developent of regularization-based paraeters estiation fro only the noisy easure by taing into account a priori nowledge about the original signal. The sae non-linear function can be applied to develop a shrinage function in a wavelets denoising context [11]. In that respect, the derivation of a SureShrin-type threshold for the wavelet based denoising approach is also thinable since the corresponding shrining function have a bounded wea derivative [12]. These issues are now being investigated. 6. REFERENCES [1] J. B. Allen and L. R. Rabiner, A unified approach to short-tie fourier analysis and synthesis, Proceedings of the IEEE, vol. 65, no. 11, pp. 1558 1564, Noveber 1977. [8] J. S. Li and A. V. Oppehei, Enhanceent and bandwidth copression of noisy speech, Proceedings of the IEEE, vol. 67, no. 12, pp. 1586 164, Deceber 1979. [9] Y. Ephrai and D. Malah, Speech enhanceent using a iniu ean-square error log-spectral aplitude estiator, IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 33, no. 2, pp. 443 445, April 1985. [1] O. Cappé, Eliination of the usical noise phenoenon with the Ephrai and Malah noise suppressor, IEEE transactions on Speech and Audio Processing, vol. 2, pp. 345 349, April 1994. [11] F. Abda, Débruitage de signaux non-stationnaires : approches teps-fréquence et teps-échelle, Tech. Rep., Master degree, Université Henri Poincaré, Nancy 1, October 25. [12] J. Bigot A. Antoniadis and T. Sapatinas, Wavelets estiators in nonparaetric regression : A coparative siulation study, Journal of Statistical Software, vol. 6, pp. 1 83, 21.