第 1 步

A = {{1., 2.}, {5., 6.}, {4.5, 6.}};
B = {5., 6., 8.};
LinearSolve[A, B]

输出

LinearSolve[{{1., 2.}, {5., 6.}, {4.5, 6.}}, {5., 6., 8.}]

第 2 步

A = {{1., 2.}, {2., 4.}, {3., 6.}};
B = {5., 6., 8.};
MatrixRank[A]

输出

第 3 步

X = PseudoInverse[A].B

输出

{0.5857142857142857, 1.1714285714285715}

第 4 步

A.X

输出

{2.928571428571429, 5.857142857142858, 8.785714285714285}

第 5 步

data = {{0.2, 4.59}, {0.4, 5.05}, {0.6, 6.2}, {0.8, 6.6}, {1., 7.4}};

无输出——这一步是在为下一步做准备。

第 6 步

plot1 = ListPlot[ data, PlotStyle -> PointSize[0.02]]

输出

-Graphics-

第 7 步

xData = data[[All, 1]];
yData = data[[All, 2]];

无输出——这一步是在为下一步做准备。

第 8 步

mat = Transpose[ {Table[1, {Length[xData]}], xData}];
mat//MatrixForm

输出

MatrixForm[{{1, 0.2}, {1, 0.4}, {1, 0.6}, {1, 0.8}, {1, 1.}}]

第 9 步

PseudoInverse[mat].yData

输出

{3.8170000000000015, 3.585}

第 10 步

sol = FindFit[data, α + β x, {α, β}, x]

输出

{α -> 3.8169999997513093, β -> 3.5849999944044773}

第 11 步

Plot[ (α + β x) /. sol, {x, 0, 1}, Epilog -> First[plot1]]

输出

-Graphics-

用到的函数

Linear algebra matrix computations least squares solutions