matlab 线性最小二乘法,求解线性方程组 - 最小二乘法

通过为 lsqr 提供用来计算 A*x 和 A'*x 的函数句柄(而非系数矩阵 A)来求解线性方程组。

创建一个非对称三对角矩阵。预览该矩阵。

A = gallery('wilk',21) + diag(ones(20,1),1)

A = 21×21

10 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

1 9 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 1 8 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 1 7 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 1 6 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 1 5 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 1 4 2 0 0 0 0 0 0 0 0 0 0 0 0 0

0 0 0 0 0 0 1 3 2 0 0 0 0 0 0 0 0 0 0 0 0

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

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

⋮

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

表达式 Ax 变为：

Ax=[1020⋯⋯01920⋮01⋱20⋮010⋱⋱⋮0⋱1⋱0⋮⋱⋱⋱20⋯⋯0110][x1x2x3⋮⋮x21]=[10x1+2x2x1+9x2+2x3⋮⋮x19+9x20+2x21x20