The effect of traditional Persian music on the cardiac functioning of young Iranian women.

In the past few decades, several studies have reported the physiological effects of listening to music. The physiological effects of different music types on different people are not similar. Therefore, in the present study, we have sought to examine the effects of traditional Persian music on the cardiac function in young women. Twenty-two healthy females participated in this study. ECG signals were recorded in two conditions: rest and music. For each of the 21 ECG signals (15 morphological and six wavelet based feature) features were extracted. SVM classifier was used for the classification of ECG signals during and before the music. The results showed that the mean of heart rate, the mean amplitude of R-wave, T-wave, and P-wave decreased in response to music. Time-frequency analysis revealed that the mean of the absolute values of the detail coefficients at higher scales increased during rest. The overall accuracy of 91.6% was achieved using polynomial kernel and RBF kernel. Using linear kernel, the best result (with the accuracy rate of 100%) was attained.

Introduction

The ECG is a biometric signal, which reflects the electrical activity of the heart muscle. An analysis of the ECG signal provides useful information about the heart’s performance. The heart’s electrical activity during normal heartbeat is characterized by five basic waves P, Q, R, S, and T waves and sometimes U wave. The P wave and T wave represent atrial depolarization and ventricle repolarization, respectively. In addition, the QRS wave represents ventricular depolarization. The intervals and amplitudes defined by ECG features provide useful clinical information.

It has been previously established that ECG signals can be affected by different factors such as age, cardiovascular diseases, and mental stress. In addition, in the past few decades, several studies in music, psychology, and medicine have reported the physiological effect of listening to music.3, 4 Loomba et al. claimed that music could decrease systolic blood pressure, diastolic blood pressure, and heart rate, significantly. This may be attributed to brain and autonomic nervous system (ANS) activity, but the relationship between music features and ANS activity is not completely understood.6, [7]

These effects have been examined in response to different kinds of music, such as sedative music and excitative music. Since the 1918s, studies were started in which the quantitative differences of cardiac autonomic function during music and rest were extracted. However, the effect of music on cardiac function have varied from research to research. Rickard et al. studied the physiological responses of sedative music. They showed that heart rate was reduced by sedative music. In other studies by Davis et al. and Vanderark et al. it has been shown that music made no changes to the heart rate. Therefore, it can be concluded that the physiological effect of different music types on different people are not similar. These different effects can be related to method, findings, music types, and gender.

Dousty et al. investigated the relationship between music types and cardiovascular function. They focused on two types of music, namely, sedative and arousal music. They claimed that the heart responded differently when different kinds of music were played. Khazaei et al. represented that music can be used as a powerful method for reducing stress and anxiety in generalized anxiety disorder patients. Other physiological effects of music have been summarized in Table 1.

Table 1

Physiologic effects of music.

Study	Number of participants	Music type	Result(s)
Salimpoor et al. 14	26	1.Self-selected intensely pleasurable music2.Neutral music	There is a strong positive correlation between ratings of pleasure and emotional arousal.
Orini et al. 6	75	1.Pleasant music2.Sequences of Shepard tones3.Unpleasant sounds	Pleasant music increased the heart and respiratory rates.
Nakahara et al. 15	13	1.Performance of Music2.Perception of Music	Musical performance would lead to a greater effect of emotion-related modulation in cardiac autonomic nerve activity than musical perception.
Kemper et al. 16	81	Designed music	The HRV parasympathetic parameter was significantly lower with music than rest.
Uggla et al. 17	24		Music significantly lowered the HRV and stress levels of children undergoing haematopoietic stem cell transplants.
Archana and Mukilan 18	30	Preferential music	Listening to pleasant music could be an effective method of relaxation, and it can shift the autonomic balance towards the parasympathetic activity.
Oh et al. 19	60		Listening to pleasant music could be an effective nursing intervention for alleviating anxiety during non-stress test.

Most of the previous studies focused on statistical measures, such as mean and variances and tried to assess the physiological effects of music. It has been previously established that frequency domain measures of ECG signals can be significantly affected by music. One appropriate tool for studying of these differences is Discrete Fourier Transform (DFT) and extraction of features from ECG signal power spectrum. An important drawback in the Fourier analysis is that the technique does not provide any information regarding the exact location of frequency components in time. Short Time Fourier (STF) analysis has been proposed for overcoming this limitation by windowing the area of interest. However, the main weakness of this technique is that its time frequency precision is not optimal. Among the different time-frequency transformations, Discrete Wavelet Transform (DWT) is used in the current study in order to investigate the physiological effect of listening to music. DWT is an effective tool for analysis of non-stationary signals like ECG. Furthermore, it has been established that applying DWT to ECG signals can be helpful in detecting clinically significant features that may be missed by other analysis techniques.22, 23, 24 In our previous study, we studied the effect of gender differences on electrical function of the heart in response to traditional Persian music. It has been shown that the mean heart rate signals in men increased during music; however, it decreased in the same protocol in the women’s group. We also reported that the heart rate signals of women's group represented maximum power spectrum, twice more than that of the men's group. In addition, it seemed the physiological effect of Persian music was higher in the women's group than the men's group. Therefore, the present study focused on other techniques for studying the ECGs of the female group during music and rest. The present work is organized as follows: In the second part, the ECG dataset, which was collected from a group of young female students before and during listening music, are briefly described. In the third part, the wavelet co-efficient based features and morphological features were developed. Next, the results and comparison between participants’ heart rate during music and rest is presented. Finally, discussion and conclusion are presented. A block diagram of the proposed method is shown in Fig. 1.

Fig. 1

Block diagram of the proposed method.

Methods

Data collection

Twenty two healthy females participated in this study. The age range of the subjects was 20–24 years. All the participants were students of Sahand University of Technology and had no previous history of neurological diseases. ECG signals were recorded under two different conditions from each subject. Initially, the participants were asked to lie in a supine position comfortably and close their eyes for five minutes. Then, five minutes of traditional Persian music was played for participants at a comfortable volume. The ECGs – lead II – were recorded by 16-channel PowerLab (ADInstruments, UK) at a sampling rate of 400 Hz. A digital notch filter was used to remove the 50 Hz Power-Line Interference from the ECG signal. The signals were recorded in MAT, which is a common MATLAB file type.

Discrete Wavelet Transform

DWT 25, 26 is a powerful technique, which can be used to analyse different kinds of biomedical signals like ECG. DWT uses a low pass filter (LPF) and a high pass filter (HPF) to decompose the signal into a number of scales related to frequency components and analyses each scale with a certain resolution. The output coefficients of the LPF and HPF are called approximations and details, respectively. Amplitude of the wavelet coefficients represents the energy density of the signal at particular frequencies and moments. The original signal can be reconstructed by using complementary filters. DWT can be defined as follows:where φ represents the mother wavelet, which is dilated by integers m and translated by integers n. DWT can be used to provide salient information about the local features of the ECG signals because it has good localization abilities in both time and frequency domain. Choosing a wavelet family, which closely matches the signal to be processed, plays an important role in wavelet applications. As shown in Fig. 2, the Daubechies wavelet family is similar in shape to the QRS complex.

Fig. 2

Daubechies wavelet family.

In the present study, Daubechies wavelet of order six (db6) with eight levels of wavelet decomposition was used to filter and analyse noisy ECG signals. Fig. 3 represents the basic multiresolution decomposition steps for signal ‘Original signal’.

Fig. 3

The algorithm of the DWT.

Preprocessing

De-noising

The ECG is often contaminated by noises of different kinds. These noises are usually generated by recording instruments or movements of subjects. In the current study, db6 is employed for removing various kinds of noises from ECG signals. The proposed de-noising algorithm can be summarized as follows:

The ECG signals decomposed at level eight using db6.

Soft thresholding technique employed for quantization of detail wavelet coefficients on each level.

ECG signal reconstructed by using approximation wavelet coefficients of last level N and the de-noised detail wavelet coefficients of all levels.

Baseline shifting removing

Baseline wandering is one of the noise artefacts that strongly affect ECG signal analysis. Electrode and respiration impedance changes due to perspiration play an important role in generating baseline wander. By removing baseline wander, we can minimize the changes in beat morphology. Normally, the frequency content of the baseline wander is concentrated in very low frequencies. In the present study, median filters were applied to minimize the influence of baseline drift. The proposed process consists of two steps. In the first step, the original ECG signal is smoothed with a moving average filter of length 150. In the second step, the filtered signal subtracted from the ‘Signal+ Baseline drift noise’ signal and a signal with baseline drift elimination is obtained.

Peak detection

Accurate detection of ECG peaks plays an important role in ECG signal analysis. This is because a lot of clinical information can be derived from the amplitudes and intervals defined by P, Q, R, S, T peaks. In the present study, the DWT method was used to detect the position of the occurrence of the P-QRS-T waves. The proposed algorithm consists of four important steps. In the first step, the ECG signals decomposed at level 8 using db6. In the second step, details between levels 3 and 5 were kept and all the rest removed. The attained signal samples were then squared to stress the signal. Then, R wave was detected using automatic thresholding techniques. In the second step, all the details down to level 5 were removed. Then the ECG signal was reconstructed. The reconstructed ECG signal and minimum signal strength about each R peak within 0.1 s registered as Q and S waves. In the third step, details between levels 1 and 5 were kept and all the rest removed and two zero crossing points of the signal were determined. Finally, for the detection of the P and T waves, only details between levels 4 and 8 were kept and all the rest were removed, Fig. 4.

Fig. 4

Detection procedure. a) R wave detection b) Q & S Waves detection. c) Zero level detection. d) P & T detection.

Feature extraction

For classifying ECG signals during rest and music, it was necessary to extract appropriate features. There are several ways to extract the feature of ECG signal. The present study focused on two types of features namely ‘Wavelet co-efficient based features’ and ‘Morphological features of ECG signal’. Applying approximation and detail coefficients to the classifier directly has negative effects on SVM operation. Therefore, statistics of wavelet coefficients are used instead of high dimensional wavelet coefficients. The following statistical features were selected to exhibit the general behaviour of time‐frequency energy distribution of the ECG waveforms:

An average of the absolute values of approximation wavelet coefficients at any sub-band.

Standard Deviation (SD) of approximation wavelet coefficients at each level.

Average of the absolute values of detail wavelet coefficients at any sub-band.

SD of detail wavelet coefficients at each level.

Besides the wavelet co-efficient based features, morphological features were also extracted from ECG signals. The selected morphological features were mean and SD values of RR interval, TT interval, PP interval, PR interval, TR interval, PT interval, and maximum and mean values of P, Q, R, S, T peaks.

Support vector machine

SVMs are supervised machine learning methods suggested by Cortes and Vapnik in 1995. SVM aims to find the best hyperplane that separates all data points of one class from those of the other classes. This is equivalent to maximizing the margin between the two classes. Consider a linearly separable sample set:where is called an input and represents the class labels. The equation of a separating line can be written as follows:where , represents the inner product of w and x, and b is bias. Therefore, the problem of best separating line is given as follows:

Then the distance from closet sample to separating line is . Thus the margin is . Maximizing the margin is equivalent to minimizing, and then the problem of finding the best separating line is converted into the following optimization problem:

The sample set can be divided into two classes by solving the above equation. While data is most likely to be not entirely separable, the slack variable and penalty parameter introduced in the formula weaken the constraint to a certain degree. Thus, the basic problem of SVM is transformed into:

C is a constant minimizing the errors by maximizing the margin. Using the Lagrangian method to solve the above equation, it is transformed to its dual problem:

Then the classifier is,where

αi represents Lagrange multipliers. In the case of nonlinear classification, by replacing the inner product term in (8) with a proper kernel; the input data is projected onto a higher dimensional space. Then a separating hyperplane is constructed to separate the data points of one class from those of the other classes.

Results

ECG signals recorded from 22 young female participants were analysed utilising the methodology described in Section 2. Following the pre-processing stage, ECG signals were decomposed into wavelet coefficients. Fig. 5 shows the mean of the absolute values of the detail and approximation wavelet coefficients of the ECG signal during rest and music.

Fig. 5

Mean of the absolute values of the detail and approximation wavelet coefficients.

Apart from wavelet co-efficient based features, the morphological feature of ECG signal is also extracted. Fig. 6 summarizes the SD and mean values of morphological features.

Fig. 6

SD and mean values of morphological features.

Good understanding of amplitudes defined by ECG peaks can provide useful clinical information. Fig. 7 shows the maximum, minimum, and mean values of R-wave amplitude, T-wave amplitude, and P-wave amplitude respectively.

Fig. 7

The maximum, minimum and mean values of P, R, T peaks.

Using all of the extracted features is not suitable for classification of ECG signals during rest and music. In the present work, the paired t-test was used for analysing the significance in extracted features during rest and music using SPSS software. P-value less than 0.05 was considered significant. The selected features are given in Table 2.

Table 2

Selected features for classification ECG signals during rest and music.

Feature type	Selected attributes
Wavelet based features	Details coefficients up to D3. (3 features) (p < 0.05, confidence interval: 95%)Approximates coefficients up to A3. (3 features) (p < 0.05, confidence interval: 95%)
Morphological features	SD values of the RR, TT, PP, RP, TR AND PT intervals. (6 features) (p < 0.05, confidence interval: 95%)
	Mean values of the RR, TT, PP, RP, TR AND PT intervals. (6 features) (p < 0.05, confidence interval: 95%)
	Mean values of P, R, T peaks. (3 features) (p < 0.05, confidence interval: 95%)

The extracted features were stored in a single feature vector. These feature vectors were used as input for the SVM classifiers. The linear kernel, polynomial kernel, and RBF kernel are used in SVMs. Seventy percent of feature vectors were selected randomly as training data set and the remaining vectors were used for validation of the classifiers. This process is repeated 10 times. The performance of the proposed algorithm is measured in terms of accuracy. We have achieved an overall accuracy of 91.6% using polynomial kernel and RBF kernel whereas 100% accuracy was achieved using linear kernel. These results confirm the physiological effects of music on cardiovascular function. The performances of different kernels have been shown in Table 3.

Table 3

Performance of different kernels.

Kernel	TP	FN	TN	FP	Accuracy(%)	Sensitivity(%)	Specificity(%)	AUC
RBF	20	2	21	1	93.18	90.90	95.45	0.9318
Linear	21	1	22	0	97.72	95.45	100	0.9773
MLP	12	10	12	10	54.54	54.54	54.54	0.5455
Polynomial	21	1	21	1	95.45	95.45	95.45	0.9545

Discussion

The ECG is a graphical representation of the heart electrical activity which is widely used for diagnosing cardiovascular diseases. ECG signal can be affected by physiological state changes and emotional factors such as listening to music. The present study aimed to examine the effects of music on cardiac function in young women. Therefore, 10 min of ECG is recorded from each of the normal participants during rest and music.

Direct visual monitoring of ECGs by human personnel is a time consuming process and usually involves a loss of information. To handle this problem, computer-based systems have been developed for detecting ECG characteristic points. Due to time-varying morphology of the QRS wave, and the presence of severe baseline drift and other high frequency noises in the ECG signal, an accurate detection of P wave, QRS complex, and R wave, in an ECG signal is a challenging issue. As ECG signal has semi-periodic patterns, an appropriate tool for detecting P-Q-R-S-T waves is DWT and extraction of wavelet co-efficient based features. As described in Section 2, in the present study the Daubechies wavelet of order 6 were used to distinguish ECG waves from noises and baseline drift.

Feature extraction is the most important step in pattern recognition. There are several ways to extract the feature of ECG signal. In this work, there are two types of extracted features of ECG waveforms. The present study focused on two types of features named ‘Wavelet co-efficient based features’ and ‘Morphological features of ECG signal’. Morphological characteristics of ECGs are well-known features that are widely used in making clinical decisions and in medical research. Therefore, these features are used in the present study. In addition, the ability to represent data in compressed parameters form (features) is one of the main application of DWT.

In order to overcome SVM over fitting, two important steps need to be taken: (1) K-fold cross validation; and (2) parameter optimization. Variations on cross-validation include leave-k-out cross-validation and k-fold cross-validation. Data sample is small in our study; therefore, leave-k-out cross-validation cannot be used in the present study. Instead, we used k-fold cross-validation. Seventy percent of feature vectors were selected randomly as training data set and the remaining vectors were used for validation of the classifiers. This process was repeated 10 times and an average of the resulting measurements were taken.35, 36, 37 In addition, it has been suggested that, to avoid the risk of SVM over fitting due to the small training sets, parameter optimization can be omitted. Therefore, we omitted parameter optimization. In addition, similar to previous studies, we used C = 1000; γ = 0.1 for all experiments and before the classification step, the features are linearly normalized to [0,1].

As shown in Fig. 7, mean R-wave amplitude, mean T-wave amplitude, and mean P-wave amplitude decreased in response to music. These results are in the line with Dousty et al. who found mean R-wave amplitude and T-wave amplitude reduced by sedative music.

The mean heart rate signals before listening to music is about 78.53, whereas this value is about 77.28 during the music. These results are in the line with Goshvarpour et al. who found mean heart rate signals in the women’s group decreased during music.

Apart from the morphological features, the wavelet coefficient based features are also extracted from ECG signals. The approximations are the low-frequency components of the signal. The results suggested that in both conditions (before and during the music) the energy spectrums are concentrated in low frequencies. In addition, the behaviours of detail and approximate coefficients at low scales are almost similar during rest and music. In contrast, the general appearance of detail and approximate coefficients at high scales are significantly different during rest and music (p < 0.001). As shown in Fig. 5, during rest the mean of the absolute values of the detail coefficients at higher scales show amplitude even more than two times that of music period.

After studying the wavelet co-efficient based features and morphological features, we tried to quantify the difference between women’s ECG signals during and before music on the basis of obtained results. In the present study, both wavelet co-efficient based features and morphological features were utilized because some studies established their capacity in the classification of ECG signals. Therefore, for each of the ECG signals, 21 morphological and wavelet co-efficient based features have been selected accordingly to provide best result. These features have provided significant distinction between participants’ ECG signals during and before music. It seems, this is the main reason for better results of linear kernel in comparison with kernel and RBF kernel. In addition, because of the SVM power in ECG Pattern classification, it is applied for classification of ECG signals during and before the music. The performance of the proposed algorithm is determined in terms of accuracy. We achieved an overall accuracy of 91.6% using polynomial kernel and RBF kernel whereas 100% accuracy is achieved using linear kernel. The results of the present assay confirm the effects of music on cardiac function in young women.

Conclusion

There is a need for more detailed studies to enrich the understanding of the physiological effects of music. However, the results of the present work confirm the effects of music on cardiac function in young women. In the other words, it seems that the heart can synchronize with traditional Persian music in young women.

Conflicts of interest

The authors have none to declare.