Implement for loop in ode45 solver
显示 更早的评论
I'm trying to vary parameter.D with a for loop so i can use the ode45 solver for different values of parameter.D
parameter.mumax = 0.42; %1/h
parameter.Km = 0.25; %g/l
parameter.Sin = 60; %g/l
parameter.Yxs = 0.25; %g/g
parameter.D = 0.1:0.1:0.5;%1/h
y0 = [ 5, 60]; %biomass substrate
tspan = [0 100]; %h
options = odeset('RelTol', 1e-6);
[t,y] = ode45(@biomassproductionmodel, tspan, y0, options, parameter);
function [dydt] = biomassproductionmodel(t,y, parameter)
mu = parameter.mumax*y(2)/(parameter.Km+y(2));
dydt(1,1) = mu*y(1) - parameter.D*y(1);
dydt(2,1) = -mu*y(1)/parameter.Yxs + parameter.D*(parameter.Sin - y(2));
end
采纳的回答
parameter.mumax = 0.42; %1/h
parameter.Km = 0.25; %g/l
parameter.Sin = 60; %g/l
parameter.Yxs = 0.25; %g/g
D = 0.1:0.1:0.5;
% parameter.D = 0.1:0.1:0.5;%1/h
y0 = [ 5, 60]; %biomass substrate
tspan = [0 100]; %h
options = odeset('RelTol', 1e-6);
for i=1:numel(D)
parameter.D = D(i);
[t,y] = ode45(@(t,y)biomassproductionmodel(t,y,parameter), tspan, y0, options);
T{i} = t;
Y{i} = y;
end
figure(1)
plot(T{1},Y{1}(:,1),T{2},Y{2}(:,1),T{3},Y{3}(:,1),T{4},Y{4}(:,1),T{5},Y{5}(:,1))
figure(2)
plot(T{1},Y{1}(:,2),T{2},Y{2}(:,2),T{3},Y{3}(:,2),T{4},Y{4}(:,2),T{5},Y{5}(:,2))
function [dydt] = biomassproductionmodel(t,y, parameter)
mu = parameter.mumax*y(2)/(parameter.Km+y(2));
dydt(1,1) = mu*y(1) - parameter.D*y(1);
dydt(2,1) = -mu*y(1)/parameter.Yxs + parameter.D*(parameter.Sin - y(2));
end
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息
Translated by 
