Why different result when sum matrix by different order
显示 更早的评论
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
采纳的回答
Chunru
2021-5-25
1 个投票
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.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Resizing and Reshaping Matrices 的更多信息
Translated by 
