matlab中已知最小二乘拟合圆求孔径精度代码？

以下是MATLAB中最小二乘拟合圆求孔径精度的代码：

% 假设已知一组圆心坐标(x,y)和半径r
% x, y, r为列向量
% 计算拟合圆的参数(a,b,R)
% a, b为圆心坐标，R为半径
% 计算孔径精度

% 拟合圆的参数
n = length(x);
A = [x, y, ones(n,1)];
b = -x.^2-y.^2;
p = A\b;
a = -p(1)/2;
b = -p(2)/2;
R = sqrt(a^2+b^2-p(3));

% 孔径精度
delta_r = sqrt(sum((sqrt(x.^2+y.^2)-R).^2)/(n-3));

其中，拟合圆的参数通过最小二乘法求解，使用矩阵求解线性方程组。孔径精度的计算使用了平均偏差的方法。