Maple 9.5/10 в математике, физике и образовании

Рис. 5.31. Пример экспоненциальной регрессии

Функция нелинейной регрессии входит в новейший пакет оптимизации Optimization, введенный в Maple 9.5, и описанный в следующей главе. Кроме того, на Интернет-сайте корпорации Waterloo Maple можно найти файлы simplenl.mws и gennlr.mws с процедурами и примерами линейной и нелинейной регрессий общего вида. Интересная реализация нелинейной регрессии для кусочной функции дается в файле nonelinearpiecewise.mws.

5.12.6. Сплайновая регрессия с помощью функции BSplineCurve

Функция BSplineCurve из пакета CurveFitting может использоваться для реализации сплайновой регрессии. Пример этого представлен на рис. 5.32. Опция order задает порядок B-сплайнов, который на 1 меньше заданного целого значения.

Рис. 5.32. Пример выполнения сплайновой регрессии В-сплайнами

Функция BsplineCurve выглядит несколько недоделанной. Так, при order=3 и 4 кривая регрессии не дотягивает до концевых точек, а при установки order=1 все точки соединяются отрезками прямых — в том числе концевые. Так что использовать эту функцию для экстраполяции нельзя.

5.13. Работа с функциями двух переменных

5.13.1. Maplet-инструмент для работы с функциями двух переменных

Для эффектной демонстрации работы с функциями многих переменных в состав пакета Student системы Maple 9.5 введен новый подпакет MultivariateCalculus. Его примеры можно запускать как с командной строки, так и из позиции Tools меню в стандартном варианте интерфейса — Tutors→Calculus→Calculus-Multi-Variables.

Approximate Integration… — открывает Maplet-окно аппроксимации двойных интегралов;

Cross Section… — открывает Maplet-окно демонстрации сечения поверхности;

Directional Derivatives… — открывает Maplet-окно вычисления производных в заданном направлении;

Gradient… — открывает Maplet-окно вычисления градиента;

Taylor Series… — открывает Maplet-окно разложения функций в ряд Тейлора.

Представленные средства носят учебный характер — не случайно они входят в пакет Student. Реально визуализация возможна только для функций двух переменных.

5.13.2. Демонстрация разложения в ряд Тейлора функции двух переменных

Команда Taylor Series… — открывает Maplet-окно разложения функции двух переменных z(х, у) в ряд Тейлора относительно заданной точки (х0, у0). Это окно представлено на рис. 5.33.

Рис. 5.33. Maplet-окно демонстрации разложения в ряд Тейлора функции двух переменных

В данном окне дан пример разложения в ряд Тейлора функции sin(x*y) в окрестности точки (0, 0) в интервале изменения х[-2, 2], у[-2, 2] и z[-1, 1]. Установки в окне совершенно очевидны. Графики в правой части представляют поверхность, описываемую исходной функцией и поверхность, представленную рядом Тейлора. Кнопка Display начинает построение графиков, кнопка Animation позволяет наблюдать анимацию разложения, а кнопка Close закрывает окно и переносит рисунок в текущий документ системы Maple 9 5.

5.13.3. Демонстрация вычисления градиента функции двух переменных

Команда Gradient… — открывает Maplet-окно демонстрации вычисления градиента функции двух переменных z(x, у) в ряд Тейлора относительно заданной точки (х0, y0). Это окно представлено на рис. 5.34.

Рис. 5.34. Maplet-окно демонстрации вычисления градиента функции двух переменных

Работа с этим окном практически не отличается от описанной для примера с рядом Тейлора. Единственное исключение — новая кнопка Gradient Field Plot. Она позволяет строить график поля градиента с помощью стрелок. Этот случай представлен на рис. 5.35.

Рис. 5.35. Maplet-окно демонстрации вычисления градиента функции двух переменных с графиком поля градиента

5.13.4. Демонстрация вычисления производной в заданном направлении

Команда Directional Derivatives… — открывает Maplet-окно демонстрации вычисления производных функции двух переменных z(х, у) в заданном направлении, указанном точкой с координатами (х, у). Это окно представлено на рис. 5.36.

Рис. 5.36. Maplet-окно демонстрации вычисления градиента функции двух переменных

Работа с этим окном практически не отличается от описанной для предшествующих примеров.

5.13.5. Демонстрация приближенного вычисления интеграла

Команда Approximate Integration… — открывает Maplet-окно демонстрации вычисления двойных интегралов с подынтегральной функцией двух переменных z(х, у). Это окно представлено на рис. 5.37.

Рис. 5.37. Maplet-окно демонстрации приближенного вычисления двойного интеграла в прямоугольной системе координат

Для вычисления интеграла нужно задать подынтегральную функцию и пределы по переменным x и у. Для построения графика можно также задать пределы по переменной z. Приближенное значение интеграла вычисляется суммированием объёмов прямоугольных столбиков, на которые разбивается пространство под поверхностью z(x, у). Число разбиений устанавливается списком Partition. Можно задать один из четырех методов расположения столбиков. В области Value отображается точное и приближенное (сумма объемов столбиков) значения интеграла.

Возможно представление интеграла и в полярной системе координат. Пример этого дан на рис. 5.38.

Рис. 5.38. Maplet-окно демонстрации приближенного вычисления двойного интеграла в полярной системе координат

5.13.6. Маплет-демонстрация сечения поверхности

Команда Cross Section… открывает Maplet-окно демонстрации сечения поверхности плоскостями. Поверхность задается функцией двух переменных z(x, у). Окно этой команды представлено на рис. 5.39.

Рис. 5.39. Maplet-окно демонстрации сечения поверхности параллельными плоскостями

Работа с этим окном вполне очевидно. На рисунке в левой части окна строится исходная поверхность, секущие плоскости и линии их пересечения.

Глава 6 Решение задач линейной алгебры, оптимизации и регрессии

Задачи линейной алгебры, оптимизации и регрессии — одни из самых массовых в науке, технике и образовании [37, 39–46]. Им и посвящена эта глава. В ней даны основные определения линейной алгебры, основы работы с массивами, векторами и матрицами, функции для работы с векторами и матрицами и для решения систем линейных уравнений. Дано описание средств оптимизации, в том числе новейших системы Maple 10.