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

House

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

3 W: @  O; k9 b$ d2 Y%% 读入离散点的坐标
' {# W; y7 @1 Q! q3 N3 ?clear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
* j2 f( `, I/ e. d3 ix=data(:,1);  % 第一列为 x 坐标
y=data(:,2);  % 第二列为 y 坐标
scatter(x,y)  % 绘制离散点的分布图
axis equal    % x 轴和 y 轴等比例输出
7 w7 K( k0 y+ y7 p4 \3 ]3 Q$ Jhold on
0 g" M& P4 c% n& W
%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
A=[x y ones(length(x),1)];  % 得到矩阵 A
B=-(x.^2+y.^2);  % 得到矩阵 B
( ^# R) g- Y. Habc=A\B;  % 得到 a b c
a=abc(1);
b=abc(2);
c=abc(3);
3 I. P% h1 G; a" ~- I
% 根据a b c 求出圆心坐标和半径
5 \1 U' `, S* v- d9 X7 X0 qx0=-0.5*a;
y0=-0.5*b;
r=sqrt(x0^2+y0^2-c);

; b) C$ z5 K6 V4 J0 d%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
k=1;
for theta=0:pi/180:2*pi      %角度从 0 到 2*pi
" G5 S. V0 e0 p9 k8 p5 M. T    X(k)=r*cos(theta)+x0;    %圆上的横坐标
& I* L7 l' U* d  m% D, d0 C( b" v    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
$ ]3 j$ V! i$ _" F2 O' b$ a# l    k=k+1;
' p! O# @) x9 l1 Q7 s5 J9 xend
" d5 _5 K" u3 y8 M: F. zplot(X,Y,'r')        %绘制圆
0 f1 k% K: c, S, Escatter(x0,y0,'r*')  %标示圆心

gaoxings

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