MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径

House

8 A6 [2 h( |" AMATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径

%% 读入离散点的坐标
1 }" T% }1 L) Y8 m& b- {clear;clc;close all
2 K9 I+ O% R: Fdata=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
# t  K( [' V; T0 X' [x=data(:,1);  % 第一列为 x 坐标
y=data(:,2);  % 第二列为 y 坐标
+ Q# Q0 c6 i1 @0 o: gscatter(x,y)  % 绘制离散点的分布图
" c( V4 I2 }$ V2 j2 K, E/ ?8 Gaxis equal    % x 轴和 y 轴等比例输出
+ B# y' b8 v2 _hold on
* q0 e1 p  m# R) D' z2 I( k
%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
B=-(x.^2+y.^2);  % 得到矩阵 B
2 E  k7 V( l( s  g& E( r. ^- jabc=A\B;  % 得到 a b c
' n1 m) J" m1 b9 Da=abc(1);
b=abc(2);
c=abc(3);

; j% E5 [" e2 k! {% 根据a b c 求出圆心坐标和半径
x0=-0.5*a;
y0=-0.5*b;
r=sqrt(x0^2+y0^2-c);

%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
k=1;
3 y1 _, {9 ~3 Y% ?6 F/ I* h' u( Sfor theta=0:pi/180:2*pi      %角度从 0 到 2*pi
+ _$ v0 C% Y+ f    X(k)=r*cos(theta)+x0;    %圆上的横坐标
    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
    k=k+1;
/ y' o7 Q% w( |$ j- c5 y4 }end
plot(X,Y,'r')        %绘制圆
scatter(x0,y0,'r*')  %标示圆心

gaoxings

MATLAB求解拟合圆的圆心和半径。