Unexpected result from linspace function

조회 수: 4 (최근 30일)
Rostislav Teryaev
Rostislav Teryaev 2018년 8월 15일
댓글: Dyuman Joshi 2023년 12월 9일
Why the result of this code
B1 = 0:0.1:1.8;
B2 = linspace(0,1.8,19)
B1(4)-0.3
B2(4)-0.3
is 5.5511e-17? I found a workaround, which is
B1 = 0:01:18;
B1 = B1/10;
But I am still wondering why the first code does not produce what it should?

채택된 답변

Rik
Rik 2018년 8월 15일
Welcome to the wondrous world of floating point numbers. Matlab is a computer program, so it works with binary. Some values are impossible to represent in a binary expansion, so they get rounded. This can be different for each algorithm, so even if the end result can be represented, rounding errors can build up. That is the use case for the eps function: it shows you what the magnitude is of possible errors.
As a decimal example, consider this: (1/3)*3. Assuming I can only store 2 decimal digits, that would be 0.33*3=0.99, even if the correct answer doesn't need 2 decimals.
  댓글 수: 9
Stephen23
Stephen23 2023년 12월 9일
"Is there any work around to deal with floating point arithmetic problems in matlab ?"
The problem is not "in MATLAB", it is simply how your computer stores floating point numbers. MATLAB runs on your computer: if there was a simple "work around" then everyone would use it, not just MATLAB.
You can certainly use symbolic or high precision classes, but then your code will be much slower.
Dyuman Joshi
Dyuman Joshi 2023년 12월 9일
"Is there any workaround to deal with floating point arithmetic problems?"
@VIGNESH BALAJI - There, I fixed the question for you.
Certainly, one option is to use a computer that runs on a decimal system. There might be a lot of complications, though.
As long as you use computers that use the binary system, you will encounter floating point errors.
Rik and Stephen have already provided suggestions as to what you can do as workarounds, I'll just reiterate them -
You can avoid using floats.
Well, what to use then? Symbolic numbers or high precision classes, which work at the cost of the code performance.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Logical에 대해 자세히 알아보기

제품


릴리스

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by