Я пытаюсь подогнать свои экспериментальные данные к модели, включающей 2 связанных ОДУ. Я создал отдельную функцию для определения двух ODE, где a
— это массив, содержащий две константы, которые я пытаюсь подогнать.
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);
Я получаю следующее:
Предупреждение: длина нижних границ >
length(x)
; игнорирование дополнительных ограничений.
Очевидно, что это не попытка подобрать правильную вещь (a
и lb
оба имеют размер 1x2). Что я делаю неправильно и как лучше всего это сделать?
Благодарю вас!