资源描述:
《数值分析讲稿8_黑白).pdf》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、§4.高斯消去法的变形二、平方根法工程实际计算中,线性方程组的系数矩阵常常具有对称正定性,即其各阶顺序主子式及全部特征值均大于零。矩阵的这一特性使它的三角分解也有更简单的形式,从而导出一些特殊的解法,如平方根法与改进的平方根法。定理:设是对称正定矩阵,则存在唯一的非A奇异下三角阵使得L,TALL且的对角元素皆为正。L定理证明(1)证:因对称正定,其各阶顺序主子式均大于零,故有AALU其中为单位下三角矩阵,为上三角阵。LU1令Ddiagu(,,),LuPDUP,则为单位上三角阵。11nnu12u1n1Luu1112Luu1n11uu1111
2、uuLMOOU2222DPOMOOunn1,uuunnnnnn1,11TTT故ALULDPAPDLTT由LU分解的唯一性PLAPDP定理证明(2)由于是对称正定的其顺序主子式均大于零。故A,uD0,uDD/0(i2,3,L,)n111iiii1T令Ddiagu(,,)LuDDD11nnTTTTT则APDPPDDP()DPDPLLT其中LD()P为非奇异下三角阵且对角元素皆为正数。,唯一性:假定存在非奇异下三角阵GL,其对角元TT素皆为正数,且使得A
3、LLGG于是有TT11TT11TT11LG()LLLG()LGGG()LGTT11因()LG为上三角阵,LG为下三角阵,故由上式得TT11LG()LGI即GL,与假设矛盾。平方根(Cholesky分解法)法ll1100LL11l21ln1llOM0lLl由ALLT212222n2MMO0MOOMllnn12LLlnn00lnn其中li0(1,2,L,).n由矩阵乘法及l0(当jk时),iijkj1122laljj()jjjk(j1,2,L,n),k1得
4、j1lalij()ijikllijjk/jj(1,L,);nk10这里规定0。计算顺序是按列进行,即k1ll(2i,3,,)LLnll(3i,,)nL。11ii1222当矩阵完成平方根分解后,求解AAxb,即求解两个三角形方程组T(1)Lyb,求y;2()Lxy,求x.i1yblii()ikyliki/i(1,2,L,n).k1nxyii()lkixlink/ii(,n1,L,1).ki13由于的对称性,平方根法的乘除运算量为An/6数量级,约是Gauss消去法的一半。上机计算时,所需存储单元也少,只要
5、存储的下三角部分和右端项,计算中存AbL放在的存储单元,Ay,x存储在的存储单元。b但这种方法在求时需作次开方运算,这样又增加了Ln计算量。为了避免开方,可使用改进的平方根方法。改进平方根法1d1ld1ALDLT212MOOllnn12L1dn1ll21Ln11ln2OM1其中llj1,0(k),由比较法得jjjk对于in1,2,L,,i12daiiildikk;k1i1lal()dldjin/(1,L,).jijijkkikik1上式虽避免了
6、开方运算,但增加了相乘因子,引进变量tld对于i1,2,L,,n有ikikki1daiiitlikik;k1i1lat()ldjin/(1,L,).jijiikjkik1TT对称正定矩阵按ALDL分解和按LL分解计算T量差不多,但LDL分解不需要开方计算。T求解LybDLx,y计算公式yb;11i1ybiilyiikk(2,,)Ln.k1xyd/;nnnnxydiii/(lxinkik1,L,2,1).ki1三、追赶法在数值计算中,如三次样条插值或用差分方法解常微分方程边值问题
7、,常常会遇到求解以下形式的方程组bc11xd11abcxd22222OOOMMabcxd简记Axd.iiiiiOOOMMabcxdnnnn1111n1abxdnnnn此系数矩阵的非零元素集中分布在主对角线及其相邻两次对角线上,称为三对角矩阵。方程组称为三对角方程组。定理:设三对角方程组系数矩阵满足下列条件:bc011bacaci0(2,3,L,n1)i