質問: a + bx + cy + dxy + ex^2 + fy^2
回答: いまちょっと検索してみましたが,見付かりませんでした. (あまりトライしてませんが・・・) 1,2と3の行列aとベクトルbを求めるところまでは, 手計算で式をぐりぐりします. 簡単な例として, 回帰曲線y(x)=a+bx+cx^2 計測点(xi,yi) の場合,手順1で, r=σ(y(xi) - yi(xi))^2 =σ(a+bxi+cxi^2 - yi)^2 これを手順2で, ∂r/∂a= 2σ(a+bxi+cxi^2 - yi)=0 ∂r/∂b= 2σ(xi×(a+bxi+cxi^2 - yi))=0 ∂r/∂c= 2σ(xi^2×(a+bxi+cxi^2 - yi))=0 と言うふうにします.これを展開して,係数a,b,c毎に分け, a,b,cの係数のかかっていない項を右辺にもってけば, 連立方程式が完成します. ここまではまず手で式を書いて,それをプログラムに反映させます. 3の連立方程式を解くものは,一番簡単なガウスの消去法とか, ガウスジョルダン法とか,lu分解法とかあります, この辺は数値計算の基本なので, 大概の本に載っていますし,ネット上にもあるような気がします. こんなものですみません.