Извлечение QT из сигнала ЭКГ

Question

По сути, я хочу извлечь QT из сигнала ЭКГ, для этого я следую некоторым шагам из статьи. Во-первых, я фильтрую сигнал полосовым фильтром, затем мне нужно применить CWT-преобразование coiflet-1 (никогда не использовал его раньше) из-за QT Pattern = (CWT Output)^2, а затем выбрать QRS Pattern(i) = {0, QT Pattern(i)

**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)**

Жду любых идей, спасибо!!