Performance Analysis of Fetal-Phonocardiogram Signal Denoising Using The Discrete Wavelet Transform

The obligation for comprehensive fetal heart rate investigation had driven to improve the passive and non-invasive diagnostic instruments despite the USG or CTG method. Fetal phonocardiography (f-PCG) utilizing the auscultation method met the above criteria, but its interpretation frequently disturbed by the presence of noise. For instance, maternal heart and body organ sounds, fetal movements noise, and ambient noise from the environment where it is recording are the noise that corrupted the f-PCG signal. In this work, the use of discrete wavelet transforms (DWT) to eliminate noise in the f-PCG signal with SNR as the performance parameters observed. It was observing the effect of changes in wavelet type and threshold type on the SNR value. The test was carried out on f-PCG data taken from physio.net. Initial SNR values ranged from -26.7 dB to -4.4 dB; after application of DWT procedure to f-PCG, SNR increased significantly. Based on the test results obtained, wavelet type coif1 with the soft threshold gave the best result with 11.69 dB in SNR value. The coif1 had a superior result than other mother wavelets that use in this work, so the fPCG signal analysis for fetal heart rate investigation suggested to use it.The obligation for comprehensive fetal heart rate investigation had driven to improve the passive and non-invasive diagnostic instruments despite the USG or CTG method. Fetal phonocardiography (f-PCG) utilizing the auscultation method met the above criteria, but its interpretation frequently disturbed by the presence of noise. For instance, maternal heart and body organ sounds, fetal movements noise, and ambient noise from the environment where it is recording are the noise that corrupted the f-PCG signal. In this work, the use of discrete wavelet transforms (DWT) to eliminate noise in the f-PCG signal with SNR as the performance parameters observed. It was observing the effect of changes in wavelet type and threshold type on the SNR value. The test was carried out on f-PCG data taken from physio.net. Initial SNR values ranged from -26.7 dB to -4.4 dB; after application of DWT procedure to f-PCG, SNR increased significantly. Based on the test results obtained, wavelet type coif1 with the soft threshold gave the best result with 11.69 dB in SNR value. The coif1 had a superior result than other mother wavelets that use in this work, so the fPCG signal analysis for fetal heart rate investigation suggested to use it

Automatic analysis and classification of cardiac acoustic signals for long term monitoring

Objective: Cardiovascular diseases are the leading cause of death worldwide resulting in over 17.9 million deaths each year. Most of these diseases are preventable and treatable, but their progression and outcomes are significantly more positive with early-stage diagnosis and proper disease management. Among the approaches available to assist with the task of early-stage diagnosis and management of cardiac conditions, automatic analysis of auscultatory recordings is one of the most promising ones, since it could be particularly suitable for ambulatory/wearable monitoring. Thus, proper investigation of abnormalities present in cardiac acoustic signals can provide vital clinical information to assist long term monitoring. Cardiac acoustic signals, however, are very susceptible to noise and artifacts, and their characteristics vary largely with the recording conditions which makes the analysis challenging. Additionally, there are challenges in the steps used for automatic analysis and classification of cardiac acoustic signals. Broadly, these steps are the segmentation, feature extraction and subsequent classification of recorded signals using selected features. This thesis presents approaches using novel features with the aim to assist the automatic early-stage detection of cardiovascular diseases with improved performance, using cardiac acoustic signals collected in real-world conditions. Methods: Cardiac auscultatory recordings were studied to identify potential features to help in the classification of recordings from subjects with and without cardiac diseases. The diseases considered in this study for the identification of the symptoms and characteristics are the valvular heart diseases due to stenosis and regurgitation, atrial fibrillation, and splitting of fundamental heart sounds leading to additional lub/dub sounds in the systole or diastole interval of a cardiac cycle. The localisation of cardiac sounds of interest was performed using an adaptive wavelet-based filtering in combination with the Shannon energy envelope and prior information of fundamental heart sounds. This is a prerequisite step for the feature extraction and subsequent classification of recordings, leading to a more precise diagnosis. Localised segments of S1 and S2 sounds, and artifacts, were used to extract a set of perceptual and statistical features using wavelet transform, homomorphic filtering, Hilbert transform and mel-scale filtering, which were then fed to train an ensemble classifier to interpret S1 and S2 sounds. Once sound peaks of interest were identified, features extracted from these peaks, together with the features used for the identification of S1 and S2 sounds, were used to develop an algorithm to classify recorded signals. Overall, 99 features were extracted and statistically analysed using neighborhood component analysis (NCA) to identify the features which showed the greatest ability in classifying recordings. Selected features were then fed to train an ensemble classifier to classify abnormal recordings, and hyperparameters were optimized to evaluate the performance of the trained classifier. Thus, a machine learning-based approach for the automatic identification and classification of S1 and S2, and normal and abnormal recordings, in real-world noisy recordings using a novel feature set is presented. The validity of the proposed algorithm was tested using acoustic signals recorded in real-world, non-controlled environments at four auscultation sites (aortic valve, tricuspid valve, mitral valve, and pulmonary valve), from the subjects with and without cardiac diseases; together with recordings from the three large public databases. The performance metrics of the methodology in relation to classification accuracy (CA), sensitivity (SE), precision (P+), and F1 score, were evaluated. Results: This thesis proposes four different algorithms to automatically classify fundamental heart sounds – S1 and S2; normal fundamental sounds and abnormal additional lub/dub sounds recordings; normal and abnormal recordings; and recordings with heart valve disorders, namely the mitral stenosis (MS), mitral regurgitation (MR), mitral valve prolapse (MVP), aortic stenosis (AS) and murmurs, using cardiac acoustic signals. The results obtained from these algorithms were as follows: • The algorithm to classify S1 and S2 sounds achieved an average SE of 91.59% and 89.78%, and F1 score of 90.65% and 89.42%, in classifying S1 and S2, respectively. 87 features were extracted and statistically studied to identify the top 14 features which showed the best capabilities in classifying S1 and S2, and artifacts. The analysis showed that the most relevant features were those extracted using Maximum Overlap Discrete Wavelet Transform (MODWT) and Hilbert transform. • The algorithm to classify normal fundamental heart sounds and abnormal additional lub/dub sounds in the systole or diastole intervals of a cardiac cycle, achieved an average SE of 89.15%, P+ of 89.71%, F1 of 89.41%, and CA of 95.11% using the test dataset from the PASCAL database. The top 10 features that achieved the highest weights in classifying these recordings were also identified. • Normal and abnormal classification of recordings using the proposed algorithm achieved a mean CA of 94.172%, and SE of 92.38%, in classifying recordings from the different databases. Among the top 10 acoustic features identified, the deterministic energy of the sound peaks of interest and the instantaneous frequency extracted using the Hilbert Huang-transform, achieved the highest weights. • The machine learning-based approach proposed to classify recordings of heart valve disorders (AS, MS, MR, and MVP) achieved an average CA of 98.26% and SE of 95.83%. 99 acoustic features were extracted and their abilities to differentiate these abnormalities were examined using weights obtained from the neighborhood component analysis (NCA). The top 10 features which showed the greatest abilities in classifying these abnormalities using recordings from the different databases were also identified. The achieved results demonstrate the ability of the algorithms to automatically identify and classify cardiac sounds. This work provides the basis for measurements of many useful clinical attributes of cardiac acoustic signals and can potentially help in monitoring the overall cardiac health for longer duration. The work presented in this thesis is the first-of-its-kind to validate the results using both, normal and pathological cardiac acoustic signals, recorded for a long continuous duration of 5 minutes at four different auscultation sites in non-controlled real-world conditions.Open Acces

Review and classification of variability analysis techniques with clinical applications

Analysis of patterns of variation of time-series, termed variability analysis, represents a rapidly evolving discipline with increasing applications in different fields of science. In medicine and in particular critical care, efforts have focussed on evaluating the clinical utility of variability. However, the growth and complexity of techniques applicable to this field have made interpretation and understanding of variability more challenging. Our objective is to provide an updated review of variability analysis techniques suitable for clinical applications. We review more than 70 variability techniques, providing for each technique a brief description of the underlying theory and assumptions, together with a summary of clinical applications. We propose a revised classification for the domains of variability techniques, which include statistical, geometric, energetic, informational, and invariant. We discuss the process of calculation, often necessitating a mathematical transform of the time-series. Our aims are to summarize a broad literature, promote a shared vocabulary that would improve the exchange of ideas, and the analyses of the results between different studies. We conclude with challenges for the evolving science of variability analysis

Biomedical Signal and Image Processing

Written for senior-level and first year graduate students in biomedical signal and image processing, this book describes fundamental signal and image processing techniques that are used to process biomedical information. The book also discusses application of these techniques in the processing of some of the main biomedical signals and images, such as EEG, ECG, MRI, and CT. New features of this edition include the technical updating of each chapter along with the addition of many more examples, the majority of which are MATLAB based

High-Performance Accelerometer Based On Asymmetric Gapped Cantilevers For Physiological Acoustic Sensing

Continuous or mobile monitoring of physiological sounds is expected to play important role in the emerging mobile healthcare field. Because of the miniature size, low cost, and easy installation, accelerometer is an excellent choice for continuous physiological acoustic signal monitoring. However, in order to capture the detailed information in the physiological signals for clinical diagnostic purpose, there are more demanding requirements on the sensitivity/noise performance of accelerometers. In this thesis, a unique piezoelectric accelerometer based on the asymmetric gapped cantilever which exhibits significantly improved sensitivity is extensively studied. A meso-scale prototype is developed for capturing the high quality cardio and respiratory sounds on healthy people as well as on heart failure patients. A cascaded gapped cantilever based accelerometer is also explored for low frequency vibration sensing applications such as ballistocardiogram monitoring. Finally, to address the power issues of wireless sensors such as wireless wearable health monitors, a wide band vibration energy harvester based on a folded gapped cantilever is developed and demonstrated on a ceiling air condition unit

Signal Processing Using Non-invasive Physiological Sensors

Non-invasive biomedical sensors for monitoring physiological parameters from the human body for potential future therapies and healthcare solutions. Today, a critical factor in providing a cost-effective healthcare system is improving patients' quality of life and mobility, which can be achieved by developing non-invasive sensor systems, which can then be deployed in point of care, used at home or integrated into wearable devices for long-term data collection. Another factor that plays an integral part in a cost-effective healthcare system is the signal processing of the data recorded with non-invasive biomedical sensors. In this book, we aimed to attract researchers who are interested in the application of signal processing methods to different biomedical signals, such as an electroencephalogram (EEG), electromyogram (EMG), functional near-infrared spectroscopy (fNIRS), electrocardiogram (ECG), galvanic skin response, pulse oximetry, photoplethysmogram (PPG), etc. We encouraged new signal processing methods or the use of existing signal processing methods for its novel application in physiological signals to help healthcare providers make better decisions

Symbolic Dynamics Analysis: a new methodology for foetal heart rate variability analysis

Cardiotocography (CTG) is a widespread foetal diagnostic methods. However, it lacks of objectivity and reproducibility since its dependence on observer's expertise. To overcome these limitations, more objective methods for CTG interpretation have been proposed. In particular, many developed techniques aim to assess the foetal heart rate variability (FHRV). Among them, some methodologies from nonlinear systems theory have been applied to the study of FHRV. All the techniques have proved to be helpful in specific cases. Nevertheless, none of them is more reliable than the others. Therefore, an in-depth study is necessary. The aim of this work is to deepen the FHRV analysis through the Symbolic Dynamics Analysis (SDA), a nonlinear technique already successfully employed for HRV analysis. Thanks to its simplicity of interpretation, it could be a useful tool for clinicians. We performed a literature study involving about 200 references on HRV and FHRV analysis; approximately 100 works were focused on non-linear techniques. Then, in order to compare linear and non-linear methods, we carried out a multiparametric study. 580 antepartum recordings of healthy fetuses were examined. Signals were processed using an updated software for CTG analysis and a new developed software for generating simulated CTG traces. Finally, statistical tests and regression analyses were carried out for estimating relationships among extracted indexes and other clinical information. Results confirm that none of the employed techniques is more reliable than the others. Moreover, in agreement with the literature, each analysis should take into account two relevant parameters, the foetal status and the week of gestation. Regarding the SDA, results show its promising capabilities in FHRV analysis. It allows recognizing foetal status, gestation week and global variability of FHR signals, even better than other methods. Nevertheless, further studies, which should involve even pathological cases, are necessary to establish its reliability.La Cardiotocografia (CTG) è una diffusa tecnica di diagnostica fetale. Nonostante ciò, la sua interpretazione soffre di forte variabilità intra- e inter- osservatore. Per superare tali limiti, sono stati proposti più oggettivi metodi di analisi. Particolare attenzione è stata rivolta alla variabilità della frequenza cardiaca fetale (FHRV). Nel presente lavoro abbiamo suddiviso le tecniche di analisi della FHRV in tradizionali, o lineari, e meno convenzionali, o non-lineari. Tutte si sono rivelate efficaci in casi specifici ma nessuna si è dimostrata più utile delle altre. Pertanto, abbiamo ritenuto necessario effettuare un’indagine più dettagliata. In particolare, scopo della tesi è stato approfondire una specifica metodologia non-lineare, la Symbolic Dynamics Analysis (SDA), data la sua notevole semplicità di interpretazione che la renderebbe un potenziale strumento di ausilio all’attività clinica. Sono stati esaminati all’incirca 200 riferimenti bibliografici sull’analisi di HRV e FHRV; di questi, circa 100 articoli specificamente incentrati sulle tecniche non-lineari. E’ stata condotta un’analisi multiparametrica su 580 tracciati CTG di feti sani per confrontare le metodologie adottate. Sono stati realizzati due software, uno per l’analisi dei segnali CTG reali e l’altro per la generazione di tracciati CTG simulati. Infine, sono state effettuate analisi statistiche e di regressione per esaminare le correlazioni tra indici calcolati e parametri di interesse clinico. I risultati dimostrano che nessuno degli indici calcolati risulta più vantaggioso rispetto agli altri. Inoltre, in accordo con la letteratura, lo stato del feto e le settimane di gestazione sono parametri di riferimento da tenere sempre in considerazione per ogni analisi effettuata. Riguardo la SDA, essa risulta utile all’analisi della FHRV, permettendo di distinguere – meglio o al pari di altre tecniche – lo stato del feto, la settimana di gestazione e la variabilità complessiva del segnale. Tuttavia, sono necessari ulteriori studi, che includano anche casi di feti patologici, per confermare queste evidenze

Advanced robust non-invasive foetal heart detection techniques during active labour using one pair of transabdominal electrodes

The thesis proposes and evaluates three state-of-the-art signal processing techniques to detect fetal heartbeats within each maternal cardiac cycle, during labour contractions, using only a pair of transabdominal electrodes. The first and second techniques are, namely, the structured third- order cumulant-slice-template matching and the bispectral-contours-template matching for fetal QRS identification, respectively. The third technique is based on the modified and appropriately weighted spectral multiple signal classification (MUSIC) with incorporated covariance matrix for uterine contraction noise-like interfering signals also contaminated with noise. Essentially, two modifications to the standard MUSIC have been developed in order to enhance the performance of the spectral estimator in our applied work. The first modification involves the introduction of an optimised weighting function to the segmented ECG covariance matrix, and is chiefly aimed at enhancing the fetal QRS major spectral peak which occurs at around 30 Hz against the mother QRS major spectral peak usually occurring around 17 Hz and all other noise contributions. Additional optional pseudo-bispectral enhancement to sharpen the maternal and fetal spectral peaks, in particular when the mother and fetal R-waves are temporally coincident, have been achieved. The second modification to the spectral MUSIC is the removal of the unjustified assumption that only white Gaussian noise is present and the incorporation of the actual measured labour uterine contraction covariance matrix in reconfigured subspace analysis. This inevitably leads to the generalised eigenvectors - eigenvalues decomposition modern signal processing. This is now coined the modified, interference incorporated pseudo-spectral MUSIC. The above mentioned first and second techniques are higher-order statistics-based (HOS) and hybrid involving both signal processing and NN classifiers. The third technique is second-order statistics-based (SOS). In all techniques, the removal of signal non-linearity with the aid of non-linear Volterra synthesisers plays a crucial part in the fetal detection integrity. Accurately assessed fetal heart classification rates as high as 95% have been achieved during labour, thus helping to provide non-invasive transparency to fetal intrapartum welfare. Performance analysis and evaluation processes involved more than 30 critical cases classified as “fetal under stress in labour” recorded in a London hospital database and used both transbadominal ECG electrodes and fetal scalp electrodes. The latter facilitates detection of the instantaneous fetal heart rate which is then used as the Reference Fetal Heart Rate in the assessment of the classification rate of each of the above mentioned techniques. It will be shown that the fetal heartbeats are completely masked by uterine activity and noise artefacts in all the recorded transabdominal maternal ECG signals. The fetal scalp electrode was, therefore, deemed necessary to provide the highest accurate measure of fetal heart functionality (from the hospital viewpoint), and in the assessment of the three non-invasive techniques presented in this thesis. The techniques may also be used during gestation and as early as 10 weeks

Novel Framework for Nonlinear HRV Analysis and its Physiological Interpretation

La inclusión de métodos no lineales aplicados a señales de variabilidad del ritmo cardiaco (HRV, del inglés Heart Rate Variability) proporciona una nueva visión en la caracterización de anomalías en el contexto de las enfermedades cardiacas o patologías como la insuficiencia cardiaca o la fibrilación auricular, por nombrar algunas. Se ha demostrado que alteraciones en el sistema nervioso autónomo (ANS, del inglés Autonomic Nervous System), el cuál modula el ritmo cardiaco, conllevan a cambios en los patrones no lineales de la HRV. Sin embargo, la incertidumbre, todavía presente, en los mecanismos que subyacen a variaciones fisiológicas o patofisiológicas en los índices no lineales de la HRV, junto con el alto tiempo que requieren los algoritmos para la estimación de estos índices, representan el cuello de botella para su aplicación en la práctica clínica.Después de una breve introducción sobre los temas abordados en esta la tesis en el capítulo 1, el segundo capítulo, el capítulo 2, está dedicado a la primera gran contribución de esta tesis, que consiste en la propuesta y desarrollo de una metodología con el fin de reducir el coste computacional asociado a la caracterización no lineal de la HRV. El esquema propuesto es muy eficaz, reduciendo el tiempo de cálculo a unos pocos segundos para el análisis no lineal de señales de HRV de corta longitud (5 minutos). Con respecto a la interpretación del análisis no lineal de la HRV, es importante señalar que hay una serie de factores que afectan a su cálculo y deben tenerse en cuenta al comparar diferentes estudios de la literatura. Las características de las series de HRV, como la frecuencia de muestreo, así como la selección de valores de parámetros en los métodos no lineales, tienen un impacto en los resultados de los índices no lineales de la HRV y, en algunas circunstancias, pueden dar lugar a interpretaciones erróneas. Uno de los principales objetivos del capítulo 3 es estudiar la influencia de la tasa de muestreo en los índices no lineales de la HRV y proponer alternativas para atenuar esta influencia. Los métodos propuestos incluyen, por una parte, la corrección de la frecuencia cardiaca de las estimaciones de la HRV mediante fórmulas de regresión individuales o basadas en la población y, por otra, el preprocesamiento de las series temporales de HRV mediante modelos de interpolación o de point-process. El capítulo 4 se centra en investigar el efecto de la selección del valor de los parámetros requeridos para el cálculo de ciertos índices no lineales de la HRV (por ejemplo, la entropía aproximada) y proponiendo un nuevo índice independiente de la definición del valor de éstos parámetros a-priori. Este novedoso índice se denomina entropía multidimensional aproximada. El análisis no lineal de la HRV, incluido el nuevo índice propuesto, se aplica al estudio de afecciones asociadas a alteraciones de la modulación cardiaca del ANS, como el envejecimiento y la insuficiencia cardiaca congestiva (CHF, del inglés Congestive Heart Failure). Por un lado, todos los índices no lineales de la HRV evaluados ven disminuidos significativamente sus valores en las personas mayores en comparación con los jóvenes ambos grupos en condiciones de reposo en posición de decubito supino. Por otro lado, los pacientes con insuficiencia cardiaca muestran valores más altos de los índices no lineales significativamente con respecto al grupo de sujetos sanos, en ambos casos analizando el período nocturno. Además, el análisis no lineal de la HRV es evaluada en respuesta a provocaciones simpáticas, inducidas por el cambio de la posición supina a la posición de pie o por la administración de atropina, donde se observa una disminución en todos los índicesno lineales estimados.El capítulo 5 está dedicado a la evaluación del rendimiento del análisis no lineal de la HRV en el triaje de la administración profiláctica con el fin de prevenir los episodios de hipotension causados por la anestesia espinal durante el parto por cesárea. El estudio se realiza en colaboración con el Servicio de Anestesia del Hospital Universitario Miguel Servet (Zaragoza, España). Debido a que la profilaxis puede producir efectos secundarios en el feto, el desafío consiste en predecir los casos normotensos para los cuales se puede prescindir del tratamiento profilactico. La hipótesis de esta tesis se basa en el hecho de que la alteración de la regulación del ANS causada por el último período de embarazo y la proximidad a la cirugía podría reflejarse en los índices no lineales de la HRV, lo que podría ayudar a predecir los casos que deriven en hipotension y normotension con mayor precisión que cuando se utilizan solamente variables demográficas. Es importante destacar que las propuestas metodológicas para el análisis no lineal de la HRV desarrolladas en la tesis se aplican en la caracterización de otras señales cardiovasculares, como la señal fotopletismografica de pulso. Las series temporales derivadas de esta señal, que incluyen información del sistema vascular periférico, se incorporan en un clasificador basado en la regresión logística junto con los índices no lineales de la HRV. El clasificador propuesto alcanza un 76,5% de sensibilidad y un 72,2% de precisión en la detección de los casos normotensos, proporcionando así información pertinente y objetiva respaldando la decisión final del equipo médico.En el capítulo 6 se presentan las principales conclusiones derivadas de la tesis y se consideran futuras ampliaciones en base a las investigaciones llevadas a cabo. Se hace hincapié en la contribución de la tesis al desarrollo de metodologías novedosas para caracterizar de manera más robusta los índices no lineales de la HRV e interpretar con fiabilidad los resultados correspondientes. Basándose en las metodologías desarrolladas, se investigan las condiciones o patologías asociadas con alteraciones en la modulación autonómica de la actividad cardiaca y se destaca la contribución del análisis no lineal de la HRV para su caracterización. En conclusión, entre los objetivos metodológicos desarrollados en esta tesis se encuentran: i) la propuesta de un esquema de trabajo para incrementar la fiabilidad de la estimación de la dimensión de correlación, usando un algoritmo que reduce la carga computacional, facilitando su aplicabilidad en la práctica clínica; ii) el desarrollo de métodos alternativos para atenuar la dependencia de los índices no lineales de la HRV con el ritmo cardiaco medio; iii) la propuesta de un índice no lineal de la HRV multidimensional independiente de la definición a priori de parámetros para su estimación. Además, los objetivos relacionados con la aplicación clínica de lascontribuciones metodológicas son: i) la caracterización del efecto del envejecimiento en los índices no lineales de la HRV; ii) la evaluación de la complejidad e irregularidad del ritmo cardiaco en pacientes que sufren de insuficiencia cardiaca comparada con sujetos sanos; iii) la mejora de la eficacia de la profilaxis para la prevención de eventos de hipotensión después de anestesia espinal durante parto programado por cesárea.<br /

Assessing Variability of EEG and ECG/HRV Time Series Signals Using a Variety of Non-Linear Methods

Time series signals, such as Electroencephalogram (EEG) and Electrocardiogram (ECG) represent the complex dynamic behaviours of biological systems. The analysis of these signals using variety of nonlinear methods is essential for understanding variability within EEG and ECG, which potentially could help unveiling hidden patterns related to underlying physiological mechanisms. EEG is a time varying signal, and electrodes for recording EEG at different positions on the scalp give different time varying signals. There might be correlation between these signals. It is important to know the correlation between EEG signals because it might tell whether or not brain activities from different areas are related. EEG and ECG might be related to each other because both of them are generated from one co-ordinately working body. Investigating this relationship is of interest because it may reveal information about the correlation between EEG and ECG signals. This thesis is about assessing variability of time series data, EEG and ECG, using variety of nonlinear measures. Although other research has looked into the correlation between EEGs using a limited number of electrodes and a limited number of combinations of electrode pairs, no research has investigated the correlation between EEG signals and distance between electrodes. Furthermore, no one has compared the correlation performance for participants with and without medical conditions. In my research, I have filled up these gaps by using a full range of electrodes and all possible combinations of electrode pairs analysed in Time Domain (TD). Cross-Correlation method is calculated on the processed EEG signals for different number unique electrode pairs from each datasets. In order to obtain the distance in centimetres (cm) between electrodes, a measuring tape was used. For most of our participants the head circumference range was 54-58cm, for which a medium-sized I have discovered that the correlation between EEG signals measured through electrodes is linearly dependent on the physical distance (straight-line) distance between them for datasets without medical condition, but not for datasets with medical conditions. Some research has investigated correlation between EEG and Heart Rate Variability (HRV) within limited brain areas and demonstrated the existence of correlation between EEG and HRV. But no research has indicated whether or not the correlation changes with brain area. Although Wavelet Transformations (WT) have been performed on time series data including EEG and HRV signals to extract certain features respectively by other research, so far correlation between WT signals of EEG and HRV has not been analysed. My research covers these gaps by conducting a thorough investigation of all electrodes on the human scalp in Frequency Domain (FD) as well as TD. For the reason of different sample rates of EEG and HRV, two different approaches (named as Method 1 and Method 2) are utilised to segment EEG signals and to calculate Pearson’s Correlation Coefficient for each of the EEG frequencies with each of the HRV frequencies in FD. I have demonstrated that EEG at the front area of the brain has a stronger correlation with HRV than that at the other area in a frequency domain. These findings are independent of both participants and brain hemispheres. Sample Entropy (SE) is used to predict complexity of time series data. Recent research has proposed new calculation methods for SE, aiming to improve the accuracy. To my knowledge, no one has attempted to reduce the computational time of SE calculation. I have developed a new calculation method for time series complexity which could improve computational time significantly in the context of calculating a correlation between EEG and HRV. The results have a parsimonious outcome of SE calculation by exploiting a new method of SE implementation. In addition, it is found that the electrical activity in the frontal lobe of the brain appears to be correlated with the HRV in a time domain. Time series analysis method has been utilised to study complex systems that appear ubiquitous in nature, but limited to certain dynamic systems (e.g. analysing variables affecting stock values). In this thesis, I have also investigated the nature of the dynamic system of HRV. I have disclosed that Embedding Dimension could unveil two variables that determined HRV