資源簡介
經驗模態分解(EMD),對imf進行分析,再采用類小波閾值去噪。
代碼片段和文件信息
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%?????強噪背景下去噪
%%%%%%%%%%??????3014.04.14.15:30
%%%%%%%%%%??功能:emd分解
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear?all;
close?all;
%%?????????讀取原始數據??????????%%%%%%%%%%%%%%%%%%%%%
load?data_203_1m.txt;
original_Data=data_203_1m;%時間數據為10s長
axis?tight;?%使坐標系的最大值和最小值和你的數據范圍一致
t=original_Data(?:?1);?%時間
y=original_Data(?:?2)?;%數據
plot(1000*ty);??
axis([0?10000?-2?2]);
title(‘原始信號圖像‘);?
f=y.‘;
%%?????????EMD?分解??????????%%%%%%%%%%%%%%%%%%%%%???
f1=f([1:10000]);
t1=t([1:10000]);
imf?=?emd(f1);??????%?調用鏡像延拓emd處理函數?[hanglie]?=?size(imf);??%?獲得矩陣imf的m*n
%?emd_visu(f1t1imf)????????%?顯示emd分解的結果
r_y?=?zeros(110000);
for?j=1:1:8??%除去最后的那個剩余分量
????y1(j:)?=?filtsoft(imf(j:));?
????r
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????421598??2014-04-07?16:00??emd\data_203_1m.txt
?????文件???????1851??2014-04-16?20:47??emd\emd_01.m
?????文件????????327??2014-04-16?20:48??emd\filtsoft.m
?????目錄??????????0??2014-04-16?20:48??emd
-----------?---------??----------?-----??----
???????????????423776????????????????????4
- 上一篇:BFGS方法解無約束優化的MATLAB程序
- 下一篇:平衡車MATLAB建模
評論
共有 條評論