cross product is wrong

8 visualizaciones (últimos 30 días)
savas yilmaz
savas yilmaz el 5 de Ag. de 2021
Editada: Paul el 6 de Ag. de 2021
>> A
A =
-2159/277
-27/2
0
>> B
B=
-579/104
-135/14
0
>> cross(A,B)
ans =
0
0
-1/70368744177664
A and B vectors are parallel. cross product must "0". But answer is wrong. Can you help me?
  3 comentarios
savas yilmaz
savas yilmaz el 6 de Ag. de 2021
A =[ -2159/277 ; -27/2 ; 0] and A=sym([-2159/277 -27/2 0]) are not equal.
Paul
Paul el 6 de Ag. de 2021
Editada: Paul el 6 de Ag. de 2021
In the first instance, A is a double and in the second instance A is a sym object. I was using the Symbolic Math Toolbox to determine if the problem you're seeing is just a rounding error. But the data provided in the question for A and B doesn't yield the result in the question using either symbolic or numeric math:
syms A B
A = sym([-2159/277 -27/2 0]);
B = sym([-579/104 -135/14 0]);
vpa(cross(A,B))
ans = 
cross(double(A),double(B))
ans = 1×3
1.0e+-4 * 0 0 -0.6695
So the symbolic math and the numeric math basically give the same result, neither of which match the result in the question:
-1/70368744177664
ans = -1.4211e-14

Iniciar sesión para comentar.

Respuestas (2)

Rik
Rik el 5 de Ag. de 2021
Editada: Rik el 5 de Ag. de 2021
1/70368744177664
ans = 1.4211e-14
eps
ans = 2.2204e-16
This number is so small that you can assume this is a rounding error.
The root cause of this rounding error is that it is not possible to store infinite decimals in finite computer memory.
Imagine you can only store decimal numbers, and only 4 digits after the decimal point.
(1/3)*3 = (0.3333)*3 = 0.9999
Is * wrong? Or /? Yes and no. You just need to be aware of the limitations of the tools you're using.
  2 comentarios
savas yilmaz
savas yilmaz el 6 de Ag. de 2021
This number is so small but I will use this for "if" function. That's why it's important to me that the result is zero.
The answer to the same question is zero in geogebra program. But I need this for matlab.
Rik
Rik el 6 de Ag. de 2021
Then compare to a tolerance:
val=1/70368744177664;
target=0;
tol=1e-10;
if abs(val-target)<=tol

Iniciar sesión para comentar.


Jan
Jan el 6 de Ag. de 2021
Geogebra treats the dta as symbolical expressions. Matlab converts -2159/277 to a numericalvalues as default. If you want to use symbolic calculations, use sym().
Do not check numerical floating point values to be exactly 0 in a condition of an if command. Calculations with numerical values include rounding effects, which cannot be neglected:
1e17 + 1 - 1e17 % 0
1e17 - 1e17 + 1 % 1
sin(2 * pi) % -2.4493e-16
...
There are no well defined fixed limits and it depends on the application what you have to "consider as 0".
A numerical algorithm requires an analysis of the rounding effects to estimate the reliability of the result.

Community Treasure Hunt

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

Start Hunting!

Translated by