- •Численные методы и универсальные математические пакеты (Номер занятия в оглавлении равен номеру лабораторной работы. Лабораторные работы состоят из практических заданий из соответствующих занятий)
- •Введение
- •Занятие 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.
Практическое задание 1
Разложить
в ряд Тейлора функцию
в окрестности точки
,
построить графики этих функций.
Для построения ряда Тейлора в Maple существует фукнция
taylor(функция, x=a, n);
где x – имя независимой переменной, от которой зависит функция, a – точка в районе, которой ищется приближение , n – количество членов ряда.
Результатом выполнения данной команды служит разложение функции в ряд Тейлора с остаточным членом. Для дальнейшей работы с этим разложение остаточный член будет только мешать, его можно отбросить командой сonvert.
В Maple необходимо сделать следующее:
задать функцию f
задать значение x0 и n=5
разложить функцию в ряд Тейлора, отбросить остаточный член и обозначить полученное разложение переменной g.
построить обе функции на графике в окрестности точки x0.
Практическое задание 2
Произвести интерполяцию функции
на промежутке [-1,1]
с 11 узлами
с помощью алгебраического полинома, полинома Лагранжа, полинома Ньютона.
Построить полином Лагранжа по 11 узлам Чебышева.
Построить функцию и полученные интерполяционные полиномы на графике.
Т.е. необходимо:
задать функцию f как функцию, зависящую от x
задать начало и конец промежутка как a, b и количество узлов n
сформировать список узлов с помощью команды seq (формула для нахождения узлов:
).используя команду map найти значения функции (последовательность значений
)
для соответствующих значениях xиспользуя пакет CurveFitting и команду PolynomialInterpolation построить три интерполяционных полинома. Обозначить их P (алгебраический), L (Лагранжа) и l (Ньютона).
сформировать список узлов Чебышева с помощью команды seq (формула для нахождения узлов:
).используя команду map найти значения функции (последовательность значений
)
для соответствующих значениях
используя пакет CurveFitting и команду PolynomialInterpolation построить интерполяционный полином Лагранжа по новым узлам, обозначить его LC.
Построить все полиномы (P, L, l, LC) саму фунцию f(x) на графике.
Занятие 3. Построение кривой по точкам
Техника численного анализа часто применяется для вычерчивания графиков экспериментальных данных.
Линия, построенная методом наименьших квадратов
Пусть
‑ совокупность
точек с различными абсциссами
.
Линия построенная, методом наименьших квадратов
коэффициетны находятся из условия минимума расстояния
Коэффициенты
линии наименьших квадратов
являются решениями следующей системы
линейных уравнений, известной под
названием нормальные уравнения:
Метод линеаризации данных для экспоненциальной кривой
Предположим, что заданы точки и требуется произвести подгонку экспоненциальной кривой
Прологорифмируем обе части равенства:
Заменим переменные и получим линейное соотношение между новыми переменными:
где
.
Исходные
точки
преобразовались в точки
,
этот процесс называется линеаризацией
данных. Тогда построенная методом
наименьших квадратов линия
является подгонкой к точкам
.
После
нахождения параметров
и
,
вычисляется параметр
:
Техника метода линеаризации данных применяется для подгонки и других кривых, например
,
.
