- •220201 - «Управление и информатика в технических системах» Коломна, 2013
- •Численные методы Методические указания к лабораторным работам
- •220201 - «Управление и информатика в технических системах» Коломна , 2013
- •Содержание
- •Введение
- •Лабораторная работа №1 Организация Windows-приложения в Delphi с использованием визуальных компонентов классов: tEdit, tLabel, tMemo, tButton. Программирование интерполяционных алгоритмов Цель работы
- •1. 1. Краткие сведения из теории
- •1.1.1. Понятия аппроксимации и интерполяции
- •1.1.2. Вычисление многочленов по схеме Горнера
- •1.2. Выполнение работы Исходные данные для выполнения работы
- •1.3. Контрольные вопросы
- •2.1. Краткие сведения из теории
- •2.1.1. Метод прямоугольников
- •2.1.2. Метод трапеций
- •2.1.3. Метод Симпсона
- •2.1.4. Алгоритм автоматического выбора шага
- •2.2. Выполнение работы Порядок выполнения работы
- •Исходные данные для выполнения работы
- •2.3. Контрольные вопросы
- •3.1. Краткие сведения из теории
- •3.1.1. Обзор методов решения систем решения линейных алгебраических уравнений
- •3.1.2. Метод Гаусса
- •Вычисление определителя
- •Прямой ход метода Гаусса
- •3.1.3. Метод Гаусса-Зейделя
- •3.2. Выполнение работы Порядок выполнения работы
- •Исходные данные для выполнения работы
- •3.3. Контрольные вопросы
- •Лабораторная работа №4 Использование графики в Delphi. Организация меню в приложении. Программирование итерационных алгоритмов решения нелинейных уравнений Цель работы
- •4.1. Краткие сведения из теории
- •4.1.1. Метод бисекции
- •4.1.2. Метод хорд
- •4.1.3. Метод касательных
- •4.1.4. Метод простой итерации
- •4.2. Выполнение работы Порядок выполнения работы
- •Исходные данные для выполнения работы
- •4.3. Контрольные вопросы
- •Лабораторная работа №5 Организация многооконного приложения в Delphi для решения о.Д.У. Одношаговыми разностными методами. Программная работа с файлами данных Цель работы:
- •5.1. Краткие сведения из теории
- •5.1.1. Метод Эйлера
- •5.1.2. Метод Эйлера с пересчётом
- •5.1.3. Метод Рунге-Кутта
- •5.1.4. Решение дифференциальных уравнений высшего порядка
- •5.2. Выполнение работы Порядок выполнения работы
- •Исходные данные для выполнения работы
- •5.3. Контрольные вопросы
- •Справка по работе с основным набором визуальных компонентов Delphi, требуемых при выполнении лабораторных работ
- •Компонент Окно выбора файла (tOpenDialog)
- •Компонент Окно сохранения файла (tSaveDialog)
- •Список использованной литературы
4.1.1. Метод бисекции
В
методе бисекции (деления отрезка пополам)
в качестве начального приближения корня
принимают середину отрезка, содержащего
корень. Затем исследуют функцию
на знак на концах двух полученных
отрезков. Отрезок, на концах которого
принимает значения разных знаков,
содержит искомый корень. Его принимают
в качестве нового отрезка исследования,
второй отрезок отбрасывают. В качестве
приближённого значения корня первой
итерации
принимают середину нового отрезка.
Снова
исследуют
на знак на концах нового отрезка и т.д.
После каждой итерации отрезок, содержащий
корень, уменьшается вдвое. Это медленная
сходимость, но важно, что метод сходится
всегда. Условие выхода из итерационного
процесса, например, может быть следующим:
,
где
– приближение корня на n-ной
итерации,
– заданная точность.
Блок-схема
метода бисекции приведена на рис. 4.1. В
данном алгоритме сужение отрезка
производится путём замены границ
или
на текущее значение корня
.
Значение
вычисляется лишь один раз, поскольку
нам нужен только знак функции
на левой границе, а он в процессе итераций
не меняется.
Да
Да
Нет
Нет
Рис. 4.1. Блок-схема метода бисекции
4.1.2. Метод хорд
В
данном методе процесс итераций состоит
в том, что в качестве приближений к корню
уравнения
принимаются значения
точек пересечения хорды с осью абсцисс.
Уравнение хорды
,
где
,
:
. (4.1)
Для
точки пересечения хорды с осью абсцисс
(
):
.
(4.2)
Далее
сравниваем знаки величин
.
Отрезок, на котором функция меняет знак
(например,
)
оставляем, вторую часть отбрасываем.
Следующая состоит в определении нового
приближения
,
как точки пересечения новой хорды с
осью абсцисс и т.д. Итерационный процесс
продолжаем до тех пор, пока значение
не станет по модулю меньше заданного
:
.
Блок-схема
метода хорд аналогична приведённой для
метода бисекции, но вместо вычисления
приближения по формуле
необходимо использовать формулу (4.2).
Так же в блок-схему необходимо ввести
операторы вычислений значений
на границах новых отрезков.
4.1.3. Метод касательных
При
уточнении корня по методу касательных
(Ньютона) в точке начального приближения
проводится касательная к функции
.
Точка
пересечения касательной с осью абсцисс
принимается за новое уточнённое значение
корня. Затем касательная к
проводится в
и т.д. Уточнение корня выполняется по
формуле:
. (4.3)
Если
процесс сходится, то
,
где
– значение корня.
Заметим, что данный метод очень быстро сходится, но, чтобы в результате очередной итерации не выйти за пределы отделения корня, необходимо выбирать начальное приближение довольно близким к решению, иначе очень вероятна расходимость.
4.1.4. Метод простой итерации
Уточнение
корня по этому методу сводится к замене
уравнения
ему равносильным:
.
Найденное начальное приближение
подставляют в правую часть уравнения
и получают новое значение корня по
формуле:
….
.
Достаточное
условие сходимости метода простой
итерации:
.
Блок-схема
метода простой итерации представлена
на рис. 4.2. На схеме
– начальное приближение корня и результат
последующей итерации,
– значение корня после каждой итерации.
Нет
Да
Рис. 4.2. Блок-схема метода простой итерации
