Calculating power in a loop has different result from the matlab power function
조회 수: 2 (최근 30일)
이전 댓글 표시
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
채택된 답변
Jan
2022년 3월 21일
편집: Jan
2022년 3월 21일
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
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Logical에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!