Storing small numbers and logarithms.

3 次查看（过去 30 天）
daniel adams
daniel adams 2021-11-30
评论： daniel adams 2021-12-1
Hi I have the following problem, in my code I need to calculate for very large. The largeness causes matlab to store and to be zero. Does anyone know a way around this? Is there a way of taking logarithms first?
  1 个评论
Walter Roberson
Walter Roberson 2021-11-30
What precision do you need for output? If a and b differ by more than about 36.04 then to numeric precision, the answer is going to be the same as -min(a,b) . If they differ by N then approximately the first (3/2)*abs(N) bits are conserved

请先登录，再进行评论。

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

回答（2 个）

Walter Roberson
Walter Roberson 2021-11-30
Ran in:
Depending on your range of values and your accuracy needs, you could try a double taylor approximation.
syms a b positive
E = log(exp(-a) + exp(-b))
E = 
Ta = taylor(E, a, 'ExpansionPoint', 10^5, 'Order', 10)
Ta = 
Tb = taylor(Ta, b, 'ExpansionPoint', 10^5, 'Order', 10)
Tb = 
Ts = simplify(expand(Tb))
Ts = 
Tsa = collect(Ts, [a])
Tsa = 
Steven Lord
Steven Lord 2021-11-30
Ran in:
How large is "very large"? Can you perform the calculation symbolically (by converting a and b into sym values before trying to pass them into exp -- if you try to perform the exp calculations in double precision then convert to sym you've already underflowed) and convert the result back to double afterwards?
a = 12345;
as = sym(a);
L = log(sym(exp(-a))) % exp(-a) underflowed in double before being converted to sym
L = 
Ls = log(exp(-as)) % No underflow
Ls = 
  3 个评论
显示 1更早的评论
Steven Lord
Steven Lord 2021-12-1
Ran in:
1) Yes, the symbolic exp and log functions operate element-wise just like the numeric exp and log functions.
a = 1:10;
as = sym(a);
e = exp(a);
es = exp(as);
format longg
double(e-es) % Small differences
ans = 1×10
1.0e+00 * 2.99524520677138e-16 1.79711394978391e-16 1.82756255255129e-16 -2.87415780158441e-15 -3.48635149004642e-15 -1.23596280244504e-14 -9.86975264043409e-14 2.71032958168736e-14 2.15308776210672e-13 1.37801347005174e-12
2) Yes, there is some additional overhead for the symbolic calculations. How large is "very large" in "very large dimensional vectors"?
daniel adams
daniel adams 2021-12-1
Hi @Steven Lord thank you greatly for your help so far this is a problem that has really been bothering me! For more context can you see my new question ( which is the full version of what I was asking on this post ) https://uk.mathworks.com/matlabcentral/answers/1600419-using-symbolic-math-to-retain-accuracy-arrays By very large the matrices will be of size .

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Assumptions 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by