ode45 + lsqcurvefit- 多个 ODE

Question

我正在尝试将我的实验数据拟合到包含 2 个耦合 ODE 的模型。我创建了一个单独的函数来定义 2 个 ODE，其中 a 是包含我要拟合的 2 个常量的数组。

function dy = CalcCoupODEs(t,y,a) 
%dy= matrix with all 2 ODEs
dy = zeros(2,1); % variable definition
dy(1) = 1.41.*a(1) .* ((y(2)*y(1)^2)^2)*y(1) - a(2) .* y(1)*t ; %ODE for R
dy(2) = -a(1)*(y(2)*y(1)^2)^2; %ODE for n
end

然后，我按以下方式调用该函数：

cycle_DMACl = [20, 25, 30, 40, 45]; % this is t variable
R_DMACl = [18.34, 11.29, 7.09, 6.51, 4.396]; %this is y(1)
n_DMACl = [92.31, 61.9, 53.82, 26.04, 11.87]; % this is y(2)
a0 = [1, 1]; %initial guess for constants
lb1 = [0, 0]; %lower bound
a = lsqcurvefit(@CalcCoupODEs, a0, cycle_DMACl(:), R_DMACl(:), n_DMACl(:), lb1);

我得到以下信息：

Warning: Length of lower bounds is > length(x); ignoring extra bounds.

显然不是在尝试拟合正确的东西（a 和 lb 都是 1x2）。我做错了什么，最好的方法是什么？

谢谢！

Answer 1

我对 lsqcurvefit(...) 经验不多，但快速浏览一下 MathWork 的函数文档，我认为您的参数不符合 MATLAB 的语法预期。你的用法看起来像

 lsqcurvefit(fun,initial_points,time,y1_data,y2_data,lower_bounds)

这是带有两个 y 数据条目的六个参数。该函数的文档显示了一个六参数语法，即

 lsqcurvefit(fun,initial_points,time,y_data,lower_bounds,upper_bounds)

我尝试通过制作一个新矩阵来快速修复，

 y_data = [y1_data', y2_data']

并插入一些高值和随机上限。我在 return 中遇到错误。

我认为您将不得不重新定义您的 ydata 数组并定义上限值。我会继续使用它并 post 任何进展，但这可能需要一些时间。祝你好运，请post一个解决方案。

Answer 2

感谢您的回复。将 y_data 放入矩阵中确实解决了问题。我注意到如果初始值关闭，lsqcurvefit 会给出一个错误，即 YDATA 和函数值的大小不同（这听起来像是与初始值无关的错误）。在使用初始值玩了很多之后，我让它收敛了。

再次感谢！

ode45 + lsqcurvefit- 多个 ODE

ode45 + lsqcurvefit- multiple ODEs

matlab

ode

non-linear-regression