当前位置: 圆柱模板 » 建筑资讯 » 工程 » 正文

一类非线性刚性有效的数值求解方法_Brown算法

放大字体  缩小字体 发布日期:2018-11-05T20:19:03  主题关键词:求解  作者:支模网  浏览次数:784
核心提示:论文导读::针对一类非线性化程度差别较大的刚性系统,改变以往采用的Newton-Raphson迭代方法,引入Brown算法解决此类非线性刚性系统的数值求解问题。数值试验结果表明,Brown算法对于解决

论文导读::针对一类非线性化程度差别较大的刚性系统,改变以往采用的Newton-Raphson迭代方法,引入Brown算法解决此类非线性刚性系统的数值求解问题。数值试验结果表明,Brown算法对于解决该类问题的有效性。
论文关键词:刚性系统,Brown算法,数值求解
 

1 引 言

刚性微分方程存在于航空、航天、热核反应、自动控制、电子网络及化学动力学等许多重要科学技术领域及实际问题中[1,2],由于方程的解中既包含有衰减十分迅速的分量,又包含有相对来说变化缓慢的分量,两者的差别可以有好几个数量级,在选定计算方法时带来很大实质困难。实际研究证明,由于数值解稳定性限制,求解刚性微分方程主要采用隐式方法,如:隐式RK方法,BDF方法,IRK方法等。而采用隐式方法将刚性方程离散化以后,其变为线性或非线性方程(组)的求解问题。目前,对线性或非线性方程(组)的求解,多采用Newton-Raphson迭代求解。但对于某些非线性方程组,由于方程之间的非线性化程度相差较大,采用Newton-Raphson迭代方法数值求解的结果并不理想。本文利用Brown算法求解此类非线性刚性系统,具有较高精度和较快迭代速度的优点,数值试验结果表明了该方法的有效性。

2 Brown算法

考虑多个实变量的非线性方程组

刚性系统(2.1)

的数值求解问题,非线性方程组可以用向量形式表示:,其中

形如:刚性系统的公式称为Newton-Raphson迭代公式。由于该方法是将同时线性化,所以它并未考虑充分利用的具体结构。如果一个非线性的向量函数,其线性精度在各个分量上的分布可能是不平衡的,有的分量是非线性函数,而有的分量是线性函数,同时非线性函数组中也有非线性程度高低的差别,在此情况下,利用Newton-Raphson迭代方法对所有分量采用完全相同的数值处理,不利于方法整体计算效率的提高。

针对以上情况,Brown于1969年提出了按分量函数方程来形成迭代过程[3],其基本思想是对各分量逐个线性化并用其中每一个线性方程消去余下非线性方程中的一个变量,最后整个方程组就简化为一个仅含单个变量的非线性方程,应用一次单步Newton-Raphson迭代并结合逐一回代,即完成一次迭代过程[4]期刊网。

Brown算法的迭代步骤如下:

第一步,设为方程组(2.1)解的第次近似,函数处近似用线性函数

刚性系统

替代,令,由此求出

定义上式右端为

第二步,对函数定义一个新函数Brown算法,且记,其中。类似地,用线性函数来近似替代。令,解出

此时,个变量的线性函数,并记此线性函数为

步,由线性函数,可得,利用Newton-Raphson迭代,求得,并由出发,利用逐一回代,即

(2.2)

从而可求出,至此完成了一次Brown迭代过程。

3 数值试验

考虑以下常微分方程组初值问题:

问题1

其中:

问题2

其中:

对于上述两问题,当时,可计算其右函数组的Jacobi矩阵的特征值,均有,其余特征值绝对值均不超过6,因此系统呈强刚性。此外,观察两问题中的右函数组,可以看出除最后一个函数是高度非线性化外,其余函数都是线性的。

对于上述两问题,采用隐式Euler方法离散方程组,并分别用Newton-Raphson迭代法与Brown迭代法求解,取步长,及相对误差界表示迭代次数)控制每步迭代,最后得到数值解的最大绝对误差界,方程真解为:问题1,;问题2,。计算结果对比分析如表1所示。

表1 数值计算结果

 

问题1

问题2

Newton-Raphson迭代次数

迭代18次收敛

不收敛,

Brown迭代次数

迭代7次收敛

迭代8次收敛

数值解的绝对误差

(Newton-Raphson迭代)

3.83e+001

溢出

数值解的绝对误差

(Brown迭代)

1.44e-002

3.82e-002

 
关键词: 求解 刚性 数值
 
[ 订阅地图 ]  [ 加入收藏 ]  [ 告诉好友 ]  [ 打印本文 ]  [ 违规举报 ]  [ 关闭窗口 ]

 
求解  相关信息
0条 [查看全部]  相关评论

 
推荐图文
浅谈不定积分运算中的灵活性_凑微分法 一类非线性刚性有效的数值求解方法_Brown算法 一类新的自适应信赖域算法_全局收敛性 一种基于改进遗传算法的车间调度问题研究_模拟退火 三角形分布的双方叫价拍卖模型在电力交易中的应用* 两个不同成本结构订货商的供应链协调[1]_订货批量
推荐建筑资讯
点击排行
资讯热搜
产品热搜
百度合作产品

本站点已加入百度原创合作,切勿转载抄袭!