Why different result when sum matrix by different order

Question

Li Fujian 2021-5-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/838710-why-different-result-when-sum-matrix-by-different-order

评论： Li Fujian 2021-5-26

采纳的回答： Chunru

When I sum a matrix by different order, the summation is different, for example:

a=rand(10)/pi;

sum(sum(a,1),2)-sum(sum(a,2),1)

the code will return a non-zeros result sometimes, which cause some problem in my calculation.

the matlab edition Im using is R2016a

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Chunru 2021-5-25

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/838710-why-different-result-when-sum-matrix-by-different-order#answer_707905

This is due to the numerical accuracy of the floating point operations.

Let's say we have a CPU with 4 significant digits. and we want to perform 0.1234+ 0.2345 + 345.6.

If we add in the first two small numbers, the result is 0.3579; Then we add to 345.6, we will get 346.0 (all results, intermediate or final, having 4 significant digits).

Now if we perform the addition in a reversed order. 345.6 + 0.2345 ~= 345.8; Then 345.8 + 0.1234 ~= 345.9. Thus we have different results by summing in different order.

In numerical calculation, one obtains more accurate result by performing summation of small numbers first. For example, sum of the series 1/n^2, n=1,..N, in reversed order, usually is a better choice (especially when N is very big).

Hope this help.