Halcon矩阵操作最小二乘

来源：互联网发布：淘宝卖家开通花呗条件编辑：程序博客网时间：2024/06/07 04:05

*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]:=1endfortuple_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)

0 0