File Exchange

image thumbnail

Methods for calculating precise logarithm of a sum and subtraction

version 1.3 (106 KB) by

Methods for operating on log's of values without calculating the original values during the process

1 Download

Updated

View License

Disclosed are methods to compute the precise value of the logarithm of a sum and the logarithm of a subtraction.

A number of practical problems can result in having too big or small values in intermediate values of a calculation. Then one tries to take logarithm of these values and operate on logarithms instead.

In the case: log (p • q) = log p + log q is very easy to compute, but the problem is then to compute (or approximate) the value of log (a + b) from the value of log a and log b.

Let us assume that log a and log b are known, and that we want to approximate log (a + b).
Most basic solution would be calculating
sum_log = log(exp(a_log) + exp(b_log),
where a_log = log a and b_log = log b and therefore sum_log = log(a + b).
But the method I propose requires calling only one exp() and one log(), instead of two exp() and one log() in the basic solution.
Additionally, the proposed method has the critical advantage of not overflowing in case of large numbers of a and b.

The method is based on the notion that
ln(a + b) = ln{exp[ln(a) - ln(b)] + 1} + ln(b).

In a similar manner we can compute the logarithm of subtraction, based on the notion that
ln(a - b) = ln{exp[ln(a) - ln(b)] - 1} + ln(b), if a > b.

Comments and Ratings (1)

I can see this being a useful utility for some applications. However the code could be improved. I would like to see you add
- a H1 line
- an example (or two) showing usage, particularly comparing your method with a method which will produce overflow
- Comments within the code explaining the procedure so a user does not need to consult the pdf file

Lastly I think the actual computation should be on just one line
i.e. R = log(exp(a_ln - b_ln) + 1) + b_ln;
rather than spread out over many lines. Let comments explain the shifting.

Darren

Updates

1.3

Seems like the updated zip file did not upload with my previous post. I will now resend it.

1.2

1) Included method for calculating logarithm of subtraction.
And according to the suggestions of Darren Rowland:
2) Added H1 line.
3) Inlined the main formula.
4) An example comparing the methods with a naive solution which will produce overflow.

MATLAB Release
MATLAB 7.4 (R2007a)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video