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方程描述。

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方程描述。
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)相拟合，可行不？有没有什么办法？

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方程描述。
调用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)
只要给出参数的具体数值和初始条件,有解.

不能！！
原因：这是一个锁死类方程。。
如果 n 也是 t 的函数，那么这个方程是无解的，有解请告诉我。
如果 n 和 t 没有关系，那么返程等价为两个ODE，分别使用ode方法就能够解出来了。。
matlab中间使用的解ode的方法有很多，比如说ode45，ode23，等等，你在matlab程序的CW里面键入 help ode45，就可以看到具体怎么使用了。。...

全部展开

不能！！
原因：这是一个锁死类方程。。
如果 n 也是 t 的函数，那么这个方程是无解的，有解请告诉我。
如果 n 和 t 没有关系，那么返程等价为两个ODE，分别使用ode方法就能够解出来了。。
matlab中间使用的解ode的方法有很多，比如说ode45，ode23，等等，你在matlab程序的CW里面键入 help ode45，就可以看到具体怎么使用了。。

收起

淘淘VS跑跑 见习魔法师 二级(332) | 我的提问 | 我的回答 | 知识掌门人 | 我的消息(0/1) | 百度首页 | 退出
新闻 网页 贴吧 知道 MP3 图片 百科
帮助


百度知道 > 理工学科 > 数学 分类上升达人排行榜
用户名 动态 上周上升
474096872 10845...

全部展开

淘淘VS跑跑 见习魔法师 二级(332) | 我的提问 | 我的回答 | 知识掌门人 | 我的消息(0/1) | 百度首页 | 退出
新闻 网页 贴吧 知道 MP3 图片 百科
帮助


百度知道 > 理工学科 > 数学 分类上升达人排行榜
用户名 动态 上周上升
474096872 10845
我不是他舅 8725
nbgauss 4195
sxzh135122 2875
jjj976431 2250
更多>>
订阅该问题
找好工作,上51job.com
中国领先的招聘网站,2000万注册会员,免费填写中英文简历,好工作立即找上门
www.51Job.com
易趣又做促销啦!谁说便宜无好货?
网上购物来易趣,底价促销抢翻天!服饰数码化妆品,超值商品任您选!
www.EachNet.com

您想在自己的网站上展示百度“知道”上的问答吗？来获取免费代码吧！
--------------------------------------------------------------------------------
如要投诉或提出意见建议，请到
百度知道投诉吧反馈。
添加到搜藏待解决
matlab高手进——该微分方程怎么解？
悬赏分：90 - 离问题结束还有 9 天 22 小时
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)相拟合，可行不？有没有什么办法？
提问者： fancyui - 试用期 一级
回答 共 2 条
不能！！
原因：这是一个锁死类方程。。
如果 n 也是 t 的函数，那么这个方程是无解的，有解请告诉我。
如果 n 和 t 没有关系，那么返程等价为两个ODE，分别使用ode方法就能够解出来了。。
matlab中间使用的解ode的方法有很多，比如说ode45，ode23，等等，你在matlab程序的CW里面键入 help ode45，就可以看到具体怎么使用了。。
回答者：sanapoint - 江湖新秀 四级 5-7 14:52
调用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)
只要给出参数的具体数值和初始条件，有解。
回答者：化学工程 - 首席运营官 十三级 5-8 00:22
我也来回答：
回答即可得2分，回答被采纳则获得悬赏分以及奖励20分。 积分规则
回答字数在10000字以内
参考资料：
如果您的回答是从其他地方引用，请表明出处。
匿名回答
©2008 Baidu

收起