[发明专利]一种利用KL散度的功率信号滤波方法和系统

权利要求书

1.一种利用KL散度的功率信号滤波方法，其特征在于，包括：

步骤101获取按时间顺序采集的信号序列S

步骤102求取信号分块块数，具体为：信号分块块数记为N^BLOCK，所用求取公式为：

其中：

为所述信号序列S的方差,

N为所述信号序列S的长度,

SNR为所述信号序列S的信噪比，

表示对下取整；

步骤103求取分块信号序列，具体为：第m个分块信号序列记为b^m，所用求取公式为：

其中：

m＝1,2,…,N^BLOCK为块序号，

为所述分块信号序列b^m的长度，

表示对上取整，

所述信号序列S的第(m-1)N₀+1个元素，

所述信号序列S的第(m-1)N₀+2个元素，

所述信号序列S的第mN₀个元素；

步骤104创建迭代控制参数并赋值，生成N^BLOCK个KL散度矢量并赋值，具体为：

迭代控制参数记为k，并赋值为0；第m个KL散度矢量记为h_m，其初始化值赋值为b^m；

步骤105更新N^BLOCK个分块信号序列的值，具体为：

其中：

为所述第m个KL散度矢量h_m的第k+1步值，

为所述第m个KL散度矢量h_m的第k步值，

为冲量因子，

为下降速度，

为相邻因子，

s_i为所述信号序列S中的第i个元素，

s_j为所述信号序列S中的第j个元素，

i＝1,2,…,N为行序号，

j＝1,2,…,N为列序号，

m₀为所述信号序列S中的均值；

步骤106求取迭代误差并结束迭代更新，具体为：迭代误差记为ε＝||b^k+1-b^k||；如果所述迭代误差ε≥0.001，则返回所述步骤105和所述步骤106重新迭代更新；否则，迭代更新过程结束，并将所述迭代控制参数k的值赋给步骤结束变量K_o，并得到第m个分块滤波后信号序列

步骤107求取重构后的信号序列，具体为：重构后的信号序列记为S_new，所用求取公式为：

2.一种利用KL散度的功率信号滤波系统，其特征在于，包括：

模块201获取按时间顺序采集的信号序列S

模块202求取信号分块块数，具体为：信号分块块数记为N^BLOCK，所用求取公式为：

其中：

为所述信号序列S的方差,

N为所述信号序列S的长度,

SNR为所述信号序列S的信噪比，

表示对下取整；

模块203求取分块信号序列，具体为：第m个分块信号序列记为b^m，所用求取公式为：

其中：

m＝1,2,…,N^BLOCK为块序号，

为所述分块信号序列b^m的长度，

表示对上取整，

所述信号序列S的第(m-1)N₀+1个元素，

所述信号序列S的第(m-1)N₀+2个元素，

所述信号序列S的第mN₀个元素；

模块204创建迭代控制参数并赋值，生成N^BLOCK个KL散度矢量并赋值，具体为：

迭代控制参数记为k，并赋值为0；第m个KL散度矢量记为h_m，其初始化值赋值为b^m；

模块205更新N^BLOCK个分块信号序列的值，具体为：

其中：

为所述第m个KL散度矢量h_m的第k+1步值，

为所述第m个KL散度矢量h_m的第k步值，

为冲量因子，

为下降速度，

为相邻因子，

s_i为所述信号序列S中的第i个元素，

s_j为所述信号序列S中的第j个元素，

i＝1,2,…,N为行序号，

j＝1,2,…,N为列序号，

m₀为所述信号序列S中的均值；

模块206求取迭代误差并结束迭代更新，具体为：迭代误差记为ε＝||b^k+1-b^k||；如果所述迭代误差ε≥0.001，则返回所述模块205和所述模块206重新迭代更新；否则，迭代更新过程结束，并将所述迭代控制参数k的值赋给步骤结束变量K_o，并得到第m个分块滤波后信号序列

模块207求取重构后的信号序列，具体为：重构后的信号序列记为S_new，所用求取公式为：