脉搏信号分析综合训练实例

第2章脉搏信号分析综合训练实例

2.1训练任务

设计视觉疲劳实验，能准确得到视觉疲劳前后两种不同状态；对受试者在实验前、实验过程和实验后的脉搏信号进行采集；基于MATLAB对脉搏信号进行滤波处理、特征提取和分析，以期找到VDT作业疲劳对脉搏信号影响最明显的特征，为视觉疲劳的评测提供依据。

2.2实验设计与数据采集

实验原理：

根据VDT作业的工效学原理，疲劳必会造成生理信号的变化。具体方案是要求受试者完成2小时的VDT作业，通过主观问卷测评受试者，对VDT试验过程中被试者的脉搏信号进行采集，对采集到的信号进行筛选，预处理，特征提取，分类。找出疲劳对生理信号的影响。通过脉搏信号的变化寻找能较快缓解疲劳的方法。

实验仪器：

相同配置电脑，显示器统一型号，无屏保，视屏显示因子各项参数保持不变，刷新频率为85Hz，MP150多导生理记录仪,采样频率设定为1000Hz。

实验步骤：

1、实验前几天用亚健康自评量表筛选被试者。提前告知被试者近期好好休息，保持良好的

生理、心理状况。

2、实验开始前通过健康状况调查表确保受试者无眼病，无眼病家族史、心理病史，无全身

疾病，右利手者，评测受试者此时感觉非常清醒，心理状态良好。

3、给被试者佩戴信号采集设备，此后让被试者调整自己的坐姿以及显示器的高度位置，尽

量使自己处于较舒适的状态。

4、VDT作业之前，采集被试者的脉搏信号数据，时间持续10分钟，并将此段数据保存为实

验前数据段。

5、被试者在放松状态下，完成120分钟的VDT作业，并对整个实验过程的脉搏信号数据进

行实时采集，将数据保存为“实验过程中数据”。

6、通过主观问卷测评受试者有无视觉疲劳症状，是否已处于视觉疲劳状态。

7、关掉电脑显示器，被试者闭眼休息，同时采集该过程的数据进行采集。30分钟以后结束

数据采集，保存为“疲劳恢复数据段”。

如图2.1，为实验过程中采集到的一段脉搏信号：

图2.1实验中采集到的脉搏信号

2.3信号预处理

人体的脉搏信号是一种微弱的信号，信噪比比较的低，正常人的脉象信号在0~20Hz的频率范内，且大约99%的能量分布在0~10Hz，在检测和采集时，由于受到仪器，人体方面的影响，所采集的信号中通常包括干扰信号,如图2.1中的脉搏信号，其中最主要的干扰信号有：

（1）基线漂移、人体呼吸等低频干扰，频率小于1Hz；

（2）工频干扰，是固定频率的干扰，频率为50Hz；

（3）其他一些干扰信号，主要是高频干扰。

所以，为了对信号做准确的分析，在分析处理之前必须做一些必要的预处理。针对信号中存在噪声的特点，基线漂移和呼吸等低频干扰在1Hz以下，而脉搏信号主要在低频范围，所以可以设计让信号先通过一个带通滤波器，借以滤除基线漂移、呼吸引起的干扰（考虑到不丢失太多的其他信息，通带截止频率设置为Wp=[0.9,50]，阻带截止频率设置为Ws=[0.1,200],通带波纹系数Rp=0.5,阻带波纹系数Rs=40），然后再通过一个陷波滤波器滤除固定的工频干扰（通带截止频率设置为Wp=[48,52],阻带截止频率Ws=[49,51],通带波纹系数Rp=0.5,阻带波纹系数Rs=40）。

在设计滤波器前，介绍一下常用数字滤波器知识：

数字滤波器可分为无限冲激响应数字滤波器（IIR）和有限冲激响应数字滤波器（FIR），其具体设计过程不再详述，就一些性能区别作简介。

从性能上说，IIR数字滤波器可用较低的阶数获得高的选择性，幅频特性好，但缺点是相位不是严格线性的。相反，FIR滤波器却可以得到严格的线性相位，但FIR只能用较高的阶数达到高的选择性；对于同样的滤波器设计指标，FIR滤波器所要求的阶数可以比IIR滤波器高5到10倍，成本较高，信号延时也较大。

综合我们选用IIR滤波器，对常用一些IIR滤波器性能做比较。

利用Butterworth滤波器、Chebysheve I型滤波器、Che-bysheveⅡ型滤波器、椭圆滤波器都可以进行低通、高通、带通、带阻滤波器的设计，但是各有特点。

其中Butterworth滤波器通带内的幅频响应曲线能得到最大限度的平滑，但牺牲了截止频率的坡度；Chebysheve I型滤波器通带内等波纹，阻带内单调；ChebysheveⅡ型滤波器通带内单调，然而阻带内等波纹；椭圆滤波器阻带和通带内都是等波纹的，但下降的坡度更大，而且可以以更低的阶数实现和其他两类滤波器一样的性能指标。

设计滤波器时要根据实际情况综合考虑截止特性和相位失真的要求等要求。截止特性好的，相位失真就严重，两者不可兼得。

Matlab信号工具箱提供直接设计IIR数字滤波器的函数，直接调用这些函数就可以很方便地设计满足我们要求的滤波器，常用数字滤波器的函数如下：

1)Butterworth数字滤波器设计：

调用格式：

[N,Wn]=buttord(Wp,Ws,Rp,Rs)；

[b,a]=butter(N,Wn,’ftype’)；

功能：

[N,Wn]=buttord(Wp,Ws,Rp,Rs)，在给定滤波器性能的情况下（通带截止频率Wp、阻带截止频率Ws、通带内最大衰减Rp和阻带内最小衰减Rs），计算Butterworth滤波器的阶数N和截止频率Wn。

[b,a]=butter(n,Wn,’ftype’)，根据阶数N和截止频率Wn计算Butterworth滤波器分子分母系数（b为分子系数的矢量形式，a为分母系数的矢量形式）,ftype为设计滤波器的类型，ftype=high时，为高通滤波器；ftype=bandpass时，为带通滤波器；ftype=stop 时，为带阻滤波器，类型不写时默认为低通滤波器。

2)ChebyshevⅠ型数字滤波器设计：

调用格式：

[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs)；

[b,a]=cheby1(N,Rp,Wn,'ftype')；

3)ChebyshevⅡ数字滤波器设计：

调用格式：

[N,Wn]=cheb2ord(Wp,Ws,Rp,Rs)；

[b,a]=cheby2(N,Rp,Wn,'ftype')；

4)椭圆滤波器设计：

[N,Wn]=ellipord(Wp,Ws,Rp,Rs)；

[b,a]=ellip(N,Rp,Rs,Wp)；

ChebyshevⅠ、ChebyshevⅡ和椭圆滤波器设计所用MATLAB函数的功能具体使用时可查看MATLAB中的帮助文件。

5)零相位数字滤波

调用格式：

y=filtfilt(b,a,x)；

功能：

y=filtfilt(b,a,x)对输入信号x做前向和反向处理，实现零相位数字滤波。b、a分别为滤波器传递函数的分子和分母系数向量（按降幂排列）。

由于Butterworth滤波器通带内有最大的平滑特性，信号经过后衰减小，因此我们选用Butterworth带通滤波器滤除基线漂移和呼吸等引起的干扰，但由于IIR滤波器本身固有的缺点，信号通过Butterworth带通滤波器后相位会失真，故我们可设计零相位Butterworth 带通滤波器去噪。50Hz固定工频在频域中是一个点，因而要求设计的带阻滤波器有好的截止特性，而切比雪夫II型滤波器有较好的截止特性，并且在其通带内单调，故而设计零相位切比雪夫II型滤波器滤除工频干扰。