Решение задачи Дирихле для уравнения Лапласа методом сеток
Страницы работы
Фрагмент текста работы
12. РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ ДЛЯ УРАВНЕНИЯ ЛАПЛАСА
12.1. Цель работы
Целью работы является закрепление знаний, полученных в лекционном курсе «Вычислительная математика» по разделу «Численные методы решения дифференциальных уравнений в частных производных», приобретение навыков использования соответствующих численных методов с применением программных средств автоматизации вычислений.
Найти приближенное решение уравнения Лапласа в заданной области с указанными граничными условиями, принять h1=h2=0.1
1. 
2. 
3. 
4. 
5. 
6. 
7. 
8. 
12.3 Теоретические сведения
Уравнение относительно неизвестной функции u(x,y) двух или более независимых переменных , которое содержит частные производные этой функции, называется уравнением в частных производных (УЧП).
В общем виде ДУ второго порядка относительно функции двух независимых переменных u(x,y) записывается так :
(12.1)
Функция F— заданная функция восьми аргументов .
Далее будем рассматривать линейные уравнения второго порядка :
Все коэффициенты и правая часть не зависят от u(x,y).
Методы классификации УЧП :
· По порядку уравнения ( наивысший порядок частных производных, входящих в уравнение)
· По числу переменных ( по числу независимых переменных )
· По критерию « линейное-нелинейное»
A,B,C,D,E,F,G – константы или заданные функции переменных x,y.
· По критерию «однородное-неоднородное»
Однородное , если G(x,y)≡0 для всех x и y.
Если G(x,y)≠0 – неоднородное .
· По виду коэффициентов
A,B,C,D,E,F,G – константы – уравнение с постоянными коэффициентами.
Основные типы УЧП :
Описывает процессы теплопроводности и диффузии : B 2 =4AC
Описывает колебательные системы и волновые движения : B 2 >4AC
Описывает установившиеся процессы : B 2 2 ).
Для уравнения Лапласа система может быть записана в виде:
В данной работе система решается методом простых итераций по формулам: 
В качестве условия окончания итерационного процесса можно принять: 
Данная разностная схема обладает свойством устойчивости и сходимости. Это означает , что выбрав достаточно малый шаг h можно сколь угодно точно численно решить исходную задачу.
12.4 Пример выполнения задания
Найти решение задачи Дирихле для уравнения Лапласа :
Нумерация элементов матриц начнется с 1
Краевые элементы матрицы искомой поверхности
Для начала расчета зададим начальное приближение , то есть вектор 
в каждом слое. Рассмотрим слой с точками (0,0.1) и (1,0.1).
Первоначально считаем, что изменение функции происходит равномерно , т.е. с некоторым постоянным шагом :
отрезок разбит на 10 частей , следовательно изменения значения функции происходит на каждом шаге на величину 
Тогда получаем, что
Только внутренние элементы матрицы
Сделав такое предположение, мы получаем полностью заполненную матрицу u , при этом надо помнить , что это только начальные (стартовые
Решение задачи Дирихле для уравнения Лапласа методом сеток — реферат по математике
 |
|
Тезисы:
- Решить численно задачу Дирихле для уравнения Лапласа.
- Описание Метода Решения Поставленной Задачи.
- Дирихле для уравнений эллептического типа.
- (x0, y0 ) = (3, 5) .
- На первом шаге область D дискретизируется.
- (xi,yj+1) ? Dh , (xi,yj-1) ? Dh .
- F (N) * [i,j] через значения F (N-1) * [i,j]) до тех пор, пока не выполнится указанное условие.
- (x,y) ?D ; u|Г=xy2=f (x,y).
- Во всех точках области Dh вычисляются значения функции f (xi,yj).
- За область Dh* принимаются все внутренние точки области Dh, т.е. удовлетворяющие требованию.
Похожие работы:
248 Кб / 15 стр / 877 слов / 6633 букв / 4 авг 2018
192 Кб / 19 стр / 2462 слов / 16433 букв / 24 июл 2013
15 Кб / 6 стр / 1115 слов / 4995 букв / 10 мая 1998
26 Кб / 14 стр / 2450 слов / 11847 букв / 6 дек 1997
8 Кб / 6 стр / 1115 слов / 4995 букв / 10 мая 1998
469 Кб / 39 стр / 4821 слов / 31103 букв / 13 мая 2010
213 Кб / 45 стр / 3504 слов / 24193 букв / 27 сен 2012
559 Кб / 39 стр / 5509 слов / 37780 букв / 25 мая 2013
283 Кб / 30 стр / 1541 слов / 9985 букв / 21 янв 2015
153 Кб / 17 стр / 1036 слов / 7343 букв / 30 дек 2013
Готовые работы
база рефератов, курсовых, дипломных работ и статей
Решение задачи Дирихле для уравнения Лапласа методом сеток
от Для Студентов Бесплатно
Решение задачи Дирихле для уравнения Лапласа методом сеток
1. ПОСТАНОВКА ЗАДАЧИ
Решить численно задачу Дирихле для уравнения Лапласа :
область D ограничена линиями: x=2 , x=4 , y=x , y=x+4 ;
Следует решить задачу сначала с шагом по x и по y : h=0.2, потом с шагом h=0.1 . Точность решения СЛАУ ?=0.01 .
2. ОПИСАНИЕ МЕТОДА РЕШЕНИЯ ПОСТАВЛЕННОЙ ЗАДАЧИ
Поставленная задача решается численно с помощью программы, реализующей метод сеток , разработанный для численного решения задачи
Дирихле для уравнений эллептического типа.
Программа написана на языке C++ , в среде Borland C++ версии 3.1. Ниже описан алгоритм работы этой программы.
1. На первом шаге область D дискретизируется. Она заменяется на область Dh путем разбиения области D параллельными прямыми по следующему правилу: yi=y0 ± ih, xj=x0 ± ih , i,j=0,1,2…. РР
Разбиение производится до тех пор, пока текущая прямая не будет лежать целиком вне области D. Получается множество точек (xi,yj).
2. За область Dh принимают те точки множества (xi,yj) , которые попали внутрь области D, а также дополняют это множество граничными точками.
3.Во всех точках области Dh вычисляются значения функции f(xi,yj)
.
4. За область Dh* принимаются все внутренние точки области Dh, т.е. удовлетворяющие требованию:
(xi,yj) ? Dh* , если (xi+1,yj) ? Dh , (xi-1,yj) ? Dh ,
(xi,yj+1) ? Dh , (xi,yj-1) ? Dh .
5. Во всех точках области Dh* вычисляется функция F(N)*[i,j] ( индекс N обозначает номер итерации, на которой производится вычисление):
F(N)*[i,j]=(f(xi+1,yj) + f(xi-1,yj) + f(xi,yj+1)
+ f(xi,yj-1))/4
6. Теперь если max | F(N+1)*[i,j] — F(N)*[i,j]| =0.005) < for(i=0;i
http://studentlib.com/referat-216257-reshenie_zadachi_dirihle_dlya_uravneniya_laplasa_metodom_setok.html
http://rectorate.ru/2019/12/reshenie-zadachi-dirixle-dlya-uravneniya-laplasa-metodom-setok/