MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   资料存档 (https://www.labfans.com/bbs/forumdisplay.php?f=72)
-   -   Extracting QT from EKG signal (https://www.labfans.com/bbs/showthread.php?t=22175)

poster 2019-11-27 05:26

Extracting QT from EKG signal
 
<p>Basically, I want to extract the QT from an EKG signal, for this I am following some steps from an article. Firstly, I am filtering the signal with a pass-band filter, then I have to apply the coiflet-1 CWT transform(never used it before) due to QT Pattern = (CWT Output)^2 and then choose
QRS Pattern(i) = {0 , QT Pattern(i) < A
QT Pattern(i), otherwise.
This is what I did until now, but i found myself in a moment when i don't know what to do next: </p>

<pre><code>**clear all
clc
close all
[filename, pathname] = uigetfile('semnal.dat');
if isequal(filename, 0) || isequal(pathname, 0)
disp('File input canceled.');
ECG_Data = [];
else
fid=fopen(filename,'r');
end;
time=10;
f=fread(fid,2*360*time,'ubit12');
semnal=f(1:2:length(f));
figure(1)
plot(semnal)

xlabel('samples');
ylabel('time')


fs = 250; % Sampling Frequency (Hz)



N=length(semnal);
fs=250;
t=[0:N-1]/fs;
fcutlow=0.5; %low cut frequency in Hz
fcuthigh=40; %high cut frequency in Hz
[b,a]=fir1(330,[fcutlow,fcuthigh]/(fs/2),'bandpass');
semnal_filtrat=filter(b,a,semnal)
figure(2)
plot(t,semnal_filtrat)

%CTW


x = [1, 2, 3, 4];
scales = [9];
wavelet_name = 'coif2';
coefs = cwt(x,scales, wavelet_name)**
</code></pre>

<p>I am waiting for any idea, thank you!!</p>



[url=https://stackoverflow.com/questions/59059438/extracting-qt-from-ekg-signal]More answer...[/url]


所有时间均为北京时间。现在的时间是 19:46

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.