[发明专利]一种利用t分布的PLC信号滤波方法和系统

权利要求书

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，具体为