
- •Численные методы и универсальные математические пакеты (Номер занятия в оглавлении равен номеру лабораторной работы. Лабораторные работы состоят из практических заданий из соответствующих занятий)
- •Введение
- •Занятие 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.
Рассмотрим
функцию
.
Используя составную формулу трапеций
с выборкой из 11 точек, вычислить
приближенное значение интеграла от
на интервале
.
План:
Задать фукнцию через оператор и обозначить ее f
Задать концы интервалов как переменные a, b
Задать количество интервалов M:=10 (если точек 11, то интервалов будет 10)
Посчитать шаг, обозначить h (
)
Создать с помощью операции $ или seq список значений x (
)
Посчитать значение интеграла по формуле (1)
Посчитать значение интеграла с помощью функции int
Практическое задание 2.
Рассмотрим функцию . Используя составную формулу Симпсона с выборкой из 11 точек, вычислить приближенное значение интеграла от на интервале .
План:
Задать фукнцию через оператор и обозначить ее f
Задать концы интервалов как переменные a, b
Задать количество интервалов M:=5 (так как число интервалов равно
)
Посчитать шаг, обозначить h (
)
Создать с помощью операции $ или seq список значений x (
)
Посчитать значение интеграла по формуле (2) (необходимо учитывать то, что нумерация в списке начинается с 1, а нумеруются с 0)
Практическое задание 3.
Выполнить задание 2 используя для вычисления интеграла по формуле (1) оператор цикла.
План:
Задать фукнцию через оператор и обозначить ее f
Задать концы интервалов как переменные a, b
Задать количество интервалов M:=10 (если точек 11, то интервалов будет 10)
Посчитать шаг, обозначить h ( )
Создать x[0]:=a;
Задать начальное значение суммы равное s:=0
Задать цикл i от 1 до M. В теле цикла посчитать x[i] (
) и s (
)
Занятие 6. Решение нелинейных уравнений
Фундаментальным принципом компьютерной науки является итерация, т.е. процесс повторяется до тех пор, пока не будет получен ответ. Техника итераций используется для нахождения корней уравнений, решения систем линейных и нелинейных уравнений и решения дифференциальных уравнений.
Метод простой итерации
Дано уравнение с одной независимой переменной
(1)
где
‑ некоторая заданная функция,
переменная
изменяется на отрезке
В методе
простой итерации члены приближающей
последовательности
находятся непосредственно из уравнения
(1)
,
(2)
Т.е. последовательность имеет такой вид:
Данный
метод будет сходиться к корню уравнения
,
когда
на интервале
Для остановки итерационного процесса используют следующий критерий
Итерационный
процесс выполняется до тех пор, пока
,
где
‑
заданная точность.
Метод бисекции (деления пополам)
Дана
функция
,
необходимо найти ее единственный корень
на промежутке
.
Если корень на этом промежутке существует, то должно выполняться неравенство
(т.е. они должны иметь разные знаки).
Так
как функции
непрерывная, то ее график пересекает
ось
в точке
‑ решение уравнения. Метод деления
пополам сдвигает крайние точки все
ближе, пока на интервале не получится
произвольно малый отрезок, содержащий
корень. Решающим шагом процесса деления
интервала пополам являетя выбор средней
точки
и анализ трех возможных ситуаций,
которые могут возникнуть.
если
, то корень содержится на интервале
(т.е. принимаестя
)
если
, то корень содержится на интервале
(т.е. принимаестя
)
если
, то корнем является точка
Данный процесс итераций прекращается, когда длина интервала станет меньше , т.е.
.