vmd参数优化 matlab,遗传算法优化VMD参数

我想用遗传算法优化两个参数，这个两个参数是求解VMD分量必要设定的参数，一个惩罚因子α，一个是模态分解分量K，这两个参数以往是经验值，现在我想以信息熵为目标函数，用遗传算法优化，求所有α和K的组合得到的所有分量中最小的信息熵，此时得到的α和K值是多少，运行MATLAB程序总是出错，请大神帮我看看代码。

% Time Domain 0 to T

N = 1000;

fs = 1/N;

t = (1:N)/N;

freqs = 2*pi*(t-0.5-1/N)/(fs);

% center frequencies of components

f_1 = 2;

f_2 = 24;

f_3 = 288;

% modes

v_1 = (cos(2*pi*f_1*t));

v_2 = 1/4*(cos(2*pi*f_2*t));

v_3 = 1/16*(cos(2*pi*f_3*t));

% for visualization purposes

wsub{1} = 2*pi*f_1;

wsub{2} = 2*pi*f_2;

wsub{3} = 2*pi*f_3;

% composite signal, including noise

x = v_1 + v_2 + v_3 + 0.1*randn(size(v_1));

for a=0:10:200

%     alpha = a

for b=1:10

%