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

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





й схемі tau - крок за часом.

Кожна функція повертає масиви, складові сітку: і і значення сіткової функції. В якості початкових наближень використовується функція rand (I, J), яка повертає матрицю випадкових чисел розмірністю, зазначеної в дужках. Елементи матриці знаходяться в діапазоні.


Розділ 3. Тестування


Приймемо наступні параметри розроблених функцій:

X=10 довжина сторони квадрата;

N=15 - кількість вузлів сітки;

fi=cos (x) + cos (y) - функція граничних умов;

E=0.001 - точність;

t=1.2- параметр для методу верхньої релаксації;

tau=0.1 - крок за часом для «методу встановлення».

Вводимо початкові дані у функцію ElipYa - схема «хрест» метод Якобі


Малюнок 1


Всього ітерацій 167.

Тепер з тими ж параметрами зробимо розрахунки за допомогою методу Гаусса-Зейделя (Elip1).


Малюнок 2

Всього ітерацій 69.

Протестуємо метод верхньої релаксації (ElipT1)


Всього ітерацій 53.


Перевіримо метод Гаусса для вирішення схеми «хрест»


Результати збігаються з попередніми схемами.

А тепер перевіримо явну схему - «метод встановлення»

Малюнок 3


Всього ітерацій 89.

Зменшимо крок по сітці до h=0.59. Тобто X=10, а N=17.


Малюнок 4

Схема проявила нестійкість. Очевидно, що при цьому розкладі умова (12) не виповнилося.

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


Висновок


Запишемо результати в таблиці

МетодКолічество ітерацій Якобі167Гаусса-Зейделя69Метод верхньої релаксаціі53Явний метод установленія89

Для методу Гаусса не доводиться враховувати кількість ітерацій.

1.Установлено що явна схема має істотний недолік: накладаються обмеження на кроки і по сітці. Чого позбавлені неявні схеми.

.Ітераціонние методи ідеально підходять для сіткової схеми «хрест», метод верхньої релаксації виявився самим бистросходящімся.

.Для методу Гауса доводиться зберігати величезну матриці в пам'яті ЕОМ, що при достатньо великих N буде накладно.


Висновок


У ході роботи була вивчена різницева схема «хрест» для знаходження чисельного рішення рівняння Лапласа (еліптичне рівняння), завдання Дирихле.

Також були засвоєні і закріплені навички:

1. Використання покажчиків в середовищі matlab.

2. Програмування в середовищі matlab.

. Розробки чисельних методів для рівнянь еліптичного типу.

Були створені чотири методи реалізують разностную схему «хрест» і один метод для «методу встановлення».


Список використаної літератури


1.Самарскій А.А. Введення в чисельні методи. Навчальний посібник для вузів. 3-е изд., Стер.-СПб .: Видавництво «Лань», 2005. - 288 с .: ил.- (Підручники для вузів. Спеціальна література).

.Бахвалов Н.С., Жидков Н. П., Кобельков Г. М. Чисельні методи - 3-е изд, доп., і перероб.- М .: БИНОМ. Лабораторія знань, 2004. - 636 с., Іл.

.Агошков В.І., Дубровський П.Б., жартуючи В.П. Методи рішення задач математичної фізики/Под ред. Г. І. Марчука: Учеб. посібник.- М: ФИЗМАТЛИТ, 2002. - 320с.

.Метьюз Джон Р., Фінкс Куртіс Д Чисельні методи використання MATLAB, 3-е видання./Пер з англ.- М. Видавничий до «Вільямс», 2001 - 720с.

.Турчак Л. І., Плотніков П. В. Основи чисельних методів: Учеб. посібник.- 2-е изд., Перераб. І доп.- М .: ФИЗМАТЛИТ, 2005. - 304 с.

. # center gt; Програми


ElipYa.mfunction ElipYa (X, N, fi, E)% сітковий метод ітераційний, процес Якобі% X - Сторона прямокутника% Y - Сторона прямокутника% N - Кількість вузлів по х% K - Кількість вузлів по y% fi - функція граничної умови% E - точність h=X/N; % Рівномірна сітка по просторовим змінним x=[0: h: h * N]; % Розбиваємо сітку по змінній х y=[0: h: h * N]; % Розбиваємо сітку по перемінної у U=zeros (N + 1, N + 1); % Матриця значень сіткової функції F=zeros (N + 1, N + 1); % Матриця значень сіткової функції f (x, y)% дискретизація початкових умов і початкове заповнення матриці значень сіткової функції U (1,:)=feval (fi, x (1), y); % Перший рядок ...


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





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

  • Реферат на тему: Рішення систем нелінійніх рівнянь. Метод ітерацій. Метод Ньютона-Канторов ...
  • Реферат на тему: Розв'язання алгебраїчніх рівнянь. Метод простий ітерацій та Ньютона
  • Реферат на тему: Мінімізація функції багатьох змінних. Наближені чисельні методи. Метод Мо ...
  • Реферат на тему: Метод Ньютона (метод дотичних). Рішення систем нелінійних алгебраїчних рів ...
  • Реферат на тему: Чисельні методи пошуку стаціонарних точок у оптимізаційних задачах: метод Н ...