基于Matlab求解拟合圆的圆心和半径（源码+数据）.rar

%% 读入离散点的坐标 clear;clc;close all data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据 x=data(:,1); % 第一列为 x 坐标 y=data(:,2); % 第二列为 y 坐标 scatter(x,y) % 绘制离散点的分布图 axis equal % x 轴和 y 轴等比例输出 hold on %% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径 A=[x y ones(length(x),1)]; % 得到矩阵 A B=-(x.^2+y.^2); % 得到矩阵 B abc=A\B; % 得到 a b c a=abc(1); b=abc(2); c=abc(3); % 根据a b c 求出圆心坐标和半径 x0=-0.5*a; y0=-0.5*b; r=sqrt(x0^2+y0^2-c); %% 根据计算得到的圆心坐标和半径, 绘制拟合圆 k=1; for theta=0:pi/180:2*pi %角度从 0 到 2*pi X(k)=r*cos(theta)+x0; %圆上的横坐标 Y(k)=r*sin(theta)+y0; %圆上的纵坐标 k=k+1; end plot(X,Y,'r') %绘制圆 scatter(x0,y0,'r*') %标示圆心