Scribd Logo
Upload

Welcome to Scribd!

  • 68 views

    ECG Wave Generation Using Fourier Series - Shreetam

    Uploaded by

    shreetam behera
    This document describes a MATLAB code that implements a Fourier series to generate the different waves of an electrocardiogram (ECG) signal. It defines the Fourier series coefficients and parameters for the P, Q, R, S, T, and U waves and combines them to create a full ECG signal. The code takes in the heart rate as input and calculates the period to use in the Fourier series. It then runs the separate functions defined to generate each wave and plots the combined ECG signal.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Download as PDF, TXT or read online from Scribd

    ECG Wave Generation Using Fourier Series - Shreetam

    Uploaded by

    shreetam behera
    68 views7 pages
    This document describes a MATLAB code that implements a Fourier series to generate the different waves of an electrocardiogram (ECG) signal. It defines the Fourier series coefficients and parameters for the P, Q, R, S, T, and U waves and combines them to create a full ECG signal. The code takes in the heart rate as input and calculates the period to use in the Fourier series. It then runs the separate functions defined to generate each wave and plots the combined ECG signal.

    Original Title

    ECG Wave Generation Using Fourier Series_shreetam

    Copyright

    © Attribution Non-Commercial (BY-NC)

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document describes a MATLAB code that implements a Fourier series to generate the different waves of an electrocardiogram (ECG) signal. It defines the Fourier series coefficients and parameters for the P, Q, R, S, T, and U waves and combines them to create a full ECG signal. The code takes in the heart rate as input and calculates the period to use in the Fourier series. It then runs the separate functions defined to generate each wave and plots the combined ECG signal.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    68 views7 pages

    ECG Wave Generation Using Fourier Series - Shreetam

    Uploaded by

    shreetam behera
    This document describes a MATLAB code that implements a Fourier series to generate the different waves of an electrocardiogram (ECG) signal. It defines the Fourier series coefficients and parameters for the P, Q, R, S, T, and U waves and combines them to create a full ECG signal. The code takes in the heart rate as input and calculates the period to use in the Fourier series. It then runs the separate functions defined to generate each wave and plots the combined ECG signal.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 7

    Shreetam 2012

    ECG Wave Generation using Fourier


    Series

    Shreetam 2012

    Shreetam 2012

    Shreetam 2012

    Shreetam 2012

    Matlab Code:Implement all the separated blocks as different .m files and store all the files in a single
    folder.
    clc;
    clear all;
    close all;
    x=0.60008:0.00008:1.4;
    rate=input('\n\nenter the heart beat rate :');
    li=30/rate;
    a_pwav=0.25;
    d_pwav=0.09;
    t_pwav=0.16;
    a_qwav=0.025;
    d_qwav=0.066;
    t_qwav=0.166;
    a_qrswav=1.6;
    d_qrswav=0.11;
    a_swav=0.25;
    d_swav=0.066;
    t_swav=0.09;
    a_twav=0.35;
    d_twav=0.142;
    t_twav=0.2;
    a_uwav=0.035;
    d_uwav=0.0476;
    t_uwav=0.433;

    Shreetam 2012
    pwav=p_wav(x,a_pwav,d_pwav,t_pwav,li);
    %qwav output
    qwav=q_wav(x,a_qwav,d_qwav,t_qwav,li);
    %qrswav output
    qrswav=qrs_wav(x,a_qrswav,d_qrswav,li);
    %swav output
    swav=s_wav(x,a_swav,d_swav,t_swav,li);
    %twav output
    twav=t_wav(x,a_twav,d_twav,t_twav,li);
    %uwav output
    uwav=u_wav(x,a_uwav,d_uwav,t_uwav,li);
    %ecg output
    ecg=pwav+qrswav+twav+swav+qwav+uwav;
    figure(1)
    %ecg
    plot(x,ecg);
    function [pwav]=p_wav(x,a_pwav,d_pwav,t_pwav,li)
    l=li;
    a=a_pwav;
    x=x+t_pwav;
    b=(2*l)/d_pwav;
    n=100;
    p1=1/l;
    p2=0;
    for i = 1:n
    harm1=(((sin((pi/(2*b))*(b-(2*i))))/(b(2*i))+(sin((pi/(2*b))*(b+(2*i))))/(b+(2*i)))*(2/pi))*cos((i*pi*x)/l);
    p2=p2+harm1;
    end
    pwav1=p1+p2;
    pwav=a*pwav1;
    function [qwav]=q_wav(x,a_qwav,d_qwav,t_qwav,li)
    l=li;
    x=x+t_qwav;
    a=a_qwav;
    b=(2*l)/d_qwav;
    n=100;
    q1=(a/(2*b))*(2-b);
    q2=0;
    for i = 1:n
    harm5=(((2*b*a)/(i*i*pi*pi))*(1-cos((i*pi)/b)))*cos((i*pi*x)/l);
    q2=q2+harm5;
    end
    qwav=-1*(q1+q2);

    function [qrswav]=qrs_wav(x,a_qrswav,d_qrswav,li)
    l=li;
    a=a_qrswav;
    b=(2*l)/d_qrswav;
    n=100;
    qrs1=(a/(2*b))*(2-b);

    Shreetam 2012
    qrs2=0;
    for i = 1:n
    harm=(((2*b*a)/(i*i*pi*pi))*(1-cos((i*pi)/b)))*cos((i*pi*x)/l);
    qrs2=qrs2+harm;
    end
    qrswav=qrs1+qrs2;

    function [swav]=s_wav(x,a_swav,d_swav,t_swav,li)
    l=li;
    x=x-t_swav;
    a=a_swav;
    b=(2*l)/d_swav;
    n=100;
    s1=(a/(2*b))*(2-b);
    s2=0;
    for i = 1:n
    harm3=(((2*b*a)/(i*i*pi*pi))*(1-cos((i*pi)/b)))*cos((i*pi*x)/l);
    s2=s2+harm3;
    end
    swav=-1*(s1+s2);
    function [twav]=t_wav(x,a_twav,d_twav,t_twav,li)
    l=li;
    a=a_twav;
    x=x-t_twav-0.045;
    b=(2*l)/d_twav;
    n=100;
    t1=1/l;
    t2=0;
    for i = 1:n
    harm2=(((sin((pi/(2*b))*(b-(2*i))))/(b(2*i))+(sin((pi/(2*b))*(b+(2*i))))/(b+(2*i)))*(2/pi))*cos((i*pi*x)/l);
    t2=t2+harm2;
    end
    twav1=t1+t2;
    twav=a*twav1;

    function [uwav]=u_wav(x,a_uwav,d_uwav,t_uwav,li)
    l=li;
    a=a_uwav
    x=x-t_uwav;
    b=(2*l)/d_uwav;
    n=100;
    u1=1/l
    u2=0
    for i = 1:n
    harm4=(((sin((pi/(2*b))*(b-(2*i))))/(b(2*i))+(sin((pi/(2*b))*(b+(2*i))))/(b+(2*i)))*(2/pi))*cos((i*pi*x)/l);
    u2=u2+harm4;
    end
    uwav1=u1+u2;
    uwav=a*uwav1;

    You might also like