3D-диаграмма рассеяния, матричные операции и обратный расчет ⇐ Python

[Anonymous]

У меня есть диаграмма рассеяния с кодом, показанным ниже.

Код: Выделить всё

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
x = np.random.rand(100)
y = np.random.rand(100)
z = np.random.rand(100)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
scatter = ax.scatter(x, y, z, c='r', marker='o')
plt.show()

Мы можем сказать

Код: Выделить всё

A = np.column_stack((x, y, z)).

Кроме того, у меня есть матрица B, которая равна матрице A, умноженной на операцию вращения, операцию сдвига + шум.

Код: Выделить всё

B = A x rot(X) + shift(Y) + Noise

Если заданы только матрицы A и B, как мне найти необходимый угол поворота и вектор сдвига, чтобы перекрыть матрицу B с матрицей A/отменить операции над матрицей B? Является ли подход с использованием центра масс лучшим или какой подход рекомендуется использовать?

Подробнее здесь: https://stackoverflow.com/questions/786 ... he-inverse