用最小二乘法求一个形如

时间:2023-11-26 02:42:33 阅读: 最新文章 文档下载
说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。
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 特征方程为

特征根均小于1Jacobi迭代法收敛。 Gauss-Seidel迭代矩阵

00.40.2

G(DL)1U00.40.7

00.040.17

32

特征方程为 |IG|0.570.0960 特征根均小于1Gauss-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.0931Jacobi迭代法不收敛;

矩阵A对称正定,故Gauss-Seidel迭代法收敛。 (b) 谱半径(B)01Jacobi迭代法收敛; 谱半径(B)21Gauss-Seidel迭代法不收敛


本文来源:https://www.wddqw.com/doc/8691af5ea45177232f60a2f7.html