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)

————————————————
版权声明：本文为CSDN博主「青箬笠绿蓑江船行」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/qq_14844633/java/article/details/71123665