- 1、本文档共17页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
毕业设计(论文)外文翻译
题 目:二阶常微分方程边值问题的数值解
系 :数理与信息系
学生姓名:施逸斌
学 号:081315123
专 业:信息与计算科学
班 级:C08信息
指导老师: 黄红英
起止日期:2012.3.17-2011.3.27
2012年 3月 27 日
11.5 Rayleigh-Ritz法
打靶法是将边值问题转化为一对初值来求边值问题的近似求解法。有限差分法将连续微分运算离散化为有限差分运算。Rayleigh-Ritz法是从第三种方法去解决这个问题的变分技术。边值问题首先被重新转述为从满足边界条件的所有充分可微函数中选择函数来使得某个积分达到最小值的问题。这样,可行性函数的集合范围被缩小,从而得到最小化问题的近似解且得到边值问题的近似解。
为了描述Rayleigh-Ritz法,我们考虑来自梁受力分析的线性两点边值问题的近似解。这个边值问题可描述为微分方程:
(11.21)
边界条件:
(11.22)
这个微分方程描述了长度为1的可变交叉截面(记为)的横梁的挠度。挠度取决于施加的压力和。
在下面的讨论中,我们假设且。进一步我们假设存在常数,使得
且,
这些假设保证了式(11.22)和式(11.23)所给出的边值问题有唯一解(见参考文献[BSW])。
如同许多描述物理现象的边值问题一样,梁方程的解满足变分特性。粱方程的变分原理对于开发Rayleigh-Ritz法是很基本的,同时它把粱方程的解表征为对中的所有函数使得某积分达到最小值的函数,其中函数集是指上那些满足性质的函数的集合。下面的定理描述了该特征。
定理11.4 设,,且
,,
函数是微分方程
, (11.23)
的唯一解,当且仅当是上唯一使得积分
(11.24)
达到最小的解。
该定理的具体证明可见参考文献[Shul,pp.88-89]。该证明步骤分三步进行。首先证明式(11.24)的任意解对所有都满足方程
(11.25)
第二步证明是式(11.25)的解当且仅当对于所有式(11.26)成立。
最后一步证明式(11.26)有唯一解。该唯一解也是式(11.25)和式(11.24)的解,故式(11.24)和式(11.25)的解是相同的。
Rayleigh-Ritz法就是通过最小化积分来逼近解,它并不是对上所有函数最小化,而是对包含某些基函数的线性组合的比较少的函数集最小化。这些基函数是线性无关的,且满足
,
对方程(11.24)的解的近似解是通过找到常数使最小化来得到的。
由方程(11.25),有
(11.26)
且知道要达到最小,将视为的函数,则必有
, (11.27)
对式(11.27)进行微分,得到
代入方程(11.28)得到
,
(11.28)
方程(11.29)所描述的法方程产生一个线性方程组Ac=b,变量为,其中A为对称矩阵,其元素为
定义b为
基函数最简单的选择是分段线性多项式。第一步是选择满足
的点来对[0,1]进行划分。设(),定义基函数为
(11.29)
其中,。(见图11.4)
图11.4
由于函数是分段线性的,所以导数不
文档评论(0)