fun5.m

来自「该算法是综合遗传算法和神经网络算法的综合性算法」· M 代码 · 共 16 行

function dx=fun5(x,F,zF)
a=1.5;L=2.706;V=3.206;H=0.5;
for i=1:41
    k(i)=a*x(i)/(1+(a-1)*x(i));
end
dx(1)=((L+F)*x(2)-(L+F-V)*x(1)-V*k(1))/H;
for i=2:20
    dx(i)=((L+F)*x(i+1)+V*k(i-1)-(L+F)*x(i)-V*k(i))/H;
end
dx(21)=(L*x(22)+V*k(20)-(L+F)*x(21)-V*k(21)+F*zF)/H;
for i=22:40
    dx(i)=(L*x(i+1)+V*k(i-1)-L*x(i)-V*k(i))/H;
end
dx(41)=(V*k(40)-L*x(41)-(V-L)*x(41))/H;
end