Problem 961. Robust alignment of coordinate system

Created by J-G van der Toorn

Solve Later

{"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","targetMode":"","relationshipId":"rId1","target":"/matlab/document.xml"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/output","targetMode":"","relationshipId":"rId2","target":"/matlab/output.xml"}],"parts":[{"partUri":"/matlab/document.xml","relationship":[],"contentType":"application/vnd.mathworks.matlab.code.document+xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n<w:document xmlns:w=\"http://schemas.openxmlformats.org/wordprocessingml/2006/main\"><w:body><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>We have to create a custom 3-D coordinate system of which one axis is defined by a vector (in the global coordinate system). The other two coordinate axes can be freely chosen. But the most intuitive would be (in this case) to align an axis of the custom coordinate system with the global one, either which one matches best.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Return the 3element normalized column vector P in a 3-D Cartesian space (regular x,y,z axes we use in general), that * is in a plane (through the origin), defined by its normal vector N * N can be a column- or row vector * N is not necessarily normalized * makes the smallest angle with any of the Cartesian axes specified A * A is a character array containing any of the character xyz and XYZ * when the smallest angle of P with multiple axes mentioned in A is equal, the order of the axes in A defines which is preferred * when A contains only one axis-name, and N is parallel to that axis, the functions return empty an empty matrix [].</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Example 1:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ align_frame([1 0 0],'zyx')]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>returns:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ 0\n 0\n 1]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Example 2:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ align_frame(1:3,'XYZ')]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>returns</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ 0.9636\n -0.1483\n -0.2224]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Example 3:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ align_frame([1;1;0],'xy')]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>returns</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ 0.7071\n -0.7071\n 0]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Example 4:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ align_frame([0 0 1],'Z')]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>returns</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ []]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Note: Solutions wrapped in (f)eval(c), inline, str2func, regexprep (dynamic regular expressions), etc, are not appreciated.</w:t></w:r></w:p></w:body></w:document>"},{"partUri":"/matlab/output.xml","contentType":"text/xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?><embeddedOutputs><metaData><evaluationState>manual</evaluationState><layoutState>code</layoutState><outputStatus>ready</outputStatus></metaData><outputArray type=\"array\"/><regionArray type=\"array\"/></embeddedOutputs>"}]}

We have to create a custom 3-D coordinate system of which one axis is defined by a vector (in the global coordinate system). The other two coordinate axes can be freely chosen. But the most intuitive would be (in this case) to align an axis of the custom coordinate system with the global one, either which one matches best.Return the 3element normalized column vector P in a 3-D Cartesian space (regular x,y,z axes we use in general), that
* is in a plane (through the origin), defined by its normal vector N
 * N can be a column- or row vector
 * N is not necessarily normalized
* makes the smallest angle with any of the Cartesian axes specified A
 * A is a character array containing any of the character xyz and XYZ 
 * when the smallest angle of P with multiple axes mentioned in A is equal, the order of the axes in A defines which is preferred
* when A contains only one axis-name, and N is parallel to that axis, the functions return empty an empty matrix [].Example 1:<pre> align_frame([1 0 0],'zyx')</pre>returns:<pre> 0
 0
 1</pre>Example 2:<pre> align_frame(1:3,'XYZ')</pre>returns<pre> 0.9636
 -0.1483
 -0.2224</pre>Example 3:<pre> align_frame([1;1;0],'xy')</pre>returns<pre> 0.7071
 -0.7071
 0</pre>Example 4:<pre> align_frame([0 0 1],'Z')</pre>returns<pre> []</pre>Note: Solutions wrapped in (f)eval(c), inline, str2func, regexprep (dynamic regular expressions), etc, are not appreciated.

We have to create a custom 3-D coordinate system of which one axis is defined by a vector (in the global coordinate system). The other two coordinate axes can be freely chosen. But the most intuitive would be (in this case) to align an axis of the custom coordinate system with the global one, either which one matches best.

Return the 3element normalized column vector P in a 3-D Cartesian space (regular x,y,z axes we use in general), that
* is in a plane (through the origin), defined by its normal vector N
  * N can be a column- or row vector
  * N is not necessarily normalized
* makes the smallest angle with any of the Cartesian axes specified A
  * A is a character array containing any of the character xyz and XYZ 
  * when the smallest angle of P with multiple axes mentioned in A is equal, the order of the axes in A defines which is preferred
* when A contains only one axis-name, and N is parallel to that axis, the functions return empty an empty matrix [].

Example 1:

align_frame([1 0 0],'zyx')

returns:

0
   0
   1

Example 2:

align_frame(1:3,'XYZ')

returns 
 
    0.9636
   -0.1483
   -0.2224
 
 
Example 3:

align_frame([1;1;0],'xy')

returns 
 
    0.7071
   -0.7071
    0
 
Example 4:

align_frame([0 0 1],'Z')

returns 
 
 []
 
Note: Solutions wrapped in (f)eval(c), inline, str2func, regexprep (dynamic regular expressions), etc, are not appreciated.

Solve

Solution Stats

15 Solutions
2 Solvers

Last Solution submitted on Jul 01, 2017

Last 200 Solutions

Solution Comments

Solution 259208

1 Comment

1 Comment

Andrew Newell on 10 Jun 2013

This solution is correct to within eps, but the test is expecting exact agreement. That is not good programming when dealing with floating point numbers.

Problem Recent Solvers2

Problem Tags

3d geometry vector

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!