Теми рефератів
> Реферати > Курсові роботи > Звіти з практики > Курсові проекти > Питання та відповіді > Ессе > Доклади > Учбові матеріали > Контрольні роботи > Методички > Лекції > Твори > Підручники > Статті Контакти
Реферати, твори, дипломи, практика » Статьи » Рішення еліптичних рівнянь декількома методами

Реферат Рішення еліптичних рівнянь декількома методами





нечномірні (розмірністю N + 1) простором сіткових функцій. Очевидно, що сіткову функцію можна розглядати як вектор [1].

Можна також провести дискретизацію і простору функцій багатьох змінних, коли - точка p-мірного евклідового простору (p gt; 1). Так на площині можна ввести сітку, як безліч точок (вузлів) перетину перпендикулярних прямих,,, де - кроки сітки за напрямами і відповідно. Сітка, очевидно рівномірна по кожному з змінних окремо. Замість функції будемо розглядати сіткову функцію


Багато стаціонарні фізично завдання (дослідження завдань, пов'язаних з магнітними, гравітаційними явищами і теплопровідності) зводяться до вирішення рівняння Лапласа виду


(1)


Рішення для цього рівняння будемо шукати для деякої обмеженої області G зміни незалежних змінних x, y. Межею області G є замкнута лінія L. Для повної формулювання крайової задачі крім рівняння Лапласа потрібно задати граничну умову на кордоні L. Приймемо його у вигляді


(*)


Різницеві схеми


В області введемо рівномірну сітку


з кроками h по x і y. Замінюючи похідну по x і y різницевими виразами



в рівнянні (1) отримуємо:


(2)


У граничних умовах замінюємо безперервну функцію дискретної на області прямокутника:



За допомогою даних рівнянь можна записати систему лінійних алгебраїчних рівнянь (2) у вигляді (схема «Хрест»)


(3)


Система (3) містить 5 невідомих. І утворює систему лінійних рівнянь, де A - блочно-трехдіагональной матриця виду


, де

а матриці - діагональні матриці з елементами на діагоналі рівними одиниці. і розмірності.


Права частина системи має вигляд.


Для вирішення такої системи можна застосувати метод Гаусса з виключенням незмістовними операцій, тобто обнулення елементів. В даному випадку загальне число операцій складе величину

Ще одним з найбільш поширених методом вирішення цієї системи є ітераційний метод. Кожне рівняння запишемо у вигляді, дозволеному щодо значення в центральному вузлі:


(4)


Дану схему можна вирішити за допомогою ітераційних методів:

1. Гаусса-Зейделя,

(5)

2. Якобі,


(6)


3. Верхньої релаксації.


(7)


В алгоритмі передбачений вибір початкових значень. Ітераційний процес контролюється максимальним відхиленням M значень сіткової функції у вузлах для послідовних ітерацій. Якщо його величина досягне деякого допустимого значення точності, ітерації припиняються.

Існують і інші різницеві схеми для вирішення рівняння Лапласа. Представивши рівняння (1) у вигляді (8), де t - мінлива часу, то вихідна задача зводиться до різницевої схемою для рівняння теплопровідності, де граничні умови збігаються з умовами .Умова (9) можна вибирати практично в довільному вигляді, погодженому з граничною умовою.

Маємо, звідси можна знайти явне вираження для значення сіткової функції на (k + 1) -му шарі:


(10)


Граничні умови приймають вид


(11)


Вищеописаний метод називається «метод встановлення».

Процес чисельного рішення являє собою ітераційний процес вирішення завдання (8) з умовою (9) і (*) полягає в переході від довільного значення (9) до шуканого стаціонарного рішенню. Рахунок ведеться до виходу рішення в стаціонарний режим. Природно, обмежуються рішенням, при деякому досить великій, якщо значення шарів збігаються із заданим ступенем точності.

Умова стійкості вищеописаної схеми визначається нерівністю (12)

Розділ 2. Бібліотека функцій


У роботі розроблені наступні функції для вирішення неоднорідного одновимірного рівняння теплопровідності:

1. ElipYa (X, N, fi, E) -реалізует схему «Хрест» за допомогою методу Якобі,

2. Elip1 (X, N, fi, E) - вирішує схему «Хрест» за допомогою методу Гаусса-Зейделя.

3. ElipR (X, N, fi, t, E) - вирішує схему «Хрест» за допомогою методу верхньої релаксації.

4. ElipGauss (X, N, fi, E) - вирішує схему «Хрест» методом Гаусса.

5. ElipT1 (X, N, tau, fi, E) - зводить решаемую задачу до задачі теплопровідності і реалізує явну схему.

де x - довжина сторони квадрата сітки, N - кількість вузлів сітки по x, fi - гранична умова, E - точність обчислень. У методі верхньої релаксації t - параметр, а в явні...


Назад | сторінка 2 з 5 | Наступна сторінка





Схожі реферати:

  • Реферат на тему: Рішення системи двох лінійних рівнянь з поданням про вирішення в числовому ...
  • Реферат на тему: Чисельне рішення рівняння теплопровідності
  • Реферат на тему: Рішення систем лінійних алгебраїчних рівнянь методом Гаусса
  • Реферат на тему: Рішення системи лінійний алгебраїчних рівнянь модифікованим методом Гаусса
  • Реферат на тему: Рішення змішаної крайової задачі для гіперболічного рівняння різницевим мет ...