matlab中cgs收敛曲线,求解线性方程组 - 共轭梯度二乘法

通过为 cgs 提供函数句柄来求解线性方程组，该函数句柄计算 A*x 而不是系数矩阵 A。

gallery 生成的 Wilkinson 测试矩阵之一是 21×21 三对角矩阵。预览该矩阵。

A = gallery('wilk',21)

A = 21×21

10 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

1 9 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 1 8 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 1 7 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 1 6 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 1 5 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 1 4 1 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 1 3 1 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 1 2 1 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0

⋮

Wilkinson 矩阵有特殊的结构，因此您可以用函数句柄来表示 A*x 运算。当 A 乘以向量时，所得向量中的大多数元素为零。结果中的非零元素对应于 A 的非零三对角元素。此外，只有主对角线具有不等于 1 的非零值。

表达式 Ax 变为：

Ax=[1010⋯⋯⋯001910001810⋮⋮0171001610⋮⋮0151001410⋮⋮013⋱000⋱⋱100⋯⋯⋯0110][x1x2x3x4x5⋮⋮x21