Чтобы установить различие в поведении физической характеристики, появляющееся при переходах от одного формализма (взаимозависимость Re ? и Im ?) к другому формализму (независимость Re ? и Im ?), вычислим и отобразим поведение ?, используя не формулу (3.15), а при тех же константах выражение:
где Т = Т0cos(kr – ?t). При малости аргументов ? тригонометрические функции f(?) можно было бы разложить по степеням ? с точностью до второго порядка. Но тогда для r получится уравнение шестого порядка, которое в радикалах не разрешается. Математические силовики Античности – Пифагор, Евклид, Евдокс, Платон, а также, возможно, философы Софокл, Прокл, Аристокл, Эмпедокл и, наверняка, Пиррон могли бы решить эту задачу. Ведь это древние придумали проблему квадратуры круга и разминку с делением угла на три равные части (трисекция угла). Но «мы пойдем своим путем», а именно: в этом случае после построения окружности Уиллера и ее фиксации для поиска зависимости r = r(?) применим метод последовательных приближений. Для этого нужно создать программу поиска функции r(?), обеспечивающую быструю сходимость и хорошую точность вычислений. Сначала, при заданных константах, на плоскости E(r, ?) определим рельеф формулы
где A = Далее нужно составить таблицу, состоящую из подтаблиц в соответствии со связностью среза. Это тоже делает программа. Теперь можно табличную функцию ri = r(?i), где i – номер «точки», преобразовать в аналитическую функцию, представленную в виде полинома по степеням ?. Этот способ называется интерполированием табличной зависимости. Высшая степень интерполяционного многочлена определяется визуально с целью сэкономить время вычислений при заданных разрешимости экрана и допустимых ошибках вычислений. На этом шаге может быть применен метод наименьших квадратов. Суть его в том, что искомые коэффициенты многочлена подбираются так, чтобы минимизировать суммарное отклонение строящейся плавной кривой от табличных «точек». Это достигается варьированием полученного (дискретного) уравнения по искомым коэффициентам (см. примеры в [23]). На этом шаге необходимо вычислить значения коэффициентов, решая систему линейных уравнений, то есть придется составлять матрицы — 24 —
|