Методические указания к курсовой работе по дисциплине «численные методы» (стр. 1 из 2)

Министерство образования и науки Российской Федерации

Федеральное агентство по образованию

Южно-Российский государственный технический университет

(Новочеркасский политехнический институт)

Шахтинский институт (филиал) ЮРГТУ (НПИ)

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

к курсовой работе по дисциплине

«ЧИСЛЕННЫЕ МЕТОДЫ»

Новочеркасск 2008

УДК 681.3(076.5)

Рецензент – доцент В.В. Луценко

Составитель Бондаренко А.И.

Методические указания к курсовой работе по дисциплине «Численные методы» / Сост. А.И. Бондаренко; Шахтинский ин-т (филиал) ЮРГТУ (НПИ). – Новочеркасск: ЮРГТУ, 2008. - 12 с.

– 50 экз.

Методические указания содержат теоретический материал, примеры выполнения и требования к оформлению курсовой работы по дисциплине «Численные методы».

Предназначены для студентов второго курса специальностей 230201«Информационные системы и технологии» и 0808001 «Прикладная информатика».

УДК 681.3(076.5)

© Шахтинский институт ЮРГТУ, 2008

© Бондаренко А.И., 2008

ВВЕДЕНИЕ

Изучение различных процессов требует наряду с глубоким пониманием физики происходящих явлений совершенного владения современными методами вычислительной математики.

Обычно математическая модель записывается в форме как угодно сложных математических структур и, как правило, получить аналитическое решение такой задачи не удаётся. Приходится использовать численные методы вычислительной математики, реализация которых на ЭВМ требует соответствующего программного обеспечения. Результаты моделирования объекта на ЭВМ позволяют “проиграть” его поведение в самых разных, подчас экстремальных условиях. Значение такого вычислительного эксперимента трудно переоценить, особенно если натурный эксперимент опасен, дорог или просто невозможен.

Большинство физических процессов можно сформулировать на языке дифференциальных уравнений с частными производными. Производные в этих уравнениях описывают важнейшие физические величины: скорость, ускорение, силу, температуру, трение, ток, потенциал и т.д.). Многие из таких уравнений не имеют аналитического решения и, чтобы их решить, приходится прибегать к численным методам.

В курсовой работе рассматривается одно из самых важных уравнений математической физики - уравнение Лапласа на примере решения задачи Дирихле в заданной плоской области. Отсутствие аналитического решения поставленной задачи требует выбора численного метода и его реализации на ЭВМ.

Курсовая работа является завершающим этапом изучения курса “Численные методы”. Цель курсовой работы:

· систематизация, закрепление и расширение теоретических и практических знаний по компьютерному моделированию типовых вычислительных алгоритмов и анализа полученной информации;

· выявление степени подготовленности студентов к самостоятельной работе в ходе решения поставленных задач.

3

1. Аналитические методы решения уравнений в частных производных

Существует целый арсенал методов для решения уравнений в частных производных. Перечислим некоторые аналитические методы решения таких уравнений.

Метод разделения переменных. Уравнение с частными производными с n независимыми переменными сводится к n обыкновенным дифференциальным уравнениям. Решение краевых задач для уравнения Лапласа может быть найдено методом разделения переменных (методом Фурье) лишь для простейших областей (круг, прямоугольник, шар цилиндр и др.).

Метод преобразования координат. Исходное уравнение с частными производными сводится к обыкновенному дифференциальному уравнению или к другому, более простому уравнению с частными производными с помощью соответствующего преобразования координат (например, поворота координатных осей и т.п.).

Введение новых переменных. Исходное уравнение с частными производными преобразуется к такому уравнению с частными производными для другой неизвестной функции, которое решается легче, чем исходное.

Метод интегральных уравнений. Уравнение с частными производными сводится к интегральному уравнению (уравнение, в котором неизвестная функция стоит под знаком интеграла).

Вариационные методы. Вместо уравнения с частными производными решается некоторая задача минимизации. Оказывается, что функция, доставляющая минимум некоторому выражению, является решением исходного уравнения.

Метод разложения по собственным функциям. Эти собственные функции находятся как решения так называемой задачи на собственные значения, которые соответствуют исходной задаче для уравнения с частными производными.

Метод функций Грина. Начальные и граничные условия заменяются системой простейших источников, и задача решается для каждого простейшего источника. Полное решение исходной задачи получается в результате суммирования решений для элементарных источников.

4

2. Численные методы решения

уравнения Лапласа

В курсовой работе требуется решить задачу Дирихле для уравнения Лапласа

(1)

одним из известных численных методов. Многие установившиеся процессы различной физической природы описываются этим уравнением.

Ставится задача о нахождении стационарного (не изменяющегося с течением времени) распределения температуры внутри многоугольника, если известно распределение температуры вдоль его сторон, т.е. необходимо решить задачу Дирихле для уравнения Лапласа.

Одна из главных трудностей, возникающих при решении этой задачи обусловлена сложной формой границы расчетной области. Аналитические решения удается получить лишь в частных случаях - для простейших областей (прямоугольник, круг, сектор, шар). Основными методами решения поставленной задачи являются численные методы.

Универсальным методом приближенного решения дифференциальных уравнений с частными производными является метод конечных разностей (или метод сеток). Суть его в следующем. Область непрерывного изменения аргументов (x,y) заменяется конечным (дискретным) множеством точек (узлов), называемым сеткой; вместо функций непрерывного аргумента рассматриваются функции дискретного аргумента, определенные в узлах сетки.

Будем считать, что заданная область покрыта прямоугольной сеткой. Узлы сетки, лежащие на границе области, называются граничными, а все остальные - внутренними. Прямоугольные сетки наиболее удобны при организации вычислительных алгоритмов (рис.1).

Уравнение (1) записывается для каждого внутреннего узла расчетной области, дифференциальное уравнение при этом заменяется системой алгебраических уравнений (разностными уравнениями).

Если полученная разностная краевая задача разрешима и ее решение при измельчении сетки приближается (сходится) к решению

5

исходной задачи для дифференциального уравнения, то оно и принимается за приближенное решение исходной задачи.

Рис.1. Прямоугольные сетки для решения задачи Дирихле

Производные, входящие в дифференциальное уравнение (1), заменяются соответствующими разностными отношениями:

;

.

Эти вычислительные шаблоны применяют там, где узлы прямоугольной сетки не являются равноотстоящими. Если a=b=h₁, а g=d=h₂, то применяют следующий вычислительный шаблон:

;
.

В случае квадратной сетки уравнение (1) для нулевого узла будет иметь вид:

Так как вычислительные шаблоны связывают лишь несколько соседних узлов сетки, то матрица коэффициентов системы линейных уравнений для определения узловых неизвестных оказывается “разряженной”, т.е. содержит много нулевых элементов.

Рассмотрим пример использования приведенных вычислительных шаблонов для решения уравнения (1).

6

Пусть на сторонах трапеции (рис. 2), заданной координатами вершин A(0,2), B(2,2), C(4,0), D(0,0) известно распределение температуры T(x,y)=xy+x+y. Необходимо найти распределение температуры внутри области.

Рис. 2. Заданная расчетная область

Покроем квадратной сеткой рассматриваемую область (рис. 2) и составим систему линейных уравнений для определения неизвестных значений T₁ и T₂:

(2)

Вычислим значение температуры в граничных узлах (Т₃ – Т₉) сетки по заданной функции на границе T(x,y)=xy+x+y и подставим в систему уравнений (2). Получим систему для определения неизвестных значений T₁ и T₂:

Искомое решение T₁ =3, T₂ = 5. Выбор числа узлов сетки зависит от требуемой точности. Зная значения температуры в узлах, можно найти температуру в любой точке области, применяя интерполяционные формулы для функций нескольких переменных.