最小二乘法求解圆方程圆形及半径

西伯利亚的长颈鹿

已于 2022-07-11 14:16:11 修改

阅读量3k

点赞数 2

分类专栏：算法文章标签：最小二乘法几何学机器学习 matlab

于 2022-06-17 13:50:05 首次发布

本文链接：https://blog.csdn.net/m0_59160272/article/details/125329324

版权

算法专栏收录该内容

1 篇文章 1 订阅

订阅专栏

ci最小二乘法定义（摘抄于百度百科）：

基本思路（摘抄于百度百科）：

简单的来说，最小二乘法为一类线性算法，将需要求解的系数当作未知数，f(x)与x当作已知数，通过多组对应关系求得系数的方法。

所以，最小二乘法仅适合 $f(x)=\alpha 1\varphi 1(x)+\alpha 2\varphi 2(x)+\alpha 3\varphi 3(x)+\alpha 4\varphi 4(x)+......+\alpha m\varphi m(x)$ 系数为一次项方程式

例如： $y=kx+b$ ，k与b作为系数是可通过最小二乘法求的

而： $y=x^{A}$ ，A作为系数，是不可通过最小二乘法求得的

理论过程：

当圆心为（0,0），半径为r时，可得圆的方程式为：

$x^{2}+y^{2}=r^{2}$

令圆心为(h,k)时，

圆方程为： $(x-h)^{2}+(y-k)^{2}=r^{2}$

这是在连续数据的情况下，我们将其代入离散数据点的情况：

圆方程为： $\left ( \sum_{i=1}^{n}X_{n}-h\right )^{2}+\left ( \sum_{i=1}^{n}Y_{n}-k\right )^{2}=r^{2}$

展开式为： $\sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}Y_{i}+h^{2}+k^{2}-r^{2}=0$

为了简化方程式，令 $h^{2}+k^{2}-r^{2}=p$

则方程为： $\sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}Y_{i}+p=0$

既误差函数为： $E=\sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}Y_{i}+p\sum_{i=1}^{n}X_{i}^{0}=0$

对误差函数的二次进行求导，已知E=0，既误差函数的二次方对任何数求导均为0:

所以： $\frac{\partial \left | E \right |^2 }{\partial h}=0;\frac{\partial \left | E \right |^2 }{\partial k}=0;\frac{\partial \left | E \right |^2 }{\partial p}=0$

推得：

$\small 2\left ( \sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}Y_{i}+p\sum_{i=1}^{n}X_{i}^{0} \right )\times (-2\sum_{i=1}^{n}X_{i})=0$ $\small 2\left ( \sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}Y_{i}+p\sum_{i=1}^{n}X_{i}^{0} \right )\times (-2\sum_{i=1}^{n}Y_{i})=0$

$\small 2\left ( \sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}Y_{i}+p\sum_{i=1}^{n}X_{i}^{0} \right )\times\left (\sum_{i=1}^{n}{X_{i}}^{0} \right )=0$

展开并化简得：

$\small \sum_{i=1}^{n}X_{i}^{3}+\sum_{i=1}^{n}X_{i}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}^{2}-2k\sum_{i=1}^{n}X_{i}Y_{i}+p\sum_{i=1}^{n}X_{i}=0$

$\small \sum_{i=1}^{n}X_{i}^{2}Y_{i}+\sum_{i=1}^{n}Y_{i}^{3}-2h\sum_{i=1}^{n}X_{i}Y_{i}-2k\sum_{i=1}^{n}{Y_{i}}^{2}+p\sum_{i=1}^{n}Y_{i}=0$

$\small \sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}-2h\sum_{i=1}^{n}X_{i}-2k\sum_{i=1}^{n}{Y_{i}}+p\sum_{i=1}^{n}Y_{i}=0$

系数放等号左边，常数放等号右边

$\small 2h\sum_{i=1}^{n}X_{i}^{2}+2k\sum_{i=1}^{n}X_{i}Y_{i}-p\sum_{i=1}^{n}X_{i}=\sum_{i=1}^{n}X_{i}^{3}+\sum_{i=1}^{n}X_{i}Y_{i}^{2}$

$\small 2h\sum_{i=1}^{n}X_{i}Y_{i}+2k\sum_{i=1}^{n}{Y_{i}}^{2}-p\sum_{i=1}^{n}Y_{i}=\sum_{i=1}^{n}X_{i}^{2}Y_{i}+\sum_{i=1}^{n}Y_{i}^{3}$

$\small 2h\sum_{i=1}^{n}X_{i}+2k\sum_{i=1}^{n}{Y_{i}}-p\sum_{i=1}^{n}Y_{i}=\sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}$

使用矩阵表示既为

$\small \begin{bmatrix} C11&C12&C13 \\ C21&C21&C23 \\ \ C31&C32&C33 \end{bmatrix}\times\begin{bmatrix} h\\ k\\ p \end{bmatrix}=\begin{bmatrix} R1\\ R2\\R3 \end{bmatrix}$

其中

$\small C11=\sum_{i=1}^{n}X_{i}^{2}$ 、 $\small C12=\sum_{i=1}^{n}X_{i}Y_{i}$ 、 $\small C13=\sum_{i=1}^{n}X_{i}$

$\small C21=\sum_{i=1}^{n}X_{i}Y_{i}$ 、 $\small C22\sum_{i=1}^{n}{Y_{i}}^{2}$ 、 $\small C23\sum_{i=1}^{n}Y_{i}$

$\small C31=\sum_{i=1}^{n}X_{i}$ 、 $\small C32=\sum_{i=1}^{n}{Y_{i}}$ 、 $\small C33=\sum_{i=1}^{n}Y_{i}$

$\small R1=\sum_{i=1}^{n}X_{i}^{3}+\sum_{i=1}^{n}X_{i}Y_{i}^{2}$

$\small R2=\sum_{i=1}^{n}X_{i}^{2}Y_{i}+\sum_{i=1}^{n}Y_{i}^{3}$

$\small R3=\sum_{i=1}^{n}X_{i}^{2}+\sum_{i=1}^{n}Y_{i}^{2}$

求解矩阵方程既可得到圆心(h,k)

$r=\sqrt{h^{2}+k^{2}-p}$

可以利用克莱姆求解矩阵方程

$\Delta=\begin{vmatrix} C11&C12&C13 \\ C21&C21&C23 \\ C31&C32&C33 \end{vmatrix}$

$\Delta h=\begin{vmatrix} R1&C12&C13 \\ R2&C21&C23 \\ R3&C32&C33 \end{vmatrix}$

$\Delta k=\begin{vmatrix} C11&R2&C13 \\ C21&R2&C23 \\ C31&R3&C33 \end{vmatrix}$

$\Delta p=\begin{vmatrix} C11&C12&R1\\ C21&C21&R2\\ C31&C32&R3\end{vmatrix}$

解得：

$h=\frac{\Delta h}{\Delta}$ 、 $k=\frac{\Delta k}{\Delta}$ 、 $p=\frac{\Delta p}{\Delta}$

此方法和推导同样适合多项式方程、椭圆方程等，这边不再做详细介绍

MATLAB代码（为MATLAB R2022版本下.mlapp回调代码）：

            [file,path] = uigetfile;
            if isempty(file)
                return;
            end 
            if ~contains(file,'.csv')
                msgbox("请确定选择文件为<.csv>文件后重试！")
                return;
            end
            filepath=strcat(path,file);
            LineData=readtable(filepath);
            LineAll=table2array(LineData);
            LineX=LineAll(:,1);
            Liney=LineAll(:,2);

            nSum=0;
            XSum=0;
            YSum=0;
            XYSum=0;
            X2Sum=0;            
            Y2Sum=0;

            X2YSum=0;
            XY2Sum=0;
            X3Sum=0;
            Y3Sum=0;


            N=size(LineX);
          
            for i=1:N
                nSum=nSum+1;
                XSum=XSum+LineX(i);
                YSum=YSum+Liney(i);
                XYSum=XYSum+LineX(i)*Liney(i);
                X2Sum=X2Sum+LineX(i)*LineX(i);
                Y2Sum=Y2Sum+Liney(i)*Liney(i);

                X2YSum=X2YSum+LineX(i)*LineX(i)*Liney(i);
                XY2Sum=XY2Sum+LineX(i)*Liney(i)*Liney(i);
                X3Sum=X3Sum+LineX(i)*LineX(i)*LineX(i);
                Y3Sum=Y3Sum+Liney(i)*Liney(i)*Liney(i);
            end 
            
            C11=2*X2Sum;
            C12=2*XYSum;
            C13=-1*XSum;
            C21=2*XYSum;
            C22=2*Y2Sum;
            C23=-1*YSum;
            C31=2*XSum;
            C32=2*YSum;
            C33=-1*nSum;

            R1=X3Sum+XY2Sum;
            R2=X2YSum+Y3Sum;
            R3=X2Sum+Y2Sum;
            
            M1=[C11,C12,C13;
                C21,C22,C23;
                C31,C32,C33;];

            M2=[R1;R2;R3];

            M3=M1\M2;

            A=M3(1);
            B=M3(2);
            R=sqrt((A^2+B^2-M3(3)));
            Resultstring=strcat("圆心(",num2str(A),",",num2str(B),")","半径：",num2str(R));
            msgbox(Resultstring);

此文章为个人知识点总结文章，转载请标注来源。

西伯利亚的长颈鹿

关注

2
点赞
踩
12

收藏

觉得还不错? 一键收藏
2
评论
最小二乘法求解圆方程圆形及半径

此文章为最小二乘法拟合圆方程，并求得圆心坐标与半径的算法推导，内附MATLAB代码。简单的来说，最小二乘法为一类线性算法，将需要求解的系数当作未知数，f(x)与x当作已知数，通过多组对应关系求得系数的方法。所以，最小二乘法仅适合系数为一次项方程式例如：，k与b作为系数是可通过最小二乘法求的而：，A作为系数，是不可通过最小二乘法求得的理论过程：当圆心为（0,0），半径为r时，可得圆的方程式为：令圆心为(h,k)时.........
复制链接

扫一扫