hht入门 碰到点问题 求指导

wzl091002

我刚接触的hht，从网上下了hht工具箱和时频工具箱，拿了论坛以前的一个人的例子进行计算结果总是不对。这是程序：
clear;
fs=1000;    %fs为采样频率;
N=1000;     %采样点数
t=1/fs:1/fs:1;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;  
end
f=(0:N-3)/N*(fs/2);
figure(2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
求大家帮忙看看 ，在你们的matlab中能不能运行，不胜感激啊

wzl091002

没有人吗 求大神指导啊

qqcd

注释掉
E=flipud(E);
能得到准确结果

wzl091002

我再顶一下