查看: 5143|回复: 7

Maple求解微分方程的结果与手算不一致 ~

字体大小: 正常放大

wangzhenyu1990

1 主题	0 听众	6 积分

升级 1.05%

该用户从未签到

自我介绍: Maple爱好者，深入学习，欢迎交流。

电梯直达

1^#

发表于 2012-10-25 23:10 |只看该作者 |正序浏览

|招呼Ta 关注Ta

我想求解二阶常系数线性齐次微分方程：未命名.jpg

即 0.125*(diff(y(x), x, x))+.75*(diff(y(x), x))+y(x) = 0

Maple求解过程为：eq:=0.125*diff(y(x),x$2)+0.75*diff(y(x),x)+y(x)=0;
dsolve(eq);
plots[odeplot](dsolve({eq,y(0)=2,D(y)(0)=4},numeric,range=0..20));
得到的结果为：y(x) = _C1*exp(-2*x)+_C2*exp(-4*x)
图像为：

但是我用手算得到的结果为：
y(x) = 6*exp(-2*x)-4*exp(-4*x)
图像为：

想问一下这是为什么 Maple求解的结果不可靠吗？如果不带数据，更能明显的看到Maple求解的结果和实际的手算不一样：
不带入数据的微分方程为：

用Maple求解为：y(x) = _C1*exp((1/2)*(-RC+sqrt(RC^2-4*LC))*x/LC)+_C2*exp(-(1/2)*(RC+sqrt(RC^2-4*LC))*x/LC)
而用手算结果为：

请问为什么？是不是我的算法有问题？谢谢！

未命名.jpg (12.76 KB, 下载次数: 238)

未命名.jpg

4.jpg (7.88 KB, 下载次数: 237)

zan

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

使用道具举报

一汪清水

1 主题	7 听众	6 积分

升级 1.05%

TA的每日心情

	无聊 2013-12-7 17:26

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 数学专业

8^#

发表于 2013-11-29 19:46 |只看该作者 |招呼Ta 关注Ta

不错的解法

我也说一句

发表

使用道具举报

Timmy金田

0 主题	0 听众	2 积分

升级 40%

该用户从未签到

自我介绍: matlab fans

7^#

发表于 2013-5-24 17:54 |只看该作者 |招呼Ta 关注Ta

请问求积分公式也可以用desolve方法吗？

我也说一句

发表

使用道具举报

lxin

0 主题	8 听众	21 积分

升级 16.84%

TA的每日心情

	开心 2016-2-27 22:46

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 自己想再学习一下数学

6^#

发表于 2013-3-8 17:42 |只看该作者 |招呼Ta 关注Ta

zsy312 发表于 2012-11-9 20:46

很不错，学习中，感谢感谢

我也说一句

发表

使用道具举报

wolegecacac

2 主题	5 听众	51 积分

升级 48.42%

TA的每日心情

	奋斗 2013-2-1 08:32

签到天数: 10 天

[LV.3]偶尔看看II

自我介绍: 好人

群组: 学术交流A

群组: 学术交流B

5^#

发表于 2013-1-28 14:12 |只看该作者 |招呼Ta 关注Ta

不错·························

我也说一句

发表

使用道具举报

久玖九王

3 主题	7 听众	46 积分

升级 43.16%

TA的每日心情

	无聊 2014-8-23 18:14

签到天数: 11 天

[LV.3]偶尔看看II

自我介绍: 因数学而存在的存在。

群组: Linux推广

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

4^#

发表于 2012-11-12 12:23 |只看该作者 |招呼Ta 关注Ta

也是呀，难道是需要付费才可以好的？

我也说一句

发表

使用道具举报

zsy312

2 主题	4 听众	64 积分

升级 62.11%

TA的每日心情

	开心 2013-9-23 18:22

签到天数: 4 天

[LV.2]偶尔看看I

3^#

发表于 2012-11-9 20:46 |只看该作者 |招呼Ta 关注Ta

eq:=0.125*(diff(y(x), x, x))+.75*(diff(y(x), x))+y(x) = 0:
+ x0 n' _8 j0 k2 x
dsolve([eq,y(0)=2,D(y)(0)=4],y(x));- P0 \, _9 e( |: ?
plot(rhs(%),x=0..5,scaling=constrained);
+ t, ^# ?+ K\" {4 C( v* F1 H! @

复制代码

eq:=L*C*(diff(u[c](t), t, t))+R*C*(diff(u[c](t), t))+u[c](t) = 0;
6 T# Z0 ?6 U+ w7 T2 b, P
dsolve([eq,u[c](0)=u[0],D(u[c])(0)=i[0]],u[c](t));+ `; g6 V+ R/ W; P( u

复制代码