【#文档大全网# 导语】以下是®文档大全网的小编为您整理的《用最小二乘法求一个形如》,欢迎阅读!
1. 用最小二乘法求一个形如yabx的经验公式,使它与下列数据拟合,并求均方误差.
19 25 31 38 44
xi
2
yi
19.0 32.3 49.0 73.3 97.8
解:
10654a14748998b738643.0010a10654b542.80ba,,
解方程得a4.00955,b0.0471846,均方误差13.0346。
2.下述矩阵能否分解为LU(其中L为单位下三角阵,U为上三角阵)?若能分解,那么分
解是否唯一?
123111126
,B221,C2515.A241
46733161546
解: 按高斯消去法,A无法进行第二次消去,换行后可以分解,B第二次消去可乘任意系数,分解不唯一,C可唯一分解。
3.设方程组
(a) 考察用雅可比迭代法,高斯-塞德尔迭代法解此方程组的收敛性;
(k1)(k)4||xx||10(b) 用雅可比迭代法,高斯-塞德尔迭代法解此方程组,要求当时迭代
5x12x2x312
x14x22x3202x3x10x3
231
终止.
解: (a) Jacobi迭代矩阵
0.40.201
BD(LU)0.2500.5
0.20.30
3
|IB|0.210.0550 特征方程为
特征根均小于1,Jacobi迭代法收敛。 Gauss-Seidel迭代矩阵
00.40.2
G(DL)1U00.40.7
00.040.17
32
特征方程为 |IG|0.570.0960 特征根均小于1,Gauss-Seidel迭代法收敛。 (b) Jacobi迭代格式为
X(k1)BX(k)f1
1T
fDb(1.250.3)1其中B如上,, 迭代18次得
X3.9999964
2.99997391.9999999
T,
Gauss-Seidel迭代格式为
X(k1)GX(k)f2 1T
其中G如上,f2(DL)b(2.42.61.53), 迭代8次得 X4.000036
5. 设方程组
2.999985
2.000003。
T
x10.4x20.4x31
0.4x1x20.8x320.4x0.8xx3
123(a) (b) x12x22x31
x1x2x312x2xx1
231
试考察解此方程组的雅可比迭代法及高斯-塞德尔迭代法的收敛性
解:. (a) 谱半径(B)1.0931,Jacobi迭代法不收敛;
矩阵A对称正定,故Gauss-Seidel迭代法收敛。 (b) 谱半径(B)01,Jacobi迭代法收敛; 谱半径(B)21,Gauss-Seidel迭代法不收敛
本文来源:https://www.wddqxz.cn/03285e77f41fb7360b4c2e3f5727a5e9846a2720.html