查看: 5140|回复: 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 微信

使用道具举报

wangzhenyu1990

1 主题	0 听众	6 积分

升级 1.05%

该用户从未签到

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

2^#

发表于 2012-10-25 23:15 |只看该作者 |招呼Ta 关注Ta

这个问题困扰着我微分方程比较简单手算也能出结果但是我想用Maple简便的算出正确的结果谢谢！

我也说一句

发表

使用道具举报

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:; s/ _# m- {8 g- t4 N
dsolve([eq,y(0)=2,D(y)(0)=4],y(x));
7 C) e# L) G/ l7 T9 S
plot(rhs(%),x=0..5,scaling=constrained);) ]1 Q# a8 x8 V

复制代码

eq:=L*C*(diff(u[c](t), t, t))+R*C*(diff(u[c](t), t))+u[c](t) = 0;
+ F0 \0 D( s( z8 x' u; |2 B; T' V M
dsolve([eq,u[c](0)=u[0],D(u[c])(0)=i[0]],u[c](t));2 x$ a$ J) Z7 N\" X# ^' X

复制代码