Numerical Precision Physics Calculations

2 次查看(过去 30 天)
I believe I may be encountering issues related to the numerical precision, but my understanding is limited so I wanted to ask the community. Hopefully I'm misunderstanding and someone can clarify!
I have a physics problem where I'm working with weak spatially varying magnetic fields, and calculating their response to conductive materials. For example, one calculation may involve a vector of doubles which range from 3e-7 to 4.9e-12 multiplied by one which has values of 0 or 5.89e7.
eps(5.89e7)
ans = 7.4506e-09
does this mean that all values below that eps value behave as zeros because there isn't sufficient precision for the calculation? Or do calculations work with say 4.9 *5.89 then e-12 *e7 seperately?

采纳的回答

Fangjun Jiang
Fangjun Jiang 2022-8-25
编辑:Fangjun Jiang 2022-8-25
eps(5.89e7)
ans = 7.4506e-09
It means that around 5.89e7 (which is a large value), the nearest values that can be represeted by double data type is 5.89e7+7.4506e-09 or 5.89e7-7.4506e-09. So you can see the minimal incremental value of 7.4506e-09 is quite small.
eps('double')
ans = 2.2204e-16
Since you are using the default double data type, there should be no problem representing your multiplication in the range of
5.89e7*4.9e-12
ans = 2.8861e-04
  1 个评论
ADSW121365
ADSW121365 2022-8-25
编辑:ADSW121365 2022-8-25
My confusion arose from that first part. If I were to add my vectors, 5.89e7 + 4.9e-12, the calculation would be below the minimal incremental value & therefore the difference between the 5.89e7 and the result would not be represented.
I'm glad to know this doesn't impact multiplications, thank you.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by