Методические указания к компьютерному практикуму по курсу «Уравнения математической физики» Москва

скачать Московский государственный технический университет имени Н.Э. Баумана К.В. Титов ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ ЗАДАЧ ДИФФУЗИИ Методические указания к компьютерному практикуму по курсу «Уравнения математической физики» Москва Издательство МГТУ им. Н.Э. Баумана 2009 УДК 518.12 ББК 22.193 Т454 Р Т454 ецензент И.О. Янов Титов К.В. Численные методы решения задач диффузии: Метод. указания к компьютерному практикуму по курсу «Уравнения математической физики». — М.: Изд-во МГТУ им. Н.Э. Баумана, 2009. — 48 с.: ил. Справочно представлены основные методы численного решения обыкновенных дифференциальных уравнений и краевых задач. Механизм и эффективность работы этих методов выявляются в процессе выполнения компьютерного практикума. Это способствует формированию у студентов необходимой теоретической и практической базы знаний для последующего решения прикладных задач диффузии. Даны все необходимые рекомендации для проведения вычислительных работ на персональных компьютерах по численным методам решения некоторых задач математической физики. Приведены теоретический материал, необходимый для работы с электронной версией методических указаний, и условия типового расчета. Для студентов старших курсов (третий и выше) специальности «Ядерные реакторы и установки». Изложенный материал может быть полезен и другим пользователям с точки зрения применения компьютерных технологий в обучении, в том числе дистанционном. УДК 518.12 ББК 22.193  МГТУ им. Н.Э. Баумана, 2009 Предисловие Поставленная перед автором задача состояла в том, чтобы за весьма ограниченное время научить студентов решать некоторые прикладные задачи по их будущей специальности. Поэтому многое из того, что составляет предмет численных методов решения обыкновенных дифференциальных уравнений (ОДУ), пришлось опустить или отослать читателя к самостоятельной проработке материала. Данные методические указания следует рассматривать лишь как введение в численные методы решения ОДУ, ориентированное в первую очередь на специальные прикладные задачи математической физики. Однако, несмотря на временные ограничения, в процессе обучения студент получает в свое распоряжение завершенный, эффективно работающий инструмент для решения рассматриваемого класса задач. Помимо «твердой копии» методические указания имеют также электронную версию, что позволяет эффективнее вести процесс обучения, в частности, проводить лабораторные работы, закрепляющие теоретические знания на практике. По данной тематике автором был прочитан курс лекций и проведен ряд лабораторных работ, в том числе зачетная лабораторная работа для студентов специальности «Ядерные реакторы и установки». В основе предлагаемого пособия лежат методические указания [1]. Автор выражает свою признательность и благодарность профессору Г.И. Богомолову за полезные советы и обсуждение и предоставленную возможность использования таблицы вариантов для проведения зачетной лабораторной работы. ^ 1. Введение в численные методы решения обыкновенных дифференциальных уравнений Численные методы и их реализация в среде систем компьютерной математики, к которым в первую очередь следует отнести MathCAD, Maple, MATLAB, Mathematica и некоторые другие, не только расширяют класс реально решаемых дифференциальных уравнений, но и существенно упрощают сам процесс их решения. В этом случае можно говорить об интегрируемости ОДУ в явном виде, когда решение может быть получено с помощью конечного числа элементарных операций, к которым относят вычисления на ЭВМ. Решение можно представить в виде графика, в том числе трехмерного, или анимационного клипа, соответствующего изменению параметров задачи в интерактивном режиме. Самый общий и хронологически более ранний подход в решении дифференциальных уравнений — представление решения рядом Тейлора, которое по разным причинам (в первую очередь из-за отсутствия мощной вычислительной техники) отвергалось и заменялось другими, более простыми методами. Этот метод интегрирования дифференциальных уравнений был исключен из широкой практики более столетия назад [2]. С появлением систем компьютерной математики вычислительные трудности были устранены, и метод разложения в ряд Тейлора можно с успехом применять при решении ОДУ. Суть метода состоит в том, что решение на отрезке [xi, xi+1] в окрестности точки xi отыскивают в виде (1) где n — целое положительное число; y(j)(xi) — производная от y(x) порядка j, вычисленная в точке xi. Такой подход был использован в [3] для решения широкого класса дифференциальных уравнений, в том числе нелинейных. В методических указаниях рассмотрены так называемые методы Рунге – Кутта: метод Эйлера, метод Адамса, уже упомянутый метод Тейлора (см. уравнение (1)) и др. Дан сравнительный анализ их точности, которая продемонстрирована на решении конкретных задач лабораторного практикума. Решение большинства задач, рассмотренных в методических указаниях, проведено в системе MathCAD как наиболее простой. Хотя с таким же успехом все они могли быть решены в любой другой из перечисленных выше систем компьютерной математики. Использование системы MathCAD приводит к необходимости достаточно специфического изложения учебного материала, что обусловлено требованиями этой системы. Имеется в виду, что в методических указаниях автор средствами редактора Word (по просьбе Издательства) пытается представить решения задач в виде, приближенном к тому, в котором эти решения записывались бы с помощью системы MathCAD. Автор предлагает читателю, изучающему материал, принять эту условность. ^ 1.1. Метод Эйлера Пусть требуется решить уравнение вида для заданных начальных условий: численно методом конечных разностей. Число точек дискретизации: Задание функции: Область определения функции: Начальные условия: Просмотр на экране дисплея численного значения функции в этой точке: Дискрет по х: . Для выбранных параметров Дискретное представление самой переменной: Таким образом, для вычисления у(х) получим выражение (2) Приведем некоторые его значения: Точное решение этого дифференциального уравнения: Значение функции на правой границе: (рис. 1; возможна анимация графика по параметру n). Рис. 1 Разностная схема (метод) не всегда работает. Если и то в силу (2) значения будут равны что не всегда так. Примером может служить процесс, описываемый уравнением или уравнением с заданными начальными условиями: . Независимо от значения q, — изменяемая во времени функция. Однако при q = 0 разностная схема дает другой результат: При картина получается достоверной. Если не ввести обозначение константы, то возникнут трудности при записи функции : Пусть переменные z и t принадлежат соответственно интервалам и . Запишем переменную t в дискретном виде (значение Т подобрано как полупериод): Рис. 2 Введем обозначение При этом (2) может быть записано в виде первых двух слагаемых ряда Тейлора: В результате получаем ошибочное решение zj, представляемое на графике прямой линией (рис. 2; возможна анимация графика). В то же время при будем иметь: Это решение также представлено на рис. 2 и является единственно правильным. Задания. 1. Снять клип по параметру q1: 2. Привести другие примеры неадекватной работы метода Эйлера. Скачать 309,22 Kb. оставить комментарий страница 1/8 Дата 02.10.2011 Размер 309,22 Kb. Тип Методические указания, Образовательные материалы