
- •Численные методы и универсальные математические пакеты (Номер занятия в оглавлении равен номеру лабораторной работы. Лабораторные работы состоят из практических заданий из соответствующих занятий)
- •Введение
- •Занятие 1. Теория погрешностей Источники и классификация погрешностей.
- •Приближенные числа.
- •2. Погрешность произведения.
- •3. Погрешность частного.
- •Введение в Maple
- •Константы
- •Переменные, неизвестные и выражения
- •Команды преобразования выражений
- •Упрощение выражения: simplify()
- •Раскрытие скобок в выражении: expand()
- •Разложение полинома на множители: factor()
- •Ограничения на неизвестные: assume()
- •Сложные типы данных
- •Последовательность выражений
- •Списки и множества
- •Внутренняя структура выражений
- •Подстановка и преобразование выражений
- •Практическое задание
- •Занятие 2. Приближение функций
- •Приближение полиномом Тейлора
- •Приближение функий заданных таблицей своих значений
- •Алгебраическое интерполирование
- •Интерполяционный многочлен Лагранжа
- •Многочлены Чебышева
- •Разделенные разности и многочлен Ньютона
- •Решение уравнений, неравенств и их систем
- •Команда solve()
- •Команда fsolve()
- •Решение неравенств
- •Дифференцирование и интегрирование
- •Int(выражение, переменная);
- •Практическое задание 1
- •Практическое задание 2
- •Занятие 3. Построение кривой по точкам
- •Линия, построенная методом наименьших квадратов
- •Метод линеаризации данных для экспоненциальной кривой
- •Нелинейный метод наименьших квадратов
- •Линейный метод наименьших квадратов
- •Интерполирование сплайнами
- •Графика в Maple
- •Команда двумерной графики plot()
- •Двумерные команды пакета plots
- •Несколько советов
- •Пространственная графика, команда plot3d()
- •Трехмерные команды пакета plots
- •Практическое задание 1.
- •Практическое задание 2.
- •Практическое задание 3.
- •Практическое задание 4.
- •Занятие 4. Численное дифференцирование Приближение производной
- •Анимация
- •Двумерная анимация
- •Трехмерная анимация
- •Практическое задание 1
- •Практическое задание 2
- •Занятие 5. Численное интегрирование
- •Квадратурные формулы
- •Основы программирования в Maple
- •If булево_выражение then последовательность_операторов
- •Практическое задание 1.
- •Практическое задание 2.
- •Практическое задание 3.
- •Занятие 6. Решение нелинейных уравнений
- •Метод простой итерации
- •Метод бисекции (деления пополам)
- •Метод Ньютона
- •Процедуры в Maple
- •Практическое задание 1.
- •Практическое задание 2.
- •Практическое задание 3.
- •Практическое задание 4.
- •Занятие 7. Решение систем Линейные системы. Метод Гаусса
- •Пакет LinearAlgebra
- •Основные типы данных
- •Элементарные операции с матрицами и векторами
- •Решение систем линейных уравнений
- •Практическое задание 1.
- •Практическое задание 2
- •Практическое задание 3
- •Занятие 8. Решение дифференциальных уравнений
- •Задача Коши
- •Метод Эйлера
- •Методы Рунге-Кутта
- •Системы дифференциальных уравнений
- •Решение обыкновенных дифференцильных уравнений в Maple
- •Практическое задание 1
- •Практическое задание 2.
- •Практическое задание 3.
Подстановка и преобразование выражений
При выполнении математических преобразований часто необходимо произвести замену переменных в выражении, функции, уравнении и т. д. Для этих целей в Maple существует команда subs(), синтаксис которой имеет следующий вид:
subs(старое_выражение=новое_выражение, выражение);
Пример. Подстановки в выражении
> ex:=x^2+x^(1/3);
> subs(x=cos(x),ex);
> subs(x=1,ex);
Иногда необходимо выполнить преобразование выражения одного типа в другой тип. Например, нельзя строить график ряда Тейлора какой-либо функции, но всегда можно построить график полинома. Следовательно, следует выражение, имеющее тип series (ряд), преобразовать в выражение, имеющее тип polynom (полином). Преобразовать выражение в другой тип можно командой convert(), первым параметром которой задается выражение, а вторым — тип, в который это выражение следует преобразовать.
Практическое задание
Дано
уравнение
Его корни определяются по формулам
(1)
Для данного примера формула для вычисления первого корня имеет большую погрешность, т.к. в числителе стоит разность близких чисел. Поэтому для вычисления этого корня необходимо избавится от этой разности:
(2)
В Maple необходимо сделать следующее:
Задать количество значащих цифр 4
Записать коэффициенты уравнения как переменные a, b, c.
Записать уравнение и обозначить его переменной eq.
Отдельно посчитать корень из дескреминанта и обозначить его переменной d (получить значение в виде числа с плавающей точкой).
Посчитать первый корень по формуле (1) и запомнить его в переменную x11 (получить его десятичное представление).
Подставить полученное значение x11 в уравнение eq и запомнить полученное значение как y11.
Посчитать первый корень по формуле (2) и запомнить его в переменную x12 (получить его десятичное представление).
Подставить полученное значение x12 в уравнение eq и запомнить полученное значение как y12.
Проанализировать полученные результаты y11 и y12.
Проделать пункты 2-9 для 8 цифр в мантиссе.
Проанализировать полученные результаты.
Занятие 2. Приближение функций
Для
вычисления стандартных функций, таких
как
,
,
и
,
в пакетах прикладных программ используются
вычислительные процедуры, включающие
в себя приближение полиномами.
Приближение полиномом Тейлора
Предположим, что функция
имеет непрерывные производные до
порядка на
и
некоторая точка из этого промежутка.
Тогда
Недостатком этого приближение является необходимость вычисления функции и ее производных в точке , которые не всегда известны.
Приближение функий заданных таблицей своих значений
Простейшая
задача, приводящая к приближению
функций, заключается в следующем. В
дискретные моменты времени
наблюдаются значения фунции
;
требуется восстановить ее значения
при других
.
Иногда из каких-то дополнительных соображений известно, что приближающую функцию следует искать в виде
Если
параметры
определяются из условия совпадения
функции
и приближающей функции
в точках
,
которые называются узлами интерполяции
то такой способ приближения называют интерполяцией или интерполированием.
Если
точка
,
в которой необходимо вычислить значение
функции
,
лежит вне отрезка
,
где
и
соответственно минимальное и максимальное
значение узлов интерполяции, то такая
операция называется экстраполяцией.
Способы
интерполирования отличаются выбором
вида аппроксимирующей функции
.