专利汇可以提供一种基于边信息先验的多导联心电信号重构估计方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于边信息先验的多 导联 心电 信号 重构估计方法,涉及计算机医学信息技术领域。本发明首先将单个导联心 电信号 的估计值作为边信息先验输入进行下一个导联心电信号的 压缩 感知 过程,然后在每个 编码器 进行编码时,对于不被用作边信息的信号,编码器选用行数较小的测量矩阵,在联合译码端,收到的测量值的数目远远小于普通的 正交 匹配追踪重构 算法 所要求的数目。使用本发明方法能够适当的将编码端的复杂度转移到译码端中,降低编码的复杂度并提高压缩过程的抗噪性和鲁棒性,节省心电采集 节点 的存储、计算能 力 和剩余 能量 ,在保持良好压缩效果的同时延长无线体域网络节点的生存周期,实现长时间伴随式监测的目的。,下面是一种基于边信息先验的多导联心电信号重构估计方法专利的具体信息内容。
1.一种基于边信息先验的多导联心电信号重构估计方法,其特征在于:将单个导联心电信号的估计值作为边信息先验输入进行下一个导联心电信号的压缩感知过程,在每个编码器进行编码时,对于不被用作边信息的信号,编码器选用行数较小的测量矩阵,再联合译码端进行译码,收到的测量值的数目小于普通的正交匹配追踪重构算法所要求的数目;包括如下步骤:
步骤1:设置导联心电信号的序号值,在译码端先通过压缩感知贪婪算法对首个导联心电信号进行信号重构,得到首个导联心电信号的重构信号
所述的重构信号,记为
其中,导联心电信号的序号值为i,第i个导联心电信号为xi,使用压缩感知贪婪算法得到的信号重构估计值为 初始化i=1,即首个导联心电信号为x1,使用压缩感知贪婪算法得到的信号重构估计值为
步骤2:将步骤1重构的心电信号 作为先验边信息,确定多导联心电信号重构估计过程的输入参数:测度矩阵和观测向量;
其中,测度矩阵,即Ai+1,通过公式(1)计算:
Ai+1=φi+1·ψ (1)
其中,Ai+1是第i+1个心电信号的Mi+1×N维测度矩阵,φi+1为Mi+1×N维的测量矩阵,ψ为N×N维的稀疏矩阵,即:
其中, EN是N×N维的单位阵;
其中,φi+1为标准正态分布的随机数矩阵,ψ是通过对离散傅里叶变换做归一化处理后得到的快速傅里叶变换基;
观测向量,即yi+1是观测所得的Mi+1×1维向量,通过如下公式(2)计算:
yi+1=φi+1·xi+1 (2)
其中,xi+1是第i+1个导联心电信号;
步骤3:初始化迭代过程的迭代次数、计数最大值、导联心电信号总数、迭代残差初值、列索引序号集合和测度矩阵中对应索引的列集合;
其中,迭代次数记为t、迭代残差初值记为r0、计数最大值Tmax、导联心电信号总数imax、第t次迭代的列索引序号集合记为It、根据It得到测度矩阵Ai+1中对应索引出的列集合记为Ωt;初始化t=0,设置 r0=yi+1-φi+1·xi|i=1,其中, 表示空集;再将t加1,跳至步骤4开始的循环;
其中,计数最大值Tmax的范围为20-200,具体为原始心电信号稀疏度的2-5倍;
步骤4:计算测度矩阵与迭代残差的内积求取列索引序号,具体为:列索引序号其中,αj表示测度矩阵Ai+1的第j列,N为测度矩阵Ai+1的总列数;μt表示在第t次迭代中,用迭代残差rt-1与αj做内积所得的最大项对应所在的矩阵列;
其中, 表示使参数|
μt表示与残差作内积运算的结果取最大值时对应的测度矩阵Ai+1中的列,将μt列所在的列序号记为jmax;则可推导出等式关系:μt=αjmax;
其中,<·,·>表示向量求内积运算;
步骤5:将步骤4求取的列索引序号jmax添加入列索引序号集合It中,根据It得到测度矩阵Ai+1中对应索引出的列μt添加到列集合Ωt中;具体为:
将步骤4中计算出的列索引序号添加到列索引序号集合It中,具体为It=It-1∪{jmax}、Ωt=Ωt-1∪{μt};
其中,符号∪表示集合并运算;
步骤6:计算第t次迭代中信号的表示系数 具体为:
通过最小二乘方法求解方程yi+1=Ωt·θt计算出第t次迭代中信号的表示系数为:
其中, 为t×1维的列向量;ΩtT表示对Ωt进行转置运算,(ΩtTΩt)-1表示对ΩtTΩt进行求逆运算;
其中,argmin||yi+1-Ωt·θt||表示求使得参数||yi+1-Ωt·θt||取最小值时所对应变量θt的值,记为 ||yi+1-Ωt·θt||表示求参数yi+1-Ωt·θt的二范数值;
步骤7:更新第t次迭代中的迭代残差rt;具体基于步骤6中计算出的信号表示系数 通过公式(4)更新迭代残差rt:
步骤8:判断迭代次数t是否已经达到计数最大值Tmax,并根据判断结果进行相应的操作,具体为:
8.1若迭代次数t大于计数最大值Tmax,则输出第t次迭代中信号的表示系数 并跳至步骤9;
8.2若迭代次数t小于等于计数最大值Tmax,则输出第t次迭代中信号的表示系数 并令迭代次数t=t+1,跳至步骤4;
步骤9:得到第i+1个导联心电信号的重构估计信号 具体为:
将步骤8输出的第t次迭代中信号表示系数 作为第i+1个心电导联信号的稀疏表示系数,通过 求出第i+1个心电导联信号得重构估计信号
步骤10:判断是否已经重构完成所有的导联心电信号,并进行相应的操作:
10.1若已重构的心电信号数量达到导联心电信号总数imax,则表示所有的导联心电信号已经完成重构估计,压缩感知过程结束;
10.2若已重构的心电信号数量未达到导联心电信号总数imax,则增加心电信号序号i=i+1,并返回步骤2重复步骤2到步骤10过程。
标题 | 发布/更新时间 | 阅读量 |
---|---|---|
基于跨通道和空间加权卷积特征的脑图像解码方法及系统 | 2020-05-08 | 415 |
一种加权多尺度字典学习框架的行星轮轴承故障辨识方法 | 2020-05-08 | 443 |
一种基于低秩和稀疏约束的地震数据重建方法 | 2020-05-12 | 531 |
一种用于互质采样星载SAR的成像处理方法 | 2020-05-12 | 839 |
一种基于物联网的钢板在线固溶的轧后快速冷却方法 | 2020-05-08 | 192 |
一种基于块稀疏贝叶斯模型的辐射源直接定位方法 | 2020-05-13 | 967 |
用于电动车辆的电池充电状态指示器 | 2020-05-11 | 212 |
一种基于信号稀疏重构的玻璃药瓶内气体浓度检测方法及装置 | 2020-05-12 | 280 |
利用宽带信号驱动负载的系统和方法 | 2020-05-14 | 224 |
一种基于混合高斯模型的多目标ISAR关键点提取方法 | 2020-05-08 | 253 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。