hht程序，帮忙翻译一下，一定要具体详细啊。新手菜鸟
本帖最后由 牛小贱 于 2015-6-14 16:12 编辑clear;
N=2048;
t=1:N;
fs=1000;
z=sin(2*pi*40*t/fs)+sin(2*pi*15*t/fs)+sin(2*pi*5*t/fs);
plot(t,z); %做频谱图
xlabel('时间');
ylabel('幅值');
title('原始信号');
imf=emd(z); %EDM分解
cemd_visu(z,1:length(z),imf)
= hhspectrum(imf); % 对IMF分量求取瞬时频率与振幅
= toimage(A,fa,tt,length(tt)); %绘图toimage是把一维的变成二维的显示出来
disp_hhs(E); % disp_hhs(E,[],fs);是非归一化频率
colormap(flipud(gray))
for k=1:size(E,1);
bjp(k) = sum(E(k,:))*1/fs;
end;
f = (0:N-3)/N*(fs/2); %进行对应的频率转换，定义横坐标频率范围的
figure(4) %创建画图窗口
plot(f,bjp); %做频谱图
xlabel('频率');
ylabel('幅值');
title('边际谱'); 包括小参数，也说明一下，感激不尽。 具体的函数可以上网查一下用法，理解的会更深刻。 这不是有工具箱的吗 程序不完整，有些参数不好写。
页:
[1]