对于数据点的HHT处理

仍小呆

现在有一些仿真测量的数据点，不知道如何进行HHT处理，希望请假大家，能不能给我帮帮忙。希望能够得到EMD分解图、HHT边际谱和HHT谱，谢谢大家了，最好能传我一份程序，让我学习下，再次感谢。。。数据点如下：

时间	力
13.60000038146973	0.08494747430086136
13.80000019073486	0.08880482614040375
14.00000095367432	0.09470564126968384
14.20000076293945	0.09777642041444778
14.40000057220459	0.09904648363590241
14.60000038146973	0.1084610298275948
14.80000019073486	0.1059361919760704
15.00000095367432	0.1063350215554237
15.20000076293945	0.103964664041996
15.40000057220459	0.1006118804216385
15.60000038146973	0.09639190882444382
15.80000114440918	0.09373065829277039
16.	0.09149065613746643
16.20000076293945	0.08934074640274048
16.40000152587891	0.0897015780210495
16.60000038146973	0.08465103060007095
16.80000114440918	0.08082866668701172
17.	0.07555326819419861
17.20000076293945	0.07031608372926712
17.40000152587891	0.06633331626653671
17.60000038146973	0.06407567113637924
17.80000114440918	0.06263621151447296
18.	0.06185589358210564
18.20000076293945	0.05966390296816826
18.40000152587891	0.05807501077651978
18.60000038146973	0.05497178062796593
18.80000114440918	0.05293424427509308
19.	0.0515085868537426
19.20000076293945	0.05191246047616005
19.40000152587891	0.05141106620430946
19.60000038146973	0.05246229097247124
19.80000114440918	0.05041955038905144
20.	0.04826271533966065
20.20000076293945	0.04873500764369965
20.40000152587891	0.04196785390377045
20.60000038146973	0.04535799100995064
20.80000114440918	0.0469869002699852
21.	0.04668370261788368
21.20000076293945	0.04857110604643822
21.40000152587891	0.04905537888407707
21.60000038146973	0.04991570115089417
21.80000114440918	0.04933634400367737
22.	0.05095796659588814
22.20000076293945	0.05239294841885567
22.40000152587891	0.05528919026255608
22.60000038146973	0.0580727644264698
22.80000114440918	0.06177578866481781
23.20000076293945	0.06557949632406235
23.39999961853027	0.06481660157442093
23.60000038146973	0.06604308634996414
23.79999923706055	0.06289757043123245
24.	0.05920350551605225
23.20000076293945	0.06557949632406235
23.39999961853027	0.06481660157442093
23.60000038146973	0.06604308634996414
23.79999923706055	0.06289757043123245
24.	0.05920350551605225
24.20000076293945	0.05572297796607018
24.39999961853027	0.05338825657963753
24.60000038146973	0.05304373800754547
24.79999923706055	0.05144522711634636
25.	0.05011733993887901
25.20000076293945	0.04753215238451958
25.39999961853027	0.04512707144021988
25.60000038146973	0.04403987526893616
25.79999923706055	0.04444540292024612
26.	0.0456012450158596
26.20000076293945	0.04656370356678963
26.39999961853027	0.04687913879752159
26.60000038146973	0.04677524045109749
26.79999923706055	0.04610222205519676
27.	0.04593367502093315
27.20000076293945	0.04659835249185562
27.39999961853027	0.04792125150561333
27.60000038146973	0.04947065934538841
27.79999923706055	0.05003156512975693
28.	0.04995182156562805
28.20000076293945	0.04851961508393288
28.39999961853027	0.04638580232858658
28.60000038146973	0.04429491981863976
28.79999923706055	0.04338868707418442
29.	0.04353157803416252
29.20000076293945	0.04445184022188187
29.39999961853027	0.044819675385952
29.60000038146973	0.04414425417780876
29.79999923706055	0.04241181537508965
30.	0.04116878286004067
30.20000076293945	0.03992815688252449
30.39999961853027	0.04145978763699532
30.60000038146973	0.04323955997824669
30.79999923706055	0.04531598091125488
31.	0.04667390882968903
31.20000076293945	0.04670368134975433
31.40000152587891	0.04582234844565392
31.60000038146973	0.04509157314896584
31.79999923706055	0.04551892727613449
32.	0.04710102453827858
32.20000076293945	0.04918133839964867
32.40000152587891	0.04964540898799896
32.59999847412109	0.04985122010111809
32.79999923706055	0.04857401922345161
33.	0.04769911989569664
33.20000076293945	0.04744653403759003
33.40000152587891	0.04861631244421005
33.59999847412109	0.04970131441950798
33.79999923706055	0.05021323263645172
34.	0.0511254221200943
34.20000076293945	0.05019605532288551
34.40000152587891	0.05007244274020195
34.59999847412109	0.05038589239120483
34.79999923706055	0.05099736154079437
35.	0.05113137513399124
35.20000076293945	0.05084780976176262
35.40000152587891	0.04914001002907753
35.59999847412109	0.04875000566244125
35.79999923706055	0.0484510101377964
36.	0.04860669746994972
36.40000152587891	0.04912392422556877
36.60000228881836	0.0492144376039505
36.79999923706055	0.04875800013542175
37.	0.04780564457178116
37.20000076293945	0.04549027606844902
37.40000152587891	0.04567569866776466
37.60000228881836	0.04448513314127922
37.79999923706055	0.04700591415166855
38.	0.04689091816544533
38.20000076293945	0.05078794062137604
38.40000152587891	0.04911050200462341
38.60000228881836	0.05094971507787705
38.79999923706055	0.04958949238061905
39.	0.04861272498965263
39.20000076293945	0.04762743413448334
39.40000152587891	0.05050629004836083
39.60000228881836	0.05171589180827141
39.79999923706055	0.04942487925291061
40.	0.05711781978607178
40.20000076293945	0.05366618558764458
40.40000152587891	0.05609699338674545
40.60000228881836	0.05420174822211266
40.79999923706055	0.05991063266992569
41.	0.05920858308672905

airazor

clear all;
close all;
clc
a=xlsread('1.xls')';
imf=emd(a(2,:));
emd_visu(a(2,:),a(1,:),imf)

你这个数据有问题，有一个时间等于24的好像是多余的。剩下的自己搞吧。emd是HHT的一部分。

仍小呆

airazor 发表于 2015-1-5 19:34
clear all;
close all;
clc

好的，谢谢你了，非常感谢

仍小呆

airazor 发表于 2015-1-5 19:34
clear all;
close all;
clc

还有一个问题，得到的IMF分量不显示采样时间啊，而且能不能帮我得出边际谱，谢谢

chybeyond

仍小呆 发表于 2015-1-6 16:53
还有一个问题，得到的IMF分量不显示采样时间啊，而且能不能帮我得出边际谱，谢谢

自己动手搜一下吧，论坛很多相关的例子http://forum.vibunion.com/forum.php?mod=viewthread&tid=83842

仍小呆

chybeyond 发表于 2015-1-7 09:38
自己动手搜一下吧，论坛很多相关的例子http://forum.vibunion.com/forum.php?mod=viewthread&tid=83842

好的。谢谢你的帮助

airazor

仍小呆 发表于 2015-1-6 16:53
还有一个问题，得到的IMF分量不显示采样时间啊，而且能不能帮我得出边际谱，谢谢

大哥，你不看文献不行啊，要对EMD有个认识的。建议你最好自己编一下EMD程序，就对EMD有个概念性的理解了。。。还有看一下matlab。。。回答：EMD的过程是对一个1×n的序列进行分解，所得到的为m行、长度为n的序列，不牵扯到横坐标的问题。至于时间轴，你在画的时候横坐标弄一下就是时间轴。例如：
t=0：0.01：10；
y=sin(t);
plot(t,y)%这个横坐标时间轴是用t加上去的，和你这个需要的时间轴类似。这句话要是写成plot(y)，横坐标就不一样了。
至于你用了emd_visu这个程序的话，将里面的东西改一下，就可以得到时间轴了。
哈哈。。。让大神见笑了

airazor

chybeyond 发表于 2015-1-7 09:38
自己动手搜一下吧，论坛很多相关的例子http://forum.vibunion.com/forum.php?mod=viewthread&tid=83842

高人啊。你给的这个边际谱是正确的。破凰的那个边际谱有问题，还被那个大神说是正确的，但是没人发现，现在网上都是那个坏程序，影响极为不好。这贴该顶。

仍小呆

airazor 发表于 2015-1-19 17:02
大哥，你不看文献不行啊，要对EMD有个认识的。建议你最好自己编一下EMD程序，就对EMD有个概念性的理解了 ...

好的，谢谢你了。。一开始学的是都是对一个仿真信号进行处理，例如什么sin、cos，结果现在一下子对数据点进行处理，就有点不知所措了，谢谢你的指导