[发明专利]动脉压信号中重搏切迹点识别方法

权利要求书

1.动脉压信号中重搏切迹点识别方法，其特征在于，将动脉压检测仪采集的动脉压信号输入计算机，由计算机进行低通滤波和采样的预处理，然后采用计算机识别动脉压信号中的重搏切迹点；采用计算机识别重搏切迹点的具体步骤包括：

A)构建模板数据库；该步骤具体包括：

a1)将重搏切迹已知的动脉压信号中一个节拍信号的采样点数归一化处理为K个，作为一个模板信号；由此建立若干个波形互不相同的模板信号；

a2)分别建立各个模板信号中已知重搏切迹所在采样点的对数极坐标分布模型，并以之构建模板数据库；

B)识别待测的动脉压信号中的重搏切迹点；该步骤具体包括：

b1)对待测的动脉压信号进行逐拍分割，将分割所得每个待测节拍信号的采样点数归一化处理为K个；

b2)识别出待测的动脉压信号经逐拍分割后得到的各个待测节拍信号中第i个待测节拍信号的重搏切迹点；该步骤的过程为：

<b21>预设曲率阈值τ₀；计算第i个待测节拍信号中各个采样点所在信号段的曲率，将所在信号段的曲率大于或等于曲率阈值τ₀的采样点作为第i个待测节拍信号中的目标采样点；

<b22>分别建立第i个待测节拍信号中各个目标采样点的对数极坐标分布模型，并分别与模板数据库中每一已知重搏切迹所在采样点的对数极坐标分布模型进行匹配运算，进而从第i个待测节拍信号的各个目标采样点中识别出第i个待测节拍信号的重搏切迹点；

b3)重复步骤b2)，逐一识别出待测的动脉压信号的各个待测节拍信号的重搏切迹点；

b4)显示并存储待测的动脉压信号中重搏切迹点的识别结果。

2.根据权利要求1所述的动脉压信号重搏切迹点识别方法，其特征在于，所述K的取值范围为100～1000。

3.根据权利要求1所述的动脉压信号重搏切迹点识别方法，其特征在于，所述采样点所在信号段，是指该采样点及其前后各E个采样点构成的信号段，E的取值范围为3～8。

4.根据权利要求1～3中任一项所述的动脉压信号重搏切迹点识别方法，其特征在于，所述建立已知重搏切迹所在采样点或目标采样点的对数极坐标分布模型的具体步骤如下：

(1)建立模型分布域，所述模型分布域是预设定对数极径半径ξ_max、极角范围为(-π，π]的对数极坐标域，并且其对数极径半径ξ_max等分为R个极径区间，极角范围(-π，π]等分为C个极角区间，从而构成模型分布域的R×C个分布区间；

(2)将已知重搏切迹或目标采样点所在的节拍信号映射到直角坐标系中，按下式求取已知重搏切迹所在采样点或目标采样点分别与其所在节拍信号中其它K-1个采样点的差向量，记为Δa(k₀→k)：

Δa(k₀→k)＝(x₀-x_k，y₀-y_k)；

其中，k₀表示已知重搏切迹所在采样点或目标采样点在其所在节拍信号中的采样点序号，(x₀，y₀)表示已知重搏切迹所在采样点或目标采样点在直角坐标系中的向量坐标，(x_k，y_k)表示已知重搏切迹或目标采样点所在节拍信号中第k个采样点在直角坐标系中的向量坐标；

由此得到已知重搏切迹所在采样点或目标采样点分别与其所在节拍信号中其它K-1个采样点的差向量集，记为ΔA(k₀→k)＝{Δa(k₀→k)|k＝1，2，...，K}；

(3)对差向量集ΔA(k₀→k)中的差向量进行对数极坐标转换，得到对应于已知重搏切迹所在采样点或目标采样点的对数极坐标集ΔB(k₀→k)＝{Δb(k₀→k)|k＝1，2，...，K}；Δb(k₀→k)表示差向量Δa(k₀→k)对应的对数极坐标，即：

Δb(k0→k)=(ξk0→k,ψk0→k);]]>其中，

ξk0→k=log(x0-xk)2+(y0-yk)2;]]>

(4)计算对应于已知重搏切迹所在采样点或目标采样点的对数极坐标集ΔB(k₀→k)中对数极坐标在模型分布域中各个分布区间的分布数量，得到已知重搏切迹所在采样点或目标采样点的对数极坐标分布模型，记为F：

F＝{f(r，c)|r＝1，2，...，R；c＝1，2，...，C}；

其中，f(r，c)表示对应于已知重搏切迹所在采样点或目标采样点的对数极坐标集ΔB(k₀→k)中分布在模型分布域第r个极径区间和第c个极角区间所构成的分布区间以内的对数极坐标的数量。

5.根据权利要求4所述的动脉压信号重搏切迹点识别方法，其特征在于，所述步骤<b22>中，第i个待测节拍信号中各个目标采样点的对数极坐标分布模型分别与模板数据库中每一已知重搏切迹所在采样点的对数极坐标分布模型进行匹配运算的具体步骤如下：

[1]按下式计算第i个待测节拍信号中一个目标采样点的对数极坐标分布模型F_Si与模板数据库中一个已知重搏切迹所在采样点的对数极坐标分布模型F_D的卡方统计检验值P[F_D，F_Si]：

P[FD,FSi]=12Σr=1r=RΣc=1c=CgD,Si,]]>

其中，f_Si(r，c)∈F_Si，f_R1(r，c)∈F_R1；

从而计算得到第i个待测节拍信号中各个目标采样点的对数极坐标分布模型分别与模板数据库中各个已知重搏切迹所在采样点的对数极坐标分布模型的卡方统计检验值；

[2]以步骤[1]中所得的各个卡方统计检验值中最小值所对应的目标采样点作为第i个待测节拍信号中的重搏切迹点。