- •Федеральное агентство по образованию
- •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. Задания
- •Библиографический список
- •Использование
8. Лабораторная работа № 4. Решение нелинейных уравнений приближенными методами
Цель работы:
Освоение методов приближенного решения нелинейных уравнений: метода деления отрезка пополам, метода Ньютона, метода простых итераций с последующей реализацией на ЭВМ.
8.1. Постановка задачи
В соответствии с выбранным вариантом (пункт 8.7) осуществить решение нелинейного f(x)=0 уравнения с точностью=0,001 тремя методами: методом деления отрезка пополам, методом Ньютона, методом простых итераций.
Произвести приближенную оценку корней.
Произвести уточнение корней тремя методами оценки.
Осуществить проверку полученного каждым методом решения.
Дать сравнительную характеристику трех методов, используемых для решения уравнения, по следующим характеристикам:
быстродействие метода (по количеству итераций);
сложность реализации метода (преимущества и недостатки).
8.2. Порядок выполнения работы
Изучить методические указания и ответить на контрольные вопросы.
Получить у преподавателя номер варианта.
В соответствии с поставленной задачей:
осуществить приближенную оценку корней заданного нелинейного уравнения, для чего построить график y=f(x)и найти область пересечения [a, b] графика функции с осью абсцисс (наиболее просто построить график в среде Mathcad);
произвести уточнение корней, представив математическую формулировку решения нелинейного уравнения по трем методам, составив схемы алгоритмов решения и написав программы в среде СИ++ или Mathcad.
Отладить программы и получить результаты расчетов.
Результаты расчетов представить графически.
Провести анализ полученных результатов (дать сравнительную характеристику использованных методов по точности, сложности реализации, быстродействию).
Оформить отчет.
8.3. Краткие теоретические сведения
8.3.1. Этапы решения нелинейного уравнения
Уравнениемназывается равенство с переменной, которое в общем виде записывается в виде:
. (59)
Значение переменной х=х*, обращающее уравнение в тождествоf(x*)=0называетсякорнем уравнения.
Решить уравнение означает найти все его корни. Уравнение, вид которого не позволяет получить формулу для расчета точного значения корня, решается приближенно, например, x=cos(x)(получить систематическое решение невозможно).
Задача приближенного решения уравнения (59) заключается в исследовании функции
с целью поиска такой точки х*на оси, в которой значение функции обращается в нуль, т. е.y*=f(x*)=0.
Численные методы решения складываются из двух этапов.
Отделение корней, т. е. нахождение такого интервала [a, b], в котором существует единственный корень. Таких интервалов может быть найдено столько, сколько существует действительных корней у решаемого уравнения.
Существует несколько методов отделения корней: аналитический, графический, графоаналитический. Чаще всего на практике пользуются комбинацией графического и аналитического методов.
Для уравнения f(x)приблизительно строится график. Отделяют интервал[a, b]предположительно содержащий корень, а затем функция в этом интервале исследуется на выполнение трех условий:
функция в интервале [a, b] должна быть непрерывна;
монотонна на [a, b], т. е. первая производная не меняет свой знак на этом интервале;
на конце интервала функция f(x)меняет знак.
Если эти условия выполняются, то интервал [a, b]содержит действительный корень, и причем единственный.
Например, требуется отделить корень уравнения
. (60)
Для этого удобно построить графики функций f(x)=sin(2x)иf(x)=ln(x)(рис. 50,а), а затем на осиOXотметить отрезки, локализующие абсциссы точек пересечения рассматриваемых кривых. Из графиков следует, что уравнение имеет корень, принадлежащий отрезку1;1,5. В другом варианте – построить график функцииf(x)=sin(2x)-ln(x). Пересечение графика с осьюОХ– определяет местонахождение корня (рис. 50,б).
а) б)
Рис. 50. Графический способ отделения корней сложной функции
Уточнение корня, т. е. нахождение его значения внутри интервала [a, b] с заданной степенью точности.
Задача уточнения корня формулируется следующим образом: пусть на интервале [а, b] имеется действительный корень и причем единственный. Необходимо найти этот корень с заданной степенью точности.
Существует большое разнообразие вычислительных методов, реализующих поставленную задачу, однако последовательность основных этапов решения задачи одинакова для всех методов и может быть представлена в виде блок-схемы (рис. 51).
Рис. 51. Этапы отделения корней
Все существующие вычислительные методы уточнения корней нелинейного уравнения условно делятся на 3 группы:
методы деления отрезка;
методы, основанные на информации о значении первой производной;
методы, использующие рекуррентные выражения.
В данной лабораторной работе рассматриваются методы, относящиеся к разным группам.