Halcon矩阵操作最小二乘

*
X0:=[1.2,1.2,3.1,4.3,5.2,6.4]
X1:=[2.2,2.44,3.2,4.15,5.00,6.1]
X2:=[1.4,2.5,3.3,4.35,5.1,6.2]
X3:=[1.45,2.37,3.12,4.44,5.2,6.33]
X4:=[1.56,2.78,3.1,4.56,5.3,6.44]
X5:=[1.56,2.5,3.23,4.78	,5.4,6.55]
X6:=[1.56,2.79,	3.42,4.79,5.55,6.66]
X7:=[3.4,5.6,3.23,4.5,5.66,6.7]
X8:=[7.8,8.9,3.45,4.67,5.77,6.8]
X:=[]
tuple_concat(X,X0,X)
tuple_concat(X,X1,X)
tuple_concat(X,X2,X)
tuple_concat(X,X3,X)
tuple_concat(X,X4,X)
tuple_concat(X,X5,X)
tuple_concat(X,X6,X)
tuple_concat(X,X7,X)
tuple_concat(X,X8,X)
Y:=[33.15,34.878,34.79,34.939,36.02,36.487,37.497,42.469,52.118]
create_matrix(9, 6, X, Xmat)
transpose_matrix(Xmat,  XmatT)
get_full_matrix(XmatT, Values1)
Ones:=[]
for i:=0 to 8  by 1  
    Ones[i]:=1
endfor
tuple_concat(Values1,Ones,Values1)
create_matrix(7, 9, Values1, Xmat)
transpose_matrix( Xmat,Xmat)
get_full_matrix(Xmat, Values2)
create_matrix(9,1,Y,YMat)

transpose_matrix(Xmat, Xt)
mult_matrix(Xt,Xmat,'AB', AtA)
invert_matrix(AtA, 'general', 0, AtAinv)
mult_matrix(AtAinv,Xt,'AB', MatrixMultID)
mult_matrix(MatrixMultID,YMat,'AB', MatrixMultID1)
get_full_matrix(MatrixMultID1, Values)