[发明专利]一种CRP伪影抑制方法、系统、电子设备和存储介质

权利要求书

1.一种CRP伪影抑制方法，其特征在于，包括以下步骤：

S1：获取待处理的包含CPR伪影的原始ECG数据；

S2：对原始ECG数据进行预处理，将ECG信号和CPR伪影进行第一次分离，得到不包含ECG信号的CPR伪影信号；

S3：计算S2中CPR伪影信号的功率谱密度；

S4：根据功率谱密度PSD找出CPR伪影信号的分量频点，再根据分量频点对CPR伪影信号进行重构，得到CPR干扰信号；

S5：采用卡尔曼自适应滤波器对CPR干扰信号进行递归求解得到CPR伪影估计量；

S6：从原始ECG数据去除CPR伪影估计量，得到无CPR伪影的ECG数据。

2.如权利要求1所述的一种CRP伪影抑制方法，其特征在于，所述S1中，原始ECG数据的类型包括至少两个波形；每个波形包括至少100个样本，每个样本至少有1000个采样点。

3.如权利要求1所述的一种CRP伪影抑制方法，其特征在于，所述S2中，预处理方法为：

采用巴特沃斯一阶滤波器设计带通频率为[1Hz,3Hz]的带通滤波器，对原始ECG数据进行低通滤波去除其中的ECG信号，实现ECG信号和CPR伪影的第一次分离。

4.如权利要求1所述的一种CRP伪影抑制方法，其特征在于，所述S3包括：

步骤S31：对CPR伪影信号进行自相关计算：

公式(1)中，R_x(τ)表示对x(t)进行自相关计算的自相关结果；t表示时间；τ表示时延；x(t)表示经过带通滤波器后时间t的时域信号；x(t+τ)表示经过带通滤波后时间t+τ的时域信号；

步骤S32：对自相关结果R_x(τ)进行傅里叶变换：

公式(2)中，S_x(f)表示傅里叶变换结果；j表示复数的虚部；f表示信号采样频率；τ表示时延；

步骤S33：对傅里叶变换结果S_x(f)进行积分得到功率谱密度：

公式(3)中，PSD表示功率谱密度；S_x(f)表示傅里叶变换结果。

5.如权利要求1所述的一种CRP伪影抑制方法，其特征在于，所述S4包括：

S41：根据功率谱密度计算CPR伪影信号的最大分量频点：

f₀＝getMaxPos(PSD) (4)

公式(4)中，f₀表示CPR伪影信号的最大分量频点，即第一主频点；getMax表示返回最大值函数；getMaxPos表示返回功率谱密度最大值的下标及功率密度最大时的频点；

S42：根据分量频点对CPR伪影信号进行重构，得到CPR干扰信号：

S_cpr[n]＝{A₀[n]cos(ω₀n+θ₀)+K*A₁[n]*cos(2ω₀n+θ₁)}；

ω₀＝2πf₀ (5)

公式(5)中，S_cpr[n]表示重构以后的CPR干扰信号；A₀[n]表示第一主频点构造的正弦函数的幅度系数；ω₀表示角速度；n表示信号的总量(例如采集了4秒的信号，采样频率为250Hz，则信号总量为1000)；θ₀表示第一主频点构造的正弦函数的初始相位；A₁[n]表示第二主频点构造的正弦函数的幅度系数；θ₁表示第二主频点构造的正弦函数的初始相位；f₀表示CPR伪影信号的分量频点；

K表示控制第二主频分量的系数：

K＝1；if f₀f_c

K＝0；if f₀≥f_c (6)

公式(6)中，f_c表示截止频率。