- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
偏微分方程之期末重要知识总复习.doc
偏微分方程数值解总复习
一、考虑一维经典的初值问题:
设函数在G=中连续,并且是关于u满足Lipschitez条件,即存在一个只依赖区域G,而与变量t,u无关的常数L(称为Lipschitez常数),使得对任意的(t,u1)和(t,u2)G,都有,这里的表示R中的任一种范数。给定等距分割:,其中步长,。
在上作:,
这一方法称为Euler方法。
如果记为微分方程在处的精确解,为差分方程在处的精确解。
1、在上,定义算子:
当时,称数值方法是相容的。
2、当时,若,,则称该数值方法是收敛的。
3、如果由初值得到精确解;由初值得到精确解,若存在常数C和充分小的步长,使得,,。则称数值方法是稳定的。
证明:Euler方法是相容的、收敛的、和稳定的。
证明
1、 将在处做Taylor展开,得
所以该数值方法是相容的。
2、设,,,
则
Euler方法是稳定的
3、设,,,,
在上,定义算子
如果是微分方程的解,则
由于
和
所以
亦即
或
从而
于是当时,若,,则Euler方法是收敛的,且有误差估计:
二、考虑一维经典的初值问题和分割均同第一题。若在上作
,
其中,,不同时为零。此类方法通称为线性k步方法。
在上,定义算子
如果,则称为该数值方法是q阶的。
试讨论线性2步方法:
当时方法是二阶的,当时方法是三阶的。
设由初值得到精确解;由得到精确解,若存在常数C和充分小的,使得,,。则称数值方法是稳定的。
试讨论Adams二步外插格式
和Milne格式的稳定性。
如果当时,若,,则称该数值方法是收敛的。
试证明:如果线性k步方法是相容的和稳定的,那么线性k步方法也是收敛的。
证明:将,在t点展开得
于是
其中 ,,…,,…
因此,线性k步方法相容的条件是:和
当时,,当时
所以,当时方法是二阶的,当时方法是三阶的。
由于线性k步方法是稳定的充分必要条件是满足根条件。
而Adams二步外插格式的,因此它是稳定的。
而Milne格式的,因此它也是稳定的。
设线性k步方法是相容的和稳定的,且
,,,
其中
记 ,,
于是可以将上式写成矩阵的形式
写出第一个分量:
即
因为满足根条件,故有界D,()有界M
于是可得
记
于是
可得
其中 ,,,
且
应用Gronwall不等式得
于是当时,有 即,因此线性k步方法是收敛的。
三、考虑一维经典的初值问题和分割均同第一题。若在上作
,
其中,且,
,,
此类方法通称为n级Runge-Kutta方法。
在上,定义算子
如果且时,则称为该数值方法是相容的。
试讨论n级Runge-Kutta相容的条件。
如果则称为该Runge-Kutta是q阶的。
试讨论3级Runge-Kutta是3阶的条件。
由于
于是
由于是微分方程的解,因此有:
因此当的充分必要条件为:
将在t点展开
其中 ,
下面考虑在点(t,u)的展开式
,
这样
由于
3级Runge-Kutta是3阶的条件为:
四、考虑一维热传导方程的初边值问题:
假设初边值问题存在唯一的、足够光滑的解。
给定等距分割:,,其中空间步长:,;时间步长,。这样以两族平行坐标轴的直线()和()将矩形域分割成矩形网格,网格节点为。记是微分方程在点处的精确解,是差分方程在点处的精确解。称为网格比,。记,()。
当时,若,,则称该差分格式是收敛的。
对于给定的古典显格式:
,,
试讨论古典显格式是收敛的。
如果由初值得到精确解;由初值得到精确解,若存在常数K和充分小的步长,使得,,。则称该差分格式是关于初值稳定的。
对于给定的古典显格式:
,,
试讨论古典显格式关于初值稳定性的条件(不使用Fourier方法)。
证明:对于给定的古典显格式:
,,
如果记,
则可以将古典显格式写成,
记,,因为,
从而
所以
于是,为使差分格式关于初值稳定,就必须要求矩阵族一致有界,即
通过计算可以求出的特征值,。
因此,可用相似变换S将C化为对角形,使得
于是有 , 其中
要使矩阵族一致有界,就必须要求特征值,即,亦即,这样只要就
文档评论(0)