デカルト空間に点の集合があります。Python と numpy を使用して、特定の 4 点 (a、b、c、d) のサブアンサンブルによって定義される二面角を計算できます。以下は私の機能です:
def getDihedral(a,b,c,d):
v1 = getNormedVector(a, b)
v2 = getNormedVector(b, c)
v3 = getNormedVector(c, d)
v1v2 = numpy.cross(v1,v2)
v2v3 = numpy.cross(v2,v3)
return getAngle(v1v2,v2v3)
def getNormedVector(a,b):
return (b-a)/numpy.linalg.norm(b-a)
def getAngle(a,b):
return numpy.rad2deg(numpy.arccos(numpy.dot(a/numpy.linalg.norm(a),b.T/numpy.linalg.norm(b))))[0,0]
二面角を 1 つだけ回転させたいのですが、python と numpy および scipy を使用して、ポイントのサブアンサンブルの新しい座標を計算するにはどうすればよいですか?