请教HHT做边际谱的问题

Legendary

我在使用HHT工具箱进行数据处理时出现了以下的问题，不知如何解决，特向各位请教！
源程序：
%emd分解
imf=emd(y_trimtrend);           % G. Rilling程序,2006年版
emd_visu(y_trimtrend,t,imf,4);    % G. Rilling程序,2006年版
                                % figure(4);figure(5);figure(6);
% 边际谱，用论坛中破凰的程序，作了小的修改 figure(7);
IMF=imf(1,:);
[row,col]=size(IMF);
N=col;
nIMF=row;
[A,fa,tt]=hhspectrum(IMF(1:nIMF,:));    %Hilbert-Huang变换，计算瞬时频率和瞬时幅值
[E,tt1]=toimage(A,fa,tt,length(tt));
NN=size(E,1);
for k=1:NN
    bjp(k)=sum(E(k,:))*1/NN;
end
f=(0:N-3)/N*(sf/2);
figure;
plot(f,bjp);
xlabel('频率(Hz)');
ylabel('幅值');

% Hilbert-Huang时频谱 figure(8)
[A,f,tt]=hhspectrum(imf(1,:));
[im,tt]=toimage(A,f);
disp_hhs(im,tt,[],sf);
ylim([0,0.3]);xlabel('时间(s)');ylabel('频率(Hz)');
运行时出现如下错误：
Error in ==> hhspectrum at 20
if  nargin < 2

??? Output argument "A" (and maybe others) not assigned during call to
"C:\Users\YIN_LAB\Desktop\hht工具箱\hht工具箱\hhspectrum.m>hhspectrum".

cwb

你好，可以解释一下toiamge得到的结果中E到底是什么意思吗？它的行和列是受toiamge函数中的参数个数控制的，我一直都没搞明白为什么默认的行是400行。。。如果加上tt和length(tt)的话又出错，，，还有为什么对E的每一行进行累加得到的就是边际谱，难道E的每一行就代表某个频率在整个时间跨度上的幅度值吗？
如能得到解答，感激不尽，谢谢！！！