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

House

MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径
! b8 t$ U" t1 |! L" _9 r3 `& ~
7 R! ~3 x$ N3 e%% 读入离散点的坐标
; s6 _, P( E4 P0 w; j: F: c- ]clear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
x=data(:,1);  % 第一列为 x 坐标
y=data(:,2);  % 第二列为 y 坐标
scatter(x,y)  % 绘制离散点的分布图
axis equal    % x 轴和 y 轴等比例输出
1 ]: Z" x3 ]! c1 Zhold on

%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
' ?5 D  t. M8 p+ l$ K" BA=[x y ones(length(x),1)];  % 得到矩阵 A
$ ~8 B0 ~' v( }! o  e1 _. b5 x3 mB=-(x.^2+y.^2);  % 得到矩阵 B
abc=A\B;  % 得到 a b c
! C/ ^# n  R3 w4 m' J1 Ua=abc(1);
' Z- z% k2 P& [) i" ub=abc(2);
9 }& e) g* _' p* b- \c=abc(3);

; N" c& j% ^! a5 c! e& ~, d% 根据a b c 求出圆心坐标和半径
2 I+ r! H2 M6 f8 n- Bx0=-0.5*a;
% }5 O% A1 ~( u* k* v7 D& yy0=-0.5*b;
; V: r. \  C0 ar=sqrt(x0^2+y0^2-c);

0 p2 r: r4 e* Q5 x%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
; n: B2 l3 F# {6 Wk=1;
for theta=0:pi/180:2*pi      %角度从 0 到 2*pi
. \! ~6 W% P  L/ k3 I# s    X(k)=r*cos(theta)+x0;    %圆上的横坐标
    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
    k=k+1;
end
& i6 t$ `; {" |2 A; u8 Xplot(X,Y,'r')        %绘制圆
scatter(x0,y0,'r*')  %标示圆心

gaoxings

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