24/04/30 08:47:53.56 VcpWQbIP.net
>>696
それを前提にして計算
変数が一つにできればNewton-Raphsonが使えるので
横軸にs,縦軸に(AB-BC)^2+(AB-CD)^2+(AB-DA)^2+(BC-CD)^2+(BC-DA)^2+(CD-DA)^2+(AC-BD)^2 をおいて
グラフ化
URLリンク(i.imgur.com)
最小値をとるsは1つだけのようなのでこれを
Newton-Raphson法(R言語ではuniroot関数)でもとめると
> opt=optimize(f,c(0,1),tol=1e-16) ; opt
$minimum
[1] 0.5921734
$objective
[1] 7.888609e-31
で 二変数でのNelder-Meadと同じ結果。
言語仕様や関数を検索しながらWolframに移植するのが次の課題。