最小二乘法拟合圆

Goodness2020

已于 2023-08-01 18:12:42 修改

阅读量713

点赞数 1

分类专栏：拟合文章标签：最小二乘法算法机器学习

于 2021-10-22 16:21:51 首次发布

本文链接：https://blog.csdn.net/Goodness2020/article/details/120908423

版权

拟合专栏收录该内容

4 篇文章 0 订阅

订阅专栏

一、说明

输入：二维点集（X1，Y1），...，（Xn，Yn）

输出：圆的方程 $X^{2}+Y^{2}+AX+BY+C=0$ 中参数A,B,C

解：

$\begin{bmatrix} \sum Xi & \sum Yi & N \\ \sum XiYi & \sum Yi^{2} & \sum Yi \\ \sum Xi^{2} & \sum XiYi & \sum Xi \end{bmatrix} \begin{bmatrix} A \\ B \\ C \end{bmatrix} = \begin{bmatrix} - \sum Xi^{2} - \sum Yi^{2} \\ - \sum Xi^{2}Yi - \sum Yi^{3} \\ - \sum Xi^{3} - \sum XiYi^{2} \end{bmatrix}$

二、代码

% 最小二乘法拟合圆 X*X+Y*Y+AX+BY+C=0  v=[A,B,C]'
clear all;
clc;
X = [1.0, -1.0, 0.0, 0.0, 0.5, -0.5];
Y = [0.0, 0.0, 1.0, -1.0, 0.866, -0.866];
N = length(X);
Xsum = sum(X);
Ysum = sum(Y);
XXsum = dot(X, X);
XYsum = dot(X, Y);
YYsum = dot(Y, Y);
XX = X .* X;
YY = Y .* Y;
XXXsum = dot(XX, X);
XXYsum = dot(XX, Y);
YYYsum = dot(YY, Y);
XYYsum = dot(X, YY);
M = [Xsum, Ysum, N; XYsum, YYsum, Ysum; XXsum, XYsum, Xsum];
u = [-XXsum-YYsum, -XXYsum-YYYsum, -XXXsum-XYYsum]';
v = M^(-1) * u;