（五） 找出波峰、剔除异搏、确定RR间期 （六） 对HRV信号进行时域分析
对前面得到的HRV信号进行时域分析，其中的技术指标有NNVGR、SDNN、RNSSD、SDSD、pNN50等。 (七) 线性内插
对提取出来的HRV信号经线性内插获得等间隔（R-R间期均值）R-R间期时间序列，然后再进行频谱分析。线性内插的公式为X(n)?d2RRidRRi?1。式中X(n)为?1d1?d2d1?d2插值，RRi和RRi?1分别为插值前、后的R-R间期序列值，d1和d2分别为X（n）所在位置与RRi和RRi?1所在位置的时间间距。
（八）根据快速傅里叶变换（FFT）得到HRV信号的频谱信号和功率谱信号。主要有TP、HF、LF、LF/HF、频谱图、功率谱图等。
六、 结果
通过MATLAB程序分别得到100、101和102三组ECG信号数据，然后对其进行相关分析，具体结果以及图像如下所述。
（一）ECG信号的提取
通过MATLAB程序分别得到100、101和102三组ECG信号数据。其时序图像具体如图6.1所示。
图 6.1 三组原始心电信号 图示 由上至下分别为来自100号、101号、102号数据的对应图像 （二）原始心电频谱图
对三组信号通过皮谱图进行频域分析。从图6.2中可以看出，信号主要集中在20Hz以下的低频端，并且在60Hz附近有一个明显的幅度值上升，可以判断其为工频干扰。
图 6.2 三组心电信号的频谱图 图示 由上至下分别为来自100号、101号、102号数据的对应图像 （三） 去除工频干扰
使用FIR滤波器法将工频噪声去除。去除工频干扰之后频谱图如图6.3所示
图 6.3 三组心电信号去除工频干扰后的频谱图 图示 由上至下分别为来自100号、101号、102号数据的对应图像