数学与信息科学学院信息与计算科学：改进的平方根法及其程序实现

1、 毕 业 论 文 题 目: 改进的平方根法及其程序实现 学 院: 数学与信息科学学院 专 业: 信息与计算科学 毕业年限: 2011年6月 学生姓名: 张cc 学 号: 200771020141 指导教师: 郭zz 目 录摘 要20 引 言31 预备知识31.1 分解理论31.2 cholesky分解法 41.3 算法描述 5 2 改进的平方根法 63 分解算法描述 74 应用举例85 程序实现 105.1 程序码源 105.2 实例计算 126 结束语 13参考文献 14致 谢 15改进的平方根法及其程序实现摘要: 针对对称正定方程组的解法, 本文先对cholesky分解法进行了分析研究,在

2、此基础上给出了改进的平方根法(即分解法), 此方法有效地避免了原平方根法开方运算所带来的误差和不便, 并通过算法描述、实例计算, 用c程序实现了分解, 进一步提高了矩阵运算的速度和精度.关键词: 对称正定矩阵, 平方根法, 分解, 算法improved methods of square root and realization of its program zhang dengke(college of mathematics and information science, northwest normal university, lanzhou, gansu, 730070)abstra

3、ct: aiming at studying solutions of symmetric positive definition matrix in linear equations. initially, the text has conducted a series of analyses and researches towards decomposition proposed by cholesky. then based on theses researches and analyses, it offers the improved methods of squareroot (

4、also called decom-position), which effectively avoid some errors and inconvenience brought by the process of extracting root. at the same time, it achieves the decomposition through the means of algorithm description, example calculation as well as applicat-ion of c program, further enhancing the sp

5、eed and accuracy in matrix operation. key words: symmetric positive definition matrix, method of square root, decomposition, algorithm0 引言很多工程中的科学计算, 例如应用有限元法解结构力学问题时, 最后往往归结为求解系数矩阵为对称正定方程组解的问题. 由于对称正定矩阵各阶顺序主子式以及全部的特征值均大于零, 这种特征也使得其三角分解具有更为简单的形式, 不同的分解也导出了一些不同的解法. 平方根法(即cholesky分解法), 就是利用对称正定矩阵的三角分解

6、而得到的求解对称正定方程组的一种有效方法, 其计算量和存储量约为普通消去法的一半, 且无需选主元就能求得较为精确的数值解, 但由于在平方根法中含有多次开方运算, 因此给计算带来了许多不便, 而在原平方根法的基础上, 给出改进的平方根法(即分解法), 成功避免了开方运算带来的的麻烦, 因此在各种工程计算中应用更为广泛.1 预备知识1.1 分解理论定理1.1（矩阵的分解定理） 设为阶矩阵, 如果的顺序主子式, 则可分解为一个单位下三角矩阵和一个上三角矩阵的乘积, 且这种分解是唯一的.定理1.2（对称正定矩阵的分解定理） 设为阶对称正定矩阵, 则必存在非奇异下的三角矩阵, 使 , (1.1) 并且当

7、的对角线元素均为正数时, 这种分解是唯一的.证明 因为对称正定, 则它的各阶顺序主子式均大于零. 由分解定理可知, 矩阵存在唯一的doolitle分解, 即 其中, 为单位下三角矩阵, 为上三角矩阵. 记 diag, 则为单位上三角矩阵, 且. 由的对称性有, 可得 由于为单位下三角矩阵, , 从而由分解的唯一性可知, , 于是 显然有 ,这样 , 故可令 diag则 其中,为非奇异的下三角矩阵. 而知限定的对角元素为正数时,这种分解是唯一的, 定理得证.1.2 cholesky分解法设为对称正定矩阵, 由上面分解定理可知, 则存在一个实的非奇异下三角矩阵, 使, 当限定的对角线元素为正数时,

8、 这时的分解是唯一的.设 (1.2)其中, , .由比较和的对应元素, 可求得的元素如下:由 , 得 , .假设的第列元素已经求得, 下面求的第列元素. 注意到 ,得计算公式, 对, 有 (1.3)由此可逐行求得的全部元素, 从而由及得到求解方程组的公式 (1.4)上述方法称为cholesky分解法. 由于计算的对角线元素要作次开平方运算, 故cholesky分解法又称平方根法.1.3 算法描述步1 输入对称正定矩阵和右端向量; 步2 cholesky分解: 对 计算: 步3 用向前消去法解下三角方程组: , 对 计算 ; 步4 解对角形方程组: 对 计算 ;步5 用回代法解上三角形方程组:

9、, 对 计算 .不难验证cholesky分解法的乘除计算总量约为, 为一般矩阵分解计算量的一半. 虽然如此, 但其增加的次开平方运算是非常不利的, 下面引出改进的平方根法分解法.2 改进的平方根法为了避免开方运算, 对作分解, 即, 则 / 提出矩阵的对角元素 / 由对称正定, 可得, 令 可证 即. (1.5)是对角元素为1的单位下三角矩阵.对矩阵作分解, 共计算个矩阵元素; 对称矩阵的分解, 只需计算个元素, 减少了近一半的工作量. 借助分解计算公式, 容易得到分解计算公式.设有分解形式其中.在分解中可套用分解公式, 只要计算下三角矩阵和的对角元素.计算中只需要保存的元素, 的行列的元素用

10、的表示. 由于对称正定矩阵的各阶主子式均大于零, 直接调用分解公式可完成分解计算.对于,有 (1.6). (1.7)3 分解算法描述步1 输入方程组阶数、系数矩阵和常数项.步2 for to ; ; for to ; ;步3 解方程组的步骤从略.由, 解方程组可分为三步完成:(1)解方程组, 其中. (1.8)(2)解方程组, 其中. (1.9)(3)解方程组. . (1.10)可以看出, 改进后的分解乘除运算量约为, 计算过程也无需开方运算, 使得其运算更加简单易行, 因此改进后的分解法相对于原平方根法具有更好的实效性和可行性.4 应用举例例 试用分解法求解方程组解 由, 可得的方程组, 令

11、, 则.计算步骤:(1)将直接分解为, 求出;(2)求解方程组;(3)求解方程组.现有, , , , , , , 即 由可得解得 由有 解得 5 程序实现 5.1 程序码源应用分解法解线性方程组, 为了便于计算, 编写如下程序, 通过计算机来实现线性方程组的求解, 简洁方便. 具体程序如下: purpose：分解法解方程组# include # include # define max_n 20 /(x_i,y_i)的最大维数/int main()int n; int i,j,k; int mi;double mx,tmp; static double amax_nmax_n,bmax_n,x

12、max_n,ymax_n,zmax_n; static double lmax_nmax_n,dmax_n; /输入ax=b的维数/ printf(n input n value(dim of ax=b):); scanf(%d,&n); if(nmax_n) printf(the input n is larger then max_n,please redefine the max_n.n);return 1; if(n=0) printf(please input a number between 1 and %d.n, max_n);return 1; /输入ax=b的a矩阵/ pri

13、ntf(now input the matrix a(i,j),i,j=0, %d:n, n-1);for(i=0;in;i+) for(j=0;jn;j+) scanf(%lf,&aij); /输入b矩阵/ printf(now input the matrix b(i),i=0, %d:n, n-1); for(i=0;in;i+) scanf(%lf,&bi); for(i=0;in;i+) lii=1; /u矩阵对角元素为1/ for(k=0;kn;k+) dk=akk; /计算d_i/ for(j=0;j=k-1;j+) dk-=(lkj *dj*lkj); for(i=k+1;in

14、;i+) /计算l矩阵/ lik=aik; for(j=0;j=k-1;j+) lik-=(lij *dj*lkj); lik/=dk; for(i=0;in;i+) /求解lz=b的z/ zi=bi; for(j=0;j=i-1;j+) zi-=lij*zj;for(i=0;i=0;i-) /求解ly=x/ xi=yi; for(j=i+1;jn;j+) xi-=lji*xj; printf(solve x_i=n); /输出/ for(i=0;in;i+) printf(%fn,xi);return 0;5.2 实例计算 用分解法求解方程组:程序输入输出input n value(dim

15、of ax=b): 3 now input the matrix a(i,j),i,j=0,2: 16 4 8 4 5 -4 8 -4 22 now input the matrix b(i),i=0,2: -4 3 10 solve x_i= 2.000000 4.000000 -2.2500006 结束语从矩阵分解角度看, 分解法与消元法本质上没有多大的区别, 其基本思想还是通过等价变换将线性方程组化为结构简单、易于求解的形式. 在实际问题中经常遇到对称正定矩阵方程组的求解问题, 对于这种具有特殊性质的系数矩阵, 采用改进的平方根法是一种行之有效的计算方法, 它是在普通消去法的基础上建立起

16、来的, 算法思想虽有些复杂, 但其计算量约为普通消去法的一半, 且有效地避免了开方运算所带来的不便, 因此改进后的分解法相对于原平方根法具有更好的实效性和可行性. 并通过c程序求解使得这类问题变得更加简单易行, 在大量工程计算中有着广泛而重要的应用.参考文献1 刘元亮. matlab在对称正定矩阵的改进平方根分解法中的应用j. 湖南: 怀化学院学报. 2004,(02).2 郭丽杰, 周硕, 秦万广. 对称矩阵的改进cholesky分解在特征值问题中的应用j. 吉林: 东北电力学院学报. 2003,(02).3 杜廷松, 沈艳君, 覃太贵. 数值分析及实验m. 北京: 科学出版社. 2007.

17、4 肖筱南, 赵来军, 党林立. 数值计算方法与上机实习指导m. 北京: 北京大学出版社. 2004.5 马昌凤, 林伟川. 现代数值计算方法m. 北京: 科学出版社. 2008.6 张韵华, 奚梅成, 陈长松. 数值计算方法与算法m. 北京: 科学出版社. 2000.7 李庆扬, 王能超, 易大义. 数值分析(4版)m. 武汉: 华中科技大学出版社. 2006.8 严克明, 欧志英, 刘树群. 数值计算基础m. 甘肃: 甘肃人民出版社. 2006.9 ferenc szidsrovzky, sidney yokowitz. principles and procedures of numer

18、ical analysism. 上海: 上海科学技术文献出版社. 1982.10 rainer kress. numerical analysism. 北京: 世界图书出版社. 2003.致 谢本论文经过将近半年的努力即将告以尾声, 从论文的选题到最后完稿, 整个过程都经过了认真地考虑、仔细地查阅和细心地修改. 在此, 我首先要感谢我的导师郭晓斌老师, 不管是我的论文选题还是论文的撰写, 以及资料的查阅等各方面, 他都给了我莫大的帮助与启发, 尤其是在论文的几次修改过程中, 郭老师以他广博的专业学识、严谨的治学精神和他耐心的指导态度, 才使我的论文能顺利完成. 谨再次向郭老师致以崇高的敬意和诚

19、挚的感谢. 也祝愿郭老师身体健康、工作顺利, 在学术研究上取得更加辉煌的成就, 为更多的学子导航.同时, 我也要对给予本论文参考文献的所有学术专家和老师致以真挚的谢意, 是他们出版的书籍与发表的学术论文给了我很大的启示与指导, 才将论文完成.其次, 在即将毕业之际, 我要感谢数信学院所有的老师对我的细心教育与培养, 让我在四年的学习生涯中不仅学到了扎实的专业知识, 而且他们的言传身教使我受益匪浅, 他们严谨的治学态度和耐心教导学生的博爱精神也是我永远学习的榜样, 并将积极影响着我今后的学习和工作.最后, 我还要再次感谢我的母校西北师范大学给了我发展的平台, 感谢各位老师四年来对我的耐心教导和栽

20、培. 2011年5月 気持今flying get！dou都斗豆逗陡抖痘兜读蚪窦篼蔸乧侸兠凟剅吺唗投斣枓梪橷毭氀浢渎渎瞗窬窦脰艔豆读逾郖酘酡钭鋀钭閗闘阧餖饾斗鬦鬪鬬鬭du读度毒渡堵独肚镀赌睹杜督都犊妒顿蠹笃嘟渎椟牍黩髑芏儥凟剢剫匵厾噣土涂妬嬻剬塅媏彖断毈瑖碫篅簖缎专腶葮褍踹躖锻鍴mr najib met search crews at pearce raaf base near perth on thursday morning, before their planes left for the day, and then later held talks with mr abbott.the

21、disappearance of mh370 has tested our collective resolve, he told a news conference.faced with so little evidence, and such a herculean task, investigators from malaysia, the us, the uk, china, australia and france have worked without pausehe thanked both search teams and the australian government f

22、or their efforts in recent weeks, and said the search would go on.continue reading the main storymh370 - facts at a glance 8 march:malaysia airlines kuala lumpur-beijing flight carrying 239 people disappears planes transponder, which communicates with ground radar, was switched off as it left malays

23、ian airspace satellite pings indicate plane was still flying seven hours after satellite contact was lost 24 march:based on new calculations, malaysian pm says beyond reasonable doubt that plane crashed in southern indian ocean with no survivors what we know the search for flight mh370i know that un

24、til we find the plane, many families cannot start to grieve. i cannot imagine what they must be going through. but i can promise them that we will not give up, he said.malaysian authorities have come in for heavy criticism over their management of the search, especially from relatives of thweather c

25、onditions were fair, with visibility of approximately 10km (6 miles), the joint agency coordination centre (jacc) - which is overseeing the search - said.the british submarine hms tireless is also in the southern indian ocean and is due to be joined by royal navy ship hms echo.the australian navy sh

26、ip ocean shield is heading to the region and has equipment for detecting the planes black-box flight recorder.experts say timing is critical as the flight recorder may only have enough battery power to send out a signal until 7 april.air chief marshall angus houston, head of the jacc, warned that th

27、e search operation faced multiple difficulties.this is one of the most demanding and challenging search and rescue operations, or search and recovery operations, that i have ever seen and i think probably one of the most complex operations of this nature that the world has ever seen, he told mr naji

28、b and mr abbott.on wednesday malaysian police chief khalid abu bakar said investigators had cleared all passengers of possible involvement in hijacking, sabotage or having personal or psychological problems that could have been connected to the disappearance.but he said that the criminal investigati

29、on could go on and on and on. we have to clear every little thing.at the end of the investigations, we may not even know the real cause. we may not even know the reason for this incident, he added.the police chief said that more than 170 interviews had been conducted with family members of the pilot

30、s and crew members, and that even cargo and food served on the plane were being investigated in case ofdui对队堆兑敦镦碓怼憝兊兑垖埻塠夺夺対对嵟憞怼捶杸濧濻瀢瀩痽磓祋綐膭薱謉譈譵追鈗锐锐錞鎚镦鐜锐陮队頧鴭dul乧dun吨顿蹲墩敦钝盾囤遁不趸沌盹镦礅炖砘伅俊吨墪壿庉忳敦憞撉撴楯橔潡炖犜獤碷腞腯蜳豚踲蹾趸逇遯钝镦鐜顿驐duo多朵夺舵剁垛跺惰堕掇哆驮度躲踱沲咄铎裰哚缍亸仛兊兑兑凙刴剟剫吋喥嚉嚲垜埵堕墯夛夺奲媠嫷尮崜嶞憜挅挆捶揣敓敚敠敪朶杂杕枤柁柂柮桗棰椯椭毲沰沱泽痥硾缍茤袳詑誃貀趓跢跥跿躱軃郸郸

31、鈬錞锗铎锗陀陊隋隋隓飿饳驮駄鬌鮵鵽点e饿哦额鹅蛾扼俄讹阿遏峨娥恶厄鄂锇谔垩锷阏萼苊轭婀莪鳄颚腭愕呃噩鹗屙亚亜亚伪佮侉偔伪伪僫匎匼卾吪呝咢咹哑唖啈啊啐哑恶囐囮垭垭垩堨堮妸妿姶娾娿媕屵岋峉峩崿庵廅悪恶戹搕搤搹擜曷枙桠櫮唉歞歹歺洝涐湂猡珴琧痷皒睋砈砐砨砵硆硪磀礘胺蒍蕚蘁蚅蝁覨讹咯誐谔譌讍豟轭軶輵迗遌遻邑鈋锇锷鑩閜阏阨阸隘頞頟额颚饿餩騀鬲魤魥鳄鰪鳄鴳鵈鹅鵞鹗齃腭齾青年（）、技术（）、恋爱（）、翡翠（）、読书（）、人（，）、幸福（）、训读词汇：青（）、术（）、恋（）、好（）、読（）、人（）、幸（）、汉字“海”可以读成“”也可以读成“”。“”是它的音读，“”是它的训读。“史宪巨”等汉字一般只使用它们的一种音读，“贝又咲（）”一般只使用它们的一种训读，但这类汉字较少。日语里的汉字通常都有两三种甚至更多的音训读发