Hilbert - Huang - guide
更新时间:2023-11-02 19:53:01 阅读量: 综合文库 文档下载
- hilbert空间推荐度:
- 相关推荐
A few examples of how to use these programs are given, with a given dataset “gsta.dat”, which is the annual mean global surface temperature anomaly. In “gsta.dat”, the first column is the time; and the second is the corresponding data value.
1) Load and display data
>> load gsta.dat;
>> plot(gsta(:,1),gsta(:,2)); >> axis([1850 2010 -0.6 0.6]);
>> title('the annual mean global surface temperature anomaly') >> xlabel('year') >> ylabel('Kelvin')
2) Using the program as a traditional EMD tool
The eemd.m can be used as an EMD decomposition tool:
>> year=gsta(:,1); >> inData=gsta(:,2); >> rslt=eemd(inData,0,1); >> plot(year,rslt(:,2)); >> hold on;
>> plot(year,rslt(:,3)-0.3); >> plot(year,rslt(:,4)-0.6); >> plot(year,rslt(:,5)-0.9);
>> plot(year,sum(rslt(:,6:8),2)-1.3,'r-'); >> hold off
>> set(gca,'yTickLabel',[]); >> axis([1850 2010 -1.8 0.3]); >> xlabel('year');
“relt” is a matrix containing the decomposition results, with the first column the original input (“gsta(:,2)”) and the last column the trend. In between are the IMFs of frequencies from high to low.
It should be noted that since in the eemd.m, the total number of IMFs m is specified as log2(N)-1, in some occasions (such as the one in this example), the components may be excessively extracted. In these cases, the sum of the latest columns may already satisfy the definition of a trend.
3) Instantaneous frequency of an IMF
The instantaneous frequency can be obtained through calling ifndq.m function:
>> omega_m3=ifndq(rslt(:,4),1); >> subplot(2,1,1); >> plot(year,rslt(:,4));
>> axis([1850 2010 -0.12 0.12]); >> title('IMF C3'); >> ylabel('Kelvin'); >> grid;
>> subplot(2,1,2);
>> plot(year, omega_m3/2/pi,'r-'); >> grid;
>> xlabel('year');
>> ylabel('cycle/year');
>> title('instantaneous frequency'); >> axis([1850 2010 0 0.12]);
It should be noted that the instantaneous frequency calculation program is not suitable for under sampled oscillations, such as the first IMF (with an averaged period about 3 data points). However, for such under sampled oscillations, the instantaneous frequency is no longer “instantaneous” any way, and any method used to obtain such a quantity will have big errors.
4) Using the program as a EEMD tool
The eemd.m can be used as an EEMD decomposition tool. In this case, the noise assed has an amplitude (standard deviation) of 0.2 of the standard deviation of the linearly detrended annual mean global surface temperature anomaly; and the number of ensemble is 100:
>> rslt=eemd(inData,0.2,100);
>> t(1)=1850; >> t(2)=2010; >> y1(1)=0; >> y1(2)=0; >> y2(1)=-0.3; >> y2(2)=-0.3; >> y3(1)=-0.6; >> y3(2)=-0.6; >> y4(1)=-0.9; >> y4(2)=-0.9; >> y5(1)=-1.2; >> y5(2)=-1.2; >> y6(1)=-1.6; >> y6(2)=-1.6; >> plot(t,y1,'k-'); >> hold on;
>> plot(t,y2,'k-'); >> plot(t,y3,'k-'); >> plot(t,y4,'k-'); >> plot(t,y5,'k-'); >> plot(t,y6,'k-'); >> plot(year,rslt(:,1)); >> plot(year,rslt(:,3)-0.3); >> plot(year,rslt(:,4)-0.6);
>> plot(year,rslt(:,5)-0.9); >> plot(year,rslt(:,6)-1.2);
>> plot(year,sum(rslt(:,7:8),2)-1.6,'r-'); >> set(gca,'yTickLabel',[]);
>> title('EEMD decomposition of GSTA (A_n=0.2; N_e_s_b=100)') >> axis([1850 2010 -2.1 0.2]); >> xlabel('year');
5) Statistical significance test
Since the annual mean global surface temperature anomaly behaves completely different from a white noise series, we use computer generated white noise to illustrate how the significance.m can be used:
>> clear; >> clf;
>> data=randn(512,1); >> rslt=eemd(data,0,1); >> imfs=rslt(:,2:8);
>> [sigline95,logep]=significance(imfs,0.05); >> [sigline99,logep]=significance(imfs,0.01); >> plot(sigline95(:,1),sigline95(:,2)); % 95 percenta line >> hold on
>> plot(sigline99(:,1),sigline99(:,2),'m-'); % 99 percenta line
>> plot(logep(:,1),logep(:,2),'r*'); >> plot(logep(1,1),logep(1,2),'k*'); >> grid;
>> xlabel('LOG2 ( Mean Period )');
>> ylabel('LOG2 ( Mean Normalized Energy )');
>> title('Significance test of IMFs of white noise'); >> axis([0 10 -7 0])
3) Trend and detrending
For example, in the previous decomposition, the sum of the last three columns satisfies the definition of trend well.
>> plot(year, rslt(:,1)); >> hold on;
>> plot(year, sum(rslt(:,7:8),2),'r-'); >> plot(year, sum(rslt(:,6:8),2),'g-'); >> plot(year, sum(rslt(:,5:8),2),'m-'); >> title('Trends of different timescales'); >> ylabel('Kelvin'); >> xlabel('year'); >> grid;
正在阅读:
精编完整版美伦河、澄江两岸带状公园景观项目工程可研报告03-16
23.带上她的眼睛 名师教案 新部编人教版七年级下册语文05-21
广东省茂名市2013届高三第一次高考模拟数学理试题及答案09-17
英文原版McGraw-Hill Construction ENR-TOP-250Rank07-19
幼儿园体育游戏10-17
DOS 概述及入门- 副本- 副本03-11
第二章 生命的基本单位——细胞习题及答案07-27
银鹭品牌策划案 - 图文12-24
资助与育人学生座谈会03-13
- 1MongoDB_Architecture_Guide
- 2第4章 希尔伯特(Hilbert)空间
- 3administration guide 11g
- 4Cisco CP Express User Guide
- 5plsql guide(基础详解)
- 6VxWorks kernal programmers guide
- 7Guide to the txtile labelling and adversing regulations
- 8PLSQL User's Guide and Reference
- 9Clean - Oil - Guide 油品清洁指南
- 10JDE9.1 for Windows Installation Guide
- 发电电气运行规程1
- 英文简历
- 最全辅导员招聘考试题库
- 4.3崇明岛的未来的样子
- 2012年上海市普通高校招生二本批次各校投档分数线
- 江苏省如皋中学2017-2018学年第一学期高三第二次阶段测试12月数
- 农业转移人口社会参与机制浅谈
- 2017-2018学年度牛津译林版8B英语初二期中试卷及答案
- 家长委员会上的讲话
- 05继电保护设备检修规程
- 组织行为学考试重点(陈春花)
- 2016年云南省公务员考试《行测》模拟试卷(十七)
- 规避“10号文”红筹系列之案例分析
- 钱寨小学学生读书活动评价方案
- 五大联赛派系
- 国际结算课件新
- 材料科学导论 - 图文
- 领导干部任前廉政法规考试模拟试题
- 汽车综合实训
- 医疗质量管理目录
- Hilbert
- Huang
- guide
- 论侦查人员的证据意识
- 学习hypermesh第一步:面板详解 - 图文
- 弥兴中学课例研修资料汇编(学校)123 - 图文
- An Analysis of Jane Eyre
- 水分析化学第三版课后习题答案
- 医学影像物理习题
- 视觉传达设计史复习题
- 云浮中学100周年校庆工作总结--钟
- 百货商场6周年店庆活动方案
- 2019年高考英语专题 阅读理解七选五-2019年高考英语抢分秘籍 Word版含解析
- 基于单片机AT89C51的无线多路灯具控制系统
- 氩弧焊简易接法及质量控制
- 舵系作业指导书 - 图文
- 上海中学顾滨
- 2016冶炼专业监理工程师继续教育考试卷(90分)
- 市海事局2012年1至5月依法行政工作情况汇报
- 流体力学讲义 第六章 流动阻力及能量损失2汇总
- 2012考试大纲答案
- 避雷器型号含义
- 2012-2013学年第一学期九年级物理期中测试卷及答案