pan_tompkin算法详细注解（便于快速理解算法代码）

function [qrs_amp_raw,qrs_i_raw,delay] = pan_tompkin(ecg,fs,gr)

%% function [qrs_amp_raw,qrs_i_raw,delay]=pan_tompkin(ecg,fs)
% Complete implementation of Pan-Tompkins algorithm
% Pan-Tompkins 算法的完整实现

%% Inputs
% ecg : raw ecg vector signal 1d signal          %原始ecg向量，1维信号
% fs : sampling frequency e.g. 200Hz, 400Hz and etc  %采样率，文献上是200Hz,这个代码支持400Hz等其他采样率
% gr : flag to plot or not plot (set it 1 to have a plot or set it zero not
% to see any plots         % 1=绘图，0=不绘图
%% Outputs
% qrs_amp_raw : amplitude of R waves amplitudes   % R波幅度(带通信号上的)
% qrs_i_raw : index of R waves                    % R波索引（带通信号上的，R波所在位置的样本序号，从1开始）
% delay : number of samples which the signal is delayed due to the
% filter