[发明专利]一种基于时频分析的信号去噪方法

摘要

本发明公开了一种基于时频分析的信号去噪方法，包括步骤：一、待处理信号同步存储；二、信号去噪：采用数据处理设备对待处理信号f(t)进行去噪，过程如下：步骤201、基于寻优算法的信号稀疏分解；采用数据处理设备对进行查找时，对的时频参数rn进行查找，过程如下：步骤C1、时频参数寻优；步骤C2、最佳时频参数确定；步骤202、信号重构。本发明方法步骤简单、设计合理且实现方便、使用效果好，采用基于寻优算法的信号稀疏分解方法搜索最佳匹配原子，同时结合适应度值与稀疏度确定最佳匹配原子，能大幅度加快信号去噪速度，并且能有效提高去噪效果，确保去噪后信号的准确性。

主权项

1.一种基于时频分析的信号去燥方法，其特征在于，该方法包括以下步骤：步骤一、待处理信号同步存储：采用数据处理设备(2)对待处理信号f(t)进行同步存储；所述待处理信号f(t)为信号采样系统(1)所采集的信号；其中，f(t)＝[f(t₁),f(t₂),...,f(t_N)]^T，t表示时间参数，t_i为信号采样系统(1)的第i个采样时刻，f(t_i)为待处理信号f(t)中第i个采样时刻的信号值，i为正整数且i＝1、2、3、…、N，N为待处理信号f(t)的信号长度；步骤二、信号去燥：采用数据处理设备(2)对步骤一中所述待处理信号f(t)进行去燥，过程如下：步骤201、基于寻优算法的信号稀疏分解：采用数据处理设备(2)且调用稀疏分解模块对步骤一中所述待处理信号f(t)进行迭代分解处理，将待处理信号f(t)转换为并获得此时的迭代分解最佳原子集合；此时的所述迭代分解最佳原子集合中包含m个最佳匹配原子，为所述迭代分解最佳原子集合中的第n个所述最佳匹配原子；式中R^m(t)为待处理信号f(t)经过m次迭代分解后的残差量，其中m为预先设定的迭代分解总次数且m为正整数，n为正整数且n＝1、2、…、m；a_n为第n次迭代分解后的最佳匹配原子与上一次迭代分解后残差量的展开系数；本步骤中进行信号稀疏分解时，采用数据处理设备(2)由先至后对所述迭代分解最佳原子集合中m个所述最佳匹配原子分别进行查找；为第n次迭代分解时采用数据处理设备(2)寻找出的最佳匹配原子，的最佳时频参数为时频参数r_n；为Gabor原子且式中函数ψ(t)为高斯窗函数且为r_n对应的Gabor原子，r_n为的时频参数，r_n＝(s_n,u_n,v_n,w_n)，s_n为尺度参数，u_n为位移参数，v_n为频率参数，w_n为相位参数；采用数据处理设备(2)对进行查找时，对的时频参数r_n进行查找，过程如下：步骤C1、时频参数寻优：调用寻优算法模块且根据预先设定的时频参数中尺度参数、位移参数、频率参数和相位参数的取值范围，找出适应度值最大的最优时频参数，并找出适应度值＞thrd_z的多个较优时频参数；本步骤中，所找出的所述较优时频参数的数量记作k_n，k_n为整数且k_n≥0；其中，c'为常数且2≤c'≤3；R_n‑1为此时采用数据处理设备(2)已查找出的前n‑1个所述最佳匹配原子的残差量，R_n‑1＝f(t)‑ψ_n‑1(t)，ψ_n‑1(t)为此时采用数据处理设备(2)已查找出的前n‑1个所述最佳匹配原子之和，||R_n‑1||₂为R_n‑1的2‑范数；当n＝1时，R₀＝f(t)；所述最优时频参数为所找出的适应度值最大的时频参数，所述较优时频参数为所找出的适应度值＞thrd_n的时频参数，每个所述时频参数的适应度值均为该时频参数所对应的Gabor原子与R_n‑1的内积；步骤C2、最佳时频参数确定：首先对步骤C1中所述较优时频参数的数量k_n进行判断：当k_n＝0或k_n＝1时，将步骤C1中找出的最优时频参数作为最佳时频参数r_n；当k_n≥2时，对步骤C1中找出的k_n个所述较优时频参数的稀疏度分别进行计算，并将稀疏度最小的所述较优时频参数作为最佳时频参数r_n；对k_n个所述较优时频参数中第d个所述较优时频参数的稀疏度XS_d进行计算时，XS_d＝||R^d||_ξ，R^d为ψ_d(t)的残差量，R^d＝f(t)‑ψ_d(t)，ψ_d(t)为k_n个所述较优时频参数中第d个所述较优时频参数对应的Gabor原子；||R^d||_ξ表示R^d的ξ‑范数，ξ为常数且0≤ξ≤1；d为正整数且d＝1、2、…、k_n；步骤202、信号重构：根据此时的所述迭代分解最佳原子集合，采用数据处理设备(2)得出待处理信号f(t)的近似信号f″(t)；其中，近似信号f″(t)为待处理信号f(t)去燥后的信号，其中为此时所述迭代分解最佳原子集合中的第n'个所述最佳匹配原子，n'为正整数且n'＝1、2、…、m；a_n'为与根据此时所述迭代分解最佳原子集合中的前n'‑1个所述最佳匹配原子对f(t)进行n'‑1次迭代分解后残差量的展开系数。