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

House

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

0 K5 N: B  k8 K5 J( l3 P%% 读入离散点的坐标
clear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
x=data(:,1);  % 第一列为 x 坐标
/ I4 j0 p; \7 m6 U5 ty=data(:,2);  % 第二列为 y 坐标
scatter(x,y)  % 绘制离散点的分布图
; P; y* d$ J5 r3 Maxis equal    % x 轴和 y 轴等比例输出
hold on
/ [7 a' {* B3 R$ D
  \% |6 P7 f, R" }1 d$ e# e%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
5 H" D1 X: h( ]% |0 i1 d  p+ cB=-(x.^2+y.^2);  % 得到矩阵 B
abc=A\B;  % 得到 a b c
, O7 Y- T% p; u+ L0 `a=abc(1);
0 t0 [, u) W( X0 }2 S8 Rb=abc(2);
c=abc(3);

. e$ D- G9 @! l1 k% 根据a b c 求出圆心坐标和半径
2 j. C& e4 g8 p) c/ lx0=-0.5*a;
y0=-0.5*b;
r=sqrt(x0^2+y0^2-c);
' O% D" U) ]: x; G  G
. ^3 I3 ?" ?$ z0 c9 N. A$ w) c%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
6 W) U$ Q, u/ n) @( p" dk=1;
for theta=0:pi/180:2*pi      %角度从 0 到 2*pi
" |9 n% m4 Z& d/ W: L  m: V+ H    X(k)=r*cos(theta)+x0;    %圆上的横坐标
( B3 d# z" Q& `& X5 l" G6 ?( x8 y    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
    k=k+1;
end
plot(X,Y,'r')        %绘制圆
) I& S- R* o  z9 d& E# |7 Zscatter(x0,y0,'r*')  %标示圆心

gaoxings

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