3.2 线性方程组求解.pdf

下载文档

0
0
约9.08千字
约 9页
2021-01-11 发布于广东
举报
版权申诉
保障服务

3.2 线性方程组求解.pdf

1、本文档共9页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

程序设计与问题求解 ——线性方程组求解线性方程组求解设线性方程组 Ax b 的系数矩阵A 可逆且主对角元素 a00 ,a11 , ,a( n1)( n1) 均不为零，求解 x • 迭代法雅克比迭代、松弛迭代等 • 解析法高斯消去法、行列式法等雅克比迭代法  数学原理 Ax b A (A D ) D 的系数矩阵 A 对角元素均不为零，将 A 分解成： a00    a 这里 D  11   ...    a  (n1)(n1)  1 1 从而原方程可写成：Dx (D A)x b x (I D A)x D b (k 1) 1 (k ) 1 于是有雅克比迭代式： x (I D A)x D b x (0) (x0(0) ,x1(0) , ,xn-1(0) )T 为初始向量 1 n1 写成分量的形式： x (k 1) [b a x (k ) ] i 0,1, , n 1 i i ij j aii j 0 j i 雅克比迭代法  关键代码分析 //前一次近似解存放在数组x2[]中 //本次近似解存放在数组x[]中设系数存放于二维数组a[][]，常 for( i=0; in; i++) { 数项存放于数组b[]，根据下面 i 对应xi tmp=0.0; 公式编程：