[发明专利]一种ECG信号联合基线校正及降噪的方法在审
申请号: | 201910018584.8 | 申请日: | 2019-01-09 |
公开(公告)号: | CN109602416A | 公开(公告)日: | 2019-04-12 |
发明(设计)人: | 舒明雷;王枭;朱清;许继勇;孔祥龙;刘照阳;高天雷;单珂 | 申请(专利权)人: | 山东省计算中心(国家超级计算济南中心);山东大学齐鲁医院 |
主分类号: | A61B5/0402 | 分类号: | A61B5/0402 |
代理公司: | 济南泉城专利商标事务所 37218 | 代理人: | 支文彬 |
地址: | 250014 山东省济*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基线校正 降噪 原始信号 恢复 滤波算法 细节信息 噪声抑制 有效地 平滑 真实性 联合 保留 保证 | ||
1.一种ECG信号联合基线校正及降噪的方法,其特征在于,包括:
a)计算机加载ECG信号数据y,提取ECG信号数据长度为N;
b)设置LTI滤波器,输入LTI滤波器的截止频率为fc、滤波阶数d、拆分矩阵的阶数K;
c)通过公式P=AT*A,Q=A*AT+B*BT计算LTI滤波器所需的稀疏带状矩阵P、Q,其中A为d阶差分矩阵,B为d-1阶差分矩阵,AT为A的转置矩阵,BT为B的转置矩阵;
d)进行LTI滤波与稀疏结合优化,引入非对称惩罚函数θ(x;r),其中式中r为非对称惩罚因子,x为原始ECG信号;
e)引入优化参数b、d,其中d=QTQP-1y-λPTb,式中λ为差分参数,P-1为矩阵P的逆矩阵,PT为矩阵P的转置矩阵,QT为矩阵Q的转置矩阵;
f)设置差分矩阵D,当K=0时,D为N阶单位矩阵,当K≥1时,N阶单位矩阵中每一行位置对应的数值减去下一行位置相对应的数值得到(N-K)*N阶矩阵,当K=1时,其为(N-1)*N阶差分矩阵;
g)设置生成的N阶矩阵为主对角线全为1,其余均为0的矩阵,引入第三方矩阵M,其中M=2*Gamma+D-1*Gamma*D;
h)设定算法迭代最大次数Nit,对最终优化x进行迭代循环计算,其中
x=P(QT*Q+PTMP)d;
i)当x前后两次的差值小于0.005或达到最大迭代循环次数Nit时停止循环,输出经过基线校正和降噪后的恢复信号x。
2.根据权利要求1所述的ECG信号联合基线校正及降噪的方法,其特征在于:步骤a)中在离散数学中通过公式y=x+f+n建立ECG信号数据y的模型,其中x为ECG原始信号,f为加入基线漂移的基线信号,n为加入稳定的方差为σ2的高斯白噪声。
3.根据权利要求1所述的ECG信号联合基线校正及降噪的方法,其特征在于:步骤a)中N为3000。
4.根据权利要求1所述的ECG信号联合基线校正及降噪的方法,其特征在于:步骤b)中LTI滤波器的截止频率为fc=0.009、滤波阶数d=1、拆分矩阵的阶数K=3。
5.根据权利要求1所述的ECG信号联合基线校正及降噪的方法,其特征在于:步骤b)中非对称惩罚因子r=1。
6.根据权利要求1所述的ECG信号联合基线校正及降噪的方法,其特征在于:步骤h)中最大次数Nit=40。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于山东省计算中心(国家超级计算济南中心);山东大学齐鲁医院,未经山东省计算中心(国家超级计算济南中心);山东大学齐鲁医院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910018584.8/1.html,转载请声明来源钻瓜专利网。