EasyChair Preprint

№ 3136

Denoising Real-Time Audio Signals Using

Matlab Filter Techniques

K. C. Jayashree, Karthik M. Bharamappa, P. Manoj Kumar,

G. Mahaa Santeep and Sugandha Saxena

EasyChair preprints are intended for rapid

dissemination of research results and are
integrated with the rest of EasyChair.

April 8, 2020
Denoising real-time audio signals using matlab filter
techniques
Jayashree K C Karthik M Bharamappa Manoj Kumar P
School of ECE School of ECE School of ECE
REVA UNIVERSITY REVA UNIVERSITY REVA UNIVERSITY
Bangalore , India Bangalore , India Bangalore , India
jayashree240899@gmail.com karthikmbfaulkner@gmail.com pkj.manojkumar@gmail.com

Mahaa Santeep G Sugandha Saxena

School of ECE School of ECE
REVA UNIVERSITY REVA UNIVERSITY
Bangalore , India Bangalore , India
santeep1315@gmail.com sugandhasxn@gmail.com

Abstract— In this paper we analyze real time audio signals and creates major problem to the audio signals. Electrostatic noise
try to reduce the noise associated with the message signal under which is generated due to the presence of the voltages during
consideration. The main drawback of noise being present in an the design implementation and other random noises which
audio signal, is that it reduces the quality of the signal that is gets added to the signal. This proposed research will solve the
being transmitted within the communication system. For drawbacks of various filtering techniques which also
analysis purpose, white gaussian noise(awgn) is concatenated
with the audio signal under consideration and the resulting
provides unique Knowledge to the reader. The project defines
noisy audio signal is subjected to the different filtering a comparison between three different filtering techniques i.e.
techniques like IIR Filter, FIR Filter, Wavelet transform IIR filtering, FIR filtering and wavelets transforms based on
techniques. The noisy audio signal is analyzed with respect to real time audio signals.
the different filter responses obtained on applying the foresaid
methods. A comparative study is done between these techniques
to arrive at a technique which would be the most efficient one
II. BACKGROUND AND RELATED WORK
for audio signal denoising.
Graps[1] came up with new analysis named Fourier transform
which could analyze the periodic function by creating
Keywords-Audio processing, Denoising, FIR filters, IIR mathematical structures that vary in scale. But the proper
filters, wavelets, Daubechies. analysis cannot be analyzed using frequency response.
I. INTRODUCTION Radhika Bhagat[2] has made an attempt of audio filtering
The importance of noise reduction in real-time audio signals using extended filters like lowpass and highpass filter using
is said to be having high significance in communication. Due FIR and IIR filtering techniques. They have designed
to this factor the noise weakens signal quality, and the different formulas and difference equations for efficient
recognition of audio signals becomes difficult and cause implementation of time varying filter applications.
serious difficulties for the users of electronic hearing aids. A Er. Harpal Singh[3] has used fast Fourier transform technique
well-established method is filtering of the signal in the for performing time domain and frequency domain analysis
frequency domain or in the simplest way is analyzing the of the signal. Mannu Singla[4] has used Butterworth filter and
signal using different filter techniques like low pass, high Chebyshev filters to reduce noise from signals with different
pass and bandpass filters. As the voice or speech signals are frequency and ripple factors.
not periodic these filters distort the signal more than they
reduce the noise. To attenuate noise, we need more advanced Seema rani[5] has proved more facts about FIR and IIR filters
methods of filtering. As it requires high advanced methods of in their paper. The paper tells FIR is more stable than IIR.
filtering, we go for the different filtering techniques like From the above research we can conclude that the error of
Fourier transforms, wavelet transforms, and other methods FIR filter is less compared to IIR filters that means the output
are used in denoising the audio signals. A Fourier transform of FIR filter is very close to the desired value and FIR filter
of a signal gives us the frequency composition of the audio is more stable than IIR filters.
signal. The disadvantage of Fourier transform is it is only C Mohan Rao[6] has presented a new algorithm that is the
valid within a certain Region of Convergence (ROC). So, we Least Mean Square (LMS) in which the awgn is added to
go for short time Fourier transform (STFT), but this method message signal and the denoising is done to reduce the noise
used the window analysis approach of defined size. IIR and with minimum or no error efficiently. But it is sensitive to the
FIR algorithms uses the Fast Fourier transform (FFT) scaling of its input.
technique for analysis of frequency spectra and signal
responses. Wavelet analysis provide more detailed analysis B. Jai Shankar[7] has proposed the use wavelet
about the signal compared to other filtering approaches. transformation technique to denoise audio signals by dividing
the signal into blocks. This technique protects each and every
This paper tells the readers about denoising of the audio unique and vital features of every individual block and
signal. As all the audio signals are continuously affected by exposes the finest detail contributed by the grouped set of
the different types of real-time noise such as electrostatic blocks. The authors Ola Ratelli, Palle, Jorgensen (2013) [8]
noise, thermal noise, channel noise, awng and etc. which in their book have proposed the about Discrete wavelet
transform, its benefits and its functionalities. Priyanka Where x(n) is the input signal, y(n) is the output signal with
Khattar[9] have published a paper in which denoising will be
filter coefficients a and b.
performed using wavelet transformation by comparing two
wavelets families, Daubechies and Haar.
III. PROPOSED WORK
The flow of the project would be as shown in figure 1. Input Noisy IIR Filter
audio signal filter Response
signal

IIR
filter

Input Noisy Filter

audio Additive
signal FIR Response
signal White
filter
Gaussian
noise
wavelets
Additive
White Figure 2.IIR Filtering method block diagram
Gaussian
noise C. FIR Filter Algorithm
The audio signal which is to be analyzed is taken as input, to
Figure.1 Block diagram of project this additive white gaussian noise is added and the resulting
noisy audio signal is obtained. White gaussian noise is
Using basics of noise theory and audio theory, analysis of the preferred as it has almost constant PSD (power spectral
audio signal is carried out. The audio signal under density) and for easy and precise analysis.
consideration is concatenated with awgn. Using FFT, this
noisy signal is converted from time domain to frequency Now the frequency domain plot of the noisy signal is obtained
domain is then analyzed in frequency domain by converting by Fast Fourier method and the signal is analyzed with
from time to frequency domain using Fast Fourier Transform respect to the peak points to decide the cut-off frequency.
algorithm (FFT). The following step would involve designing Once the cut-off frequency is determined based on the
of the different filters taking into consideration the different analysis, the filter design process is initiated by passing the
parameters involved in filter design like order of the system, essential parameters like the order, normalized cut-off
cutoff frequency. The filtering approach consists of frequency, filter type (low pass filter) and appropriate
normalization of signal, decomposition technique and windowing technique is selected. The filter coefficients (the
reconstruction technique which we use in the filtering process numerator and denominator coefficients) are obtained the
to achieve our objectives. magnitude plot of the filter is plotted and analyzed. Now that
the filter is designed the normalized noisy signal is passed to
A. IIR Filter Algorithm low pass filter to obtain the filtered signal. The filtered signal
The audio signal which is to be analyzed is taken as input, to is finally plotted and used for analysis and comparison with
this awgn is concatenated with audio signal and the resulting the other results obtained.
noisy audio signal is obtained. White gaussian noise(awgn) is
preferred as it has almost constant PSD (power spectral
D. Equations
density) and for easy analysis of the audio signal.
Difference Equation of FIR filter:
Now the frequency domain plot of the noisy signal is obtained
by Fast Fourier Transform method and the signal is analyzed 𝑦(𝑛) = ∑𝑁−1
𝑘=0 ℎ(𝑘)𝑥(𝑛 − 𝑘) (2)
to decide the cut-off frequency. Once the cut-off frequency is Hamming Window function:
determined based on the analysis, the filter design process is
initiated by passing the essential parameters and low pass
2𝜋𝑛
filter is built. 𝑤(𝑛) = 0.54 − 0.46cos { } 0≤ 𝑛 ≤ 𝑁−1 (3)
𝑁−1
The filter coefficients are obtained and the magnitude plot of
the filter is plotted and analyzed. Now that the filter is
designed the normalized noisy signal is passed to low pass Input FIR Filter
Noisy
filter to obtain the filtered signal. The filtered signal is finally audio filter Response
signal
signal
plotted and used for analysis and comparison with the other
results obtained.
Additive
B. Equations White
Gaussian
Difference Equation of IIR Filter: noise
𝑦(𝑛) = − ∑𝑁 𝑀
𝑘=1 𝑎𝑘 (𝑛 − 𝑘) + ∑𝑘=1 𝑏𝑘 𝑥(𝑛 − 𝑘) (1)
Figure 3. FIR Filtering method block diagram
E. Wavelets
Wavelet transform method can be employed for analysis of T Approximate
H coefficients
the audio signal with respect to approximation and detailed R A0, A1…. An
coefficients. Wavelet transforms can be classified as, E
continuous wavelet transform and discrete wavelet transform. S
We are using discrete wavelet transform technique in our H Concatenation Reconstructed
O
project as it is more suitable for denoising of audio signals. L audio signal
The audio signal which is to be analyzed is taken as input, to D Detailed
this awgn is added and the resulting noisy audio signal is coefficients
obtained. White gaussian noise is preferred as it has almost S
/ D0, D1… Dn
constant PSD (power spectral density) and for easy and H
precise analysis. This signal is passed through Level filter.
The noisy signal is decomposed into two parts, detailed Figure.4b Block diagram of DWT reconstruction technique
coefficients and approximation coefficients. The number of
levels required for decomposition generally depends upon
IV. RESULTS AND WAVEFORMS
nature of the signal.
The audio signal is imported to MATLAB and plotted as
Multilevel decomposition is done to repeat the process of shown in figure.5
decomposition so that many lower resolution components of
the signal can be obtained through wavelet decomposition
trees. Wavelet thresholding, the final step is to reconstruct the
original audio signal without much loss of information. A
construction process which involves using the wavelet
coefficients and considering the levels of iteration, successful
reconstruction of the original audio signal is obtained.
DWT has two functions wavelet and scaling function

Scaling function ∅(𝑡) = ∑𝑁−1

𝑛=0 ℎ[𝑛]√2∅(2𝑡 − 𝑛) (4)

Wavelet function 𝜑(𝑡) = ∑𝑁−1

𝑛=0 𝑔[𝑛]√2∅(2𝑡 − 𝑛 ) (5)

Approximation coefficients:

1
𝑊∅ [𝑗0 , 𝑘] = ∑𝑛 𝑓[𝑛]∅𝑗0,𝑘 [𝑛] (6)
√𝑀
Figure.5 – Audio signal.
Detailed coefficients:
This audio signal is added with additive white gaussian noise
1 and plotted as shown in figure.6
𝑊𝜑 [𝑗, 𝑘] = ∑𝑛 𝑓[𝑛]𝜑𝑗,𝑘 [𝑛] j≥j0 (7)
√𝑀

T
Noisy Approx. H
Input Signal Coeffici
Low R
Audio ents
Pass E
Signal A0, A1, S
L-Filter …An H
Additive
White O
Gaussian High L
Detailed D
Noise Pass coefficie
nts S
D0, /
D1… Dn H
Figure.4b Block diagram of DWT decomposition technique

Figure.6 Noisy audio signa

The noisy signal is passed through IIR low pass filter having
magnitude & phase response is as shown in figure.7a.1 &
figure.7a.2 respectively and pole & zeroes plot is as shown in The output of IIR lowpass filter is as shown in figure 8.
figure.7b.

Figure.7a.1 Magnitude and phase response of IIR filter

Figure8.IIR low pass filter output (denoised output)

The same noisy signal is passed through the FIR filter having
magnitude & phase response as shown in figure.9a.1 & 9a.2
respectively and pole & zeroes plot as shown in figure.9b.

Figure.7a.2 Phase response of IIR filter

Figure.9a.1 Magnitude response of FIR lowpass filter

Figure.7b pole/zero plot of IIR filter

Figure.9a.2 Phase response of FIR lowpass filter

 Figure.9b pole/zero plot of FIR filter Figure.12 Scaling function-phi and Wavelet function psi of db3 wavelet

The output of FIR lowpass filter is as shown in figure.10

V. DISCUSSION
An analog system is said to be stable if all its poles lie in the
left half of the ‘s’ plane, but IIR filter has a both poles and
zeros as we can see in figure 7b. So IIR is said to be unstable.
IIR filters are difficult to implement as it has delays and
distortions due to large number of poles. IIR filters are better
in lower orders as IIR filters may become unstable with
higher orders. As IIR filters are unstable they cannot have a
linear phase as shown in the phase plot in figure.7a.
Theoretically we know that IIR filters depends on present
input and previous output. By the waveforms in figure.8, we
observe that noise introduced is denoised to certain limit but
the output audio signal will also get damped.
Due to these drawbacks we would prefer FIR over IIR as FIR
filters do not have delays and distortions. Therefore, FIR
filters are stable than IIR as FIR filters have only zeros on the
unit circle in the s plane and one pole at origin as shown in
Figure.10 Output of FIR lowpass filter figure.9b. FIR filters are better for higher orders and they
maintain stability. FIR filters have linear phase
When the same noisy signal is denoised using wavelets characteristics as shown in phase plot of figure.9a and it’s
transform technique, the output is as plotted in figure.11 output depends only on the present inputs. In order to avoid
damping of the output audio signal we have tried various FIR
filters using different windowing techniques like hamming,
Kaiser and rectangular. These windows were compared to get
the Filtered audio output but Hamming window was preferred
as it had linear phase. FIR filters had denoised to limit much
greater than IIR but the audio signal was also damped.
Due to the drawbacks we would prefer wavelets over FIR to
remove noise and get the filtered signal. The wavelets used
here is Daubechies wavelets have highest number N of
vanishing moments with the support width 2N-1. db wave
solves the problem of signal discontinuities and is applicable
for continuous and discrete wavelet transforms.
db wave belongs to orthogonal family and it has daughter
wavelets like db1, db2 till db45. db wavelet removes noise to
get filtered output. The figure.12 shows the Scaling function
phi and Wavelet function psi of db3 wavelet. Low pass filters
are used to avoid aliasing effect and is applicable in
communication circuits as anti-aliasing filters. Bandpass
Figure.11 Denoised signal from wavelet transform method filters are used to pass only certain range of frequencies
which is applicable in wireless transmitters and receivers in
order to avoid noise.
VI. FUTURE SCOPE
We have analyzed the audio signal by three filtering
techniques that is IIR, FIR and wavelets which suppresses the
noise using MATLAB functions. In some cases, the sound
signal which we consider as noise may also contain some
important information. For example, in military, satellite
applications minute information present in the sound signal
from the base stations plays a major role. So, the future scope
of this project is to not suppress the noise but segregate it
from the audio signal and play it as store and play it as another
signal. This project can also be used to separate the signal
which have undergone stereo mixing. For example, in a song
we can separate the male voice, female voice and background
music and play all the three as separate audio files. In this
way this algorithm can be used for generating karaoke from
any given song.

REFERENCES
[1] A. GRAPS, “AN INTRODUCTION TO WAVELETS” ,22724 MAJESTIC OAK
WAY, CUPERTINO, CA, USA,1995.
[2] Radhika Bhagat, Ramandeep Kaur, “Improved Audio Filtering Using
Extended High Pass Filters”,2013
[3] Er. Mannu Singla, Er. Harpal Singh, “Review Paper on Frequency Based
Audio Noise Reduction Using Different Filters”, International Journal of
Science, 2015.
[4]Er Mannu Singla, Er Harpal Singh, “Paper on frequency based audio noise
reduction using Butterworth, Chebyshev and Elliptical filters”, International
Journal of Science, 2015.
[5] Seema rani, Amanpreet Kaur, J S Ubhi, “Comparative study of FIR and
IIR filters for the removal of Baseline noises from ECG signal”, International
Journal of Computer Science and Information Technologies, Vol. 2 (3),
2011, 1105-1108.
[6] C Mohan Rao, Dr. B Stephen Charles, Dr. M N Giri Prasad, “A Variation
of LMS Algorithm for Noise Cancellation”, International Journal of
Advanced Research in Computer and Communication Engineering Vol. 2,
Issue 7, July 2013.
[7] Shanar B. J & Duraiswamy K. (2010). “Wavelet-Based Block Matching
Process: n Efficient Audio Denoising Technique” , European Journal of
Scientific Research, Vol.48 No.1, pp.16-28. ISSN 1450-216X.
[8] Ola Ratelli , Palle Jorgensen , ” Wavelets through a looking glass :
The world of the spectrum ” , Springer science & Business media ,
Illustrator : B. Treadway , 2013.
[9] Priya Khattar1, Dr. Amrita Rai2, Mr. Subodh Tripathi, “Audio De-
noising using Wavelet Transform”, 2016
[10] Adri E. Villanueva- Luna, Alberto Jaramillo-Nunez, Daniel Sanchez-
Lucero, Carlos M. Ortiz-Lima, “De-Noising Audio Signals Using MATLAB
Wavelets Toolbox”, Mexico.