Calculating power in a loop has different result from the matlab power function
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
kuku hello
am 21 Mär. 2022
Kommentiert: kuku hello
am 21 Mär. 2022
Hello, I tried to calculate 13^(15) using a for loop, and using the power function of matlab and I get different results.
Can someone explain to me how come? For example:
num = 13;
pow = 15;
for i = 2:pow
num = num*13;
end
result = 13^15;
num - result % The result here is 8, which means the difference between the two numbers is 8?
result ~= num % It comes out to be true
0 Kommentare
Akzeptierte Antwort
Jan
am 21 Mär. 2022
Bearbeitet: Jan
am 21 Mär. 2022
Neither the power operator nor the loop produce the exact result:
x = sym(13);
x^15
fprintf('%.0f\n', 13^15)
y = 13;
for k = 2:15, y = y * 13; end
fprintf('%.0f\n', y)
This effect is expected, because the value is larger than the maximum value, which has an exact representation in the IEEE754 format. This format is used in Matlab, other numerical software and in modern CPUs. It uses 52 bits for the mantissa. Then numbers, which need a 53'th bit or more cannot be store accurately anymore.
This happens in the 15th iteration of the loop, when the value exceeds 2^53, which is about 9e15. See: doc flintmax
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Logical finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!