Численные методы Часть 2
..pdfМинистерство общего и профессионального образования
Российской Федерации
Пермский государственный технический университет
М. Г. Боярш инов
Чи с л е н н ы е м е т о д ы
Часть 2
Учебное пособие для студентов направления
“Прикладная математика и информатика”
Пермь 1999
УДК 681.3
Б86
Численные методы. Часть 2: Учебное пособие для студентов направления «Прикладная математика и информатика» / М. Г. Бояршинов; Перм. гос. техн. ун-т. Пермь, 1999. 200 с.
Учебное пособие написано на основе курса, читаемого студентам специальности “Прикладная математика” (специализация “Математическое моделирование’') в Пермском государственном техническом университете.
Введены основные понятия теории разностных схем. Рассмотрены методы решения обыкновенных дифференциальных уравнений (с начальными и граничными условиями). Решения дифференциальных уравнений в частных производных строятся с привлечением методов Галеркина, Ритца, наименьших квадратов; рассматриваются различные разностные схемы, аппроксимирующие дифференциальные операторы. Основное внимание уделяется оценкам порядка аппроксимации дифференциальных уравнений разностными операторами, вопросам устойчивости и сходимости численных решений. Рассматриваются методы сведения решения пространственных краевых задач к последовательности одномерных.
Пособие предназначено для студентов и аспирантов вузов, специалистов, занимающихся вопросами построения моделей систем и процессов. Может быть использовано как пособие для учителей средних учебных заведений при проведении факультативных занятий по компьютерному моделированию.
Печатается по решению редакционно-издательского совета Пермского государственного технического университета
Т*абл. 7. Ил. 58. Библиогр.: 27 назв.
Рецензенты: д-р. физ.-мат. наук Е. К. Хеннер, зав. кафедрой информатики Перм. гос. пед. ун-та; канд. техн. наук А. Н. Аношкин
ISBN 5-88151- 191-3 |
С Пермский государственный |
|
технический университет, 1999 |
Введение. |
6 |
1. Обыкновенные дифференциальные уравнении |
7 |
2. Задача Коши |
9 |
Устойчивость решения задачи Коши . |
9 |
Метод Пикара. |
11 |
Метод Эйлера. |
14 |
Метод Рунге - Купы . |
20 |
Методы Рунге - Кутгы третьего и четвертого порядков |
.25 |
Метод Адамса. |
27 |
Неявные схемы интегрирования . |
29 |
Дифференциальные уравнения высших порядков. Системы |
|
дифференциальных уравнений. |
32 |
Метод Эйлера для системы дифференциальных уравнений |
33 |
Метод Рунге - Купы для системы дифференциальных уравнений |
. 36 |
Разностные схемы интегрирования дифференциальных уравнений |
|
второго порядка. |
37 |
Контрольные вопросы и задания . |
39 |
3. Граничные задачи для обыкновенных дифференциальных уравнений |
41 |
Метод пристрелки. |
42 |
Метод дифференциальной прогонки . |
45 |
Метод моментов. |
54 |
Метод Галеркина. |
58 |
Разрешимость системы алгебраических уравнений метода |
|
Галеркина. |
64 |
Метод наименьших квадратов . |
66 |
Разрешимость системы уравнений метода наименьших квадратов |
. 69 |
Сходимость метода наименьших квадратов . |
70 |
Метод Ритца. |
72 |
Сходимость метода Ритца. |
74 |
Сеточный метод решения линейной граничной задачи . |
77 |
Разрешимость системы алгебраических уравнений метода сеток. |
78 |
Оценка порядка аппроксимации . |
79 |
Метод прогонки для решения сеточной задачи . |
80 |
Контрольные вопросы и задания . |
. 83 |
4. Дифференциальные уравнения в частных производных . |
84 |
Некоторые дифференциальные уравнения второго порядка . |
. 85 |
Уравнение теплопроводности. |
85 |
Уравнение свободных поперечных колебаний струны |
86 |
Уравнение стационарной диффузии . |
87 |
Дифференциальные уравнения для функций нескольких |
|
переменных. |
87 |
Метод Фурье разделения переменных . |
88 |
Основные понятия и определения теории разностных схем |
.91 |
Аппроксимация уравнения разностной схемой. Порядок |
|
аппроксимации. |
93 |
Устойчивость разностных схем |
96 |
Принцип максимума. |
100 |
Оценка устойчивости разностных схем методом Неймана . |
105 |
Сходимость решения разностной схемы . |
108 |
Контрольные вопросы и задания . |
110 |
5. Сеточные схемы для уравнений в частных производных . |
111 |
Уравнения первого порядка . |
] 11 |
Схемы бегущего счета. |
111 |
Явно-неявная схема . |
115 |
Многомерное уравнение переноса. |
115 |
Контрольные вопросы и задания . |
118 |
Уравнения параболического типа |
.119 |
Схема с «весами» . |
] 19 |
Трехслойная схема Ричардсона . |
122 |
Схема Дюфорта и Френкеле. |
. 123 |
Схема бегущего счета. |
124 |
Многомерные уравнения. |
125 |
Схема переменных направлений (продольно-поперечная схема). |
127 |
Метод расщепления. |
128 |
Контрольные вопросы и задания . |
130 |
Уравнения гиперболического типа . |
131 |
Схема «крест». |
131 |
Разностная схема с «весами» . |
134 |
Многомерные уравнения. |
135 |
Факторизация разностной схемы с «весами» . |
.136 |
Контрольные вопросы и задания . |
138 |
Уравнения эллиптического типа. |
.139 |
Контрольные вопросы и задания . |
144 |
6. Прикладные задачи математического моделирования . |
145 |
Аккумуляция загрязняющих веществ лесным массивом и их |
|
вторичный вынос. |
145 |
Движение потока жидкости по системе каналов городского |
|
коллектора. |
166 |
Моделирование состояния проволоки при знакопеременном изгибе |
|
с натяжением . |
176 |
Моделирование поведения тяжелой полосы металла в процессе |
|
свободного петлеобразования . |
187 |
Контрольные вопросы и задания |
195 |
Библиографический список. |
196 |
Предметный указатель. |
.198 |
Методы решения дифференциальных уравнений условно
классифицируются по следующему признаку:
точное решение, позволяющее представить искомую функцию в элементарных функциях;
приближенные решения, в которых точное решение получается как предел некоторой последовательности; в этом случае, как правило, используются разложения искомой функции в ряды Тэйлора, Фурье и так далее; - численные решения, когда искомая функция определяется для конечного
числа значений аргумента в узлах разностной сетки.
Благодаря развитию современных средств вычислительной техники методы решения дифференциальных уравнений с использованием компьютеров, в том числе персональных, получили широкое распространение. Задачи теплопроводности, механики жидкостей и газов, механики твердого деформируемого тела и многие другие были решены в основном благодаря широкому использованию сеточных методов. Вместе с тем следует иметь в виду, что неквалифицированное применение разностных схем к решению дифференциальных уравнений приводит к получению решений, далеких от истинных. Поэтому понятен интерес к теории разностных схем, позволяющей еще на стадии разработки алгоритма решения сложной инженерной проблемы выяснить условия успешной реализации вычислительной модели.
Одним из важнейших вопросов является оценка порядка аппроксимации разностной схемой исходного дифференциального уравнения, позволяющая судить о степени адекватности используемой сеточной модели исходной дифференциальной задаче. Основная идея решения дифференциальных уравнений численными методами на ЭВМ заключается, как правило, в сведении исходной задачи к решению систем алгебраических (линейных или нелинейных) уравнений. При этом, естественно, возникает вопрос о разреишмости этой получаемой системы. Численное решение сходится к точному, если при неограниченном увеличении числа алгебраических уравнений решение дискретизированной задачи стремится к решению исходной задачи. Поскольку решить систему с бесконечным числом алгебраических уравнений невозможно, весьма аюуальным представляется вопрос об оценке погрешности получаемых численных решений исходной дифференциальной задачи.
1. ОБЫКНОВЕННЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ
Решением, интегралом или интегралыюй кривой обыкновенного дифференциального уравнения первого порядка
jj^ = f(x,y) или y' = f(x,y)
называется дифференцируемая функция у = у(х), удовлетворяющая этому уравнению, то есть такая, что у'(х) ■ f(x,y(x)) на некотором участке изменения аргумента х.
Различают три типа задач:
- задача Коши, когда для одного из значений аргумента (например, х«),
принимаемого за начальное, известно значение функции уо, то есть |
у(х 0) = уп ; |
||||
- граничная задача1, когда условия на искомую функцию задаются для |
|||||
нескольких |
различных |
значений |
аргумента, |
у(х0) = Уо» |
У(х|) = Уи |
y (x j)= y 2.
- задачи на собственные значения, когда в формулировку задачи входят неопределенные параметры, определяемые из самой задачи (нахождение частот собственных колебаний многомассовых систем).
Численные методы следует применять лишь к корректно поставленным задачам, то есть таким, для которых решение существует, единственно и непрерывно зависит от входных данных, то есть от начальных (граничных) условий и коэффициентов уравнения.
Численное решение задачи считается корректным, если исходная задача поставлена корректно и ее дискретный аналог сохраняет свойства корректности, то есть получаемая система алгебраических уравнений однозначно разрешима и устойчива1по входным данным.
Пример 1.1. На материальную точку с массой т , находящуюся на гладкой горизонтальной поверхности, действует горизонтальная гармоническая сила F(t) = A sin(t). Найти закон движения точки, если ее начальная скорость v().
Уравнение, описывающее движение точки, имеет вид:
d2x а • М m dtг = А я п (0
Очевидно, что граничная задача может быть сформулирована лишь для дифференциальных уравнений второго и более высокого порядков.
2 Под устойчивостью в данном случае понимается непрерывная зависимость решения от входных данных [2].
х(0) ■ О, v(0) = |
|
= v0. |
|
Интегрирование исходного уравнения с заданными начальными условиями |
|||
дает решение |
|
|
|
4 0 =v0+ —О - cos(t)), x(t)=fv0+—\ t - —sin(t), |
|||
ш |
v |
niJ |
m |
непрерывно зависящее от начальных |
|
условий |
х(0) = 0, v(0) = v0 и |
коэффициентов А, ш.
д
Положим v0 = -----.В этом частном случае решение исходной задачи
ш
x(0 = ~ sin(0
ш
соответствует случаю гармонического колебания точки вблизи ее начального
положения. Иными словами, Vt > 0 |x(t]| £ — . m
Очевидно, что при малом отклонении величины v0 от выбранного значения, например за счет погрешности округления данных в ЭВМ, при численном решении задачи получаем (рис. 1.1)
Таким образом, даже при корректно поставленной задаче ее численное решение может бьггь связано со значительными трудностями.
x(t)
Рис. 1.1. Действительное решение (нижняя кривая) и возмущенное решение (верхняя кривая) задачи из примера 1.1
2. ЗАДАЧА КОШИ
Требуется найти функцию у(х)э удовлетворяющую дифференциальному уравнению
! = f(*,y(x)). 0< х< ;т, у(о)=у0. |
(2.1) |
Здесь f(x, у(х)) - заданная непрерывная функция двух аргументов. Условия существования и единственности решения задачи (2.1) устанавливает
Теорема 2.1 (Пеано1 [3]). Пусть функция f(x, у(х)) непрерывна в открытой области D. Тогда через каждую точку (u, v) этой области проходит хотя бы одна интегральная кривая. Каждая интегральная кривая может быть продолжена в обе стороны вплоть до границы любой замкнутой области G, целиком содержащейся в D и содержащей точку (u, v) внутри себя.
Кроме того, если функция f(x, у(х)) имеет в области D непрерывные частные производные первого порядка или удовлетворяет условию Липшица
|f(x.yi)-f(x,y,)|£K |yJ - y I|, V(x,y,),(x,y2)e G . K = const> 0,
то задача (2.1) имеет единственное решение.
В дальнейшем будем предполагать, что правая часть f(x, у) уравнения дифференцируема достаточное число раз.
Устойчивость решения задачи Коши
Наряду с уравнением (2.1) рассмотрим задачу
£ = f(x,?)+8(x). У(0) = У„+у, |
(2.2) |
где 5(х), у - возмущения, вносимые в правую часть уравнения (2.1) и начальные
условия.
Задача Коши (2.1) абсолютно устойчива [4], если Ve>0 3т|>0 такое, что Vx > 0 из условия ||5|| < е, |у| < е следует ||у- у\ й т]•
Иначе говоря, решение устойчиво, если «возмущенное» решение «не слишком сильно» отличается от исходного.
Рассмотрим условия устойчивости решения по начальному условию, то есть положим б(х) г 0.
1 Пеано Джузеппе [27.8.1858 - 20.9.1930] - итальянский математик. С 1890 года - профессор Туринского университета.
Оценим погрешность решения |
|
z(x) = y(x)-y(x). |
(2.3) |
Для этого, вычитая уравнение (2.1) из выражения (2.2), получим:
7 ^= f { \ j ) ~ % У ) = f(x.y + z) - f(x,y) ax
Воспользовавшись теоремой Лагранжа о среднем, для правой части последнего выражения получим:
f(x ,y+ z)-f(x,y)= z f;(x,y+0z), е е(0,1).
Решаем полученное дифференциальное уравнение
g - z . f ; ( x . y + ez),
разделяя переменные:
Y= f»'(x>y+0z) dx-
}f = }f;(t.y+ ez).dt,
*• |
О |
|
Ц 2)!. = J f»'(t»y+0z) d t* |
|
|
|
о |
|
z = z 0 ee |
|
|
Очевидно, что погрешность z(t) не будет возрастать, если |
Vx z(x) < z0. |
|
|
х |
|
Но это возможно в том случае, |
когда Jfy'(t,y + 0z) dt £ 0 . |
Поскольку это |
|
о |
|
требование должно выполняться для |
любого значения х, необходимо |
потребовать, чтобы выполнялось условие |
Ц й 0. |
Пример 2.1. Расмотрим задачу: — = к • у, у(0) = у0.
ах |
|
Точное решение этого уравнения имеет вид |
у(х) = у0екж Если к < 0, |
условие устойчивости выполняется: |
|
f; = k < 0 , |
|
К х)| = И х) - у(х)1 = |Уо - Уо| • е1> |
> 0 |