- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
B3: Maple中的偏微分方程求解
西希安工程模拟软件 (上海)有限公司,2008
11.0 Maple中的微分方程求解器介绍
Maple中微分方程求解器使用领先的算法求解 以下问题:
常微分方程 (ODEs):dsolve命令用于求解线性和非线性ODEs,初始值问题 (IVP),以及边界
值问题 (BVP),可以通过参数项选择求符号解 (解析解)或数值解。
微分方程分析器助手提供一个交互式用户界面方便用户求解 ODE 以及显示结果的图形。了
解更多信息,参考帮助系统中的 dsolve,dsolve/numeric,和 ODE Analyzer.
偏微分方程 (PDEs):pdsolve命令用于求 PDEs和含边界值问题的 PDEs的符号解或数值
解。使用Maple的PDE 工具可以完成对PDE 系统的结构分析和指数降阶处理。了解更多信
息,参考帮助系统中的 pdsolve andpdsolve/numeric.
微分-代数方程 (DAEs):dsolve/numeric命令是符号-数值混合求解器,使用符号预处理和降
阶技术,让Maple能够求解高指数的DAE问题。Maple内置三个求解器用于处理DAEs:1)
修正的 Runge-Kutta Fehlberg方法,2)Rosenbrock方法,以及 3)修正的拓展后向差分隐
式方法。
11.1求解偏微分方程PDE 问题 (BVP和IVP)
Maple求解经典力学难题的能力是非常著名的,它的数值和符号偏微分方程求解器是其中的
重要工具。
例子:在不同的边界条件下,求波动方程的数值解、解析解、和图形解。
11.1.1初始化
下面的Maple代码定义了一个名为PX的程序,生成函数的周期展开。
local L, D, var;
var := op(var); end if;
end if;
end:
11.1.2数值解和图形解
一个空间变量的波动方程是:
(1)(1)
假设初始形状 由下面的函数给出:
(2)(2)
对应的图形如下:
0
0 2 4 6 8 10
x
假设波动方程的初始条件和边界条件如下:
(3)(3)
Maple中的命令pdsolve将求解单变量演化方程 (双 曲和抛物)的数值解 (有限差分)。
(4)(4)
这个命令创建了一个模块,可以看到模块的输出函数是 plot,plot3d,animate,和 value。
时间 的图形:
图 2
0
0 2 4 6 8 10
x
解的表面图:
Figure 3
动画:
time = 0.000000
1
0
2 4 6 8 10
x
给出计算值 的函数:
文档评论(0)