查看: 1484|回复: 3

[问题求助] 求解常微分方程遇到的问题

字体大小: 正常放大

4 主题	11 听众	89 积分

升级 88.42%

TA的每日心情

	开心 2015-1-31 19:30

签到天数: 33 天

[LV.5]常住居民I

电梯直达

1^#

发表于 2014-8-14 10:23 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

clc,clear
dsolve('m*D2s=m*g-rou*g*V-k*Ds','Ds(0)=0','t')
m=239.46;g=9.8;rou=1035.71;V= 0.2058;k=0.6;
运行结果
ans =

(g*m*(m - V*rou))/(k^2*exp((k*t)/m)) - ((C4*exp((k*t)/m))/k^2 - (g*t*exp((k*t)/m)*(m - V*rou))/k)/exp((k*t)/m)
怎样将m、g等的值代入？

zan

转播1 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

madio

3万主题	1310 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

2^#

发表于 2014-8-14 15:52 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

clc,clear
syms m g rou V k s
s=dsolve('m*D2s=m*g-rou*g*V-k*Ds','Ds(0)=0','t');
m=239.46;g=9.8;rou=1035.71;V= 0.2058;k=0.6;
subs(s)

复制代码

你应该要这样去写！

我也说一句

发表

数学建模社会化

使用道具举报

MCR董事长

5 主题	10 听众	421 积分

升级 40.33%

TA的每日心情

	开心 2015-10-26 08:32

签到天数: 143 天

[LV.7]常住居民III

自我介绍: 新手

群组: 国赛讨论

3^#

发表于 2014-8-14 15:52 |只看该作者 |招呼Ta 关注Ta

>> dsolve('m*D2s=m*g-rou*g*V-k*Ds','Ds(0)=0','t')
ans =
(m*exp(-(k*t)/m)*(g*m - V*g*rou))/k^2 - exp(-(k*t)/m)*((C4*exp((k*t)/m))/k^2 - (t*exp((k*t)/m)*(g*m - V*g*rou))/k)
>> m=239.46;g=9.8;rou=1035.71;V= 0.2058;k=0.6;
>> b=subs(ans)
b =
(188578330473244179*exp(-(10*t)/3991))/1099511627776 + exp(-(10*t)/3991)*((236254485684345*t*exp((10*t)/3991))/549755813888 - (25*C4*exp((10*t)/3991))/9)

我也说一句

发表

使用道具举报

MCR董事长

5 主题	10 听众	421 积分

升级 40.33%

TA的每日心情

	开心 2015-10-26 08:32

签到天数: 143 天

[LV.7]常住居民III

自我介绍: 新手

群组: 国赛讨论

4^#

发表于 2014-8-14 15:52 |只看该作者 |招呼Ta 关注Ta

>> dsolve('m*D2s=m*g-rou*g*V-k*Ds','Ds(0)=0','t')
ans =
(m*exp(-(k*t)/m)*(g*m - V*g*rou))/k^2 - exp(-(k*t)/m)*((C4*exp((k*t)/m))/k^2 - (t*exp((k*t)/m)*(g*m - V*g*rou))/k)
>> m=239.46;g=9.8;rou=1035.71;V= 0.2058;k=0.6;
>> b=subs(ans)
b =
(188578330473244179*exp(-(10*t)/3991))/1099511627776 + exp(-(10*t)/3991)*((236254485684345*t*exp((10*t)/3991))/549755813888 - (25*C4*exp((10*t)/3991))/9)