QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5141|回复: 7
打印 上一主题 下一主题

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

[复制链接]
字体大小: 正常 放大

1

主题

0

听众

6

积分

升级  1.05%

该用户从未签到

自我介绍
Maple爱好者,深入学习,欢迎交流。
跳转到指定楼层
1#
发表于 2012-10-25 23:10 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
我想求解二阶常系数线性齐次微分方程: 未命名.jpg " n; j/ N4 _& q( n3 h
" E# \1 i7 C3 C9 F1 @$ m, q  v
即  0.125*(diff(y(x), x, x))+.75*(diff(y(x), x))+y(x) = 01 y' u4 @! |: I, k3 o. y" p
2 }0 c' H/ `( _8 w5 E  g, m
Maple求解过程为:eq:=0.125*diff(y(x),x$2)+0.75*diff(y(x),x)+y(x)=0;
5 |; X1 h! l! T( b. g# P! L dsolve(eq);
4 S& G' s; y- ]# D/ t7 ~/ h/ l plots[odeplot](dsolve({eq,y(0)=2,D(y)(0)=4},numeric,range=0..20));
: Q% P3 `& F2 Q/ e3 M/ }  F( ]4 l 得到的结果为:y(x) = _C1*exp(-2*x)+_C2*exp(-4*x); \8 Y* s* O: i- l$ }4 W7 y
图像为: 2.gif ; Q3 e4 ?! B3 |. o  z
但是我用手算得到的结果为:6 V8 U) e0 _5 R& h' Z( ~
y(x) = 6*exp(-2*x)-4*exp(-4*x)1 a% C3 r3 E1 z) v
图像为: 3.jpg 8 U* W3 r# N% B/ t" Y" Z  m
想问一下这是为什么 Maple求解的结果不可靠吗?如果不带数据,更能明显的看到Maple求解的结果和实际的手算不一样:
* V4 I( w4 a6 w/ V 不带入数据的微分方程为: 1.gif 2 G0 h/ q- e) ~' o* o" n0 f
用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)! ^2 C& D1 q5 c
而用手算结果为: 5.jpg
* A$ X- W& i: o请问为什么 ?是不是我的算法有问题?谢谢!

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

未命名.jpg

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

4.jpg

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

1

主题

0

听众

6

积分

升级  1.05%

该用户从未签到

自我介绍
Maple爱好者,深入学习,欢迎交流。
这个问题困扰着我  微分方程比较简单  手算也能出结果  但是我想用Maple简便的算出正确的结果  谢谢!
回复

使用道具 举报

zsy312        

2

主题

4

听众

64

积分

升级  62.11%

  • TA的每日心情
    开心
    2013-9-23 18:22
  • 签到天数: 4 天

    [LV.2]偶尔看看I

    新人进步奖

    1. eq:=0.125*(diff(y(x), x, x))+.75*(diff(y(x), x))+y(x) = 0:* L. C1 y/ n0 w! C
    2. dsolve([eq,y(0)=2,D(y)(0)=4],y(x));
      ' t9 M/ n6 `9 B
    3. plot(rhs(%),x=0..5,scaling=constrained);
      / b3 `' N, c2 T$ a3 R. u# O. x
    复制代码
    1. eq:=L*C*(diff(u[c](t), t, t))+R*C*(diff(u[c](t), t))+u[c](t) = 0;' Y! |5 K+ |0 c5 |8 l% B$ y# k6 H
    2. dsolve([eq,u[c](0)=u[0],D(u[c])(0)=i[0]],u[c](t));
      , j1 @9 N\" k# n$ p! {
    复制代码

    点评

    lxin  很不错,学习中,感谢感谢  详情 回复 发表于 2013-3-8 17:42
    回复

    使用道具 举报

    3

    主题

    7

    听众

    46

    积分

    升级  43.16%

  • TA的每日心情
    无聊
    2014-8-23 18:14
  • 签到天数: 11 天

    [LV.3]偶尔看看II

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

    群组Linux推广

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

    回复

    使用道具 举报

    2

    主题

    5

    听众

    51

    积分

    升级  48.42%

  • TA的每日心情
    奋斗
    2013-2-1 08:32
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    自我介绍
    好人

    群组学术交流A

    群组学术交流B

    回复

    使用道具 举报

    lxin        

    0

    主题

    8

    听众

    21

    积分

    升级  16.84%

  • TA的每日心情
    开心
    2016-2-27 22:46
  • 签到天数: 1 天

    [LV.1]初来乍到

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

    使用道具 举报

    0

    主题

    0

    听众

    2

    积分

    升级  40%

    该用户从未签到

    自我介绍
    matlab fans
    回复

    使用道具 举报

    1

    主题

    7

    听众

    6

    积分

    升级  1.05%

  • TA的每日心情
    无聊
    2013-12-7 17:26
  • 签到天数: 1 天

    [LV.1]初来乍到

    自我介绍
    数学专业
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-10 21:09 , Processed in 0.360316 second(s), 95 queries .

    回顶部