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

House

MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径
/ n; {% G- w) q- B
%% 读入离散点的坐标
clear;clc;close all
4 m& v, C. E: b$ \9 T- g$ kdata=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
: U) O" Y* O( G) V; Bx=data(:,1);  % 第一列为 x 坐标
y=data(:,2);  % 第二列为 y 坐标
! j- b( a- y5 L' dscatter(x,y)  % 绘制离散点的分布图
5 a" x1 I2 x7 D7 G' [axis equal    % x 轴和 y 轴等比例输出
! b+ Z$ D9 P* Q. x" v4 Uhold on
9 [* ~+ J7 Z# i9 e3 ?- K" [
  }+ v% y! |' Q3 `/ z% |, A%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
B=-(x.^2+y.^2);  % 得到矩阵 B
8 q7 ^) L, v( O/ |abc=A\B;  % 得到 a b c
a=abc(1);
b=abc(2);
1 `5 O: T" u( s- ]- Kc=abc(3);
6 n4 e0 u- S3 C+ }* I8 c' y
% 根据a b c 求出圆心坐标和半径
x0=-0.5*a;
/ a" Y$ `5 e! fy0=-0.5*b;
9 b4 g% E9 i' E6 y, m( V* Wr=sqrt(x0^2+y0^2-c);
0 g( X- O1 a# M
%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
k=1;
; e% m1 V6 B) N3 |& C& Ufor theta=0:pi/180:2*pi      %角度从 0 到 2*pi
1 C9 o: ?& N; u/ P/ G: T    X(k)=r*cos(theta)+x0;    %圆上的横坐标
    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
; P5 ~8 F  `  E8 u  K1 e; q' X    k=k+1;
$ c+ b5 g9 ~# @) A6 r& c8 qend
plot(X,Y,'r')        %绘制圆
scatter(x0,y0,'r*')  %标示圆心

gaoxings

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