Дьяконов Владимир Павлович
Шрифт:
Однобуквенные имена этих функций отождествляются с первой буквой в наименовании ортогональных полиномов. Вопреки принятым в Maple правилам, большие буквы в названиях этих полиномов не указывают на инертность данных функций — все они являются немедленно вычисляемыми. В данном разделе функции этого пакета будут полностью описаны.
Отметим определения указанных функций:
G(n,a,x) — полином Гегенбауэра (из семейства ультрасферических полиномов);
H(n,x) — полином Эрмита;
L(n,x) — полином Лагерра;
L(n,a,x) — обобщенный полином Лагерра;
P(n,x) — полином Лежандра;
P(n,a,b,x) — полином Якоби;
T(n,x) — обобщенный полином Чебышева первого рода;
U(n,x) — обобщенный полином Чебышева второго рода.
Свойства ортогональных многочленов хорошо известны. Все они характеризуются целочисленным порядком n, аргументом х и иногда дополнительными параметрами а и b. Существуют простые рекуррентные формулы, позволяющие найти полином n-го порядка по значению полинома (n-1)-го порядка. Эти формулы и используются для вычисления полиномов высшего порядка.
5.4.2. Вычисление ортогональных полиномов
Ниже представлены примеры вычисления ортогональных полиномов (файл orthpol):
В отличие от ряда элементарных функций, ортогональные многочлены определены только для действительного аргумента х. При комплексном аргументе ранее результат просто повторял исходное выражение с многочленом:
Но уже в Maple 9 ортогональные полиномы с комплексными аргументами могут вычисляться:
Ортогональные многочлены не определены и для дробного показателя n. Впрочем, надо отметить, что такие многочлены на практике используются крайне редко.
5.4.3. Построение графиков ортогональных полиномов
Представляет интерес построение графиков ортогональных многочленов. На рис. 5.6 построены графики ряда многочленов Гегенбауэра и Эрмита. На рис. 5.7 построены графики ортогональных многочленов Лагерра и Лежандра. Наконец на рис. 5.8 даны графики ортогональных многочленов Чебышева T(n, х) и U(n, x).
Рис. 5.6. Графики ортогональных многочленов Гегенбауэра и Эрмита
Рис. 5.7. Графики ортогональных многочленов Лагерра и Лежандра
Рис. 5.8. Графики ортогональных многочленов Чебышева
Приведенные графики дают начальное представление о поведении ортогональных многочленов. К примеру, многочлены Чебышева имеют минимальное отклонение от оси абсцисс в заданном интервале изменения х. Это их свойство объясняет полезное применение таких многочленов при решении задач аппроксимации функций, которые рассматриваются в этой главе далее. Можно порекомендовать читателю по их образцу и подобию построить графики ортогональных многочленов при других значения параметра n и диапазонах изменения аргумента х.
5.4.4. Работа с рядами ортогональных многочленов
Для работы с рядами ортогональных многочленов имеется пакет OrthogonalSeries для работы с рядами ортогональных многочленов. Он имеет довольно представительный набор функций: