Solving for Euler Angles












3












$begingroup$


I would like to determine Euler angles according to the following example.



My example:



I have the three vectors in an original set of axes:



r1e = {-0.517853, 0., -0.759239}
r2e = {-0.517853, 0., 0.759239}
r3e = {0.0647316, 0., 0.}


And after expressing them in a new reference frame they obtain the following components:



rt1e={0.310733, -0.358839, -0.786917}
rt2e={0.690333, 0.298661, 0.527983}
rt3e={-0.0625667, 0.00376111, 0.0161833}


In reality, I know the Euler angles to be $(30,60,120)$ degrees. How can I get Mathetmatica to give me this?










share|improve this question











$endgroup$












  • $begingroup$
    Have you tried EulerAngles?
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:09










  • $begingroup$
    That command won't work, as the rotation matrix itself is not known. Only the two vectors.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:10






  • 1




    $begingroup$
    In general, a rotation matrix is not uniquely defined by the action on a single vector...
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:12










  • $begingroup$
    Fair enough, but in this case the problem has been constructed such that the three Euler angles are known to exist. The question specifically relates to why NSolve is not working.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:14






  • 1




    $begingroup$
    What I tried to say: If you prescribe a pair $u$ and $v$ of same length $neq 0$ in $mathbb{R}^3$, then there will be a one-parameter family of rotations (and thus a one-parameter family of Euler angles) that map $u$ to $v$. So, it is as it is: Your problem is underdetermined.
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:18


















3












$begingroup$


I would like to determine Euler angles according to the following example.



My example:



I have the three vectors in an original set of axes:



r1e = {-0.517853, 0., -0.759239}
r2e = {-0.517853, 0., 0.759239}
r3e = {0.0647316, 0., 0.}


And after expressing them in a new reference frame they obtain the following components:



rt1e={0.310733, -0.358839, -0.786917}
rt2e={0.690333, 0.298661, 0.527983}
rt3e={-0.0625667, 0.00376111, 0.0161833}


In reality, I know the Euler angles to be $(30,60,120)$ degrees. How can I get Mathetmatica to give me this?










share|improve this question











$endgroup$












  • $begingroup$
    Have you tried EulerAngles?
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:09










  • $begingroup$
    That command won't work, as the rotation matrix itself is not known. Only the two vectors.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:10






  • 1




    $begingroup$
    In general, a rotation matrix is not uniquely defined by the action on a single vector...
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:12










  • $begingroup$
    Fair enough, but in this case the problem has been constructed such that the three Euler angles are known to exist. The question specifically relates to why NSolve is not working.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:14






  • 1




    $begingroup$
    What I tried to say: If you prescribe a pair $u$ and $v$ of same length $neq 0$ in $mathbb{R}^3$, then there will be a one-parameter family of rotations (and thus a one-parameter family of Euler angles) that map $u$ to $v$. So, it is as it is: Your problem is underdetermined.
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:18
















3












3








3


1



$begingroup$


I would like to determine Euler angles according to the following example.



My example:



I have the three vectors in an original set of axes:



r1e = {-0.517853, 0., -0.759239}
r2e = {-0.517853, 0., 0.759239}
r3e = {0.0647316, 0., 0.}


And after expressing them in a new reference frame they obtain the following components:



rt1e={0.310733, -0.358839, -0.786917}
rt2e={0.690333, 0.298661, 0.527983}
rt3e={-0.0625667, 0.00376111, 0.0161833}


In reality, I know the Euler angles to be $(30,60,120)$ degrees. How can I get Mathetmatica to give me this?










share|improve this question











$endgroup$




I would like to determine Euler angles according to the following example.



My example:



I have the three vectors in an original set of axes:



r1e = {-0.517853, 0., -0.759239}
r2e = {-0.517853, 0., 0.759239}
r3e = {0.0647316, 0., 0.}


And after expressing them in a new reference frame they obtain the following components:



rt1e={0.310733, -0.358839, -0.786917}
rt2e={0.690333, 0.298661, 0.527983}
rt3e={-0.0625667, 0.00376111, 0.0161833}


In reality, I know the Euler angles to be $(30,60,120)$ degrees. How can I get Mathetmatica to give me this?







equation-solving rotation






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited Jan 6 at 18:29







Spherical Cow

















asked Jan 6 at 17:59









Spherical CowSpherical Cow

162




162












  • $begingroup$
    Have you tried EulerAngles?
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:09










  • $begingroup$
    That command won't work, as the rotation matrix itself is not known. Only the two vectors.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:10






  • 1




    $begingroup$
    In general, a rotation matrix is not uniquely defined by the action on a single vector...
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:12










  • $begingroup$
    Fair enough, but in this case the problem has been constructed such that the three Euler angles are known to exist. The question specifically relates to why NSolve is not working.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:14






  • 1




    $begingroup$
    What I tried to say: If you prescribe a pair $u$ and $v$ of same length $neq 0$ in $mathbb{R}^3$, then there will be a one-parameter family of rotations (and thus a one-parameter family of Euler angles) that map $u$ to $v$. So, it is as it is: Your problem is underdetermined.
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:18




















  • $begingroup$
    Have you tried EulerAngles?
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:09










  • $begingroup$
    That command won't work, as the rotation matrix itself is not known. Only the two vectors.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:10






  • 1




    $begingroup$
    In general, a rotation matrix is not uniquely defined by the action on a single vector...
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:12










  • $begingroup$
    Fair enough, but in this case the problem has been constructed such that the three Euler angles are known to exist. The question specifically relates to why NSolve is not working.
    $endgroup$
    – Spherical Cow
    Jan 6 at 18:14






  • 1




    $begingroup$
    What I tried to say: If you prescribe a pair $u$ and $v$ of same length $neq 0$ in $mathbb{R}^3$, then there will be a one-parameter family of rotations (and thus a one-parameter family of Euler angles) that map $u$ to $v$. So, it is as it is: Your problem is underdetermined.
    $endgroup$
    – Henrik Schumacher
    Jan 6 at 18:18


















$begingroup$
Have you tried EulerAngles?
$endgroup$
– Henrik Schumacher
Jan 6 at 18:09




$begingroup$
Have you tried EulerAngles?
$endgroup$
– Henrik Schumacher
Jan 6 at 18:09












$begingroup$
That command won't work, as the rotation matrix itself is not known. Only the two vectors.
$endgroup$
– Spherical Cow
Jan 6 at 18:10




$begingroup$
That command won't work, as the rotation matrix itself is not known. Only the two vectors.
$endgroup$
– Spherical Cow
Jan 6 at 18:10




1




1




$begingroup$
In general, a rotation matrix is not uniquely defined by the action on a single vector...
$endgroup$
– Henrik Schumacher
Jan 6 at 18:12




$begingroup$
In general, a rotation matrix is not uniquely defined by the action on a single vector...
$endgroup$
– Henrik Schumacher
Jan 6 at 18:12












$begingroup$
Fair enough, but in this case the problem has been constructed such that the three Euler angles are known to exist. The question specifically relates to why NSolve is not working.
$endgroup$
– Spherical Cow
Jan 6 at 18:14




$begingroup$
Fair enough, but in this case the problem has been constructed such that the three Euler angles are known to exist. The question specifically relates to why NSolve is not working.
$endgroup$
– Spherical Cow
Jan 6 at 18:14




1




1




$begingroup$
What I tried to say: If you prescribe a pair $u$ and $v$ of same length $neq 0$ in $mathbb{R}^3$, then there will be a one-parameter family of rotations (and thus a one-parameter family of Euler angles) that map $u$ to $v$. So, it is as it is: Your problem is underdetermined.
$endgroup$
– Henrik Schumacher
Jan 6 at 18:18






$begingroup$
What I tried to say: If you prescribe a pair $u$ and $v$ of same length $neq 0$ in $mathbb{R}^3$, then there will be a one-parameter family of rotations (and thus a one-parameter family of Euler angles) that map $u$ to $v$. So, it is as it is: Your problem is underdetermined.
$endgroup$
– Henrik Schumacher
Jan 6 at 18:18












1 Answer
1






active

oldest

votes


















7












$begingroup$

As noted, you can use FindGeometricTransform in tandem with EulerAngles:



r = {{-0.517853, 0., -0.759239}, {-0.517853, 0., 0.759239}, {0.0647316, 0., 0.}};
rt = {{0.310733, -0.358839, -0.786917}, {0.690333, 0.298661, 0.527983},
{-0.0625667, 0.00376111, 0.0161833}};

fg = FindGeometricTransform[r, rt];

EulerAngles[Drop[TransformationMatrix[Last[fg]], -1, -1]]/°
{60.0048, 30.0019, 120.}





share|improve this answer









$endgroup$













    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: "387"
    };
    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',
    autoActivateHeartbeat: false,
    convertImagesToLinks: false,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: null,
    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
    },
    onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    });


    }
    });














    draft saved

    draft discarded


















    StackExchange.ready(
    function () {
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f188949%2fsolving-for-euler-angles%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









    7












    $begingroup$

    As noted, you can use FindGeometricTransform in tandem with EulerAngles:



    r = {{-0.517853, 0., -0.759239}, {-0.517853, 0., 0.759239}, {0.0647316, 0., 0.}};
    rt = {{0.310733, -0.358839, -0.786917}, {0.690333, 0.298661, 0.527983},
    {-0.0625667, 0.00376111, 0.0161833}};

    fg = FindGeometricTransform[r, rt];

    EulerAngles[Drop[TransformationMatrix[Last[fg]], -1, -1]]/°
    {60.0048, 30.0019, 120.}





    share|improve this answer









    $endgroup$


















      7












      $begingroup$

      As noted, you can use FindGeometricTransform in tandem with EulerAngles:



      r = {{-0.517853, 0., -0.759239}, {-0.517853, 0., 0.759239}, {0.0647316, 0., 0.}};
      rt = {{0.310733, -0.358839, -0.786917}, {0.690333, 0.298661, 0.527983},
      {-0.0625667, 0.00376111, 0.0161833}};

      fg = FindGeometricTransform[r, rt];

      EulerAngles[Drop[TransformationMatrix[Last[fg]], -1, -1]]/°
      {60.0048, 30.0019, 120.}





      share|improve this answer









      $endgroup$
















        7












        7








        7





        $begingroup$

        As noted, you can use FindGeometricTransform in tandem with EulerAngles:



        r = {{-0.517853, 0., -0.759239}, {-0.517853, 0., 0.759239}, {0.0647316, 0., 0.}};
        rt = {{0.310733, -0.358839, -0.786917}, {0.690333, 0.298661, 0.527983},
        {-0.0625667, 0.00376111, 0.0161833}};

        fg = FindGeometricTransform[r, rt];

        EulerAngles[Drop[TransformationMatrix[Last[fg]], -1, -1]]/°
        {60.0048, 30.0019, 120.}





        share|improve this answer









        $endgroup$



        As noted, you can use FindGeometricTransform in tandem with EulerAngles:



        r = {{-0.517853, 0., -0.759239}, {-0.517853, 0., 0.759239}, {0.0647316, 0., 0.}};
        rt = {{0.310733, -0.358839, -0.786917}, {0.690333, 0.298661, 0.527983},
        {-0.0625667, 0.00376111, 0.0161833}};

        fg = FindGeometricTransform[r, rt];

        EulerAngles[Drop[TransformationMatrix[Last[fg]], -1, -1]]/°
        {60.0048, 30.0019, 120.}






        share|improve this answer












        share|improve this answer



        share|improve this answer










        answered Jan 6 at 18:38









        J. M. is computer-lessJ. M. is computer-less

        96.3k10301461




        96.3k10301461






























            draft saved

            draft discarded




















































            Thanks for contributing an answer to Mathematica Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid



            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.


            Use MathJax to format equations. MathJax reference.


            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f188949%2fsolving-for-euler-angles%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

            Human spaceflight

            Can not write log (Is /dev/pts mounted?) - openpty in Ubuntu-on-Windows?

            張江高科駅