我将一项功能另存为' calculateResults.m'其中的功能' ode45'嵌套:
function y = calculateResults (concV, ksV, xTrace)
options = odeset('NonNegative', [1:size(concV,2)], 'RelTol', 1e-6, 'AbsTol', 1e-12);
[x,y] = ode45(@(t, cY)odeSet(t, cY, ksV), xTrace, concV, options);
我想将此功能用作fittype:
ft = fittype('calculateResults( concV, ksV, x )','independent',{'x'});
但这样做会导致错误:
使用fittype / testCustomModelEvaluation时遇到错误(第12行) 表达式calculateResults(concV,ksV,x)不是有效的MATLAB表达式,有非标量系数,或者不能 被评估: fittype表达式中的错误==> calculateResults(concV,ksV,x) ???指数超过矩阵维度。
我的目标是改变' concV'为了适应功能' calculateResults.m'一些数据(x和y):
mdl = fit( x, y, ft);
最后返回最适合数据的concV值:
concVCalc = mdl.Coefficients.Estimate;
请帮忙!谢谢!