判定h矩阵个迭代算法修正翻译

判定h矩阵个迭代算法修正翻译

ID:23020274

大小:302.50 KB

页数:7页

时间:2018-11-02

判定h矩阵个迭代算法修正翻译_第1页
判定h矩阵个迭代算法修正翻译_第2页
判定h矩阵个迭代算法修正翻译_第3页
判定h矩阵个迭代算法修正翻译_第4页
判定h矩阵个迭代算法修正翻译_第5页
资源描述:

《判定h矩阵个迭代算法修正翻译》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、判定H-矩阵的一个迭代算法的修正*李耀堂,张讯(云南大学数学系,云南昆明650091)摘要:随着H一矩阵在科学与工程计算中的广泛应用,如何判定一个给定矩阵是否为H一矩阵引起了许多研究者的兴趣.本文对一个现有判定H一矩阵的迭代算法进行了修正,得到了一个新的迭代算法.数值算例表明该算法是有效的.关键词:H-矩阵;对角占优矩阵;迭代算法令表示所有矩阵,.的比较矩阵定义如下:如果是M-矩阵,那么矩阵是H-矩阵.众所周知H-矩阵通常称作广义对角占优矩阵(GDDM).那就是,存在一正对角矩阵使得为严格对角占优矩阵(SDDM).一类H-矩阵是在工程和科学计算中具有广泛应用的一类

2、矩阵.许多迭代算法求解线性方程是收敛的如果系数矩阵是H-矩阵.因此,如何判断一个给定的矩阵是H-矩阵在以上研究领域起着重要作用.在文章[2]—[6],一些研究者给出了一些方法来确定H-矩阵.尽管这类方法有很多种,但是它们都来源于一个事实:试着去寻找一个正对角矩阵使得为严格对角占优矩阵.本文的方法也是基于这一点.做为改性的方法,我们从两对角占优行和非对角占优行入手.这在之前的学习研究中被忽略.1、H-矩阵的一些引理为了不失一般性,令为一个不可约矩阵,在本文中,我们使用以下符号:;;;;;为了不产生歧义,我们用表示,表示,表示,表示.显然,,,,.引理1令为H-矩阵,

3、.引理2令为H-矩阵,.引理3令为H-矩阵,.引理4令为H-矩阵,是一个正对角矩阵,当且仅当为H-矩阵时为H-矩阵.引理4指出了和具有同一属性的对角优势.我们的任务是确定一个正对角矩阵使得为严格对角占优矩阵.2、迭代标准矩阵算法A(1)输入矩阵,如果,此时不是H-矩阵,停止并且输出“不是H-矩阵!”;(2)如果,此时不是H-矩阵,停止并且输出“不是H-矩阵!”;如果,此时是H-矩阵,停止并且输出“是H-矩阵!”.(1)令显然,,,(当时可约),现取,这时,使,这里.(4)计算,回到(2)算法A是从对角占优行开始,找到一个正对角矩阵使得为严格对角占优矩阵.显然,如果

4、这个算法终止在有限迭代次数,这时我们可以得到一个明确的结论:时A是一个H-矩阵,或者说时A不是H-矩阵.用表示k重正对角矩阵,表示k重生成矩阵,通常的,.引理5在算法A中,.证明一般情况下,当为一重正对角矩阵时,我们只需证明.令,.,,.那就是,即.引理5指出,在算法A中,的对角占优行数随迭代次数的增加而增加.所以经过有限次迭代,我们可以得到一个正对角矩阵,当为H-矩阵时,为严格对角占优矩阵.但是这里有两个问题需要解答:首先,如果不是H-矩阵会怎么样;第二,算法什么时候终止?这些问题会在我们给出我们的改进算法后得到解答.算法B(1)输入矩阵,如果或者,此时不是H-

5、矩阵,停止并且输出“不是H-矩阵!”;(2)令,,;(3)如果,此时是H-矩阵,停止并且输出“是H-矩阵!”;(4)令选取,这时,规定,这里(5)使,回到(3).注意:很显然.这说明对所有,输入,减少速度比其他的要快.尽管的对角元素的非对角占优行在减少,但输入的其他的在同一行的减少比也同样如此甚至比这更多.3、算法的定理证明定理1算法B中(i);(ii).证明为了完成证明,我们列举三个例子.1)令,这时,并且,.所以=那么,那就是说.结论(i)得证.2)令,这时且.那么有.因为且,因此因为且,因此,并且.1)令,此时.那么有.因为且.此时.从2)和3)中,我们得出

6、了(ii)的结论综上所述,我们可以得到明确的结论;算法A和B都能保持矩阵的对角占优性.但是,作为一个改性方法,算法B能同时处理对角占优和非对角占优列,于是我们可以我们可以说如果矩阵没有对角占优行,那么算法B的收敛速度比算法A更快.定理2算法B在有限的迭代终止或产生一种独特的无限序列对所有都有.证明通过定理1我们知道在算法B中

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。