matlab使用经验模式分解emd对信号进行去噪
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
matlab使⽤经验模式分解emd对信号进⾏去噪
:
对于这个例⼦，考虑由具有明显频率变化的正弦波组成的⾮平稳连续信号。
⼿提钻的振动或烟花声是⾮平稳连续信号的例⼦。
以采样频率加载⾮平稳信号数据fs，并可视化混合正弦信号。
1.
load（'sinusoidalSignalExampleData.mat'，'X'，'fs'）;
2.
3.
xlabel（'Time（s）'）;
观察到混合信号包含具有不同幅度和频率值的正弦波。
为了创建希尔伯特谱图，您需要信号的IMF。
执⾏经验模式分解以计算信号的固有模式函数和残差。
由于信号不平滑，请指定' pchip'作为Interpolation⽅法。
[imf，residual，info] = emd（X，'Interpolation'，'pchip'）;
1.
⽬前的IMF | #Sift Iter | 相对Tol | 停⽌标准命中
2.
1 |
2 | 0.026352 | SiftMaxRelativeTolerance
3.
2 | 2 | 0.003957
3 | SiftMaxRelativeTolerance
4.
3 | 1 | 0.024838 | SiftMaxRelativeTolerance
5.
4 | 2 | 0.05929 | SiftMaxRelativeTolerance
6.
5 | 2 | 0.11317 | SiftMaxRelativeTolerance
7.
6 | 2 | 0.12599 | SiftMaxRelativeTolerance
8.
7 | 2 | 0.13802 | SiftMaxRelativeTolerance
9.
8 | 3 | 0.15937 | SiftMaxRelativeTolerance
10.
9 | 2 | 0.15923 | SiftMaxRelativeTolerance
11.
分解停⽌是因为残留信号的极值数⼩于'MaxNumExtrema'。
在命令窗⼝中⽣成的表指⽰每个⽣成的IMF的筛选迭代次数，相对容差和筛选停⽌标准。
此信息也包含在info。
您可以通过指定Display为隐藏表0。
HHT（IMF，FS）;
频率对时间图是⼀个稀疏图，其中垂直颜⾊条表⽰IMF中每个点的瞬时能量。
该图表⽰从原始混合信号分解的每个分量的瞬时频谱。
从该图中可以观察到三个IMF，其频率在1s处有明显变化。
可视化信号的残余和内在模式功能
对于这个例⼦，考虑由具有明显频率变化的正弦波组成的⾮平稳连续信号。
⼿提钻的振动或烟花声是⾮平稳连续信号的例⼦。
加载⾮平稳信号数据，并可视化混合正弦信号。
load（'sinusoidalSignalExampleData.mat'，'X'，'fs'）;
观察到混合信号包含具有不同幅度和频率值的正弦波。
执⾏经验模式分解以绘制固有模式函数和信号残差。
由于信号不平滑，请指定' pchip'作为Interpolation⽅法。
emd（X，'Interpolation'，'pchip'）;
1.
⽬前的IMF | #Sift Iter | 相对Tol | 停⽌标准命中
2.
1 |
2 | 0.026352 | SiftMaxRelativeTolerance
3.
2 | 2 | 0.003957
3 | SiftMaxRelativeTolerance
4.
3 | 1 | 0.024838 | SiftMaxRelativeTolerance
5.
4 | 2 | 0.05929 | SiftMaxRelativeTolerance
6.
5 | 2 | 0.11317 | SiftMaxRelativeTolerance
7.
6 | 2 | 0.12599 | SiftMaxRelativeTolerance
8.
7 | 2 | 0.13802 | SiftMaxRelativeTolerance
9.
8 | 3 | 0.15937 | SiftMaxRelativeTolerance
10.
9 | 2 | 0.15923 | SiftMaxRelativeTolerance
11.
分解停⽌是因为残留信号的极值数⼩于'MaxNumExtrema'。
⽣成具有原始信号，前3个IMF和残差的交互式图。
在命令窗⼝中⽣成的表指⽰每个⽣成的IMF的筛选迭代次数，相对容差和筛选停⽌标准。
您可以通过指定Display为隐藏表0。
右键单击图中的空⽩区域以打开IMF选择器窗⼝。
使⽤IMF选择器有选择地查看⽣成的IMF，原始信号和残差。
从列表中选择要显⽰的IMF。
选择是否在图上显⽰原始信号和残差。
选定的IMF现在显⽰在图上。
使⽤该图可视化从原始信号中分解的各个组件以及残差。
请注意，残差是根据IMF总数计算的，并且不会根据IMF选择器窗⼝中选择的IMF进⾏更改。
⾮常感谢您阅读本⽂，有任何问题请在下⽅留⾔！。