[發(fā)明專利]一種基于稀疏成分分析的房顫信號提取方法

有效

• 申請?zhí)枺?/td>: 201410722832.4
• 申請日: 2014-12-02
• 公開（公告）號: CN104323772A
• 公開（公告）日: 2015-02-04
• 發(fā)明（設(shè)計）人: 李燈熬;陳園;趙菊敏
• 申請（專利權(quán)）人: 太原理工大學(xué)
• 主分類號: A61B5/046
• 分類號: A61B5/046
• 代理公司: 太原科衛(wèi)專利事務(wù)所(普通合伙) 14100
• 代理人: 朱源

1.一種基于稀疏成分分析的房顫信號提取方法，其特征在于包括以下步驟：

將含有房顫信號和噪聲信號的心電信號作為混合信號X；

稀疏化：用Haar母小波對混合信號X進行小波變換，將混合信號X從空域變換到頻域，并從頻域中選擇對角矩陣為(j＝1,2,...,T/4)，T指混合信號X的周期；

去除零列，同一方向：對于對角矩陣(j＝1,2,...,T/4)中的每個列向量如果則如果則去掉該列得到中間矩陣X’；

線性聚類:中間矩陣X’中若任意兩列向量X'_p和X'_q滿足公式(1)，(1)，則這兩列共線，當(dāng)且僅當(dāng)X'_p與X'_q共線，則X'_p∈θ(k)，X'_q∈θ(k)，通過將中間矩陣X’中所有的列向量線性聚類，可以得到θ(k)k＝1,2,...,t(t≥m)，m為中間矩陣X’行向量的個數(shù)；

估計混合矩陣A_SCA：對每個線性聚類θ(k)，如果該類樣本點個數(shù)少于閾值，則直接將該類刪除；否則，計算該類聚類中心，θ(k)中包含n個元素X′₁,X'₂,...,X'_n，則聚類中心C由公式(2)獲得；(2)，由聚類中心C得到聚類中心矩陣，聚類中心矩陣則為混合矩陣A_SCA；

計算源信號：根據(jù)S'＝A_SCA^-1X，即可恢復(fù)出源信號S′==S1′S2′S3′;]]>

選取歸一化相關(guān)系數(shù)NCC對源信號S'進行反饋，當(dāng)源信號S'的min∑₂NCC(S',X)最小時，則認為它是最純凈的信號，輸出該心電信號S₁'，

構(gòu)造新混合信號X_new：令輸出的心電信號S₁'為零并通過X_new＝A_SCAS'反作用于混合信號X，構(gòu)成新混合信號X_new，

對新混合信號X_new重復(fù)前述步驟，即可在第二次迭代中得到房顫信號，剩余部分為噪聲信號。