ie参数不足咋么办

15646

function [N,IE]=ie(x,n,dt,fs)
%x 原始信号
%n采样点数
%dt采样周期
%fs=1/dt;%采样频率
imf=emd(x);
%[S,freq]=hhte(imf,2) %Y.Lzhang
[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
for k=1:size(E,2)
    IE(k)=sum(E(k,:))^2*2*pi*fs;  %瞬时能量谱
end
N=(0:n-3)/n


ie输入参数的数目不足。