- •Федеральное агентство по образованию
- •2404000000-35 Удк 681.142:519.6
- •Оглавление
- •1. Основы работы в Mathcad 10
- •2. Роль численных методов 36
- •3. Методы аппроксимации и интерполирования 42
- •4. Лабораторная работа № 1. Интерполирование степенными многочленами 44
- •5. Лабораторная работа № 2. Параметрическая идентификация математических моделей методами аппроксимации 65
- •6. Решение систем линейных алгебраических уравнений 93
- •7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами 109
- •8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами 142
- •9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами 174
- •10. Лабораторная работа № 6. Решение обыкновенных дифференциальных уравнений первого порядка численными методами 197
- •11. Лабораторная работа № 7. Численное интегрирование 212
- •12. Лабораторная работа № 8. Моделирование реактора идеального вытеснения для многостадийной химической реакции с линейной кинетикой 226
- •13. Лабораторная работа № 9. Расчет моделей процессов диффузии и теплопроводности с помощью явной разностной схемы 249
- •Предисловие
- •1. Основы работы в Mathcad
- •1.1. Панели инструментов
- •1.2. Ввод и вывод данных
- •1.3. Осуществление несложных вычислений
- •1.4. Построение и настройка графиков
- •1.5. Программирование в Mathcad
- •1.5.1. Программирование без программирования
- •1.5.2. Язык программирования Mathcad
- •1.5.3. Создание программы (Add Line)
- •1.5.4. Редактирование программы
- •1.5.5. Локальное присваивание ()
- •1.5.6. Условные операторы (if, otherwise)
- •1.5.7. Операторы цикла (for, while, break, continue)
- •1.5.8. Возврат значения (return)
- •1.5.9. Перехват ошибок (on error)
- •1.5.10. Примеры программирования
- •2. Роль численных методов
- •2.1. Этапы решения задачи на компьютере
- •2.2. Математические модели
- •2.3. Численные методы
- •3. Методы аппроксимации и интерполирования
- •4. Лабораторная работа № 1. Интерполирование степенными многочленами
- •4.1. Постановка задачи
- •4.2. Порядок выполнения работы
- •4.3. Краткие теоретические сведения
- •4.3.1. Метод неопределенных коэффициентов
- •4.3.2. Интерполяционный многочлен Лагранжа
- •4.3.3. Интерполяционные формулы Ньютона для равностоящих узлов
- •4.4. Примеры выполнения
- •4.4.1. Интерполирование степенными многочленами с использованием метода неопределенных коэффициентов
- •4.4.2. Интерполирование степенными многочленами с использованием второй интерполяционной формулы Ньютона
- •4.5. Требования к отчету
- •4.6. Контрольные вопросы и задания
- •4.7. Задания
- •5. Лабораторная работа № 2. Параметрическая идентификация математических моделей методами аппроксимации
- •5.1. Постановка задачи
- •5.2. Порядок выполнения работы
- •5.3. Краткие теоретические сведения
- •5.3.1. Метод выбранных точек
- •5.3.2. Метод средних
- •5.3.3. Метод наименьших квадратов
- •5.4. Примеры выполнения
- •5.4.1. Аппроксимация с использованием метода выбранных точек
- •5.4.2. Аппроксимация с использованием метода средних
- •5.4.3. Аппроксимация с использованием метода наименьших квадратов
- •5.4.4. Сравнительный анализ методов аппроксимации
- •5.5. Требования к отчету
- •5.6. Контрольные вопросы и задания
- •5.7. Задания
- •6. Решение систем линейных алгебраических уравнений
- •6.1. Общие положения
- •6.2. Точные методы решения систем линейных уравнений
- •6.2.1. Метод Крамера
- •6.2.2. Метод Гаусса
- •6.2.3. Метод обращения матриц
- •7. Лабораторная работа № 3. Решение систем линейных уравнений приближенными методами
- •7.1. Постановка задачи
- •7.2. Порядок выполнения работы
- •7.3. Краткие теоретические сведения
- •7.3.1. Математическое описание реактора идеального смешения непрерывного действия
- •7.3.2. Математическое описание кинетических закономерностей химических превращений
- •7.3.4. Приближенные методы решения систем линейных уравнений
- •7.3.4.1. Метод простых итераций
- •7.3.4.2. Метод Зейделя
- •7.4. Примеры выполнения
- •7.4.1. Пример выполнения задания точным методом
- •7.4.2. Пример выполнения задания методом итераций и методом Зейделя
- •7.5. Требования к отчету
- •7.6. Контрольные вопросы и задания
- •7.7. Задания
- •8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами
- •8.1. Постановка задачи
- •8.2. Порядок выполнения работы
- •8.3. Краткие теоретические сведения
- •8.3.1. Этапы решения нелинейного уравнения
- •8.3.4. Метод деления отрезка пополам (вилки, дихотомии)
- •8.3.5. Метод Ньютона (метод касательных)
- •8.3.6. Метод простых итераций
- •8.4. Пример выполнения задания методом итераций
- •8.5. Требования к отчету
- •8.6. Контрольные вопросы и задания
- •8.7. Задания
- •9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами
- •9.1. Постановка задачи
- •9.2. Порядок выполнения работы
- •9.3. Краткие теоретические сведения
- •9.3.1. Метод Ньютона
- •9.3.2. Метод итераций
- •9.4. Примеры выполнения
- •9.4.1. Метод Ньютона
- •9.4.2. Метод итераций
- •9.5. Требования к отчету
- •9.6. Контрольные вопросы и задания
- •9.7. Задания
- •10. Лабораторная работа № 6. Решение обыкновенных дифференциальных уравнений первого порядка численными методами
- •10. 1. Постановка задачи
- •10.2. Порядок выполнения работы
- •10.3. Краткие теоретические сведения
- •10.3.1. Метод Эйлера
- •10.3.2. Модифицированный метод Эйлера
- •10.3.3. Метод Эйлера-Коши
- •10.3.4. Метод Рунге-Кутта 4-го порядка
- •10.4. Примеры выполнения
- •10.4.1. Реализация метода Эйлера в математическом редактореMathcad
- •10.4.2. Решение обыкновенного дифференциального уравнения 1-го порядка с помощью функции rkfixed
- •10.5. Требования к отчету
- •10.6. Контрольные вопросы и задания
- •10.7. Задания
- •11. Лабораторная работа № 7. Численное интегрирование
- •11. 1. Постановка задачи
- •11.2. Порядок выполнения работы
- •11.3. Краткие теоретические сведения
- •11.3. 1. Метод прямоугольников
- •11.3.2. Метод трапеций
- •11.4. Пример выполнения
- •11.5. Требования к отчету
- •11.6. Контрольные задания
- •1 Таблица 101.7. Задания
- •12. Лабораторная работа № 8. Моделирование реактора идеального вытеснения для многостадийной химической реакции с линейной кинетикой
- •12. 1. Постановка задачи
- •12.2. Порядок выполнения работы
- •12.3. Краткие теоретические сведения
- •12.3.1. Математическая модель реактора идеального вытеснения
- •12.3.2. Численное решение систем дифференциальных уравнений
- •12.4. Пример выполнения
- •12.5. Проверка расчета с помощью функции rkfixed
- •12.6. Требования к отчету
- •12.7. Контрольные вопросы и задания
- •1 Таблица 112.8. Задания
- •13. Лабораторная работа № 9. Расчет моделей процессов диффузии и теплопроводности с помощью явной разностной схемы
- •13.3.2. Решение уравнений в частных производных
- •13.3.3.Метод сетки
- •13.3.4. Явная разностная схема
- •13.3.5. Условия устойчивости явной разностной схемы
- •13.4. Пример выполнения
- •13.5. Требования к отчету
- •13.6. Контрольные вопросы и задания
- •13.7. Задания
- •Библиографический список
- •Использование
9. Лабораторная работа № 5. Решение систем нелинейных уравнений приближенными методами
Цель работы:
Освоение методов приближенного решения систем нелинейных уравнений: метода Ньютона, метода итераций с последующей реализацией на ЭВМ.
9.1. Постановка задачи
В соответствии с выбранным вариантом (пункт 9.7, табл. 8) осуществить решение системы нелинейных уравнений f(x,y)=0 и(x,y)=0 с точностью=0.001 двумя методами – методом Ньютона и методом итераций:
произвести приближенную оценку корней, для чего построить графики y=f*(x) иy=*(x) и определить область их пересечения, на осиX[a1,b1], на осиY [a2,b2].
произвести уточнение корней двумя методами оценки, для чего составить алгоритм и программу решения системы нелинейных уравнений.
Дать сравнительную характеристику методов, используемых для решения уравнения, по следующим характеристикам:
быстродействие метода (оценивается по количеству итераций);
сложность реализации метода (указать преимущества и недостатки).
9.2. Порядок выполнения работы
Изучить методические указания и ответить на контрольные вопросы.
Получить у преподавателя номер варианта.
В соответствии с поставленной задачей:
осуществить приближенную оценку корней заданной системы нелинейных уравнений;
произвести уточнение корней, с этой целью представить математическую формулировку решения системы нелинейных уравнений по двум методам, составить схемы алгоритмов решения и написать программы в среде СИ++ или Mathcad.
Отладить программы и получить результаты расчетов.
Результаты расчетов представить графически.
Провести анализ полученных результатов (дать сравнительную характеристику двух использованных методов по сложности реализации, быстродействию).
Оформить отчет.
9.3. Краткие теоретические сведения
Решение систем нелинейных уравнений состоит из двух этапов:
отделение корней;
уточнение корней в отдельных интервалах.
Задача отделения корней даже для нелинейных уравнений может вызвать значительные трудности. А в случае систем нелинейных уравнений эта задача еще более усложняется.
(89)
Отделение корней будет заключаться в приближенном определении точек пересечения двух кривых f(х,у), (х,у). Для этого уравненияf(x,y)=0 и(x,y)=0 следует представить в видеyy=f*(x) иxx=*(y) и построить графики. Точка пересечения определит корни уравнения.
Например, имеется система нелинейных уравнений:
(90)
Преобразуем её в
(91)
и построим графики (рис. 79).
Из рис. 79 видно, что корень уравнения находится в области 6 x 9, –7y –4.
Рассмотрим два метода уточнения корней, которые используются для решения систем нелинейных уравнений.
9.3.1. Метод Ньютона
Рассмотрим систему нелинейных уравнений 2-го порядка (89) в общем виде. Однако необходимо учитывать, что метод Ньютона может использоваться и для систем уравнений более высокого порядка.
Выберем х0иу0в качестве грубого приближения к решению,х0иу0могут выбираться либо на этапе отделения корней, либо из физических соображений или при постановке задачи.
Рис. 79. Отделение корней
Для реализации последовательного приближения грубого приближения к решению (точному) х и у необходимо записать в алгоритм поиска наиболее точного решения в следующем виде:
(92)
где i=0, 1,, при условии, чтох0иу0известны.
То есть для решения системы уравнений необходимо указать правила нахождения добавок hi,ki.
Для этого введем следующее обозначение
. (93)
Подставим (93) в систему (89):
. (94)
Разложим уравнения системы (94) в ряд Тейлора в окрестностях точки (х0,у0), при этом ограничимся линейными членами разложения:
,(95)
где 1,2– нелинейные относительно(x-x0)и(y-y0)члены разложения.
Если пренебречь 1и2, обозначить,,
то получим систему линейных уравнений относительно hиk:
. (96)
Решение системы (96) даст неизвестные hиk, которые приближаютх0,y0к точному решению, но точного значения корня не обеспечивают (вследствие пренебрежения остаточными членами разложения1и2).
Т. к. f(x0, y0), (x0, y0), f’x(x0, y0), f’y(x0, y0), x’(x0, y0), y’(x0, y0)– константы, то систему уравнений (96) можно представить в привычном виде:
. (97)
Решение системы (97) может быть получено с использованием методов решения систем линейных уравнений, например, правила Крамера:
;, где
Определитель матрицы, составленный из первых производных системы уравнений, называется Якобианом. После того как hиk найдено, необходимо повторять процесс поиска новых значенийhиkи до тех пор, пока решение не достигнет заданной степени точности. При этом в качестве начального приближения выбирается всякий раз очередное приближение к решению, т. е. итерационный процесс поиска можно представить в следующем виде:
. (98)
Условием прекращения поиска решения является выполнение следующего условия:
, (99)
причем эти условия должны выполняться одновременно.
Таким образом, поиск решения выполняется при реализации следующей последовательности действий:
выбираются начальные значения х0иу0;
для этих значений рассчитываются значения функций f, , f'x, f'y, ’x, ’y;
решается система линейных уравнений (96), т. е. находятся значения hиk;
по формулам уравнений (98) находятся значения х, у.Процесс поиска, т. е. действия (2)–(4) выполняются до тех пор, пока не выполнится условие достижения заданной степени точности (99).
Блок-схема решения системы нелинейных уравнений методом Ньютона представлена на рис. 80.