QT间期分析算法代码 参考论文：New algorithm for QT interval analysis in 24-hour Holter ECG: performance and

% 测试QT间期分析算法
% 参考论文：New algorithm for QT interval analysis in 24-hour Holter ECG:
% performance and applications
% P. Laguna N.V.Thakor etc

clc; clear; close all;
load ECG_sample_noisy.mat
curpath = pwd;
%cd ('..');
%pan_tompkin(ECG1,ECGinfo.Fs_ECG1,1);
%ECG2 = ECG2(1:2500);
%ECG2=ECG1;
% 预处理
[d,f]=Preprocessing(ECG2);
% plot(ECG2); hold on;
% plot(d,'g');
% plot(f,'r'); grid

% QRS检测
[PK,RRav]=QRSdetection(f,ECGinfo.Fs_ECG2);
fs=ECGinfo.Fs_ECG2;

%plot(PK,f(PK),'ro'); hold on;
Rp=r_position(f,PK);
Qb=qrs_onset(d,Rp,PK,fs);
[Tpeak,Te]=t_end(f,Rp,fs,RRav);

figure;
ax(1)=subplot(211);
plot(f,'g'); hold on;
plot(PK,f(PK),'ro');
plot(Rp,f(Rp),'*'); grid;
title('qrs detection and R position');

% plot(PK,H,'g');plot(PK,-H,'g');
% plot(PK,Hs,'k');
ax(2)=subplot(212);
%plot(PK,RR);
plot(d,'g'); hold
plot(Qb,d(Qb