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

House

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

/ ~+ s" Z* T8 y, e( |; B%% 读入离散点的坐标
! `7 u! d% J1 `4 Fclear;clc;close all
- T2 x( B1 u& K5 ~& Ydata=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
x=data(:,1);  % 第一列为 x 坐标
y=data(:,2);  % 第二列为 y 坐标
5 n8 B! h- z/ w- Q% `scatter(x,y)  % 绘制离散点的分布图
' D' U- ~$ G4 T3 |- eaxis equal    % x 轴和 y 轴等比例输出
hold on

%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
B=-(x.^2+y.^2);  % 得到矩阵 B
abc=A\B;  % 得到 a b c
a=abc(1);
" |$ U4 ^8 K' X7 _7 h( p4 ~) Db=abc(2);
c=abc(3);

% 根据a b c 求出圆心坐标和半径
* [. m5 Z" _  x, Ox0=-0.5*a;
y0=-0.5*b;
r=sqrt(x0^2+y0^2-c);

%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
k=1;
* J1 O" p3 L( n- c) nfor theta=0:pi/180:2*pi      %角度从 0 到 2*pi
    X(k)=r*cos(theta)+x0;    %圆上的横坐标
! |, ~; a0 R$ h7 q8 [" [    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
    k=k+1;
8 }6 A7 P! H9 k; a! oend
plot(X,Y,'r')        %绘制圆
scatter(x0,y0,'r*')  %标示圆心

gaoxings

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