![在这里插入图片描述]()
增长率为k=r(M-P)其中,r是常数,M为环境容纳量的上限(比如人口/或某些生物在某个环境下的数量上限)
>> dsolve('DP=r*(M-P)*P')
ans =
M/(1+exp(-r*M*t)*C1*M)
>> dsolve('DP=r*(M-P)*P','P(t0)=P0')
ans =
M/(1+exp(-r*M*t)*exp(t0*r*M)*(M-P0)/P0)
使用dsolve求解微分程注意
方程中的乘法号’*'必须要显示的写出来不然matlab以为那是某个标识符,导致出错
另外,matlab求得的结果往往不是最简的.