fconc12.m

来自「Incorporating Prior Knowledge in Cubic S」· M 代码 · 共 19 行

function y = fconc12(x)

global b1 b2 b3 b4
global c12
global cnew

b1 = x(1);
b2 = x(2);
b3 = x(3);
b4 = x(4);
[t,cnew] = ode45('freact1',[0:0.02:2],[1 0 0]);
c12new = cnew(round(linspace(1,101,12)),:);

y = sum(sum((c12new-c12).*(c12new-c12)));