企业级网站开发平台wordpress king模板

当前位置： 首页 > news >正文

企业级网站开发平台,wordpress king模板,怎么设计logo用什么软件,详细介绍网站建设的整个流程今天来介绍一种效果颇为不错的降噪方法。#xff08;针对高频白噪声#xff09; 上一篇文章我们讲到了FastICA方法。在现实世界的许多情况下#xff0c;噪声往往接近高斯分布#xff0c;而有用的信号#xff08;如语音、图像特征等#xff09;往往表现出非高斯的特性。F…今天来介绍一种效果颇为不错的降噪方法。针对高频白噪声 上一篇文章我们讲到了FastICA方法。在现实世界的许多情况下噪声往往接近高斯分布而有用的信号如语音、图像特征等往往表现出非高斯的特性。FastICA通过最大化输出信号的非高斯性来恢复这些有用的信号从而有效地从噪声中分离出信号。 将类经验模态分解EMD方法与FastICA结合使用可以创建一种巧妙的信号处理策略这种结合利用了两种方法的互补优势通过EMD分解得到的IMFs可以简化信号的结构但单独使用EMD可能不足以有效分离信号中的噪声和有用信息自适应得到的imf数量也可能存在冗余FastICA则能够从imf信号中进一步提取独立成分更有效地识别噪声成分。 下面将详细解释这种结合的算法流程、优势以及MATLAB代码实现。 一、算法流程 我们直接开门见山通过案例来讲解算法流程。 首先生成一段待降噪信号。 rng(123456) t (0:0.001:(1-0.001)); x1 0.6*sin(15*pi*tpi/5); x2 cos(60*pi*tsin(10*pi*t)); x3 (10.3*cos(10*pi*t)).*sin(200*pi*t); x4 wgn(1000,1,-10);fs 1000; %采样频率 ps x1x2x3; %未加噪声的纯净信号 s psx4; figure(color,w) subplot(2,1,1);plot(t,ps,k);title(无噪声的仿真信号) subplot(2,1,2);plot(t,s,k);title(加入噪声的仿真信号) 这段信号加入噪声前后的图像如下我们的任务就是将下边的含噪信号尽可能回复成上图中的无噪声信号。 1.类EMD分解 关注我的专栏的老读者们都知道我将EMD、EEMD、CEEMD、CEEMDAN、ICEEMDAN、VMD等一系列模态方法统称为“类EMD”方法在本篇文章中该步骤也可以是上述分解方法中的任意一种毕竟这一步的目的是将复杂信号分解为一定数量的模态分量。 此案例中我们使用VMD分解算法。 alpha1000; % alpha - 惩罚因子 tol1e-6; % tol - 收敛容差是优化的停止准则之一可以取 1e-6~5e-6 K5; % K - 指定分解模态数 imf pVMD(s,fs, alpha, K, tol); 上边代码中使用了pVMD函数这是笔者封装的一个易用的VMD画图函数其介绍可以看之前的文章。分解结果如下 也可以在调用pVMDandFFT函数绘制imf分量与频谱对比图。 imf pVMDandFFT(s,fs, alpha, K, tol); %函数见此处https://zhuanlan.zhihu.com/p/396775790 当然使用MATLAB内置的vmd函数也是可以的但是要注意输入FastICA的imf分量的行列方向需要保证imf是每一行为一个分量。 2.FastICA盲源分离 以类EMD分解得到的imf分量作为输入进行ICA混解。 在这里一定要算出混合矩阵A忘记概念的同学可以看上一篇。 %% 3.ICA混解并绘制混解后的图像 numOfIC 0; % 需要提取的独立成分数目如果不指定数目则输入0 g pow3; % 使用的非线性函数类型可选pow3, tanh, gauss [icasig, A, W] pFastICA(imf, numOfIC, g); 此时我们将会得到如下分解结果 从FastICA盲源分离结果可以推断独立成分5是噪声分量但是在自动化实现降噪的代码中我们自然希望对噪声分量的判断是不需要人为介入的。此时我们引入一个功率谱熵阈值判断方法来实现。 3.功率谱熵阈值判断 在之前的文章Mr.看海【熵与特征提取】基于“信息熵”的特征指标及其MATLAB代码实现功率谱熵、奇异谱熵、能量熵中讲到信息熵越大代表不确定性越大信号中包含的信息量越少信号越无秩序/越接近于白噪声则信息熵越大。 由于盲源分离本身的特性分离出来的独立成分中最多包含一个高斯分量也就是噪声分量。 所以我们只需要判断几个独立成分的功率谱熵中是否有数值相较于其他分量大得多的成分就是噪声分量。 也就是找到功率谱熵中的离群值。 这里我们阈值的计算标准是阈值功率谱熵均值功率谱熵的标准差 我们将大于阈值的独立成分直接赋值为0。 计算结果如下图可以看到噪声分量被识别了出来。 4.重构滤波信号 不知道大家是否还记得FastICA方法的重要性质之一输出信号幅度的不确定性。FastICA分解得到的独立成分是不能像类EMD分解得到的imf那样直接相加进行重构的。 但是FastICA也有自己的重构方法就是使用上边提到的混合矩阵A。重构公式为 XA∗icasig 上式中的icasig就是经过处理过后的独立成分对于大于阈值的成分我们已经将其幅值全部置0X就是重构结果也就是该方法中的滤波结果。 此时得到的滤波结果为 可以看出滤波效果还是颇为不错的经计算此时的信噪比为15.9dB。 二、算法流程图 将上述算法梳理成流程图如下当然vmd分解这步是可以替换成其他分解算法的另外功率谱熵也可以换成其他熵值指标。 三、MATLAB实现 为了方便大家使用按照惯例笔者对上述流程进行了封装。封装的流程就是上述流程图中红框内的部分之所以没有将vmd也一并封装进去是为了大家方便替换其他的模态分解算法。 封装后的函数只需要将分解好的imf变量导入即可一行代码实现滤波 reSig filEMDsICA(imf); % 调用filEMDsICA函数实现滤波 经实测这个方法的滤波效果还是很不错的。 我为大家提供了完整的演示案例可以一键实现下述图像和滤波结果的输出 要你做的基本就只需替换数据啦 需要上边这个函数文件以及测试代码的同学可以在公众号 khscience看海的城堡中回复“ICA滤波”获取。 扩展阅读 Mr.看海【滤波专题-第1篇】数字滤波器15分钟入门——这可能是最简单的FIR有限冲激响应滤波讲解 Mr.看海【滤波专题-第2篇】数字滤波器15分钟入门——这可能是最简单的IIR无限冲激响应滤波讲解 Mr.看海【滤波专题-第3篇】IIR无限冲激响应和FIR有限冲激响应数字滤波器有什么区别 Mr.看海【滤波专题-第4篇】滤波器滤波效果的评价指标信噪比SNR、均方误差MSE、波形相似参数NCC Mr.看海【滤波专题-第5篇】FIR、IIR滤波器设计及MATLAB实现 Mr.看海【滤波专题-第6篇】小波阈值去噪方法看这一篇就明白了~附MATLAB实现 Mr.看海【滤波专题-第7篇】“类EMD”算法分解后要怎样使用3——EMD降噪方法及MATLAB代码实现 Mr.看海【盲源分离】快速理解FastICA算法附MATLAB绘图程序