Как пользоваться таблицами Брадиса
Каждая таблица Брадиса дает значения какой-либо величины (функции) в зависимости от значения некоторой другой величины (аргумента). Например, таблица 3 дает значения квадрата в зависимости от значений возводимого в квадрат числа (функция у = х2 аргумента х), таблица 7 — значения площади круга в зависимости от значений его диаметра (функция К = πd2/4 аргумента d) и т. д. Ради экономии места все таблицы сборника расположены «в два хода»: каждое табличное значение функции находится на пересечении строки, имеющей в заголовке (слева) некоторые первые цифры соответствующего назначения аргумента, и столбца, имеющего в заголовке (сверху) остальные его цифры. Например, для квадрата числа 5,67 находим в таблице 3 на строке 5,6 в столбце 7 значение 32,15, представляющее собой результат округления до 4 значащих цифр точного квадрата 5,672 = 32,1489.
Все табличные значения функций, приведенные в сборнике, получены путем округления до 4 или 5 значащих цифр соответствующих точных значений, а потому отличаются от точных не более как на половину единицы разряда последней цифры. Например, найдя из таблицы 7, что площадь круга диаметра 2,16 линейных единиц равна 3,664 соответствующих квадратных единиц (строка 2,1, столбец 6), мы можем быть уверены, что точное значение этой площади отличается от этой табличной не больше чем на половину тысячной, то есть что 3,6635
Значения аргумента в каждой таблице равномерно растут (по крайней мере в некотором интервале), и постоянное значение разности двух соседних значений аргумента носит название «ступени» таблицы. Так, в таблице 3 ступень везде 0,01, а в таблице 4 сперва 0,01, потом 0,1. Значения функции в большинстве таблиц тоже растут, но равномерным их рост оказывается только для линейных функций, т. е. функций вида у = ах + Ь, где а и Ь — постоянные. Увеличение х на ступень h дает у таких функций увеличение функции на постоянное число ah. Например, при увеличении диаметра на 0,01 длина окружности С = πd увеличивается на 0,01π = 0,0314…. Просматривая табличные значения длины окружности, замечаем, что при возрастании d на 0,01 они возрастают то на 31, то на 32 тысячных. Это небольшое колебание вызвано приближенным характером табличных значений.
Разность двух соседних табличных значений функции называется «табличной разностью». Имея дело с таблицей функции, изменяющейся неравномерно, следует различать два случая: случай «почти равномерного» изменения функции, когда табличные разности изменяются очень медленно, и случай «резко неравномерного» ее изменения, когда уже соседние табличные разности отличаются друг от друга на несколько единиц последнего разряда. Так, в таблице кубов 13 = 1, 23 = 8, 33 — 27, 43 = 64,… мы имеем пример таблицы с резко неравномерным изменением функции, но если ту же таблицу кубов взять со ступенью не в 1, а в 0,001 и округлять кубы до 4 значащих цифр, то получится таблица с почти равномерным изменением функции, которую мы и имеем, где на протяжении всей строки 1,00 табличные разности равны 3 (тысячным), а на нескольких следующих — то 3, то 4. Различие между таблицами с равномерным, почти равномерным и резко неравномерным изменением функции проявляется особенно наглядно при изображении этих функций посредством графиков в прямоугольных координатах: в первом случае получается график в виде прямой, во втором — в виде кривой, небольшие участки которой искривлены едва заметно, в третьем — в виде кривой с заметной кривизной уже на каждом малом участке. Одна и та же таблица может быть таблицей с почти равномерным изменением функции на одном интервале и с резко неравномерным ее изменением на другом. Такова, например, таблица 9, где на передних строках изменение функции резко неравномерно. Имея таблицу с резко неравномерным изменением функции, можно превратить ее в таблицу с почти равномерным изменением двумя способами: уменьшением ступени таблицы, то есть заменой ее другой, более подробной, что делается не так просто (надо либо иметь такую более подробную таблицу, либо заново ее составить), или округлением табличных значений, что делается очень просто, но связано с потерей точности. Например, тангенсы углов, указанные в таблице 10 на строке 89°20′ с точностью до сотых и изменяющиеся резко неравномерно, после округления до десятых становятся изменяющимися почти равномерно.
Каждая таблица содержит значения функции не для всех, а лишь для некоторых значений аргумента. Возникает вопрос: как получить значения функции для промежуточных значений аргумента? Операция получения таких значений носит название «интерполяции». Ее иногда образно называют «чтением между строками таблицы».
В случае таблицы с равномерным или почти равномерным изменением функции применяется так называемая «линейная интерполяция», состоящая в следующем. Если при увеличении значения аргумента на И единиц какого-либо разряда функция увеличивается на d единиц некоторого разряда, то в силу равномерности изменения функции увеличение аргумента на 1 вызывает увеличение функции на d/h единиц, а увеличение аргумента на u — увеличение функции на v = du/h единиц. Очевидно, что для получения искомого значения функции надо взять ближайшее меньшее табличное ее значение и прибавить эту «поправку» v. Например, чтобы узнать, чему равен квадрат числа 8,053, берем в таблице 3 8,052 = 64,80, 8,062 = 64,96, 8,072 = 65,12 и убеждаемся, что изменение функции здесь почти равномерно: при ступени h = 0,01 или 10 тысячным табличная разность составляет здесь 16 сотых. Данное значение аргумента 8,053 превосходит ближайшее меньшее табличное его значение 8,05 на u = 3 (тысячным), а потому поправка v равна 16·3/10 = 4,8 = 5 (сотым). Прибавив ее к ближайшему меньшему табличному значению функции 8,052 = 64,80, получим 8,0532 = = 64,80 + 0,05 = 64,85 (непосредственное умножение дает точно 8,0532 = 64,850809).
Вместо того чтобы брать поправку на «избыток» данного значения аргумента над ближайшим меньшим табличным его значением, как мы это только что делали, можно дать поправку на его «недостаток» по сравнению с ближайшим большим табличным его значением и вычитать поправку из ближайшего большего значения функции. Например, для получения квадрата числа 8,057 берем 8,062 = 8,0602 = 64,96 и вычитаем поправку на 3 тысячных, равную 5 сотым, получая 8,0572 = 64,91 (при точном значении, равном 64,915249). Поправка на избыток выгоднее, если избыток не превосходит половины ступени; в противном случае выгоднее брать поправку на недостаток.
Операцию линейной интерполяции можно объяснять, исходя не из равномерности изменения функции, как мы это сейчас делали, а из пропорциональности приращений аргумента и функции, т. е. из пропорциональности избытка аргумента и поправки для функции, прекрасно иллюстрируемой на графике, где получаются два подобных прямоугольных треугольника, один с катетами h и d, другой с катетами u и v. По существу, оба способа, конечно, тождественны, так как оба основаны на одной и той же пропорции u_v=h:d.
Какова точность результатов, получаемых посредством линейной интерполяции? Здесь имеются три источника погрешностей: неточность взятого ближайшего табличного значения функции, не превосходящая половины единицы разряда последней его цифры; неточность поправки, обусловленная неточностями табличных значений и округлением поправки, и, наконец, неточность поправки, вызванная неполной равномерностью изменения функции. Более глубокое рассмотрение вопроса показывает, что при разнице двух соседних значений табличной разности, не превосходящей 4 единиц, третий источник погрешности сколько-нибудь заметного влияния не имеет, и общая погрешность результата линейной интерполяции лишь в исключительно редких случаях может немного превзойти единицу разряда последней цифры. Это заключение легко проверяется на опыте. Например, пользуясь таблицей квадратов, находим, применяя линейную интерполяцию, квадраты чисел, приведенных ниже в первой строке, и пишем их во второй строке, а в третьей строке помещаем соответствующие точные квадраты, округленные до четвертого десятичного знака, в четвертой же — разности чисел второй и третьей строк, выраженные в сотых долях.
x |
8,051 | 8,052 | 8,053 | 8,054 | 8,055 | 8,056 | 8,057 | 8,058 | 8,059 |
у1≈х2 | 64,82 | 64,83 | 64,85 | 64,86 | 64,88 | 64,90 | 64,91 | 64,93 | 64,94 |
y1≈х2 | 64,8186 | 64,8374 | 64,8508 | 64,8669 | 64,8830 | 64,8991 | 64,9152 | 64,9314 | 64,9475 |
(y1 — y2) · 100 | 0,14 | -0,47 | -0,08 | -0,69 | -0,30 | -0,09 | -0,52 | -0,14 | -0,75 |
Как видим, погрешности интерполированных значений нигде не превосходят единицы разряда последней цифры.
Чтобы облегчить выполнение линейной интерполяции, в большинстве таблиц настоящего сборника даны «готовые поправки» в столбцах справа, набранные курсивом. Если табличные разности мало меняются на протяжении целой строки, то поправки по формуле v = du/h можно вычислить для всех чисел строки. Например, для строки 8,0 таблицы квадратов №3 поправка на 0,001 в начале строки равна (8,012 — 8,002):10 = 0,01601, или 1,601 (сотых), а на конце ее (8,102 — 8,092): 10 = 0,01619, или 1,619, а в среднем 1,610 (сотых). Умножая эту среднюю поправку на числа от 1 до 9, получаем 1,61; 3,22; 4,83; 6,44; 8,05; 9,66; 11,27; 12,88; 14,49 или после округления до целых 2; 3; 5; 6; 8; 10; 11; 13; 14.
Именно эти числа и приведены на строке 8,0 таблицы квадратов справа (набраны курсивом). Как показывает опыт, применение этих готовых поправок сберегает до 50% времени, затрачиваемого на работу с таблицами.
Если табличные разности на протяжении строки меняются более заметно, готовые поправки приходится вычислять для частей строки, как это сделано, например, в таблице 9 для строк 73°, 74°, 75° или для нескольких первых строк таблицы 13 мантисс логарифмов. Если изменение табличных разностей на протяжении строки выражено еще резче, от готовых поправок приходится отказаться. В таких случаях операцию линейной интерполяции приходится проводить полностью, находя h, d, u, v = du/h, как, например, в таблице 15 и нескольких других.
При большой табличной разности поправку следует вводить не только на первую цифру избытка, но и на вторую, если она имеется, уменьшая приведенные в таблице готовые поправки в 10 раз. Так, чтобы найти 2,93452, в таблице 3 берется 2,932 = 8,585 и прибавляется поправка на 4 тысячных, равная 24 (тысячным), а затем поправка на 5 десятитысячных, равная 29:10 = 3 (тысячным), и получается окончательно 8,612 (непосредственное умножение дает 8,61129…).
Как мы уже видели, если избыток данного значения аргумента больше половины ступени, выгоднее пользоваться ближайшим большим значением функции, отнимая от него поправку на недостаток данного значения аргумента по сравнению с ближайшим большим его значением. Поэтому во всех таблицах, где аргументом служит угол и где ступень равна 6′, готовые поправки даны только на 1′, 2′, 3′. Если избыток составляет 4′ или 5′, надо брать поправку на 2′ или 1′, вычитая ее из ближайшего большего значения функции. Кроме экономии места, занимаемого таблицей, это дает некоторый выигрыш в точности получаемых результатов, так как малые поправки точнее больших.
Необходимо решительно предостеречь от применения линейной интерполяции в случае резко неравномерного изменения функции. Всякий раз, когда готовые поправки не даны, а нужно интерполировать, следует выяснить, насколько равномерен ход функции, и применять линейную интерполяцию лишь в том случае, когда соседние табличные разности мало отличаются друг от друга (не больше чем на 4 единицы), а в противном случае искать других путей. Так, например, желая найти lg sin 1°04’36», берем в таблице 15, где готовых поправок нет, lg sin 1°04′ — 2,2699, lg sin 1°05′ = 2,2766, lg sin 1°06′ = 2,2832 и убеждаемся, что линейная интерполяция здесь допустима, так как табличные разности равны 67 и 66. Вычисляя v = 67·36/60 = 40,2 = 40 и прибавляя эту поправку к табличному логарифму 2,2699, получаем lg sin 1°04’36» = 2,2739 (по семизначным таблицам получается 2,2739331). Но если надо получить lg sin 0°05’30» и мы применим тот же способ линейной интерполяции, то получим lg sin 0°05′ = 3,1627, d = 792, h =60″, u = 30″, v = 792/60 · 30 = 396, lg sin 0°05’30» =60= 3,2023, в то время как более точное значение этого логарифма, найденное по семизначным таблицам, есть 3,2040886. Недопустимо большая погрешность нашего результата обусловлена резко неравномерным изменением функции: рядом с использованной нами табличной разностью 792 находится разность 669, линейная интерполяция здесь недопустима. Здесь можно воспользоваться тем обстоятельством, что при очень малых углах синус весьма мало отличается от радианной меры (меньше чем на шестую часть куба этой радианной меры). В таблице 11 берем радианную меру угла в 5′, равную 0,0014544, а также угла в 30″, равную 0,00014544, и, складывая, получаем число 0,0015998, представляющее собой приближенное значение sin 0°05’30» с 7 точными десятичными знаками. Найдя по таблице 13 его логарифм, получаем 3,2041, т. е. как раз то, что надо.
Во многих случаях таблицы дают непосредственно значения функции лишь в одном ограниченном интервале значений аргумента, но путем несложных дополнительных расчетов, производимых обычно в уме, можно существенно расширить этот интервал. Так обстоит дело с таблицами квадратов, кубов, обратных значений и ряда других. Возьмем, например, таблицу 7, дающую непосредственно значения площади круга с диаметром от d = 1 до d = 10, замечая, что при увеличении диаметра круга в 10 раз его площадь увеличивается в 102 = 100 раз, мы можем по этой же таблице находить площадь круга любого диаметра. Например, желая найти площадь круга диаметра d = 49,52, находим по таблице сперва площадь круга диаметра 4,952 (строка 49, столбец 5, поправка на 2), равную 19,26, а затем увеличиваем этот результат в 100 раз и получаем окончательно 1926. Чтобы найти площадь круга с диаметром d = 0,04567, получаем сперва площадь круга диаметра 4,567 (строка 45, столбец 7, вычитается поправка на 3), равную 16,38, потом уменьшаем ее в 1002 = 10000 раз и получаем 0,001638.
Разобрав во всех деталях вопрос о разыскании посредством таблиц значения функции по данному значению ее аргумента, то есть так называемый «прямой вопрос», переходим к «обратному вопросу», когда по данному значению той функции, для которой таблица составлена, надо найти соответствующее значение аргумента.
Если данное значение функции имеется в таблице, все дело сводится к выписыванию соответствующего значения аргумента. Если же данного значения функции в таблице нет, то пользуются той же операцией линейной интерполяции, внеся в нее надлежащие изменения и предварительно убедившись в ее допустимости. Берут ближайшее меньшее табличное значение функции и находят, сколько надо добавить к соответствующему значению аргумента, чтобы довести это ближайшее меньшее значение функции до данного. Здесь используется та же пропорция u:v = h:d, что и раньше, с той лишь разницей, что теперь у дано, и ищем а по формуле u = hv/d. Так, чтобы найти с помощью таблицы квадратов число, квадрат которого равен 4,235, т. е. квадратный корень из числа 4,235, берут ближайший меньший и ближайший больший табличные квадраты 4,203 = 2,052 и 4,244 = 2,062. Здесь ступень h = 10 (тысячным), табличная разность d = 41 (тысячной), следующая табличная разность тоже 41, линейная интерполяция допустима. Чтобы довести ближайшее меньшее табличное значение до данного, надо увеличить это ближайшее на 4,235 — 4,203 = 0,032, откуда v = 32 (тысячным). Поэтому u = 10·32/41 =8 и искомый корень равен 2,050 + 0,008 = 2,058. Можно взять не ближайшее меньшее, а ближайшее большее значение функции и уменьшать его до данного, выясняя, каково соответствующее уменьшение ближайшего большего значения аргумента. В данном примере соответственно этому берем 4,244 — 4,235 = 0,009, т. е. и = 9 (тысячным), и находим u = 10·9/41 ≈ 2, а затем искомый корень 2,060 — 0,002 = 2,058. Вообще, лучше пользоваться тем из ближайших табличных значений функции, какое ближе к искомому.
Применение готовых поправок и здесь существенно облегчает работу: найдя разность между данным значением функции и ближайшим табличным ее значением (меньшим или большим), смотрим, какая поправка из напечатанных курсивом на той же строке ближе всего к этой разности, и берем цифру, находящуюся в заголовке соответствующего столбца. Для получения квадратного корня из числа 4,235 достаточно заметить, что это число отличается от ближайшего меньшего табличного квадрата на 32 (тысячных) и что среди поправок, напечатанных на этой же строке, ближайшей к этому числу 32 является 33. Прибавив к соответствующему табличному значению аргумента 2,05 число 8 (тысячных), взятое из заголовка этого столбца поправок, получаем окончательно 2,05 + 0,008 = 2,058. Если взять ближайшее большее значение функции (4,244), то получается разность 4,244 — 4,235 = 0,009. В столбцах поправок находим ближайшую цифру 8 в столбце 2 и выполняем вычитание 2,06 — 0,002, приводящее к тому же результату 2,058.
Вопрос о точности, с какой обратная линейная интерполяция дает искомое значение функции, довольно сложен. Оказывается, что здесь возможны самые различные случаи и что результат здесь тем более точен, чем больше табличная разность (предполагается, что линейная интерполяция допустима). Например, если дано приближенное значение sin А = 0,9997 с 4 точными десятичными знаками, то в таблице 8 мы находим целых три угла с таким синусом (88°30′, 88°36′, 88°42′). Полагая А = 88°36′, надо иметь в виду, что это значение искомого угла весьма неточно: оно может отличаться от точного до 9′. Если же sin A = 0,1070, то находимое по таблице 8 с помощью готовых поправок значение 6°08′ отличается от точного, как можно показать, не больше чем на 1′: применение способа границ приводит к заключению, что 6°08′
Итак, каждая таблица служит не только для получения значений той функции, для которой она составлена, но и для получения значений аргумента, т. е. для получения значений обратной функции: по таблице квадратов можно находить и квадратные корни, по таблице логарифмов — антилогарифмы и т. д. Однако опыт показывает, что решение обратного вопроса требует несколько большей затраты труда, чем решение прямого, а потому в настоящем сборнике наряду с таблицей логарифмов помещена таблица антилогарифмов, наряду с таблицей квадратов — таблица квадратных корней, хотя можно было бы обойтись и без них.
До сих пор речь шла только о таблицах возрастающих функций. Легко видеть, как изменяется способ пользования таблицей, если функция убывает, как например, в таблице 2, дающей значения дробей вида 1/и, или в таблице 8 при разыскании косинусов. При работе с таблицей возрастающей функции ошибки от недостаточного внимания случаются реже, а потому можно рекомендовать заменять подыскание косинусов подысканием синусов дополнительных углов, подыскание котангенсов — подысканием тангенсов дополнительных углов.
Таблицы настоящего сборника, вообще говоря, обеспечивают получение искомых значений с 4, иногда с 5 значащими цифрами. Но бывают особо неблагоприятные случаи вычислений (вычитание из приближенного числа другого приближенного числа, близкого к первому, возведение приближенного числа в степень с большим показателем и т. д.), когда окончательный результат получается с меньшей точностью. Если точность результата требуется большая, надо либо обратиться к более подробной таблице (пятизначной, семизначной и т.д.), либо проводить вычисление непосредственно, что не представляет непреодолимых трудностей при возведении в степень, извлечении корня и некоторых других операциях. Ниже приведены некоторые «ряды», позволяющие находить с произвольно высокой точностью значения логарифмов, антилогарифмов, синусов, косинусов, тангенсов, корней квадратных и кубических.