返回信息流为这三个微分方程的解,已经泪奔好几天了,希望大家都帮帮忙,感谢感谢!
这是一条镜像帖。来源:北邮人论坛 / math-model / #8505同步于 2011/12/2
该镜像源已超过 30 天没有更新,可能在源站已被删除。
MathModel机器人发帖
急求帮忙:解微分方程,已经泪奔了。重谢!
james527227
2011/12/2镜像同步7 回复
订阅后,新回复会通过你的通知中心匿名送达。
7 条回复
matlab程序
clc
clear
y1=dsolve('Dy=p*exp(-d*t)*(m-y)+q/m*y*(m-y)','y(0)=0','t')%D2用d表示
y2=dsolve('Dy=p*(m-y)+q/m*y*(m-y)*exp(-d*t)','y(0)=0','t')
y3=dsolve('Dy=p*exp(-d1*t)*(m-y)+q/m*y*(m-y)*exp(-d2*t)','y(0)=0','t')
大牛,这个跑出来的结果,我还是没办法继续整理出结果,不知你能否帮忙,能有类似这样的结果吗?
其中N(t)即是y
【 在 jellyray 的大作中提到: 】
: matlab程序
: clc
: clear
: ...................
恩,我大概知道什么意思了
你这个图应该是不加e的微分方程的解,想得到相同形式加上e的d次方的解
我试了一下y1的化简
y1=simple(y1)
pretty(y1)
得不到所要的形式
又用math试了一下
a = DSolve[{y'[t] == p*(m - y[t]) + q/m*y[t]*(m - y[t]), y[0] == 0}, y[t], t]
Simplify[a]
结果更不理想
无能为力了,期待高手解答。。。
ps:还有个问题,我得到不加e也就是math下
a = DSolve[{y'[t] == p*(m - y[t]) + q/m*y[t]*(m - y[t]), y[0] == 0}, y[t], t]
Simplify[a]
结果
忽略符号问题,这个式子的分母和你的那个图还是有点出入(p,q),是条件问题还是我理解错了?
【 在 james527227 的大作中提到: 】
: 大牛,这个跑出来的结果,我还是没办法继续整理出结果,不知你能否帮忙,能有类似这样的结果吗?
: 其中N(t)即是y
: 【 在 jellyray 的大作中提到: 】
: ...................
恩,出入就是正负号的出入。这个可能是一点小问题。恩,那这样的意思,是否这三个微分方程没有解析解?
【 在 jellyray 的大作中提到: 】
: 恩,我大概知道什么意思了
: 你这个图应该是不加e的微分方程的解,想得到相同形式加上e的d次方的解
: 我试了一下y1的化简
: ...................
恩,还真是,呵呵,没仔细看还真没看出来。辛苦了!哎,求大牛啊求大牛。泪奔了啊
【 在 jellyray 的大作中提到: 】
: 不是正负号!是p,q有个反了,化简下就看出来了,求解释= =
: --