function [fitness] = fun(x,f,tau, DC, init, tol)
x= round(x);
alpha = x(1); % moderate bandwidth constraint
K = x(2);
%% 包络熵
[u, ~,~] = VMD(f, alpha, tau, K, DC, init, tol);
p=zeros(size(u));
E = zeros(size(u,1),1);
for i = 1:K
a = abs(hilbert(u(i,:))); % hilbert解调,abs() 函数返回数字的绝对值。
智能算法及其模型预测