EDA365电子论坛网

标题: MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径 [打印本页]

---

作者: House    时间: 2020-3-18 09:41
标题: MATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径

( Z# ~& a, a$ CMATLAB源程序代码分享：MATLAB求解拟合圆的圆心和半径
- I7 u8 M1 `3 B4 d( E5 a
%% 读入离散点的坐标
- Z& k% W; O$ h( c. u/ Y" gclear;clc;close all
data=xlsread('xy.xlsx'); % 读入 excel 文件中的数据
( Y4 @2 q- v; Q8 dx=data(:,1);  % 第一列为 x 坐标
9 j4 j# _- _6 n: ny=data(:,2);  % 第二列为 y 坐标
scatter(x,y)  % 绘制离散点的分布图
( s6 H8 d1 b2 g- w$ haxis equal    % x 轴和 y 轴等比例输出
7 a$ G+ ^- C2 b4 S; Y2 jhold on
- E$ y7 R+ A+ N' F* R# I
* d0 m" @8 M8 }7 q4 z9 k9 h%% 使用线性拟合的方法, 计算拟合圆的圆心坐标和半径
7 w4 Y" n6 E' ?A=[x y ones(length(x),1)];  % 得到矩阵 A
B=-(x.^2+y.^2);  % 得到矩阵 B
& p) _3 ^& k! U. \abc=A\B;  % 得到 a b c
a=abc(1);
, |: x# u5 b5 m! J8 ab=abc(2);
c=abc(3);
, f1 b! d6 ?7 z2 e5 v& a
! [6 {& i$ i# _% V% 根据a b c 求出圆心坐标和半径
x0=-0.5*a;
y0=-0.5*b;
r=sqrt(x0^2+y0^2-c);

5 ~6 ~6 F4 X; y%% 根据计算得到的圆心坐标和半径, 绘制拟合圆
k=1;
6 B" g( }( _2 J2 Hfor theta=0:pi/180:2*pi      %角度从 0 到 2*pi
3 a$ S& d. K; m/ y6 b; q    X(k)=r*cos(theta)+x0;    %圆上的横坐标
    Y(k)=r*sin(theta)+y0;    %圆上的纵坐标
    k=k+1;
, j' j$ s  I/ k6 e% c0 ?end
8 R& z  s  h7 y- D( Kplot(X,Y,'r')        %绘制圆
scatter(x0,y0,'r*')  %标示圆心

---

作者: gaoxings    时间: 2020-3-19 18:21
MATLAB求解拟合圆的圆心和半径。

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2