想使用fmincon函数优化求微分方程的参数,写完程序之后提示这个错误,大佬们能不能帮我看看哪里错了?
也帮忙看一下我写的程序对不对,有没有有需要改正的地方。
k0=[0.023 0.000613 0.1098 0.00038];
lb=[0 0 0 0];
% ub=[+inf +inf +inf +inf];
x0=[4000000 62166 1044 38364 53793 3311];
load Expdata;
t=Expdata(:,1);
yexp=Expdata(:,[2,5]);
%使用fmincon函数进行参数估计
[k,fval,flag] = fmincon('ObjFunc4Fmincon',k0,[],[],[],[],lb,[],[],[],x0,yexp);
fprintf('\n使用函数fmincon()估计得到的参数值为:\n')
fprintf('\tk1 = %.4f\n',k(1))
fprintf('\tk2 = %.4f\n',k(2))
fprintf('\tk3 = %.4f\n',k(3))
fprintf('\tk4 = %.4f\n',k(4))
fprintf(' The sum of the squares is: %.1e\n\n',fval)
k_fmincon = k;
%---------------------------
function f=ObjFunc4Fmincon(k,x0,ye