- 1、本文档共19页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析大作业牛顿下山法,拉格朗日法,切比雪夫法及Matlab程序
课 程 设 计
课程名称: 数值分析 设计题目: 学 号: 姓 名:
完成时间: 2014.11.18
题目一: 解线性方程组的直接法
设方程组,其中
,
矩阵中,由相应的矩阵元素计算,使解向量。
不变,对的元素加一个扰动,求解方程组;
不变,对的元素和分别加一个扰动,求解方程组;
对上述两种扰动方程组的解做误差分析。
一.数学原理:
本计算采用直接法中的列主元高斯消元法,高斯列主元消元法原理如下:
1、设有n元线性方程组如下:
=
2、
第一步:如果a11!=0, 令
li1= ai1/a11, I= 2,3,……,n
用(-li1)乘第一个方程加到第i个方程上,得同解方程组:
a(1)11 a(1)12 . . . a(1)1n x1 b(1)1
a(1)21 a(1)22 . . . a(1)2n x2 b(1)2
. . . . . . . = .
a(1)n-11 a(1)n-12 . . a(1)n-1n xn-1 b(1)n-1
a(1)n1 a(1)n2 . . . a(1)nn xn b(1)n
简记为:
A(2) x = b(2)
其中
a(2)ij = a(1)ij – li1 * a(1)1j , I ,j = 2,3,..,n
b(2)I = b(1)I – li1 * b(1)1 , I = 2,3,...,n
第二步:如果a(2)22 != 0,令
li2= a(2)i2/a(2)22, I= 3,……,n
依据同样的原理,对矩阵进行化间(省略),依次下去,直到完成!
最后,得到上三角方程组:
a(1)11 a(1)12 . . . a(1)1n x1 b(1)1
0 a(1)22 . . . a(1)2n x2 b(1)2
. . . . . . . = .
0 0 . . a(n-1)n-1n xn-1 b(n-1)n-1
0 0 . . . a(n)nn xn b(n)n
简记为:
A(n) x = b(n)
最后从方程组的最后一个方程进行回代求解为:
Xn = b(n) / a(n)nn
Xi = ( b(k)k - a(k)kjxj ) / a(k)kk
二.解题过程:
1.由题中所给条件可求出b。
B =
6.0000
7.7156
9.9299
12.7560
16.3238
20.7813
(1) 不变,对的元素加一个扰动,求解方程组。
B =[6.0000 7.7156 9.9299 12.7560 16.3238 20.7813+0.0001]
解得x =[0.5997 2.6920 -1.8500 3.3917 0.0000 1.1667]’
(2)不变,对的元素和分别加一个的扰动,求解方程组。
A =
1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000
1.000000000000000 1.100001000000000 1.210000000000000 1.331000000000000 1.464100000000001 1.610510
文档评论(0)