Надеюсь правильно сформулировал название темы.
по мотивам следущих тем(оказалось, что нужно решить не матричное уравнение, но смысл тот же)
topic61117.htmlнеобходимо решить систему уравнений
(в форме "операторов" (не уверен в строгой математической записи, но в целом должно быть понятно))
![$M[Rot[(x,y)]]=L[M[(x,y)]]$ $M[Rot[(x,y)]]=L[M[(x,y)]]$](https://dxdy-01.korotkov.co.uk/f/c/d/5/cd5b20c1f4f398cfa11ad76c6274102a82.png)
где каждый "оператор" (L,M,Rot)представляет собой перспективное искажение
по формуле(если быть точным то для Rot формула проще
![$m_{31}=0,m_{32}=0,m_{33}=1$ $m_{31}=0,m_{32}=0,m_{33}=1$](https://dxdy-01.korotkov.co.uk/f/c/f/e/cfef2ba00e16ccbe60e70bff792c3e5182.png)
и вообще для всех справедливо
![$m_{33}=1$ $m_{33}=1$](https://dxdy-04.korotkov.co.uk/f/b/1/2/b12220b1c4260099f7b95635ca2c699582.png)
)
![Изображение](http://dxdy-02.korotkov.co.uk/f/5/e/8/5e8e8984a34d8aed41e5b5020de5adf082.png)
попробовал решить задачу в пакете математика, причём подставляя численные значение
![$a_{i,j}$ $a_{i,j}$](https://dxdy-03.korotkov.co.uk/f/2/7/1/271c219e1b5625e22b2b455d15fc50b582.png)
,
![$b_{i,j}$ $b_{i,j}$](https://dxdy-01.korotkov.co.uk/f/0/2/1/02159efc56c78f7b61a81e70d35609f482.png)
, математика долго думает и не видно ни конца не края, хотелось бы хотя бы понять правильно ли я всё делаю или всё просто зависло.
в коде
![$a_{i,j}$ $a_{i,j}$](https://dxdy-03.korotkov.co.uk/f/2/7/1/271c219e1b5625e22b2b455d15fc50b582.png)
для матрицы Rot,
![$b_{i,j}$ $b_{i,j}$](https://dxdy-01.korotkov.co.uk/f/0/2/1/02159efc56c78f7b61a81e70d35609f482.png)
для матрицы L
параметры
![$b_{i,j}$ $b_{i,j}$](https://dxdy-01.korotkov.co.uk/f/0/2/1/02159efc56c78f7b61a81e70d35609f482.png)
были определены для задачи примерно, т.е. есть некоторая погрешность, не знаю как это может повлиять на решение.
получаю 8 уравнений 8 неизвестных для 4 точек(точки могут быть любыми или лучше всё таки взять некоторые определенные\более простые?)
код для математики:
(в коде
![$a_{i,j}$ $a_{i,j}$](https://dxdy-03.korotkov.co.uk/f/2/7/1/271c219e1b5625e22b2b455d15fc50b582.png)
для матрицы Rot,
![$b_{i,j}$ $b_{i,j}$](https://dxdy-01.korotkov.co.uk/f/0/2/1/02159efc56c78f7b61a81e70d35609f482.png)
для матрицы L)
Цитата:
Solve[(b13+(b11 (m13+m11 x1+m12 y1))/(m33+m31 x1+m32 y1)+(b12 (m23+m21 x1+m22 y1))/(m33+m31 x1+m32 y1))/(b33+(b31 (m13+m11 x1+m12 y1))/(m33+m31 x1+m32 y1)+(b32 (m23+m21 x1+m22 y1))/(m33+m31 x1+m32 y1))==
(m13+(m11 (a13+a11 x1+a12 y1))/(a33+a31 x1+a32 y1)+(m12 (a23+a21 x1+a22 y1))/(a33+a31 x1+a32 y1))/(m33+(m31 (a13+a11 x1+a12 y1))/(a33+a31 x1+a32 y1)+(m32 (a23+a21 x1+a22 y1))/(a33+a31 x1+a32 y1))&&
(b23+(b21 (m13+m11 x1+m12 y1))/(m33+m31 x1+m32 y1)+(b22 (m23+m21 x1+m22 y1))/(m33+m31 x1+m32 y1))/(b33+(b31 (m13+m11 x1+m12 y1))/(m33+m31 x1+m32 y1)+(b32 (m23+m21 x1+m22 y1))/(m33+m31 x1+m32 y1))==
(m23+(m21 (a13+a11 x1+a12 y1))/(a33+a31 x1+a32 y1)+(m22 (a23+a21 x1+a22 y1))/(a33+a31 x1+a32 y1))/(m33+(m31 (a13+a11 x1+a12 y1))/(a33+a31 x1+a32 y1)+(m32 (a23+a21 x1+a22 y1))/(a33+a31 x1+a32 y1))&&
(b13+(b11 (m13+m11 x2+m12 y2))/(m33+m31 x2+m32 y2)+(b12 (m23+m21 x2+m22 y2))/(m33+m31 x2+m32 y2))/(b33+(b31 (m13+m11 x2+m12 y2))/(m33+m31 x2+m32 y2)+(b32 (m23+m21 x2+m22 y2))/(m33+m31 x2+m32 y2))==
(m13+(m11 (a13+a11 x2+a12 y2))/(a33+a31 x2+a32 y2)+(m12 (a23+a21 x2+a22 y2))/(a33+a31 x2+a32 y2))/(m33+(m31 (a13+a11 x2+a12 y2))/(a33+a31 x2+a32 y2)+(m32 (a23+a21 x2+a22 y2))/(a33+a31 x2+a32 y2))&&
(b23+(b21 (m13+m11 x2+m12 y2))/(m33+m31 x2+m32 y2)+(b22 (m23+m21 x2+m22 y2))/(m33+m31 x2+m32 y2))/(b33+(b31 (m13+m11 x2+m12 y2))/(m33+m31 x2+m32 y2)+(b32 (m23+m21 x2+m22 y2))/(m33+m31 x2+m32 y2))==
(m23+(m21 (a13+a11 x2+a12 y2))/(a33+a31 x2+a32 y2)+(m22 (a23+a21 x2+a22 y2))/(a33+a31 x2+a32 y2))/(m33+(m31 (a13+a11 x2+a12 y2))/(a33+a31 x2+a32 y2)+(m32 (a23+a21 x2+a22 y2))/(a33+a31 x2+a32 y2))&&
(b13+(b11 (m13+m11 x3+m12 y3))/(m33+m31 x3+m32 y3)+(b12 (m23+m21 x3+m22 y3))/(m33+m31 x3+m32 y3))/(b33+(b31 (m13+m11 x3+m12 y3))/(m33+m31 x3+m32 y3)+(b32 (m23+m21 x3+m22 y3))/(m33+m31 x3+m32 y3))==
(m13+(m11 (a13+a11 x3+a12 y3))/(a33+a31 x3+a32 y3)+(m12 (a23+a21 x3+a22 y3))/(a33+a31 x3+a32 y3))/(m33+(m31 (a13+a11 x3+a12 y3))/(a33+a31 x3+a32 y3)+(m32 (a23+a21 x3+a22 y3))/(a33+a31 x3+a32 y3))&&
(b23+(b21 (m13+m11 x3+m12 y3))/(m33+m31 x3+m32 y3)+(b22 (m23+m21 x3+m22 y3))/(m33+m31 x3+m32 y3))/(b33+(b31 (m13+m11 x3+m12 y3))/(m33+m31 x3+m32 y3)+(b32 (m23+m21 x3+m22 y3))/(m33+m31 x3+m32 y3))==
(m23+(m21 (a13+a11 x3+a12 y3))/(a33+a31 x3+a32 y3)+(m22 (a23+a21 x3+a22 y3))/(a33+a31 x3+a32 y3))/(m33+(m31 (a13+a11 x3+a12 y3))/(a33+a31 x3+a32 y3)+(m32 (a23+a21 x3+a22 y3))/(a33+a31 x3+a32 y3))&&
(b13+(b11 (m13+m11 x4+m12 y4))/(m33+m31 x4+m32 y4)+(b12 (m23+m21 x4+m22 y4))/(m33+m31 x4+m32 y4))/(b33+(b31 (m13+m11 x4+m12 y4))/(m33+m31 x4+m32 y4)+(b32 (m23+m21 x4+m22 y4))/(m33+m31 x4+m32 y4))==
(m13+(m11 (a13+a11 x4+a12 y4))/(a33+a31 x4+a32 y4)+(m12 (a23+a21 x4+a22 y4))/(a33+a31 x4+a32 y4))/(m33+(m31 (a13+a11 x4+a12 y4))/(a33+a31 x4+a32 y4)+(m32 (a23+a21 x4+a22 y4))/(a33+a31 x4+a32 y4))&&
(b23+(b21 (m13+m11 x4+m12 y4))/(m33+m31 x4+m32 y4)+(b22 (m23+m21 x4+m22 y4))/(m33+m31 x4+m32 y4))/(b33+(b31 (m13+m11 x4+m12 y4))/(m33+m31 x4+m32 y4)+(b32 (m23+m21 x4+m22 y4))/(m33+m31 x4+m32 y4))==
(m23+(m21 (a13+a11 x4+a12 y4))/(a33+a31 x4+a32 y4)+(m22 (a23+a21 x4+a22 y4))/(a33+a31 x4+a32 y4))/(m33+(m31 (a13+a11 x4+a12 y4))/(a33+a31 x4+a32 y4)+(m32 (a23+a21 x4+a22 y4))/(a33+a31 x4+a32 y4))&&m33==1,{m11,m12,m13,m21,m22,m23,m31,m32}]
собственно вопрос как можно решить, хотелось бы услышать названия методов и как эт оможно потом вставить в код на с++.