Help with an ECG signal

You are provided with .wav file which contains a simulated electrocardiogram signal (ECG) with a low frequency (drift) signal, an interfering tone and white Gaussian noise added. Using MATLAB you will analyse this signal in the time and frequency domain to quantify the interference/noise. You will then design and implement digital filters to remove as much of the interference/noise as possible while preserving the underlying ECG signal.

Link to the noised ECG signal and desired output is given below:-

https://drive.google.com/drive/folders/12rrvzTg9nBdeMY-dp5IylfMwsmYgsztb?usp=sharing

Problem 1:-

First read the provided signal (ECG.wav) into MATLAB then plot and inspect the signal in the time domain, recording your observations on the nature of the interference. Next, view the signal in the frequency domain and extract the key parameters of interest e.g. the frequency of tonal interference and the level and spectral distribution of noise. These parameters are required for the filter design tasks to follow.

Problem 2:-

Design an infinite impulse response (IIR) highpass filter based on a 2nd order Butterworth prototype using the bilinear z-transform method. This should have a cut off frequency of 0.1 Hz to remove the low frequency drift on the signal. Verify the response of this filter in MATLAB and then apply to the ECG signal, observing the result in the time domain.

Problem 3:-

Using the impulse response truncation method, design a finite impulse response (FIR) lowpass filter with a passband edge at 10 Hz and at least 50 dB attenuation at the frequency of the tonal interference identified in problem 1. Verify the response of this filter, then apply this to the highpass filtered ECG signal and inspect the output signal in the time and frequency domain.

Problem 4:-

Increase the cut-off frequency of the FIR lowpass filter applied in problem 3 to 50 Hz to preserve more of the ECG signal spectrum. Now remove the tonal interference signal using an adaptive filter in the noise cancelling configuration (direct model) and compare the output signal with that produced in problem 3. (The reference interference signal may be obtained by generating a sine wave at the interference frequency identified in problem 1).

Problem 5:-

The ECG signal is significantly oversampled and hence the above operations could be carried with fewer operations by using a multirate DSP approach. Analyse the computations required to achieve the output of problem 3 and propose an alternative approach using a two stage filter to reduce computation. Verify that this produces a similar output with fewer computations and faster execution time in MATLAB.

[Abdullah Javed]

Not assignment but a bonus task

[Abdullah Javed]

i have code for the first three ...need help in 4 and 5 only

[Walter Roberson]

You did not ask a MATLAB question.

[Abdullah Javed]

it is a matlab question sir..

[Walter Roberson]

It looks like an algorithm question to me, not a MATLAB question. A MATLAB question would be like, "Could someone explain logical indexing?", or "What is the meaning of the 'same' option for conv2()?" or "could someone help me figure out why this code is telling me index out of range on line 17?" or "Is there a function for Debroglis Wavelet 7?"

[Abdullah Javed]

Is there any function for adaptive filter in matlab??

[Walter Roberson]

https://www.mathworks.com/help/dsp/ug/overview-of-adaptive-filters-and-applications.html

[Sumit]

can u please provide that code?

[Walter Roberson]

https://www.mathworks.com/matlabcentral/answers/724308-help-with-an-ecg-signal#comment_2735789

You would be after Muhammed in line.

Answer 1

[Constantino Carlos Reyes-Aldasoro]

Please read the following live script, it describes several of the issues you raised:

https://uk.mathworks.com/matlabcentral/fileexchange/68246-visualization-and-analysis-of-an-electrocardiogram-signal?s_tid=srchtitle

https://github.com/reyesaldasoro/Electrocardiogram

[Abdullah Javed]

thanks sir

[Muhammad Bilal]

can you provide me the answers of these tasks.

[Walter Roberson]

@Muhammad Bilal

Yes, yes I could provide you with the answers of the tasks posted in the Question.

Shall I add it to the end of my "Homework to do for other people" list ? At the moment you would be 13973'rd in line.