IEEE Arithematic
조회 수: 1 (최근 30일)
이전 댓글 표시
I have a row vector [3 (6-3e)/(2-2e) 1] and a column vector [3 14 -5-((6-3e)/(2-2e))*14]' where e = 1+0.5*10^-15. When I multiply these two vectors I get a 0 and I don't understand why? The second and third product is equal in magnitude but opposite in sign, so they should cancel out and the answer should be close to 9.This is how matlab computes these two vectors:
row vector = [3.000000000000000e+000 -3.377699720527871e+015 1.000000000000000e+000]
column vector = [ 3.000000000000000e+000 1.400000000000000e+001 4.728779608739018e+016]
댓글 수: 0
채택된 답변
James Tursa
2012년 1월 29일
Limitations of floating point arithmetic, again. Particularly with matrix and vector operations, you are not guaranteed any particular order of operations in MATLAB. If you want a particular order then you would have to manually write the expression using parentheses. So in general, relatively "small" values can get swamped by relatively "large" values in a computation and lose significance. When you subtract two relatively "large" values the relatively "small" value significance may not be there to recover, and you end up with a lot of garbage bits. I could go through your calculation in detail to explain it, but I think I will leave that up to you instead. Just write out the expression in detail and use the num2strexact utility to examine each intermediate result and then you should see what is going on with the lost precision.
댓글 수: 6
James Tursa
2012년 1월 29일
If you download a supported compiler, it should show up in the list when you type mex -setup.
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Downloads에 대해 자세히 알아보기
제품
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!