solution of $B = XAX^{-1}$ in terms of $AX = BX$ system of linear equations











up vote
0
down vote

favorite












I am looking for general solution of the form $B = XAX^{-1}$
with the following constraints



$X^{-1}=X'$



where $X$ is unknown matrix and $A,,B,, X$, are $3times 3$ matrices.



The problem appears peculiar and thus I wish to simplify it further to more practical problem
X is rotation matrix ($R$), and $A$ is plane equation matrix, B is resulted homography in the sense of
$B =XAX^{-1} = XAX^{T}= R(I-tn^T/d)R^{T}$ ,



where $t$ is translation vector, $n$ is plane normal vector, and $d$ is signed distance, all three are known



My dirty solution so far is as follows :
Finding the Jordan Canonical form:



$MBM^{-1}=J=NAN^{-1}$,



and then $X$ is obviously



$X = N^{-1}M$



but this works only if the canonical forms are the same, otherwise cannot be utilized for solution.



The second dirty trick is to solve a system of Linear equations
in terms of



$AX = BX$



but i cannot figure out how to program this in the software (e.g. Matlab linsolve).
Could you help me to clarify how to write system of linear equations
in the shape of 9 rows?



The trick is how to notate them on the "right side" to use in ready to use functions for solving SLE.



Thank, you










share|cite|improve this question
























  • well, let's simplify: X is always invertible, orthogonal and X' = X^{1} and is also 3x3 matrix
    – Misha Bolgarskiy
    Nov 14 at 13:32












  • Are $A$ and $B$ symmetric then ?
    – nicomezi
    Nov 14 at 13:32










  • This will not work. Take $A=0$ and $B=I$, the identity.
    – Dietrich Burde
    Nov 14 at 13:33










  • No, A and B are not symmetric. For the symmetric case, there is already a valuable thread. Even I can generalize a solution of the kind MAM' = J = NBN' (finding the jordan canonical form) - if the canonical form is same, then there is solution, if not - no solution, then X = N^{-1}M
    – Misha Bolgarskiy
    Nov 14 at 13:38















up vote
0
down vote

favorite












I am looking for general solution of the form $B = XAX^{-1}$
with the following constraints



$X^{-1}=X'$



where $X$ is unknown matrix and $A,,B,, X$, are $3times 3$ matrices.



The problem appears peculiar and thus I wish to simplify it further to more practical problem
X is rotation matrix ($R$), and $A$ is plane equation matrix, B is resulted homography in the sense of
$B =XAX^{-1} = XAX^{T}= R(I-tn^T/d)R^{T}$ ,



where $t$ is translation vector, $n$ is plane normal vector, and $d$ is signed distance, all three are known



My dirty solution so far is as follows :
Finding the Jordan Canonical form:



$MBM^{-1}=J=NAN^{-1}$,



and then $X$ is obviously



$X = N^{-1}M$



but this works only if the canonical forms are the same, otherwise cannot be utilized for solution.



The second dirty trick is to solve a system of Linear equations
in terms of



$AX = BX$



but i cannot figure out how to program this in the software (e.g. Matlab linsolve).
Could you help me to clarify how to write system of linear equations
in the shape of 9 rows?



The trick is how to notate them on the "right side" to use in ready to use functions for solving SLE.



Thank, you










share|cite|improve this question
























  • well, let's simplify: X is always invertible, orthogonal and X' = X^{1} and is also 3x3 matrix
    – Misha Bolgarskiy
    Nov 14 at 13:32












  • Are $A$ and $B$ symmetric then ?
    – nicomezi
    Nov 14 at 13:32










  • This will not work. Take $A=0$ and $B=I$, the identity.
    – Dietrich Burde
    Nov 14 at 13:33










  • No, A and B are not symmetric. For the symmetric case, there is already a valuable thread. Even I can generalize a solution of the kind MAM' = J = NBN' (finding the jordan canonical form) - if the canonical form is same, then there is solution, if not - no solution, then X = N^{-1}M
    – Misha Bolgarskiy
    Nov 14 at 13:38













up vote
0
down vote

favorite









up vote
0
down vote

favorite











I am looking for general solution of the form $B = XAX^{-1}$
with the following constraints



$X^{-1}=X'$



where $X$ is unknown matrix and $A,,B,, X$, are $3times 3$ matrices.



The problem appears peculiar and thus I wish to simplify it further to more practical problem
X is rotation matrix ($R$), and $A$ is plane equation matrix, B is resulted homography in the sense of
$B =XAX^{-1} = XAX^{T}= R(I-tn^T/d)R^{T}$ ,



where $t$ is translation vector, $n$ is plane normal vector, and $d$ is signed distance, all three are known



My dirty solution so far is as follows :
Finding the Jordan Canonical form:



$MBM^{-1}=J=NAN^{-1}$,



and then $X$ is obviously



$X = N^{-1}M$



but this works only if the canonical forms are the same, otherwise cannot be utilized for solution.



The second dirty trick is to solve a system of Linear equations
in terms of



$AX = BX$



but i cannot figure out how to program this in the software (e.g. Matlab linsolve).
Could you help me to clarify how to write system of linear equations
in the shape of 9 rows?



The trick is how to notate them on the "right side" to use in ready to use functions for solving SLE.



Thank, you










share|cite|improve this question















I am looking for general solution of the form $B = XAX^{-1}$
with the following constraints



$X^{-1}=X'$



where $X$ is unknown matrix and $A,,B,, X$, are $3times 3$ matrices.



The problem appears peculiar and thus I wish to simplify it further to more practical problem
X is rotation matrix ($R$), and $A$ is plane equation matrix, B is resulted homography in the sense of
$B =XAX^{-1} = XAX^{T}= R(I-tn^T/d)R^{T}$ ,



where $t$ is translation vector, $n$ is plane normal vector, and $d$ is signed distance, all three are known



My dirty solution so far is as follows :
Finding the Jordan Canonical form:



$MBM^{-1}=J=NAN^{-1}$,



and then $X$ is obviously



$X = N^{-1}M$



but this works only if the canonical forms are the same, otherwise cannot be utilized for solution.



The second dirty trick is to solve a system of Linear equations
in terms of



$AX = BX$



but i cannot figure out how to program this in the software (e.g. Matlab linsolve).
Could you help me to clarify how to write system of linear equations
in the shape of 9 rows?



The trick is how to notate them on the "right side" to use in ready to use functions for solving SLE.



Thank, you







linear-algebra matrices






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Nov 15 at 1:46

























asked Nov 14 at 13:23









Misha Bolgarskiy

41




41












  • well, let's simplify: X is always invertible, orthogonal and X' = X^{1} and is also 3x3 matrix
    – Misha Bolgarskiy
    Nov 14 at 13:32












  • Are $A$ and $B$ symmetric then ?
    – nicomezi
    Nov 14 at 13:32










  • This will not work. Take $A=0$ and $B=I$, the identity.
    – Dietrich Burde
    Nov 14 at 13:33










  • No, A and B are not symmetric. For the symmetric case, there is already a valuable thread. Even I can generalize a solution of the kind MAM' = J = NBN' (finding the jordan canonical form) - if the canonical form is same, then there is solution, if not - no solution, then X = N^{-1}M
    – Misha Bolgarskiy
    Nov 14 at 13:38


















  • well, let's simplify: X is always invertible, orthogonal and X' = X^{1} and is also 3x3 matrix
    – Misha Bolgarskiy
    Nov 14 at 13:32












  • Are $A$ and $B$ symmetric then ?
    – nicomezi
    Nov 14 at 13:32










  • This will not work. Take $A=0$ and $B=I$, the identity.
    – Dietrich Burde
    Nov 14 at 13:33










  • No, A and B are not symmetric. For the symmetric case, there is already a valuable thread. Even I can generalize a solution of the kind MAM' = J = NBN' (finding the jordan canonical form) - if the canonical form is same, then there is solution, if not - no solution, then X = N^{-1}M
    – Misha Bolgarskiy
    Nov 14 at 13:38
















well, let's simplify: X is always invertible, orthogonal and X' = X^{1} and is also 3x3 matrix
– Misha Bolgarskiy
Nov 14 at 13:32






well, let's simplify: X is always invertible, orthogonal and X' = X^{1} and is also 3x3 matrix
– Misha Bolgarskiy
Nov 14 at 13:32














Are $A$ and $B$ symmetric then ?
– nicomezi
Nov 14 at 13:32




Are $A$ and $B$ symmetric then ?
– nicomezi
Nov 14 at 13:32












This will not work. Take $A=0$ and $B=I$, the identity.
– Dietrich Burde
Nov 14 at 13:33




This will not work. Take $A=0$ and $B=I$, the identity.
– Dietrich Burde
Nov 14 at 13:33












No, A and B are not symmetric. For the symmetric case, there is already a valuable thread. Even I can generalize a solution of the kind MAM' = J = NBN' (finding the jordan canonical form) - if the canonical form is same, then there is solution, if not - no solution, then X = N^{-1}M
– Misha Bolgarskiy
Nov 14 at 13:38




No, A and B are not symmetric. For the symmetric case, there is already a valuable thread. Even I can generalize a solution of the kind MAM' = J = NBN' (finding the jordan canonical form) - if the canonical form is same, then there is solution, if not - no solution, then X = N^{-1}M
– Misha Bolgarskiy
Nov 14 at 13:38










1 Answer
1






active

oldest

votes

















up vote
0
down vote













Case 1. $A,B$ are not orthogonally similar; no solutions.



Case 2. $A,B$ are orthogonally similar; $A=PBP^T$, where $Pin O(n)$. Let $C(A)={Xin M_n(mathbb{R});AX=XA}$. Then, it is not difficult to show



$textbf{Proposition}$ ${Xin O(n);A=XBX^T}=(C(A)cap O(n))P$.



It remains to obtain $C(A)cap O(n)$. When $n=3$, $C(A)$ is a vector space of dimension $3,5$ or $9$ (when $A$ is a scalar matrix).



Generically (randomly choose $A$), $A$ is cyclic



https://en.wikipedia.org/wiki/Cyclic_subspace



and $C(A)={aI+bA+cA^2;a,b,cin mathbb{R}}$. Generically, the eigenspaces of $A$ are not orthogonal and, consequently, $C(A)cap O(n)=pm I_n$ and $X=pm P$.



EDIT. $A,B$ are orthogonally similar iff the couples $(A,A^T)$ and $(B,B^T)$ are similar. Then, practically, one studies the linear system



$(S)$ $AX=XB,A^TX=XB^T$.



If $A,B$ are not orthog. similar, then the sole solution is $X=0$.



If $A,B$ are orthog. similar, then, for a generic $A$, $(S)$ admits a vector space of solutions of dimension $1$: $X=uX_0$. It remains to calculate $u$ so that $||uX_0||_2=1$.






share|cite|improve this answer























  • thank you, I need time to digest this..
    – Misha Bolgarskiy
    Nov 14 at 19:03










  • you are somehow, point towards similarity transformations..?
    – Misha Bolgarskiy
    Nov 14 at 19:04










  • Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
    – Misha Bolgarskiy
    Nov 15 at 1:38










  • @Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
    – loup blanc
    Nov 15 at 11:08










  • I am really sorry Mr. Blanc, i am pretty newbie here
    – Misha Bolgarskiy
    Nov 15 at 17:03











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














 

draft saved


draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2998265%2fsolution-of-b-xax-1-in-terms-of-ax-bx-system-of-linear-equations%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























1 Answer
1






active

oldest

votes








1 Answer
1






active

oldest

votes









active

oldest

votes






active

oldest

votes








up vote
0
down vote













Case 1. $A,B$ are not orthogonally similar; no solutions.



Case 2. $A,B$ are orthogonally similar; $A=PBP^T$, where $Pin O(n)$. Let $C(A)={Xin M_n(mathbb{R});AX=XA}$. Then, it is not difficult to show



$textbf{Proposition}$ ${Xin O(n);A=XBX^T}=(C(A)cap O(n))P$.



It remains to obtain $C(A)cap O(n)$. When $n=3$, $C(A)$ is a vector space of dimension $3,5$ or $9$ (when $A$ is a scalar matrix).



Generically (randomly choose $A$), $A$ is cyclic



https://en.wikipedia.org/wiki/Cyclic_subspace



and $C(A)={aI+bA+cA^2;a,b,cin mathbb{R}}$. Generically, the eigenspaces of $A$ are not orthogonal and, consequently, $C(A)cap O(n)=pm I_n$ and $X=pm P$.



EDIT. $A,B$ are orthogonally similar iff the couples $(A,A^T)$ and $(B,B^T)$ are similar. Then, practically, one studies the linear system



$(S)$ $AX=XB,A^TX=XB^T$.



If $A,B$ are not orthog. similar, then the sole solution is $X=0$.



If $A,B$ are orthog. similar, then, for a generic $A$, $(S)$ admits a vector space of solutions of dimension $1$: $X=uX_0$. It remains to calculate $u$ so that $||uX_0||_2=1$.






share|cite|improve this answer























  • thank you, I need time to digest this..
    – Misha Bolgarskiy
    Nov 14 at 19:03










  • you are somehow, point towards similarity transformations..?
    – Misha Bolgarskiy
    Nov 14 at 19:04










  • Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
    – Misha Bolgarskiy
    Nov 15 at 1:38










  • @Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
    – loup blanc
    Nov 15 at 11:08










  • I am really sorry Mr. Blanc, i am pretty newbie here
    – Misha Bolgarskiy
    Nov 15 at 17:03















up vote
0
down vote













Case 1. $A,B$ are not orthogonally similar; no solutions.



Case 2. $A,B$ are orthogonally similar; $A=PBP^T$, where $Pin O(n)$. Let $C(A)={Xin M_n(mathbb{R});AX=XA}$. Then, it is not difficult to show



$textbf{Proposition}$ ${Xin O(n);A=XBX^T}=(C(A)cap O(n))P$.



It remains to obtain $C(A)cap O(n)$. When $n=3$, $C(A)$ is a vector space of dimension $3,5$ or $9$ (when $A$ is a scalar matrix).



Generically (randomly choose $A$), $A$ is cyclic



https://en.wikipedia.org/wiki/Cyclic_subspace



and $C(A)={aI+bA+cA^2;a,b,cin mathbb{R}}$. Generically, the eigenspaces of $A$ are not orthogonal and, consequently, $C(A)cap O(n)=pm I_n$ and $X=pm P$.



EDIT. $A,B$ are orthogonally similar iff the couples $(A,A^T)$ and $(B,B^T)$ are similar. Then, practically, one studies the linear system



$(S)$ $AX=XB,A^TX=XB^T$.



If $A,B$ are not orthog. similar, then the sole solution is $X=0$.



If $A,B$ are orthog. similar, then, for a generic $A$, $(S)$ admits a vector space of solutions of dimension $1$: $X=uX_0$. It remains to calculate $u$ so that $||uX_0||_2=1$.






share|cite|improve this answer























  • thank you, I need time to digest this..
    – Misha Bolgarskiy
    Nov 14 at 19:03










  • you are somehow, point towards similarity transformations..?
    – Misha Bolgarskiy
    Nov 14 at 19:04










  • Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
    – Misha Bolgarskiy
    Nov 15 at 1:38










  • @Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
    – loup blanc
    Nov 15 at 11:08










  • I am really sorry Mr. Blanc, i am pretty newbie here
    – Misha Bolgarskiy
    Nov 15 at 17:03













up vote
0
down vote










up vote
0
down vote









Case 1. $A,B$ are not orthogonally similar; no solutions.



Case 2. $A,B$ are orthogonally similar; $A=PBP^T$, where $Pin O(n)$. Let $C(A)={Xin M_n(mathbb{R});AX=XA}$. Then, it is not difficult to show



$textbf{Proposition}$ ${Xin O(n);A=XBX^T}=(C(A)cap O(n))P$.



It remains to obtain $C(A)cap O(n)$. When $n=3$, $C(A)$ is a vector space of dimension $3,5$ or $9$ (when $A$ is a scalar matrix).



Generically (randomly choose $A$), $A$ is cyclic



https://en.wikipedia.org/wiki/Cyclic_subspace



and $C(A)={aI+bA+cA^2;a,b,cin mathbb{R}}$. Generically, the eigenspaces of $A$ are not orthogonal and, consequently, $C(A)cap O(n)=pm I_n$ and $X=pm P$.



EDIT. $A,B$ are orthogonally similar iff the couples $(A,A^T)$ and $(B,B^T)$ are similar. Then, practically, one studies the linear system



$(S)$ $AX=XB,A^TX=XB^T$.



If $A,B$ are not orthog. similar, then the sole solution is $X=0$.



If $A,B$ are orthog. similar, then, for a generic $A$, $(S)$ admits a vector space of solutions of dimension $1$: $X=uX_0$. It remains to calculate $u$ so that $||uX_0||_2=1$.






share|cite|improve this answer














Case 1. $A,B$ are not orthogonally similar; no solutions.



Case 2. $A,B$ are orthogonally similar; $A=PBP^T$, where $Pin O(n)$. Let $C(A)={Xin M_n(mathbb{R});AX=XA}$. Then, it is not difficult to show



$textbf{Proposition}$ ${Xin O(n);A=XBX^T}=(C(A)cap O(n))P$.



It remains to obtain $C(A)cap O(n)$. When $n=3$, $C(A)$ is a vector space of dimension $3,5$ or $9$ (when $A$ is a scalar matrix).



Generically (randomly choose $A$), $A$ is cyclic



https://en.wikipedia.org/wiki/Cyclic_subspace



and $C(A)={aI+bA+cA^2;a,b,cin mathbb{R}}$. Generically, the eigenspaces of $A$ are not orthogonal and, consequently, $C(A)cap O(n)=pm I_n$ and $X=pm P$.



EDIT. $A,B$ are orthogonally similar iff the couples $(A,A^T)$ and $(B,B^T)$ are similar. Then, practically, one studies the linear system



$(S)$ $AX=XB,A^TX=XB^T$.



If $A,B$ are not orthog. similar, then the sole solution is $X=0$.



If $A,B$ are orthog. similar, then, for a generic $A$, $(S)$ admits a vector space of solutions of dimension $1$: $X=uX_0$. It remains to calculate $u$ so that $||uX_0||_2=1$.







share|cite|improve this answer














share|cite|improve this answer



share|cite|improve this answer








edited Nov 15 at 0:27

























answered Nov 14 at 15:54









loup blanc

22.1k21749




22.1k21749












  • thank you, I need time to digest this..
    – Misha Bolgarskiy
    Nov 14 at 19:03










  • you are somehow, point towards similarity transformations..?
    – Misha Bolgarskiy
    Nov 14 at 19:04










  • Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
    – Misha Bolgarskiy
    Nov 15 at 1:38










  • @Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
    – loup blanc
    Nov 15 at 11:08










  • I am really sorry Mr. Blanc, i am pretty newbie here
    – Misha Bolgarskiy
    Nov 15 at 17:03


















  • thank you, I need time to digest this..
    – Misha Bolgarskiy
    Nov 14 at 19:03










  • you are somehow, point towards similarity transformations..?
    – Misha Bolgarskiy
    Nov 14 at 19:04










  • Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
    – Misha Bolgarskiy
    Nov 15 at 1:38










  • @Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
    – loup blanc
    Nov 15 at 11:08










  • I am really sorry Mr. Blanc, i am pretty newbie here
    – Misha Bolgarskiy
    Nov 15 at 17:03
















thank you, I need time to digest this..
– Misha Bolgarskiy
Nov 14 at 19:03




thank you, I need time to digest this..
– Misha Bolgarskiy
Nov 14 at 19:03












you are somehow, point towards similarity transformations..?
– Misha Bolgarskiy
Nov 14 at 19:04




you are somehow, point towards similarity transformations..?
– Misha Bolgarskiy
Nov 14 at 19:04












Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
– Misha Bolgarskiy
Nov 15 at 1:38




Maybe I could simplify further the problem, A is plane equation matrix in the shape of $i-pn/h$ and $X$ is a rotation matrix.
– Misha Bolgarskiy
Nov 15 at 1:38












@Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
– loup blanc
Nov 15 at 11:08




@Misha Bolgarskiy , If you want to change the problem, then open a new file. Anyway, I did the job; now it's your business.
– loup blanc
Nov 15 at 11:08












I am really sorry Mr. Blanc, i am pretty newbie here
– Misha Bolgarskiy
Nov 15 at 17:03




I am really sorry Mr. Blanc, i am pretty newbie here
– Misha Bolgarskiy
Nov 15 at 17:03


















 

draft saved


draft discarded



















































 


draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2998265%2fsolution-of-b-xax-1-in-terms-of-ax-bx-system-of-linear-equations%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

How to change which sound is reproduced for terminal bell?

Title Spacing in Bjornstrup Chapter, Removing Chapter Number From Contents

Can I use Tabulator js library in my java Spring + Thymeleaf project?