Центроискатель

[Jakob Brontfeyn]

Товарищи знатоки, может кто знает, где есть уже готовый
суб VI, для нахождения координат центра окружности и ее радиуса,
если окружность проходит через 3 точки
с координатами М1(Х1,Y1), М2(Х2,Y3), М3(Х3,Y3)
Имеем систему уравнений:

(Х-Х1)^2+(У-У1)^2=R^2
(Х-Х2)^2+(У-У2)^2=R^2
(Х-Х3)^2+(У-У3)^2=R^2

надо получить на выходе Х, У и R

[mzu2006]

Уравнения известны (ближе к концу топика):
Coordinates of a triangle's circumcenter

осталось лишь реализовать несложные вычисления

[AndreyDmitriev]

IMAQ Fit Circle:

http://zone.ni.com/reference/en-XX/help ... _circle_2/

Если, конечно, Vision пакет под рукой имеется.

[IvanLis]

пусть координаты точек x₁,y₁,x₂,y₂,x₃,y₃, координаты центра x,y
тогда (x-x₁)²+(y-y₁)²=(x-x₂)²+(y-y₂)² ⇒ 2x(x₂-x₁)+2y(y₂-y₁)=x₂²-x₁²+y₂²-y₁²
аналогично 2x(x₃-x₁)+2y(y₃-y₁)=x₃²-x₁²+y₃²-y₁²
проще всего не решать эту систему в общем виде а запрограммировать решение линейной системы по коэффициентам
a₁₁=2(x₂-x₁), a₁₂=2(y₂-y₁), b₁=x₂²-x₁²+y₂²-y₁²
a₂₁=2(x₃-x₁), a₂₂=2(y₃-y₁), b₂=x₃²-x₁²+y₃²-y₁²
x=(a₂₂*b₁-a₁₂*b₂)/(a₁₁*a₂₂-a₁₂*a₂₁)
y=(a₁₁*b₂-a₂₁*b₁)/(a₁₁*a₂₂-a₁₂*a₂₁)
радиус можно найти как √((x-x₁)²+(y-y₁)²)

[AndreyDmitriev]

Похоже на правду, за исключением случая, когда точки стоят на одной прямой.

[IvanLis]

Похоже на правду, за исключением случая, когда точки стоят на одной прямой.

случай для окружности исключительный
проверку можно сделать на всякий случай...

[Jakob Brontfeyn]

Большое спасибо за помощь,
важно выиграть время и быть всегда на шаг впереди.
Вот публикую, может кому пригодится