Bước 1

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

Đầu ra

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

Bước 2

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

Đầu ra

Bước 3

X = PseudoInverse[A].B

Đầu ra

{0.5857142857142857, 1.1714285714285715}

Bước 4

A.X

Đầu ra

{2.928571428571429, 5.857142857142858, 8.785714285714285}

Bước 5

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

Không có đầu ra - bước này chuẩn bị cho bước kế tiếp.

Bước 6

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

Đầu ra

-Graphics-

Bước 7

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

Không có đầu ra - bước này chuẩn bị cho bước kế tiếp.

Bước 8

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

Đầu ra

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

Bước 9

PseudoInverse[mat].yData

Đầu ra

{3.8170000000000015, 3.585}

Bước 10

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

Đầu ra

{α -> 3.8169999997513093, β -> 3.5849999944044773}

Bước 11

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

Đầu ra

-Graphics-

Các hàm được dùng

Linear algebra matrix computations least squares solutions