Friday, March 5, 2010

Pounds, shillings, pence as a system of logarithms

To the memory of Orson Pratt, 1811-1881. [Professor of Mathematics in the University of Nauvoo, and LDS (Mormon) apostle. Pratt presided over the LDS church's British Mission at a time when there were twice as many LDS in Britain as in the United States, so he would certainly have been familiar with the pounds, shillings, pence system of currency, in which 1 pound is worth 240 pence. A couple of decades later, Pratt, apparently at his own request, was assigned by Brigham Young to redivide the Book of Mormon into chapters and verses, a division which stands today: a canonical and somewhat lengthy colophon known as the Title Page, and 239 chapters, so 240 chapter-length divisions in all. At this remove, it is not possible to determine whether this coincidence of numbers was intentional, unconscious, or happenstance. Nevertheless, this system would not have been devised by this author without Pratt's subdivision, because it arose originally in connection with an effort to devise a mnemotechnical scheme for the chapter-by-chapter content of the Book of Mormon.]

The purpose of this posting is to describe a mixed system of logarithmic units.

By the end of the posting, with just the use of a table of 19 (easily memorized) values, you will be able to find arbitrary powers and roots of arbitrary numbers, generally to within one percent accuracy.

This system could, in principle, have been invented a couple of hundred years ago.  It is contructed mostly by stitching together a large variety of preexisting pieces: decades, decibels, R-series prefered numbers, pounds-shillings-pence money, order-of-magnitude arguments, percentage tweaking calculations.

My particular contribution here is putting together a mixed system of logarithmic units.  I have never seen that done anywhere else.  The system particularly exploits the numerical coincidence that

            101/2401.01,

where
  • 10 is the base of the decimal number system we use, 
  • 1.01 is the factor corresponding to a tweak of one percent, (percent being the usual unit of tweaking arguments,) and 
  • 240 is the number of pence in a pound, a highly composite number.

An earlier posting defined the (well-known) logarithmic units decade, neper, and octave as lg(10), lg(e), and lg(2), respectively.  All logarithmic units in common use are at most simple rational multiples of one of these three units.  So, there are decadic units, neperic units, and octavic (or binary) units.

First, let's define three decadic units, and then show why they are particularly useful as a system.

These logarithmic units are the decade (dec), the decibel (dB), and the point (pt). Let's defer discussion of the merits of these particular names.

The old British Imperial currency, its lineage dating back past Charlemagne, had a pound equal to twenty shillings, and a shilling equal to twelve pence.

These logarithmic units share the same arithmetic.  The decade is twenty decibels, and the decibel is twelve points.  Thus
  • 1 decade=lg(10),
  • 1 decibel=120lg(10)=lg(101/20)lg(1.12), and
  • 1 point=1240lg(10)=lg(101/240)lg(1.01).
The decade is the multiplying strength of (a factor of) ten.

Arguments that use round powers of ten are a commonplace in technical disciplines.  They are known as order-of-magnitude arguments, and people use expressions like "three orders of magnitude" to mean "very roughly a factor of a thousand".  The decade is a precise order of magnitude.  Order of magnitude arguments are, at root, logarithmic arguments.

The common logarithm, the one found by the 'log' key on a calculator, is the number of decades in the multiplying strength of a number.

The decibel is the multiplying strength of the twentieth root of ten, which is about 1.12.  Integer powers of this, i.e. numbers of the form 10n/20 for n an integer, constitute the R20 series of prefered numbers, as now maintained in an international standard by the ISO.  The 'R' in R20 honours the French engineer Charles Renard, the inventor of this system.  The 20 tells us that we go from one power of ten to the next in twenty steps.  These steps---this is Renard's invention---go up in constant proportion.

All of the values in stripy multiplication table

            [11.783.165.62101.783.165.621017.83.165.621017.831.65.621017.831.656.21017.831.656.2100].

come from what could reasonably be called the R4 series of prefered numbers.  Interposing rows and columns, we can build a similar table for the R20 series.  The R20 values between 1 and 10, rounded to three significant figures, are the values in the second column.  Their logarithms, expressed in decibels, appear in the first column ":

             0 dB    11 dB    1.122 dB    1.263 dB    1.414 dB    1.585 dB    1.786 dB    2.007 dB    2.248 dB    2.519 dB    2.8210 dB    3.1611 dB    3.5512 dB    3.9813 dB    4.4714 dB    5.0115 dB    5.6216 dB    6.3117 dB    7.0818 dB    7.9419 dB    8.9120 dB    10

Put another way, if n is the numerical part of the logarithm in the first column, then the number in the second column is 10n/20.

In the seventh edition of the SI Handbook, published by the BIPM, the bel was defined as half of what we have been calling a decade.  The decibel is a tenth of this, and so it is a twentieth of a decade.  The eighth edition of the SI Handbook removed this definition of the bel, but did not supersede it.

The R20 values form a bridge between the two other kinds of units.  Also, since, to three significant digits, the number 2 is an R20 value, the octave, lg(2), is approximately a round number of decibels.  This gives us a tie between octavic units and decadic units.

The point is the multiplying strength of the 240th root of ten, which is about 1.01.  A change of one point is roughly equivalent to a change of one percent.  It is also roughly equivalent to a change of one centineper, which gives a tie to neperic units.

Percentage tweaking arguments are widely used in technical disciplines.  If a quantity goes up by three percent, say, then its square goes up by about six percent.  If the length of a rectangle goes up five percent, while its width goes up by two percent, than its area increases by about seven percent.  Addition of small percentage tweaks is used to do the work of multiplication---this is the hallmark of logarithms.  Small percentage tweaks used in this way are being used as approximate logarithms.

This mixed system of logarithmic units is particularly well suited to calculation of approximate logarithms.

As an example, let's calculate the logarithm of 72900.  First, we rewrite this in reverse scientific notation, i.e. we factor it as a whole power of ten multiplied by a number between 1 and 10:

            72900=104×7.29.

Next, we break the second factor into two factors.  The first factor is selected from the R20 values in the table.  Since 7.29 is between 7.08 and 7.94, we will pick one of these values.  We usually pick the lower table value, unless our value is within about two percent of the higher value.  Then

            72900=104×7.08×7.297.08.

We now calculate the ratio 7.29/7.08.  This could be done by long division, but we prefer instead the folowing :

        7.297.08=729708=708+21708=1+21708.

To get three or four digits of accuracy in the sum 1+21/708, we need only one or two digits of accuracy in 21/708.  We can therefore afford to round the denominator, so that

            2170821700=3100.

Then 1+21/7081.03.

The same conclusion can be reached by making the proportion table :

             708100%7110%71%21?.

The first row always consists of the target divisor or denominator, here 708, and 100%.  One creates the second row by dividing through by ten or another simple divisor, until the value in the first column is comparable to or smaller than the target dividend or numerator, here 21.  Numbers in the first column are usually rounded to the same number of decimal places (here zero) as in the first row.  If the target numerator is a small multiple of the target numerator---here 21=3×7, this gives us what we are looking for, i.e. the desired percentage is 3×1%=3%.

(We don't want a high multiple, especially if rounding has been severe.  If, for instance, the numerator had been 51 instead of 21, then we would created yet another row by taking three times the [7,1%] row from the [71,10%] row and we would have had [50,7%], so that 51/710.07.)

Either way, here, we get 

            72900=104×7.08×1.03.

From the table, we can rewrite 7.08=1.1217.  The binomial theorem gives us that (1+x)n1+nx, and so 1.031.013.  So

            72900104×1.1217×1.013.

Then

            lg(72900)lg(104×1.1217×1.013),
                =lg(104)+lg(1.1217)+lg(1.013),
                =4lg(10)+17lg(1.12)+3lg(1.01),
                =4 dec+17 dB+3 pt.

Compound arithmetic with mixed units can be used for adding or subtracting the logarithms of different numbers, and even for multiplying or dividing by simple integers.  These are the processes one needs when one is using logarithms to multiply, divide, or take simple roots or powers.

For more complicated powers, and often in other calculations, it more convenient to work in a consistent units.  We can rewrite this mixed form using wiggle notation, and use that to calculate values in one or other of the given units, so that

            lg(72900)4;2017123 dec,
                =42017;123 dB,
                =42017123; pt.

For instance,

            4;2017123=4+17×12+320×12,

                =4+207240,
                =4+207×1.04240×1.04,
                4+215250,
                =4+8601000,
                =4.86.

i.e. lg(72900)4.86 dec, or log10(72900)4.86.  This compares well with the more accurate value of 4.8627275.

It is usually more convenient to work in points.  The conversion doesn't require further division, at least here, and each point corresponds to a one percent difference in the answer, so it is easier to make decisions about acceptable error.

            lg(72900)42017123; pt,
            =(80+17)123; pt,
            =97×12+3 pt,
            =1164+3 pt,
            =1167 pt.

Suppose we wanted to find 729003/7.

The logarithm of 729003/7 is 3/7 of the logarithm of 72900, i.e. about 3/7 of 1167 points.  1167×3=35013501÷7500.1.  So the logarithm of our answer is aboout 500.1 pt.

Reverting this into mixed units, we use divmod division, first by 12, and then by 20.

            500.1÷12=41;128.1, and


            41÷20=2;201,

so that

             500.1=2201128.1,

and

            lg(7290037)2 dec+1 dB+8.1 pt.

Then

            7290037102×1.12×1.081,
                =121.,

which compares well with the more exact value 121.346, especially considering that table values such as 1.12 have already been rounded to three significant figures.