Implementación de programación en Matlab de ajuste de círculos por mínimos cuadrados para encontrar el radio
función [R, A, B] = circ (x, y, N) define la función circ
Método de mínimos cuadrados para ajustar curva circular
x1 = 0;
x2 = 0;
x3 = 0;
y1 = 0;
y2 = 0;
y3 = 0;
x1y1 = 0;
x1y2 = 0;
x2y1 = 0;
para i = 1: N
x1 = x1 x(i);
x2 = x2 x(i)*x(i);
x3 = x3 x( i)*x(i)*x(i);
y1 = y1 y(i);
y2 = y2 y(i)*y(i); p>
p>
y3 = y3 y(i)*y(i)*y(i);
x1y1 = x1y1 x(i)*y(i); p>
x1y2 = x1y2 x(i)*y(i)*y(i);
x2y1 = x2y1 x(i)*x(i)*y(i); p>
end
Se puede resolver:
C = N * x2 - x1 * x1;
D = N * x1y1 - x1 * y1 ;
E = N * x3 N * x1y2 - (x2 y2) * x1;
G = N * y2 - y1 * y1;
H = N * x2y1 N * y3 - (x2 y2) * y1;
a = (H * D - E * G)/(C * G - D * D);
b = (H * C - E * D)/(D * D - G * C);
c = -(a * x1 b * y1 x2 y2)/N;
Obtener A , B, valor de ajuste estimado de R
A = a/(-2);
coordenada x
B = b/(- 2);
coordenada y
R = sqrt(a * a b * b - 4 * c)/2;