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

House

MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径
$ \& W+ ^3 ?, P6 [  m0 s  N
%% 读入离散点的坐标
clear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
x=data(:,1);  % 第一列为 x 坐标
y=data(:,2);  % 第二列为 y 坐标
/ ?$ c7 T; |, Vscatter(x,y)  % 绘制离散点的分布图
axis equal    % x 轴和 y 轴等比例输出
hold on

4 M$ m0 {" i5 o: m! z" @1 U" H%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
2 i) x( g) H/ G. R2 a0 V! G. y$ cB=-(x.^2+y.^2);  % 得到矩阵 B
6 }% O1 e  a7 ^0 U0 B8 d& W) ~" `& {abc=A\B;  % 得到 a b c
a=abc(1);
8 O8 A5 X3 T+ K/ B/ f8 qb=abc(2);
: ?* s, ~- w* n) \# b6 bc=abc(3);
! O- m* ]) e) @. [7 n; [1 w2 z# p
% 根据a b c 求出圆心坐标和半径
1 e' S. }% X) r1 z  ox0=-0.5*a;
! d. ^/ S9 z& W; j+ Fy0=-0.5*b;
4 W4 p( s: s0 g& ?r=sqrt(x0^2+y0^2-c);

%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
+ I  y' y/ I* V4 ]  O4 F+ Pk=1;
for theta=0:pi/180:2*pi      %角度从 0 到 2*pi
    X(k)=r*cos(theta)+x0;    %圆上的横坐标
2 O7 m, _/ F  o+ W1 v  b* T6 N    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
8 A$ s: J* {3 t0 U2 s8 i2 u2 k. N    k=k+1;
- u8 {* j* E  w8 s/ S: uend
9 p, |' g8 E. N7 r4 ^* Kplot(X,Y,'r')        %绘制圆
. R6 K8 N2 E  Y" Bscatter(x0,y0,'r*')  %标示圆心

gaoxings

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