Next Article in Journal
Dimension Reduction of Digital Image Descriptors in Neural Identification of Damaged Malting Barley Grains
Next Article in Special Issue
Healthcare Application of In-Shoe Motion Sensor for Older Adults: Frailty Assessment Using Foot Motion during Gait
Previous Article in Journal
Path Loss Investigation in Hall Environment at Centimeter and Millimeter-Wave Bands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coherence Function and Adaptive Noise Cancellation Performance of an Acoustic Sensor System for Use in Detecting Coronary Artery Disease

School of Electrical Engineering, Computing and Mathematical Sciences (EECMS), Faculty of Science and Engineering, Curtin University, Bentley, WA 6102, Australia
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(17), 6591; https://doi.org/10.3390/s22176591
Submission received: 26 July 2022 / Revised: 26 August 2022 / Accepted: 28 August 2022 / Published: 31 August 2022
(This article belongs to the Special Issue Wearable Sensors and Internet of Things for Biomedical Monitoring)

Abstract

:
Adaptive noise cancellation is a useful linear technique to attenuate unwanted background noise that cannot be removed using traditional frequency-selective filters. Usually, this is due to the signal and noise co-existing in the same frequency band. This paper tests a weighted least mean squares (WLMS) algorithm on a stethoscope system for use in detecting coronary artery disease in the presence of background noise. Each stethoscope is equipped with two microphones: one used to detect heart signals and one used to detect background noise. The WLMS method was used for four different sources of background noise whilst measuring a heartbeat, including a single tone, multiple tones, hospital/clinic noise, and breathing noise. The magnitude-squared coherence between both microphones was unity for the tone scenarios, resulting in complete attenuation. For the other background noise sources, a less-than-unity magnitude-squared coherence resulted in minor and no attenuation. Thus, the coherence function is a tool that can be used to predict the amount of attenuation achievable by linear adaptive noise-cancellation techniques, such as WLMS, as presented in this article.

1. Introduction

Heart auscultation—the interpretation of heart sounds by a physician—is a common method utilised for cardiac diagnosis [1]. Many physicians are documented to have poor auscultatory skills due to their difficulty to acquire, yet it remains the primary method used for diagnosis and screening in primary healthcare [2,3,4]. The leading cause of mortality and morbidity in the world is cardiovascular disease, contributing 31% towards all global deaths in 2019 according to the World Health Organisation [5]. This is particularly worrisome for low-to-middle-income countries where high-quality diagnostics are difficult to obtain [6]. To address this large and complex problem, there is a global need for the affordable pre-screening of coronary artery disease (CAD).
Data quality will be a limiting factor in the achievable accuracy of any classification technique of CAD. Signal corruption from unwanted noise will limit data quality. In the available public datasets of normal and abnormal heartbeats, noise limits the accuracy achievable by machine learning and classifiers [7], with frequency-selective filtering being the only technique used to reduce noise. Thus, cancelling noise in the frequency band of interest is a crucial stage in acquiring data and has been implemented in the field of biomedical signals [8,9,10]. A comprehensive review of adaptive filters for noise cancellation in [11] identified existing methods derived from the least mean squares (LMS) method, highlighting it as a preferred method for noise cancellation due to its robustness and simplicity. This study shows that these techniques can be implemented on a digital stethoscope with background-noise-sensing capability and is preliminary to future classification studies. A weighted least mean squares (WLMS) adaptive filter algorithm is implemented in this study.
A system has been developed that incorporates six digital stethoscopes onto a wearable vest that simultaneously measures heartbeat signals and uses machine learning methodologies to pre-screen for CAD. Each stethoscope, shown in Figure 1, has two microphones. The first microphone is the heart-sensor microphone (HM) and is located behind the diaphragm. The second microphone is the background-noise microphone (BNM) located at the other end of the stethoscope. The assumption is that the HM will pick up the heart signal plus part of the background noise, and that the BNM will pick up just the background noise. Using the two-microphone configuration, it is hypothesised that the background noise can be reduced from the heart-sensor microphone. This will contribute to successful diagnosis techniques, as acquired signals from the system will be cleaner with higher signal integrity. The coherence function is the measure of the causal relationship between two signals with the presence of other signals [12] and is utilised in this paper to predict the attenuation of the background noise in the HM signal. Using a single tone whilst measuring a heartbeat, we tuned a WLMS adaptive filter algorithm to achieve complete attenuation from the HM. The parameters from this exercise were then used in this study and tested on other noise sources.

2. Background

Background noise can couple into the HM and corrupt acoustic heartbeat measurements. This can decrease signal integrity, as it cannot be filtered out using conventional frequency-selective filters if the noise lies within the frequency band of interest of 25–250 Hz [13]. To attenuate unwanted background noise or enhance the desired signal, adaptive filtering techniques and other non-linear techniques can be utilised.

2.1. Adaptive Noise Cancellation

Noise cancellation via a Wiener filter is a linear process, where the data are assumed to be linear-time-invariant (LTI). As explained in [14], a (p − 1)st-order Wiener filter, given by Equation (1), can be derived, which produces the minimum mean-square estimate (MMSE) of a given process d ( n ) by filtering a statically related process x ( n ) that has unwanted noise. The output of the filter d ^ ( n ) is the convolution of w ( n ) and x ( n ) given by Equation (2).
W ( z ) = n = 0 p 1 w ( n ) z n
d ^ ( n ) = l = 0 p 1 w ( l ) x ( n 1 )
Here, we assume that the auto-correlation functions r x ( k ) = E [ x ( n ) x ( n + k ) ] and r d ( k ) = E [ d ( n ) d ( n + k ) ] and the cross-correlation r d x ( k ) = E [ d ( n ) x ( n + k ) ] are known. The Wiener coefficients are calculated to minimize the mean-square error given by Equation (3).
E { | e ( n ) | 2 } = E { | d ( n ) d ^ ( n ) | 2 }
Minimizing Equation (3) derives the Winer–Hopf equations, shown in Equation (4):
R x w = r d x ,
where w represents the filter coefficients, and R x is a p × p Hermitian Toeplitz matrix of auto-correlations [14]. A signal- and noise-cancellation model is shown in Figure 2. The desired heartbeat signal to be measured, d ( n ) , is corrupted by various background noise sources, v ( n ) , to produce the signal measured by the HM, x ( n ) = d ( n ) + v 1 ( n ) . Without any information about v 1 ( n ) , it is not possible to remove it from d ( n ) . The BNM measures the noise sources directly, depicted as v 2 ( n ) ; however, the two sensors do not detect the noise the same way v 1 ( n ) v 2 ( n ) , meaning the desired signal cannot be yielded via a direct subtraction [9]. Instead, a Wiener filter produces an estimate of v 1 ( n ) , denoted as v ^ 1 ( n ) , via the observational measurements of v 2 ( n ) . This estimate is subtracted from x ( n ) to attenuate the background noise. Here, the Wiener–Hopf equations are given by Equation (5).
R v 2 w = r x v 2
In environments where the noise source is constantly changing, adaptive filters are needed that update the filter coefficient w in Equation (5) in real time. An adaptive filter model is shown in Figure 3.
Here, the output of the adaptive filter aims to minimise the MSE of estimating v 1 ( n ) ; thus, e ( n ) is the MSE of the desired signal d ( n ) . The optimum filter coefficients that minimise this error are found by solving Equation (5), which can be difficult, as they constantly change for each value of n. To simplify this requirement, the form of Equation (6) can be utilised to update the filter coefficients at each time n.
w n + 1 = w n + Δ w n
Here, the filter coefficient w at time n is corrected by Δ w n to form a new set of coefficients, w n + 1 , at time n + 1 . Defining how this correction is to be applied is at the heart of adaptive algorithms.

2.2. Coherence Function

The coherence between the HM and BNM channels indicates how much noise can be attenuated at particular frequencies and is a measure of the accuracy between the assumed linear channels [12]. As derived in [15], the coherence function, sometimes referred to as the coherence-squared function, may be defined by Equation (7):
γ v x 2 ( f ) = | S v x ( f ) | 2 S v v ( f ) S x x ( f ) ,
where S v x ( f ) is the cross-spectral density, and S v v ( f ) , S x x ( f ) are the auto spectral densities, which are assumed to exist. This is valid for all frequencies f and is bounded by the condition:
0 γ v x 2 ( f ) 1
If we apply this ideology to a noise-free LTI system governed by x ( t ) = v ( t ) h ( t ) , where h ( t ) is the system impulse response, and ∗ denotes the convolution, or X ( f ) = V ( f ) H ( f ) , where H ( f ) is the transfer function, the auto- and cross-spectral densities are governed by Equations (9) and (10):
S x x ( f ) = | H ( f ) | 2 S v v ( f )
S v x = H ( f ) S v v ( f )
Here, γ v x 2 ( f ) = 1 and S v v ( f ) and S x x ( f ) 0 for any frequency. If a heartbeat signal is introduced in x ( t ) depicted as x ( t ) = v ( t ) h ( t ) + d ( t ) , or if the system is non-linear, then γ v x 2 ( f ) 1 . We can model this system as a linear system, shown in Figure 4, to deduce if there is a linear relationship between v ( t ) and x ( t ) .
The true system is modelled by Equation (11), and the model is obtained by minimising the error given by Equation (12).
x ^ ( t ) = v ( t ) h ^ ( t )
e ( t ) = x ( t ) x ^ ( t )
The MMSE method is a standard way to achieve this in the time or frequency domain. The frequency domain MMSE method is as follows.
min H ^ R e e ( f ) = lim T E ^ T X T ( f ) X ^ T ( f ) = R x x ( f ) R x x ^ ( f ) R x ^ x ( f ) + R x ^ x ^ ( f ) = R v v ( f ) | H ( f ) | 2 + R d d ( f ) H ( f ) R v v ( f ) H ^ ( f ) H ^ ( f ) R v v ( f ) H ( f ) + | H ^ ( f ) | 2 R v v ( f )
Here, X T ( f ) and X ^ T ( f ) are the frequency–domain representation of x ( t ) and x ^ ( t ) , respectively. Taking the derivative of Equation (13) with respect to H ^ yields:
H ( f ) R v v + H ^ ( f ) R v v ( f ) = 0 ,
which leads to:
H ^ ( f ) = H ( f ) ,
implying:
R e e ( f ) = R d d ( f ) .
Here, we assumed that E ( d v ) = 0 , resulting in R d v ( f ) = 0 . We see that, for this scenario, the coherence function is given as:
γ v x 2 ( f ) = | H ( f ) | 2 R v v 2 ( f ) R v v ( f ) ( | H ( f ) | 2 R v v ( f ) + R d d ( f ) ) = 1 1 + R d d ( f ) | H ( f ) | 2 R v v ( f ) = 1 1 + R d d ( f ) R x x ( f ) R d d ( f ) .
thus,
1 γ v x 2 ( f ) = R d d R x x ,
which connects the LTI noise-cancellation approach to the coherence function. One of the major problems encountered when applying adaptive noise-cancellation techniques in real acoustic environments was the low coherence between the noise signal corrupting the desired signal and the noise sensor [16]. For the scenario depicted in Figure 2, a coherence value less than unity indicates that the system relating x and v 2 is not linear [15], thus the noise at those particular frequencies cannot be completely removed. The amount of possible noise cancellation for this scenario as a function of frequency, R ( ω ) , is given by Equation (18) and plotted in Figure 5 [16].
R ( ω ) = 1 1 | γ v 2 x 2 ( ω ) |

3. Literature Survey

An adaptive filter has four main aspects, as stated in [17]: the signal, the structure, the parameters, and the adaptive algorithm. The signal is processed by the filter, the structure relates the output signal computed from the input signal, the parameters alter the input–output relationship iteratively, and the adaptive algorithm describes how these parameters are changed over time. Adaptive noise cancellation is an application of adaptive filtering that has been studied vastly in the literature [11], with the theory being explained throughout Section 2. Research in this area dealing with bio-signals has grown over the recent years, in particular due to the filtering of EEG and ECG signals that are corrupted by line interference [18,19]. The following part of this research discusses studies that directly apply to enhancing bio-generated sounds in electronic stethoscopes.
Cancelling environmental noise from cardiac sounds was investigated by Giustina et al. using a multi-channel adaptive algorithm [9]. Sound acquired from two microphones was used to reconstruct the transfer function of a stethoscope head allowing distortionless noise reduction. The first microphone was placed in a pipe directly connected to the stethoscope head, and the second microphone was placed in an adjacent pipe. It was experimentally verified that the inner microphone detected cardiac sounds and noise and the outer microphone detected only noise. A weighted difference (WD) adaptive algorithm was implemented and resulted in significant noise-reduction ratios (NRR) for different scenarios. In particular, a noisy office with background voices, a synthetic chirp and a high-frequency noise obtained NRRs of −14.04, −15.5 and −16.11 dB, respectively.
A similar study was conducted by Patel et al. in [10], where the aim was to attenuate simulated USAF C-130 aircraft noise whilst measuring and processing lung sounds. A two-channel stethoscope was held to a subject’s right anterior upper chest while they breathed through a pneumotachograph. The two microphones, labelled x (background noise pickup microphone) and d (subject microphone), offer the same functionality as seen in [9] and the proposed stethoscope shown in Figure 1. The signals from x and d were amplified and high-pass filtered at a 100 Hz cut-off frequency, then digitized at 5120 samples per second. Some 20-s data segments were recorded during breathing for aircraft noise levels of 80, 90 and 100 dB sound pressure levels (SPL). These epochs were processed with the adaptive-noise-reduction least mean squares (LMS) and normalised least mean squares (NLMS) algorithms using MATLAB. Both of these algorithms apply Equation (6) in different ways and are explained comprehensively in [10,14]. The LMS approach achieved approximately 15 dB of noise reduction in the frequency range 100–600 Hz at 100 dB SPL. The NLMS algorithm provided faster convergence with frequencies above 450 Hz attenuated by an additional 5 dB.

4. System Model

This paper uses a weighted least mean squares (WLMS) algorithm to update the filter coefficients shown in Equation (6). Let the HM signal corrupted with noise be denoted as x and the BNM signal be denoted as y . For each sample from the HM, x ( i ) , the error history is given by Equation (19):
e ( i ) = x ( i ) w T y ( i )
Here, y ( i ) = ( y [ i ] , y [ i + 1 ] , , y [ i + N ] ) T , where N is the filter length (FL). The weights of the filter were updated for each sample as per Equation (20):
w = ( 1 0.001 ) w + μ y ( i ) e ( i ) | | y ( i ) | | 2
Here, the FL determines the size of vectors y and w , and the step size, μ , determines the size at which the filter coefficients are updated. This algorithm is a “variable leaky LMS” algorithm, which has the potential to significantly outperform the standard LMS algorithm [20].

5. Materials and Methods

Four different sources of background noise were tested on heartbeat measurements: a single 300 Hz tone; multiple tones consisting of 200 Hz, 300 Hz and 500 Hz; hospital/clinic noise; and breathing noise. The tones were generated through Audacity, and the hospital/clinic noise was used from YouTube. The heartbeat measurements were of the first author and offer no diagnostic insight; thus, no ethical approval was required for this research. A FireFace UCX was used with a MATLAB interface to allow the simultaneous playback of the background noise through a Fostex 6301B speaker whilst recording from the stethoscope shown in Figure 1. The stethoscope was taped to the first author’s chest, making sure that the BNM was exposed. Figure 6 displays this configuration.
Heartbeat measurements with a duration of 15 s were taken for each background noise scenario. The first author was under breath-held conditions, except for the breathing noise scenario, where the speaker was turned off. The FireFace UCX collected data at 44.1 kHz sampling frequency, which was re-sampled down to 2 kHz. Version R2022a of MATLAB was used to complete this, as well as implement the algorithm from Section 4 on a laptop (HP laptop; Processor—Intel(R) Core(TM) i7-10510U; RAM—16 GB; Operating system—Windows 10; Processor Speed—2.3 GHz).
The FL and step size μ were first tuned on the single 300 Hz tone. The combination that achieved the best attenuated signal was then used for the other background noise scenarios after comparing the performance to a conventional LMS algorithm reviewed in [11]. For each case, the coherence function between the HM and the BNM was plotted via a Welch estimation in MATLAB to indicate the expected noise cancellation performance. Power spectral density (PSD) plots and spectrograms generated through MATLAB allowed visualisation of the performance.

6. Results and Discussion

This section consists of four subsections each containing the results from the four tested background noise sources; a single tone, multiple tones, hospital/clinic noise, and breathing noise. As the HM and BNM were re-sampled down to 2 kHz, the x-axis and y-axis range from 0–1 kHz for coherence function and spectrogram plots, respectively.

6.1. Single 300 Hz Tone

When a single 300 Hz tone was played through the Fostex speaker, the coherence function between the HM and BNM, shown in Figure 7, displayed unity at 300 Hz. Thus, we can expect complete attenuation of the tone using linear techniques.
The WLMS algorithm was tested for FL = 256, 512 and 1024. The PSD of the filtered HM signals plotted with the unfiltered HM signal is displayed in Figure 8.
For each scenario, the amount of attenuation remained constant at 27 dB. A FL of 512 was selected for tuning the step size value, μ , of the WLMS algorithm. Values of 0.01, 0.05, 0.1 and 0.5 were tested, with the PSDs compared once again in Figure 9.
The attenuation increased as μ increased, where almost complete attenuation (approx. 35 dB) was achieved at μ = 0.5 . However, we noted that the low frequency content was affected, evident in Figure 9a, which must be investigated further. Nevertheless, this agrees with the coherence function that suggests that complete attenuation is possible. A 300 Hz tone is a periodic signal, suggesting why the measured attenuation was high. This may not be the case when other non-periodic noise sources are present.
The performance using a conventional LMS algorithm is shown in Figure 10 based on the highest-performing parameters achieved in the WLMS algorithm. The filter length was kept at 512, and the algorithm tested for μ values of 0.1 and 0.5. The maximum attenuation achieved was 32 dB, when μ was 0.5, and was not fully attenuated.
The LMS algorithm under-performed our proposed WLMS algorithm, resulting in lower signal integrity that will negatively affect the classification of CAD. Thus, the following results in this section use the WLMS algorithm with FL = 512 and μ = 0.1 , which was tuned in this subsection to achieve high attenuation without affecting the lower-frequency content.

6.2. Multiple Tones

The 200 Hz, 300 Hz and 500 Hz tones were played simultaneously through the speaker, and the coherence between the HM and BNM is displayed in Figure 11.
The three tones had unity magnitude-squared coherence, as well as the first harmonic of each tone. This suggests that the first two harmonics can be fully attenuated. The WLMS algorithm was tested on the scenario where each tone was played for 5 s while a heartbeat was recorded, so as to test the ability of the algorithm to adapt to different tones. Near-complete attenuation was observed in the PSD comparison of the filtered and unfiltered HM signal shown in Figure 12a, which is visually seen in the spectrogram shown in Figure 12b. Here, the BNM, HM and filtered HM spectrograms are displayed side by side, where the 200 Hz, 300 Hz and 500 Hz tones were approximately attenuated by 24.5 dB, 21.4 dB and 20.3 dB, respectively.
The 200 Hz tone lies within the signal of interest. It did not affect the heartbeat signal when attenuated, highlighting the strength of the adaptive filtering. The BNM and HM also support the hypothesis of the stethoscope design, where the BNM detects the background noise, and the HM detects the heartbeat plus the background noise. This is essential for adaptive noise cancellation, so as to not affect the signal of interest. The HM did not detect the higher-order harmonics to the same degree as the tones, suggesting the reason unity coherence was not seen in Figure 11.

6.3. Hospital/Clinic Background Noise

This type of noise can be expected when measuring heartbeats in practical situations, such as in hospitals or clinics. Thus, it is of high interest to attenuate this type of non-stationary noise within the frequency of interest. The magnitude-squared coherence between the HM and BNM, shown in Figure 13a, is not in unity for any frequency band. Between 200–500 Hz, it varies from 0.3 to 0.7. Thus, adaptive filtering will not result in the complete attenuation of the background noise. Around 3–5 dB will be possible at a magnitude-squared coherence value of 0.5. There is visual evidence of attenuation displayed in the spectrogram comparison in Figure 13b. The energy at 500 Hz and the band between 200 and 300 Hz is slightly suppressed. To increase the attenuation, applying the WLMS algorithm across separate frequency bands can be explored. A re-plot of the coherence function for each frequency band will provide a better indication of achievable suppression.

6.4. Breathing Noise

Many studies in the literature on the detection of CAD ensures that patients hold their breath whilst heartbeat measurements are taken [21,22,23]. Although this is ideal for CAD detection, extending the stethoscope system to detect lung disease could be performed in unison, where the breathing noise can contain critical information. Thus, it is of great interest to separate the heartbeat sound from the breathing sound, so they can be studied and processed separately. For this test, the first author breathed with no external background noise. Across the whole frequency band, the magnitude-squared coherence function was less than 0.3, as seen in Figure 14a. This suggests that less than 1 dB of breathing noise suppression is achievable. This was confirmed in Figure 14b, where the HM signal before and after filtering appear identical. Non-linear methods may be explored to suppress the breathing noise and enhance the heart signal. These other methods should be carefully designed, so as to not suppress murmurs that are present in CAD patients, as the diagnostic process will be negatively affected.
The spectrogram comparison between the BNM and HM is worth noting. The BNM has observed power in the phase of the recording corresponding to exhaling. This is not seen in the HM, indicating that both microphones detect the breathing sounds differently. From the HM perspective, the sound generated by the air moves through body tissue and bone from inside the body and is coupled through the diaphragm. In contrast, breathing sounds from one’s mouth and nose travel through the air before being detected by the BNM. Both of these pathways have different acoustic properties, suggesting why the breathing appears differently on each spectrogram, indicating its difficulty to suppress or separate.

6.5. Summary of Results

A summary of our results and previously reported results from Section 3 is presented in Table 1.

7. Conclusions

In this paper, the performance of the WLMS algorithm to suppress background-noise-corrupting heartbeat measurements was predicted using the coherence function and tested. When periodic signals, such as the single and multiple tones, were played as background noise, the magnitude-squared coherence between the HM and BNM was unity. Thus, we expected that full suppression was attainable and saw this by adjusting the FL and step size, μ , of the WLMS algorithm. Confidence in using this algorithm increased after it outperformed a conventional LMS algorithm on the single-tone background-noise scenario. These parameters were used for the following background noise sources and breathing noise. The magnitude-squared coherence in the hospital/clinic background noise scenario suggested that 3–5 dB noise suppression was achievable using these linear techniques, and this was evident in the spectrogram comparison. For the breathing-noise scenario, the low magnitude-squared coherence resulted in no visible suppression on the spectrogram comparison. To achieve higher attenuation in scenarios where there is low magnitude-squared coherence between the HM and BNM, linear techniques, such as the WLMS algorithm presented in this papers system model, could be applied to separate frequency bands. If this does not achieve desirable suppression, non-linear techniques can be explored.
Our tests produced comparable performance to previously reported data presented in Section 3 and Table 1. Our highest achieved attenuation arose from the scenarios when tones were the background noise source under test. The synthetic chirp and high-frequency noise produced the highest attenuation reported in [9], and the NLMS algorithm led to 20 dB attenuation for 100 dB SPL aircraft noise reported in [10]. Similar to the periodic tones, these sources are stationary, and, although coherence function analysis was not included in these studies, we can assume that there was a high magnitude-squared coherence between both microphones in each case to explain the results. The scenario in [9], where background voices were attenuated by 14.04 dB, is encouraging towards our future studies, as this was considerably higher than the hospital/clinic background noise attenuation that we measured. These are non-stationary sources, indicating our algorithm must be tuned to achieve comparable results.

Author Contributions

Conceptualization, M.F.; data curation, M.F.; formal analysis, M.F.; investigation, M.F.; methodology, M.F.; software, S.N.; supervision, S.N. and Y.R.; writing—original draft, M.F.; writing—review and editing, Y.R. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

We would like to thank Ticking Heart Pty Ltd. for use of their stethoscope design and Curtin University for approving this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Karnath, B.; Thornton, W. Auscultation of the heart. Hosp. Physician 2002, 38, 39–45. [Google Scholar]
  2. Mahnke, C.B.; Nowalk, A.; Hofkosh, D.; Zuberbuhler, J.R.; Law, Y.M. Comparison of two educational interventions on pediatric resident auscultation skills. Pediatrics 2004, 113, 1331–1335. [Google Scholar] [CrossRef]
  3. Roy, D.; Sargeant, J.; Gray, J.; Hoyt, B.; Allen, M.; Fleming, M. Helping family physicians improve their cardiac auscultation skills with an interactive CD-ROM. J. Contin. Educ. Health Prof. 2002, 22, 152–159. [Google Scholar] [CrossRef] [PubMed]
  4. Roy, D.L. The paediatrician and cardiac auscultation. Paediatr. Child Health 2003, 8, 561–563. [Google Scholar] [CrossRef]
  5. WHO. Cardiovascular Diseases (CVDs); WHO: Geneva, Switzerland, 2021.
  6. Raghu, A.; Praveen, D.; Peiris, D.; Tarassenko, L.; Clifford, G. Engineering a mobile health tool for resource-poor settings to assess and manage cardiovascular disease risk: SMARThealth study. BMC Med. Inform. Decis. Mak. 2015, 15, 36. [Google Scholar] [CrossRef]
  7. Clifford, G.D.; Liu, C.; Moody, B.; Springer, D.; Silva, I.; Li, Q.; Mark, R.G. Classification of normal/abnormal heart sound recordings: The PhysioNet/Computing in Cardiology Challenge 2016. In Proceedings of the 2016 Computing in Cardiology Conference (CinC), Vancouver, BC, Canada, 11–14 September 2016; IEEE: Piscataway, NJ, USA, 2016; pp. 609–612. [Google Scholar]
  8. Kim, P.U.; Lee, Y.; Cho, J.H.; Kim, M.N. Modified adaptive noise canceller with an electrocardiogram to enhance heart sounds in the auscultation sounds. Biomed. Eng. Lett. 2011, 1, 194–198. [Google Scholar] [CrossRef]
  9. Della Giustina, D.; Riva, M.; Belloni, F.; Malcangi, M. Embedding a multichannel environmental noise cancellation algorithm into an electronic stethoscope. Int. J. Circuits Syst. Signal Process. 2011, 5, 184–191. [Google Scholar]
  10. Patel, S.B.; Callahan, T.F.; Callahan, M.G.; Jones, J.T.; Graber, G.P.; Foster, K.S.; Glifort, K.; Wodicka, G.R. An adaptive noise reduction stethoscope for auscultation in high noise environments. J. Acoust. Soc. Am. 1998, 103, 2483–2491. [Google Scholar] [CrossRef]
  11. Dixit, S.; Nagaria, D. LMS adaptive filters for noise cancellation: A review. Int. J. Electr. Comput. Eng. 2017, 7, 2520. [Google Scholar] [CrossRef]
  12. Fu, Z.F.; He, J. Modal Analysis; Elsevier: Amsterdam, The Netherlands, 2001. [Google Scholar]
  13. Schmidt, S.E.; Holst-Hansen, C.; Hansen, J.; Toft, E.; Struijk, J.J. Acoustic features for the identification of coronary artery disease. IEEE Trans. Biomed. Eng. 2015, 62, 2611–2619. [Google Scholar] [CrossRef] [PubMed]
  14. Hayes, M.H. Statistical Digital Signal Processing and Modeling; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  15. Bendat, J.S.; Piersol, A.G. Random Data: Analysis and Measurement Procedures; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  16. Elko, G.W. Adaptive noise cancellation with directional microphones. In Proceedings of the 1997 Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, NY, USA, 19–22 October 1997; IEEE: Piscataway, NJ, USA, 1997; p. 4. [Google Scholar]
  17. Douglas, S.C. Introduction to adaptive filters. In Digital Signal Processing Handbook; Madisetti, V.K., Williams, D.B., Eds.; CRC Press LLC: Boca Raton, FL, USA, 1999. [Google Scholar]
  18. Correa, A.G.; Laciar, E.; Patiño, H.; Valentinuzzi, M. Artifact removal from EEG signals using adaptive filters in cascade. J. Phys. Conf. Ser. 2007, 90, 012081. [Google Scholar] [CrossRef]
  19. Tudosa, I.; Adochiei, N. LMS algorithm derivatives used in real-time filtering of ECG signals: A study case on performance evaluation. In Proceedings of the 2012 International Conference and Exposition on Electrical and Power Engineering, Iasi, Romania, 25–27 October 2012; IEEE: Piscataway, NJ, USA, 2012; pp. 565–570. [Google Scholar]
  20. Kamenetsky, M.; Widrow, B. A variable leaky LMS adaptive algorithm. In Proceedings of the Conference Record of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 7–10 November 2004; IEEE: Piscataway, NJ, USA, 2004; Volume 1, pp. 125–128. [Google Scholar]
  21. Samanta, P.; Pathak, A.; Mandana, K.; Saha, G. Classification of coronary artery diseased and normal subjects using multi-channel phonocardiogram signal. Biocybern. Biomed. Eng. 2019, 39, 426–443. [Google Scholar] [CrossRef]
  22. Pathak, A.; Samanta, P.; Mandana, K.; Saha, G. An improved method to detect coronary artery disease using phonocardiogram signals in noisy environment. Appl. Acoust. 2020, 164, 107242. [Google Scholar] [CrossRef]
  23. Pathak, A.; Samanta, P.; Mandana, K.; Saha, G. Detection of coronary artery atherosclerotic disease using novel features from synchrosqueezing transform of phonocardiogram. Biomed. Signal Process. Control 2020, 62, 102055. [Google Scholar]
Figure 1. Single stethoscope. HM facing upward (left) and BNM facing upward (right).
Figure 1. Single stethoscope. HM facing upward (left) and BNM facing upward (right).
Sensors 22 06591 g001
Figure 2. Using a secondary sensor for noise cancellation. Adapted from [14].
Figure 2. Using a secondary sensor for noise cancellation. Adapted from [14].
Sensors 22 06591 g002
Figure 3. Adaptive noise canceller. Adapted from [14].
Figure 3. Adaptive noise canceller. Adapted from [14].
Sensors 22 06591 g003
Figure 4. Linear system model.
Figure 4. Linear system model.
Sensors 22 06591 g004
Figure 5. Adaptive cancellation vs. the squared coherence between v 2 and x.
Figure 5. Adaptive cancellation vs. the squared coherence between v 2 and x.
Sensors 22 06591 g005
Figure 6. (a) Stethoscope taped to chest recording heartbeat with BNM exposed. (b) Background noise played through Fostex 6301B Speaker.
Figure 6. (a) Stethoscope taped to chest recording heartbeat with BNM exposed. (b) Background noise played through Fostex 6301B Speaker.
Sensors 22 06591 g006
Figure 7. The 300 Hz tone background noise coherence function via Welch.
Figure 7. The 300 Hz tone background noise coherence function via Welch.
Sensors 22 06591 g007
Figure 8. PSD comparison of different filter lengths with 300 Hz tone background noise and μ = 0.5 . (a) Comparison from 0–400 Hz. (b) Comparison from 260–340 Hz.
Figure 8. PSD comparison of different filter lengths with 300 Hz tone background noise and μ = 0.5 . (a) Comparison from 0–400 Hz. (b) Comparison from 260–340 Hz.
Sensors 22 06591 g008
Figure 9. PSD comparison of different μ values with 300 Hz tone background noise and FL = 512 (a) Comparison from 0–400 Hz. (b) Comparison from 260–340 Hz.
Figure 9. PSD comparison of different μ values with 300 Hz tone background noise and FL = 512 (a) Comparison from 0–400 Hz. (b) Comparison from 260–340 Hz.
Sensors 22 06591 g009
Figure 10. PSD comparison for different μ vales with 300 Hz tone background noise using a conventional LMS algorithm. FL = 512.
Figure 10. PSD comparison for different μ vales with 300 Hz tone background noise using a conventional LMS algorithm. FL = 512.
Sensors 22 06591 g010
Figure 11. The 200, 300, and 500 Hz tone background noise coherence function via Welch.
Figure 11. The 200, 300, and 500 Hz tone background noise coherence function via Welch.
Sensors 22 06591 g011
Figure 12. WLMS algorithm on changing tones. FL = 512, μ = 0.1 (a) PSD Comparison. (b) Spectrogram comparison.
Figure 12. WLMS algorithm on changing tones. FL = 512, μ = 0.1 (a) PSD Comparison. (b) Spectrogram comparison.
Sensors 22 06591 g012
Figure 13. WLMS algorithm on hospital/clinic background noise. (a) Coherence function via Welch (b) Spectrogram comparison.
Figure 13. WLMS algorithm on hospital/clinic background noise. (a) Coherence function via Welch (b) Spectrogram comparison.
Sensors 22 06591 g013
Figure 14. (a) Breathing noise coherence function via Welch; (b) Spectrogram comparison.
Figure 14. (a) Breathing noise coherence function via Welch; (b) Spectrogram comparison.
Sensors 22 06591 g014
Table 1. Comparison of our results to other studies.
Table 1. Comparison of our results to other studies.
Noise under TestStudyAlgorithmMeasurementMagnitude Squared CoherenceAttenuation [dB]
Single 300 Hz ToneUsWLMSHeartbeat135 (full)
Single 300 Hz ToneUsLMSHeartbeat132
Multiple Tones:
200 Hz
UsWLMSHeartbeat124.5
Multiple Tones:
300 Hz
UsWLMSHeartbeat121.4
Multiple Tones:
500 Hz
UsWLMSHeartbeat120.3
Hospital/Clinic Background Noise:
200–500 Hz
UsWLMSHeartbeat0.3–0.7<2
Breathing NoiseUsWLMSHeartbeat<0.30
Background voices[9]WDHeartbeatN/A14.04
Synthetic chirp[9]WDHeartbeatN/A15.5
High Frequency Noise[9]WDHeartbeatN/A16.11
100 dB SPL aircraft noise: 100–600 Hz[10]LMSLung SoundsN/A15
100 dB SPL aircraft noise: 450–600 Hz[10]NLMSLung SoundsN/A20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fynn, M.; Nordholm, S.; Rong, Y. Coherence Function and Adaptive Noise Cancellation Performance of an Acoustic Sensor System for Use in Detecting Coronary Artery Disease. Sensors 2022, 22, 6591. https://doi.org/10.3390/s22176591

AMA Style

Fynn M, Nordholm S, Rong Y. Coherence Function and Adaptive Noise Cancellation Performance of an Acoustic Sensor System for Use in Detecting Coronary Artery Disease. Sensors. 2022; 22(17):6591. https://doi.org/10.3390/s22176591

Chicago/Turabian Style

Fynn, Matthew, Sven Nordholm, and Yue Rong. 2022. "Coherence Function and Adaptive Noise Cancellation Performance of an Acoustic Sensor System for Use in Detecting Coronary Artery Disease" Sensors 22, no. 17: 6591. https://doi.org/10.3390/s22176591

APA Style

Fynn, M., Nordholm, S., & Rong, Y. (2022). Coherence Function and Adaptive Noise Cancellation Performance of an Acoustic Sensor System for Use in Detecting Coronary Artery Disease. Sensors, 22(17), 6591. https://doi.org/10.3390/s22176591

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop