[發明專利]一種利用t分布的PLC信號濾波方法和系統

在審

• 申請號: 202010463465.6
• 申請日: 2020-05-27
• 公開（公告）號: CN111541635A
• 公開（公告）日: 2020-08-14
• 發明（設計）人: 翟明岳
• 申請（專利權）人: 廣東石油化工學院
• 主分類號: H04L25/03
• 分類號: H04L25/03;H04B1/10
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.所述一種利用t分布的PLC信號濾波方法，其特征在于，包括：

步驟101獲取按時間順序采集的信號序列S；

步驟102求取分塊塊數N^BLOCK，求取公式為其中，SNR為所述信號序列S的信噪比；為所述信號序列S的方差；表示下取整；N為所述信號序列S的長度；

步驟103生成第m個分塊信號序列s_m，求取公式為其中，m為分塊信號序列的序號，其取值范圍為m＝1,2,…,N^BLOCK；為所述信號序列S的第(m-1)N_s+1個元素；為所述信號序列S的第(m-1)N_s+2個元素；為所述信號序列S的第mN_s個元素；N_s為所述第m個分塊信號序列s_m的長度，其計算公式為表示上取整；

步驟104創建迭代控制參數k并賦值為0；求取第m個t分布序列h的初始值h₀，求取公式為h₀＝s_m；

步驟105所述迭代控制參數k的值加1；求取沖量因子α_k，求取公式為α_k＝[1/ln(SNR+1)]^k；

步驟106求取所述第m個t分布序列h的第k步值h_k，求取公式為其中，h_k-1為所述第m個t分布序列h的第k-1步值；η為下降速度，其計算公式為p_ij為第一相鄰矢量，其計算公式為q_ij為第二相鄰矢量，其計算公式為為所述第m個分塊信號序列s_m的第i個元素，i為元素第一序號，取值范圍為i＝1,2,…,N_s；為所述第m個分塊信號序列s_m的第j個元素，j為元素第二序號，取值范圍為j＝1,2,…,N_s；

步驟107求取逼近誤差e，求取公式為e＝|h_k-h_k-1|；

步驟108判斷所述逼近誤差e是否大于或等于預設閾值ε₀，得到第一判斷結果。如果所述第一判斷結果顯示所述逼近誤差e大于或等于所述預設閾值ε₀，則返回所述步驟105、所述步驟106、所述步驟107和所述步驟108；直至所述第一判斷結果顯示所述逼近誤差e小于所述預設閾值ε₀；并得到第m個分塊濾波信號序列B_m，其計算公式為B_m＝h_k。其中，所述預設閾值∈₀為ε₀＝0.001；

步驟109記錄濾除了噪聲的信號序列S_NEW，具體為

2.所述一種利用t分布的PLC信號濾波系統，其特征在于，包括：

模塊201獲取按時間順序采集的信號序列S；

模塊202求取分塊塊數N^BLOCK，求取公式為其中，SNR為所述信號序列S的信噪比；為所述信號序列S的方差；表示下取整；N為所述信號序列S的長度；

模塊203生成第m個分塊信號序列s_m，求取公式為其中，m為分塊信號序列的序號，其取值范圍為m＝1,2,…,N^BLOCK；為所述信號序列S的第(m-1)N_s+1個元素；為所述信號序列S的第(m-1)N_s+2個元素；為所述信號序列S的第mN_s個元素；N_s為所述第m個分塊信號序列s_m的長度，其計算公式為表示上取整；

模塊204創建迭代控制參數k并賦值為0；求取第m個t分布序列h的初始值h₀，求取公式為h₀＝s_m；

模塊205所述迭代控制參數k的值加1；求取沖量因子α_k，求取公式為α_k＝[1/ln(SNR+1)]^k；

模塊206求取所述第m個t分布序列h的第k步值h_k，求取公式為其中，h_k-1為所述第m個t分布序列h的第k-1步值；η為下降速度，其計算公式為p_ij為第一相鄰矢量，其計算公式為q_ij為第二相鄰矢量，其計算公式為為所述第m個分塊信號序列s_m的第i個元素，i為元素第一序號，取值范圍為i＝1,2,…,N_s；為所述第m個分塊信號序列s_m的第j個元素，j為元素第二序號，取值范圍為j＝1,2,…,N_s；

模塊207求取逼近誤差e，求取公式為e＝|h_k-h_k-1|；

模塊208判斷所述逼近誤差e是否大于或等于預設閾值ε₀，得到第一判斷結果。如果所述第一判斷結果顯示所述逼近誤差e大于或等于所述預設閾值ε₀，則返回所述模塊205、所述模塊206、所述模塊207和所述模塊208；直至所述第一判斷結果顯示所述逼近誤差e小于所述預設閾值ε₀；并得到第m個分塊濾波信號序列B_m，其計算公式為B_m＝h_k。其中，所述預設閾值∈₀為ε₀＝0.001；

模塊209記錄濾除了噪聲的信號序列S_NEW，具體為