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

House

4 H0 @7 P7 N1 @: i  L& Z/ _MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径
# R; B+ y5 l: `% H
%% 读入离散点的坐标
, h+ ?8 s/ C* Dclear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
7 a" L* U+ k! |. Dx=data(:,1);  % 第一列为 x 坐标
1 A9 |5 f9 c% {" W; }y=data(:,2);  % 第二列为 y 坐标
scatter(x,y)  % 绘制离散点的分布图
' t/ \# T9 o/ \/ C: U7 `axis equal    % x 轴和 y 轴等比例输出
( `8 f2 b) H! U' ]' ~8 _6 ihold on
- ~5 A8 w  d4 k' N. @9 D- p8 q8 F
4 L9 F$ @% U& H) s' [- ~# P( w%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
B=-(x.^2+y.^2);  % 得到矩阵 B
abc=A\B;  % 得到 a b c
  r1 R2 A, w& ^, o% I4 ca=abc(1);
6 w) Y# l, [: K4 P5 U+ }: x$ Pb=abc(2);
c=abc(3);

! _0 _) W1 x3 J) Q2 [1 |% 根据a b c 求出圆心坐标和半径
x0=-0.5*a;
' v. V# O) Q* a" q; {# Ty0=-0.5*b;
r=sqrt(x0^2+y0^2-c);

%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
8 b; A1 P$ L! M" sk=1;
& X! \4 ^* P  Ofor theta=0:pi/180:2*pi      %角度从 0 到 2*pi
/ z' z) r# ]5 j    X(k)=r*cos(theta)+x0;    %圆上的横坐标
    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
3 p1 W4 L# h/ `: t" c    k=k+1;
end
, W4 |. p7 ^; \( e3 z  w2 @" x0 [9 I0 @plot(X,Y,'r')        %绘制圆
8 E) ?# ?0 R$ q8 ?scatter(x0,y0,'r*')  %标示圆心

gaoxings

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