- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
线性代数计算方法的
第3章 线性代数计算方法 ; 在自然科学和工程技术中很多问题的解决常常归结为解线性代数方程组。
例如:
电学中的网络问题,
船体数学放样中建立三次样条函数问题,
用最小二乘法求实验数据的曲线拟合问题等;Evaluation only.
Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0.
Copyright 2004-2011 Aspose Pty Ltd.;
直接法:
经过有限次运算后可求得方程组精确解的方法(不计舍入误差!)
迭代法:
从解的某个近似值出发,通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解);§1 高斯消去法 ;且aii≠0,i=1,2,…,n;下三角阵?;化成上(下)三角!;高斯消去法;②-①×2;第2行:计算比例因子;-1;第 i 行 - li1 ? 第1行,得到:;第k步:消去;消元过程总体流程:;(2)回代过程:;二、 选主元消去法;例如:用高斯消去法求解下列方程组(用四位有效数字计算):;化简可得
x2=0.6000
回代求得
x1=105(0.6-0.6000)=0
而方程组的解应为
x1=0.4000
x2=0.6000
; 显然用上述方法求出的解x1与方程组的实际解相差很大。若改变两个方程的顺序,即
x1+x2=1 ①
10-5 x1+x2=0.6 ②
②-①×10-5得
(1.000-1.000×10-5)x2=0.6-1.000×10-5
0.99999x2=0.59999
化简得
x2=0.6000
回代求得
x1=(1-0.6000)=0.4000; 高斯主元素消去法是顺序消去法的一种改进。它的基本思想是在逐次消元时总是选绝对值最大的元素(称之为主元)做除数,按顺序消去法的步骤消元。
这里主要介绍求解线性方程组最常用的列主元素消去法和全主元素消去法。; 所谓列主元素消去法就是在每一步消元过程中取系数子矩阵的第一列元素中绝对值最大者作主元。对线性方程组进行n-1次消元后,可得到上三角形方程组; 取四位有效数字计算。
解 ②中-18为主元,交换②和①得 ; ②+①×12/18,③+①×1/18得 ; ③+②×1/1167得 ;开始;找列最大值;/*本算法用高斯列主元素消去法求解矩阵方程AX=B。
其中:A是N×N矩阵;B是N×1矩阵。
输入:n—A的行数;
a—二维矩阵A
b—矩阵B
算法结束后,函数返回值为ERROR_CODE时,表示A是奇异的或病态的;否则,A代表行列式的值。
a— A消元后的上三角矩阵
b—矩阵方程的解X
*/;Evaluation only.
Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0.
Copyright 2004-2011 Aspose Pty Ltd.;double Gaussian_elimination(int n, double a[n][n], double b[n])
{
int i, j, k, mk;
double mm, f;
for (k=0; kn-1; k++)
{
mm=a[k][k]; mk=k;
for (i=k+1; in; i++ )
if (fabs(mm)fabs(a[i][k]))
{mm=a[i][k]; mk=i;}
if (fabs(mm)=0) return (ERR
文档评论(0)