基于经验模态分解改进的乘性噪声去除技术制造技术

本发明专利技术公开了一种基于经验模态分解改进的乘性噪声去除技术，其特征在于，该技术包括如下步骤：1)谐波信号及其同态变型的噪声能量分布特性分析；2)本征模函数幅值滤波阈值准则的调整；3)“基‑2”同态变换预处理；4)基于经验模态分解的信号去噪处理；5)基于去噪本征模函数分量的信号重构；6)基于“基‑2”同态逆变换的去噪观测信号后处理。本发明专利技术的有益效果是，方法设计合理，使用简明，去噪性能高。

Improved multiplicative noise removal technology based on empirical mode decomposition

The invention discloses a multiplicative noise removal technique based on empirical mode decomposition, which is characterized by the following steps: 1) the analysis of the noise energy distribution of the harmonic signal and its homomorphic variant; 2) the adjustment of the threshold criterion for eigenmode function amplitude filtering; 3) \base 2\ homomorphic transformation; 4) Signal de-noising in empirical mode decomposition; 5) signal reconstruction based on denoising eigenmode function component; 6) post-processing of de-noised observation signal based on \base 2\ homomorphic inversion. The beneficial effect of the invention is that the method is reasonable in design, concise in use and high in denoising performance.

【技术实现步骤摘要】
基于经验模态分解改进的乘性噪声去除技术
本专利技术基于信号处理理论，由于利用同态变换可以将信号的乘法混叠模型转化为加法模型，因此在此基础上进一步应用加性噪声消除技术即可达到消除乘法性噪声并恢复真实源信号的目的。该消噪技术为解决图像处理、雷达、遥感、激光学等许多学科领域中面临的乘性噪声消除问题奠定了理论基础。
技术介绍
乘(法)性噪声去除，在图像处理、雷达、遥感、激光学等许多学科领域中都会涉及到[1]。乘性噪声往往由信道不理想引起，而乘性随机性可以看成系统的时变性或者非线性所造成的。目前常用的乘性噪声去除方法有滤波技术和盲源分离等[2-4]，像自适应维纳滤波、非线性滤波及同态滤波等技术相对成熟，特别是同态滤波应用较多[3-4]。但是，该方法要求对真实信号的频谱特征有先验的了解，否则难以取得满意的去噪效果。至于盲源分离去噪，由于算法本身的幅值、相位及排序不确定性容易造成去噪结果的不稳定，其实际应用也受到限制[5]。20世纪末，Huang等提出了经验模态分解(EMD)技术，用于解决非线性、非稳态数据的分析问题[6]。将EMD应用于“纯”白噪声信号分解的研究结果显示：随着分解深度的增加，所得的本征模函数(IMF)能量以对数规律递减[7-8]，这表明受白噪声污染的观测信号其低阶的IMFs主要是高频的噪声成分，高阶的IMFs则主要包含低频的源信号信息。基于此，多种EMD基信号去噪方法相继被提出[9-11]。其中，KOPSINIS与MCLAUGHLIN受小波去噪原理的启发，建立了多种适应EMD特性的幅值滤波规则，进而提出了改进的EMD基去噪方法。与已有的EMD基去噪方法相比，改进方法获得了更好的去噪效果[9]。近年来，多位学者还根据不同应用领域中面临的信号去噪需求，对已有的EMD基去噪方法从不同角度进行了改进[10-12]。不过，现有的EMD基去噪方法主要针对加(法)性噪声设计，无法直接用来解决乘(法)性噪声的去除问题，有必要加以改进。
技术实现思路
本专利技术的目的是为了解决上述问题，开发了一种基于经验模态分解改进的乘性噪声去除技术。实现上述目的本专利技术的技术方案为，一种基于经验模态分解改进的乘性噪声去除技术，其特征在于，该方法包括如下步骤：1)谐波信号及其同态变型的噪声能量分布特性分析；2)本征模函数幅值滤波阈值准则的调整；3)“基-2”同态变换预处理；4)基于经验模态分解的信号去噪处理；5)基于去噪本征模函数分量的信号重构；6)基于“基-2”同态逆变换的去噪观测信号后处理。所述谐波信号及其同态变型的噪声能量分布特性分析计算式为：式中x为带噪观测信号向量，xi为第i个信号分量，E(·)为数学期望估计函数，σ2|S(·)为标准方差估计函数，median(·)为中值函数，I1为第1个本征模函数(IMF)分量，σ2|R(·)为鲁棒方差估计函数，N为样本数。所述本征模函数幅值滤波阈值准则的调整计算式为：式中E1为原始乘性带噪信号经过经验模态分解(EMD)抽取的第1阶IMF分量的噪声能量的估计，E1′为乘性带噪信号经过同态变换解耦后由EMD抽取的第1阶IMF分量的噪声能量的估计，PK(·)为基于方差散点数据学习得到的K阶回归多项式，ak,k＝0,1,...,K为多项式系数，Ei′为乘性带噪信号经过同态变换解耦后由EMD抽取的第i阶IMF分量的噪声能量估计，Ti′为幅值滤波阈值，N为信号样本点数，C为常数通常取为C＝0.7，β与ρ均为能量估算参数。所述“基-2”同态变换预处理计算式为：式中x、s与u分别为带噪观测信号、真实源信号与噪声信号向量，x'、s'与u'分别为经过同态变换预处理后所获得的新的虚拟带噪观测信号、虚拟源信号与虚拟噪声信号向量。所述基于经验模态分解的信号去噪处理计算式为：刚性取阈：柔性取阈：式中为取阈运算前、后第i个期望带噪IMF分量，Ti为作用于I(i)(t)上的阈值，zj(i)＝[wj(i)wj+1(i)]，i＝1,2,...,M，j＝1,2,...,Nz(i)为第i个期望带噪IMF分量的过零点区间，Nz(i)为过零区间zj(i)的个数，wj(i)、wj+1(i)为过零区间端点，I(i)(rj(i))为过零区间信号的极值，rj(i)为过零区间极值点。所述基于去噪本征模函数分量的信号重构计算式为：式中L为EMD处理所获得的所有IMF(包括残余分量)的个数，M2-M1+1为期望的带噪IMF分量个数，s'ˉ为经过消噪处理所得的虚拟源信号估计。所述基于“基-2”同态逆变换的去噪观测信号后处理计算式为：式中s～为经过消噪处理所得的真实源信号估计，为信号后处理规整平滑函数。附图说明图1是本专利技术所述基于经验模态分解改进的乘性噪声去除技术的流程示意图；图2是方差散点数据的高阶多项式回归图；图3是方差散点数据的低阶多项式回归图；图4仿真源及其虚拟带噪观测信号波形图；图5基于SVE回归阈值调整(回归阶次：10)的不同去噪算法对乘性混叠噪声的去噪效果；图6基于RVE回归阈值调整(回归阶次：9)的不同去噪算法对乘性混叠噪声的去噪效果；图7基于SVE回归阈值调整(回归阶次：2)的不同去噪算法对乘性混叠噪声的去噪效果；图8基于RVE回归阈值调整(回归阶次：2)的不同去噪算法对乘性混叠噪声的去噪效果；图9不同EMD基去噪算法对IMF1噪声能量的估计。具体实施方式下面结合附图对本专利技术进行具体描述，如图1是本专利技术所述基于经验模态分解改进的乘性噪声去除技术的流程示意图，采用基于谐波信号及其同态变型的噪声能量分布特性分析方法，对本征模函数幅值滤波阈值准则进行调整，并基于“基-2”同态变换预处理去除乘性噪声混叠模型中源信号与噪声信号的相倚性，以基于经验模态分解的信号去噪处理作为主去噪算法进行乘性噪声的消除。本技术方案以添加不同程度的均匀分布乘性白噪声的“撞击(Bumps)”混叠带噪观测信号的消噪为例子阐述基于经验模态分解改进的乘性噪声去除的过程，其基本消噪原理为：分别对乘性混叠带噪观测集合{x}及其同态变型集合log2({x})进行EMD分解，按照标准方差估计(SVE)和鲁棒方差估计(RVE)两种方法估计所抽取的第一阶IMF分量的方差，并分别进行低阶和高阶多项式回归分析以求解方差模型的参数，如图2和图3。即式中σ2|S(Ii)、σ2|R(Ii)分别为EMD抽取的第i个IMF分量Ii的SVE和RVE估计值。实施例1在一个“撞击(Bumps)”源信号s(t)中，通过添加不同程度的均匀分布乘性白噪声获得具有不同信噪比(SNR)的乘性混叠带噪观测信号集合{x(t)}，如图4所示。其中，噪声分量集合{u}＝(12{σu})1/2(u0-μu)，其中均值μu＝0，标准差集合{σu}＝{0.5i}，i＝1,2,...,30。实验数据样本的长度集合{L}＝2N，N＝9,10,…,14。采样频率Fs＝2048Hz。分别应用不同去噪算法对含有不同水平乘性噪声的带噪“Bumps”信号进行去噪处理效果，结果如图5到图8所示。数据样本长度L＝211＝2048。除了改进的EMD基乘性噪声去除算法以外，还对比了小波基去噪算法以及K-M算法，分别以HEMD、WT和K-M来表示。其中的HEMD算法又分为两类，一类是采用调整后的IMF幅值滤波阈值准则，以HEMD-A表示；另一类本文档来自技高网...

【技术保护点】
1.一种基于经验模态分解改进的乘性噪声去除技术，其特征在于，该方法包括如下步骤：1)谐波信号及其同态变型的噪声能量分布特性分析；2)本征模函数幅值滤波阈值准则的调整；3)“基‑2”同态变换预处理；4)基于经验模态分解的信号去噪处理；5)基于去噪本征模函数分量的信号重构；6)基于“基‑2”同态逆变换的去噪观测信号后处理。

【技术特征摘要】
1.一种基于经验模态分解改进的乘性噪声去除技术，其特征在于，该方法包括如下步骤：1)谐波信号及其同态变型的噪声能量分布特性分析；2)本征模函数幅值滤波阈值准则的调整；3)“基-2”同态变换预处理；4)基于经验模态分解的信号去噪处理；5)基于去噪本征模函数分量的信号重构；6)基于“基-2”同态逆变换的去噪观测信号后处理。2.根据权利要求1所述的基于经验模态分解改进的乘性噪声去除技术，其特征在于，所述谐波信号及其同态变型的噪声能量分布特性分析计算式为：式中x为带噪观测信号向量，xi为第i个信号分量，E(·)为数学期望估计函数，σ2|S(·)为标准方差估计函数，median(·)为中值函数，I1为第1个本征模函数(IMF)分量，σ2|R(·)为鲁棒方差估计函数，N为样本数。3.根据权利要求1所述的基于经验模态分解改进的乘性噪声去除技术，其特征在于，所述本征模函数幅值滤波阈值准则的调整计算式为：式中E1为原始乘性带噪信号经过经验模态分解(EMD)抽取的第1阶IMF分量的噪声能量的估计，E1′为乘性带噪信号经过同态变换解耦后由EMD抽取的第1阶IMF分量的噪声能量的估计，PK(·)为基于方差散点数据学习得到的K阶回归多项式，ak,k＝0,1,...,K为多项式系数，Ei′为乘性带噪信号经过同态变换解耦后由EMD抽取的第i阶IMF分量的噪声能量估计，Ti′为幅值滤波阈值，N为信号样本点数，C为常数通常取为C＝0.7，β与ρ均为能量估算参数。4.根据权...

【专利技术属性】
技术研发人员：焦卫东，杨志强，
申请(专利权)人：浙江师范大学，
类型：发明
国别省市：浙江,33