原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第 32 卷第 1 期 数学的实践与认识 V o l32 N o 1
2002 年 1 月 M A TH EM A T ICS IN PRA CT ICE AND TH EO R Y Jan. , 2002
多约束二阶非线性常微分方程极值点的数值求解法
杜志明, 冯长根
(北京理工大学机电工程学院, 北京 100081)
摘要: 本文对科学研究和工程应用问题中常见的多约束二阶非线性常微分方程的极值点进行了讨论, 提
出了一种数值解法. 以带有化学反应项的热传导方程为例, 给出了相应的计算结果.
关键词: 非线性; 常微分方程; 极值; 数值方法; 热传导方程
1 引 言
大量的科学研究问题和工程应用问题最终可以归结为微分方程的求解问题. 在数学上
已经发展了多种求解微分方程分析解和近似解的方法, 但仍有一些复杂类型的微分方程需
进行数值计算. 数值计算有时不仅需求解微分方程的一般解, 而且还需要求解带约束参数
( )
的微分方程的极值解. 极值解往往对应于实际问题的一些敏感性质 如临界性 , 因而受到
研究工作者的特别关注.
计算机技术的飞速发展为一些复杂科学技术问题的解决提供了前所未有的手段, 本文
将给出一种求解复杂的多约束二阶非线性常微分方程极值点的数值方法, 并以在热爆炸研
[ 1 ]
究中常见的带有化学放热反应项的热传导方程 为例, 给出其极值点的计算结果, 这些数值
正是实际问题发生热爆炸与否的关键数值, 称为临界值或爆炸判据.
2 定态极值条件
带有约束参数的定态方程组一般可表示为:
f i (x 1 , x 2 , …, x n ; z 1 , z 2 , …, z m ) = 0
i = 1, 2, …, n ( 1)
其中x j (j = 1, 2, …, n ) 是状态变量, z s (s= 1, 2, …,m ) 是约束参数. f i 表示某种函数关系. 显
( )
然, 当各约束参数都有确定值时, 方程 1 才有确定解.
( )
将式 1 两边微分, 可得:
n f i m f i
df i = ∑x j d x j + ∑z s d z s = 0 (2)
j = 1 s= 1
( ) ( )
若只让一个约束参数变化 不妨记为 p 而固定其余约束
文档评论（0）