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

House

# t+ u! t$ f9 Y. N1 |MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径
& V3 j5 n5 _0 T  {$ _- J
%% 读入离散点的坐标
clear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
x=data(:,1);  % 第一列为 x 坐标
$ i1 w( ?9 J" |3 S9 q4 m8 {y=data(:,2);  % 第二列为 y 坐标
! J- P1 D5 [4 v+ v8 kscatter(x,y)  % 绘制离散点的分布图
1 m. p8 ?+ e, E5 R8 eaxis equal    % x 轴和 y 轴等比例输出
. r* W0 n" {& k$ V, Hhold on

# X1 \: I9 A' U/ p, p+ K( q% Y%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
" N6 ]3 H" i8 ^$ T$ WB=-(x.^2+y.^2);  % 得到矩阵 B
abc=A\B;  % 得到 a b c
4 W' x2 W4 B( L6 ~a=abc(1);
b=abc(2);
$ W) U: _+ C) c: i0 b9 qc=abc(3);

$ ^$ }( l$ B" c% R  Z/ L% 根据a b c 求出圆心坐标和半径
+ F1 p( z% x$ _& s5 Ax0=-0.5*a;
& B. t8 W( z+ G+ B* E; a5 z4 W' Ly0=-0.5*b;
+ q! L$ @$ @5 w# I- Qr=sqrt(x0^2+y0^2-c);

%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
k=1;
6 y) X+ j  e& \6 ~8 H$ Ufor theta=0:pi/180:2*pi      %角度从 0 到 2*pi
    X(k)=r*cos(theta)+x0;    %圆上的横坐标
& q& L) u  B& c9 ?    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
, n3 V$ M( }! b- a+ u5 l    k=k+1;
/ H8 a7 H2 ^' r2 X- @- y: gend
6 g0 \7 U2 b9 ?" Y  Q9 }plot(X,Y,'r')        %绘制圆
scatter(x0,y0,'r*')  %标示圆心

gaoxings

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