matlab高手进——该微分方程怎么解?
dN/dt=n*r*(1-N/M)
dn/dt=n*s*(1-n/K)
其中r,s,M,K为常数,能否求出 N=f(t)?
如果把n视为常量呢?如何用matlab解?
n(t)是已知的,为观测到的数据,但可以用Logistic方程描述。
我想把另一组观测到的数据与N=f(t)相拟合,可行不?有没有什么办法?
人气:229 ℃ 时间:2020-05-01 01:10:07
解答
调用maple
>> maple('dsolve','{diff(N(t),t)=n(t)*r*(1-N(t)/M),diff(n(t),t)=n(t)*s*(1-n(t)/K)}')
ans =
[{n(t) = K/(1+exp(-s*t)*_C2*K)},{N(t) = (exp(r/M*Int(n(t),t))*M+_C1)*exp(Int(-n(t)*r/M,t))}]
可以验证n(t)
>> n=dsolve('Dn=n*s*(1-n/K)')
n =
K/(1+exp(-s*t)*C1*K)
只要给出参数的具体数值和初始条件,有解.
推荐
猜你喜欢
- 在一个乘法算式里,已知积是其中一个因数的18倍,则另一个因数是多少
- 一项工程,甲乙合作8天可以完成,甲做6天,乙做了5天,共完成27/40.求甲单独要多少天能完成!
- 1、下面的加法算式是按一定规律排列的:1+2,2+5,3+8,1+11,2+14,3+17,...第40个加法算式是:
- 某种细菌每分由1个分裂成2个 (1)经过5min,1个细菌分裂成多少个?(2)这些细菌再继续分裂
- 才饮长沙水,又食武昌鱼.万里长江横渡,极目楚天舒.不管风吹浪打,胜似闲庭信步,今日得宽余.子
- 一个初中文言文题目~
- y=2arccos乘x减2分之1求函数定义域和值域 急
- 著名的雯波那契数列是这样的:1、1、2、3、5、8、13、21……这串数列当中第2010个数除以3所得的余数为多少