Understanding Precision in Matlab

3 次查看（过去 30 天）
Ahmad Gad
Ahmad Gad 2020-7-1
评论： Walter Roberson 2020-7-1
I have a simple question. How to split the different between the two numbers into 1e6 increment (unisg linspace)?
The numbers: 5.047639467288342e+4 and 5.047639467288343e+4.
Thanks
  2 个评论
Stephen23
Stephen23 2020-7-1
"How to split the different between the two numbers into 1e6 increment (unisg linspace)?"
You can't. There are no binary floating point numbers inbetween the two values that you gave.
We can confirm this very easily by noting that their hex representations differ only by 1:
>> A = 5.047639467288342e+4;
>> B = 5.047639467288343e+4;
>> num2hex(A)
ans = 40e8a58ca12906dd
>> num2hex(B)
ans = 40e8a58ca12906de
Ahmad Gad
Ahmad Gad 2020-7-1
Thank you! Then, I have to consider changing the algorithm then.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2020-7-1
You cannot do that with double precision.
If you have the symbolic toolbox then:
r = linspace(sym('5.047639467288342e+4'), sym('5.047639467288343e+4'), 1e6+1);
Note that the results will be symbolic numbers with higher precision than double(), so if you were to double() them, it would be with loss of precision
  2 个评论
Ahmad Gad
Ahmad Gad 2020-7-1
The code is numeric and those values will be inserted a matrix determinant function. I think it will be time consuming.
Walter Roberson
Walter Roberson 2020-7-1
Determinants can have a lot of problems with cancellation. If you were to use symbolic numbers it would be better to use
r = linspace(sym(5.047639467288342e+4,'f'), sym(5.047639467288343e+4,'f'), 1e6+1);
However, symbolic determinants are slow. To put 1e6 numbers into an array you would need an array at least 1e3 x 1e3, and a proper accurate determinent of that would involve over 4 * 10^2567 terms.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by