医学信号处理原理和方法

1、第八讲信号的时域测量n时域测量常用的算法信号微分（差分）n二点公式f (a) = (-f (a)+f (a+h) / hn三点公式f (a) = (-f (a-h)+f (a+h) / (2*h)微分（差分）计算图解微分（差分）计算图解 信号积分n梯形公式s = (f (a)+f (b) * (b-a) / 2n辛卜生(simpson)公式s = (f (a)+4f(c)+f(b) * (b-a) / 6n柯特斯(cotes)公式s = (f (a)+32(d)+12f(c)+32f(e)+f(b) * (b-a) / 90插值的常用方法（线性插值和拉格朗日插值公式拉格朗日插值多项式图示手工测

2、量和自动测量1、手工测量2、自动识别 周期和波群的判断方法周期和波群的判断方法 幅度单阈值法。 幅度双阈值法。根据信号的最大值和最小值，确立上、下二个判断阈值g1和g2 g1=pmin+(pmax-pmin)/4g2=pmax- (pmax-pmin)/4 综合判别法（包括：幅度、变化率、时程等）周期判断图解周期判断图解测量实例 （一）血压波测量 （二）心电图测量典型的正常心电模式图五种实验动物的心电图间期均值猴 兔豚鼠大鼠小鼠p波(s)0.0370.0310.0200.0170.070p-r间期(s)0.0780.0680.0630.0480.120qrs波(s)0.0370.0420.01

3、30.0130.060q -t间期(s)0.2000.1400.1300.0740.210s -t间期(s)0.078t波(s)0.0370.0650.0440.0630.056心率(次/min) 215 247 261 385心电图自动识别方法举例caceres自动识别qrs波各顶点的方法 采样率为625次秒，测量导联为、3导联。 r波顶点：根据微分值来确定r波的出现。如果在最大的负微分值前0.128秒（80个采样点）内存在微分值大于9.375mv/s的点，则认为有r波存在，然后在该点后的40个样点内找到最大点定为r波顶点。（三）脉搏波分析脉像研究a. d/h 值 反映血管弹性b. 谱分析

4、中心频率(fc)和er值er=e1/(e2+e3+e4+e5)（四）心肌动作电位四）心肌动作电位 n振幅c点与a点对应值的差值n静息电位a点对应值n超射幅度c点对应值napd10复极到10%的间隔时间napd20复极到20%的间隔时间napd50复极到50%的间隔时间napd90复极到90%的间隔时间napd50/apd10n最大上升速度a、c间上升速度的最大值n最大复极速度 c、b间复极速度的最小值n复极化平均斜率c、b间复极速度的平均值tmax ：收缩最大张力:a、b区域间的最大值tmin ：舒张最小张力:a、b区域间的最小值t：张力增量:a点对应值与tmax的差值sti：收缩间期:a点到

5、tmax对应点的间隔时间dti：舒张间期:tmax对应点到tmin对应点的间隔时间dti50：舒张50%间期:肌肉开始舒张到舒张50%的间隔时间dti90：舒张90%间期:肌肉开始舒张到舒张90%的间隔时间+dt/dt max:肌肉收缩时张力最大变化速率-dt/dt max:肌肉舒张时张力最大变化速率t-dt/dt max:肌肉开始收缩至发生dt/dt max的间隔时间n最大呼气峰压a、c间的最大值n最小吸气谷压a、c间的最小值n+dp/dt max从谷压点到峰压点上升速度最大值n-dp/dtmax从峰压点到c点下降速度最小值n呼气间期b、c间隔时间n吸气间期a、b间隔时间n呼吸时比呼气/吸气间期n平均呼气压b、c间的平均值 n平均吸气压a、b间的平均值n呼吸频率 由a、c为一个周期计算出的频率（