MATLAB Function Reference

spaugment

最小二乗拡大システムの作成

表示

S = spaugment(A,c)

詳細

S = spaugment(A,c) は、スパース、正方、対称半正定行列 S = [c*I A; A' 0] を作成します。行列 S は、最小二乗問題に関連しています。

min norm(b - A*x)

by

r = b - A*x
S * [r/c; x] = [b; 0]

残差スケーリングファクタの最適値 c は、min(svd(A)) と norm(r) を含んでいて、通常、計算に要するコストを非常に必要とします。

S = spaugment(A) において、c を指定しない場合、max(max(abs(A)))/1000 を使用します。

注意   MATLAB の前のバージョンで、拡大行列は、スパース線形方程式ソルバ、 \ と非正方問題に対して / で使われていました。現バージョンでは、A の qr 分解を使った最小二乗ソルバを使います。

参考

spparms

sparse

spconvert