关于线性最小二乘法:https://blog.linjincan.cn/archives/362

Mx:=[100:10:500]
tuple_length(Mx,len)
tuple_gen_const(len,5,r)
Ma:=2
Mb:=40
tuple_rand(len , noise)
My:= Ma *Mx + Mb*noise
* (My,Mx):数据点集
gen_circle(ContCircle, My, Mx, r) 
*创建一个My矩阵y
create_matrix(len,1,My,y)
*创建一个两列的Mx矩阵x,第一列为Mx(Mx肯定为满秩),第二列为1,y=kx+b,b的系数为1
create_matrix(len,2,1,x)
set_value_matrix(x, [0:len-1], gen_tuple_const(len, 0),Mx)

* XT 代表X的转置 inv(*)代表*的逆,beta:[k,b]
* x beta = y
* xT x beta = xT y
* beta = inv( xT x) xT y
mult_matrix(x,x,'ATB',xtx)
mult_matrix(x,y,'ATB',xty)
invert_matrix(xtx,'general', 0, invxtx)
mult_matrix(invxtx,xty,'AB', beta)
*Values[0]:k,Values[1]:b
get_full_matrix(beta, Values)
* 画一条直线
Newy:=Values[0] * [10,800] + Values[1]
gen_contour_polygon_xld(Contour, Newy, [10,800])

运行效果:

线性最小二乘法的halcon实现-萤火

参考链接:
https://www.cnblogs.com/xiaomaLV2/p/6697167.html