请教各位前辈。。急。。

youshen

1》请问EMD分解得到IMF分量以后，是不是要安装个HHT谱变换工具箱才能出HHT谱图的？？？
2》关于论坛里面有位前辈的程序举例：
t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
plot(t,z)

imf=emd(z);
emd_visu(z,1:length(z),imf)

[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))


t=1:500;
t=t*1/2000;      
这个两句程序是什么意思？请哪位前辈帮解释下啊。

youshen

先谢谢各位了额。。

ChaChing

逐步执行不就知道了? (将;改为,)

blackiris

2》关于论坛里面有位前辈的程序举例：
t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
plot(t,z)

imf=emd(z);
emd_visu(z,1:length(z),imf)

[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))


t=1:500;
t=t*1/2000;      
这个两句程序是什么意思？请哪位前辈帮解释下啊。

这两句应该是抽样的，抽500个点，时间是从0到1/4

blackiris

1》请问EMD分解得到IMF分量以后，是不是要安装个HHT谱变换工具箱才能出HHT谱图的？

你要是安的是Rilling的工具箱，里面本身就有计算谱和画谱的函数。[A,f,tt]=hhspectrum(imf);
是计算谱的。disp_hhs(im);是画的，这句后面的那句你要是注释掉就能看到彩色谱图

youshen

明白了。。我少了instfreq函数。。
下了个时频工具箱就出谱图了。
谢谢前辈。。

youshen

我还有一个问题：请问前辈 这句   emd_visu(z,1:length(z),imf)    能降解一下吗？谢谢啊。。

blackiris

  emd_visu(z,1:length(z),imf)
z是原始信号，1:length(z)是x轴坐标（这句话就是以1为步长取了length(z)个点），imf是个矩阵，里面有imf分量及余量，这条语句是显示emd分解结果的，应该会出现三个图，第一个是分解结果，后两个分别是由高频到低频，由低频到高频的重建信号的过程图。
我也是新学这块内容的，彼此交流吧~

youshen

回复 blackiris 的帖子

谢谢你啊。。

youshen

回复 blackiris 的帖子

load('data.txt');

tspan=30;
%时长

fs=50;
%采样频率

t=(1/fs:1/fs:tspan);

N=length(t);
%采样点数

x=data(1:N,2);

figure(10);

plot(t,x);

[imf,ort,nbits]=emd(x);

%emd_visu(x,1:length(x),imf)

%%% 求边际谱

[A,fa,tt]=hhspectrum(imf(1:end-1,:));

[E,tt1,ff]=toimage(A,fa,tt,length(tt));

disp_hhs(E);colorbar;

% colormap(flipud(gray))

for k=1:size(E,1)

bjp(k)=sum(E(k,:))*1/fs*1/tspan;

end

f=(0:N-3)/N*(fs/2);

figure(60)

plot(f,bjp);

xlabel('频率 / Hz');

ylabel('·幅值');

请问这句 “ x=data(1:N,2); 是什么意思啊？

youshen

这是一个数据处理的HHT程序。

blackiris

回复 youshen 的帖子

请问这句 “ x=data(1:N,2); 是什么意思啊？

data应该是个矩阵型的数据，1:N是指在这个矩阵中取行的第一到第N个数据，2是指取矩阵的第2列。你可以看看MATLAB的基础，这些是MATLAB的基础语句。

youshen

回复 blackiris 的帖子

t=(1/fs:1/fs:tspan);

请问这一句是什么意思？ 应该是控制时间的吧？

blackiris

回复 youshen 的帖子

对，是时间轴，指的是从1/fs开始，以1/fs为步长，取到tspan为止。针对这个程序，得到的是从0.02开始，以0.02为步长，一直取到30，共1500个点的向量

youshen

回复 blackiris 的帖子

非常感谢。。今天跟导师讨论了下。。初步做出了想要的图。。谢谢你啊，帮了我很多。。以后多讨论讨论啊。。