正如标题中指定的那样,我需要强制一个函数始终为正(如果它运行为负,则它必须等于0)。这个函数在ODE中:
这是我调用ODE的脚本:
clear
t=[0,276];
Tie=0.3;
mumax=2;
Qmin=1;
X0=[4,2];
[t,X]=ode45(@(t,X) odeset(mumax,Qmin,Tie,X),t,X0);
这是我希望始终保持正面的功能
function [ func2 ] = func2 (mumax,Qmin,Q)
func2=mumax*(1-Qmin/Q);
end
这是ODE
function [ dXdt ] = odeset(mumax, Qmin, Tie,X)
dXdt=zeros(2,1);
dXdt(1)=func2(mumax,Qmin,X(1))-X(2)*Tie;
dXdt(2)=func2(mumax,Qmin,X(1))-X(2)*Tie;
end
答案 0 :(得分:1)
将func2
的输出更改为
function [ out ] = func2 (mumax,Qmin,Q)
out=max(mumax*(1-Qmin/Q),0);
end
这将强制输出为正或零。