41620

Решение задачи Дирихле для уравнения Пуассона методом Чебышева

Лабораторная работа

Информатика, кибернетика и программирование

Разностную задачу 5 будем решать явным итерационным методом с чебышевским набором параметров который выражается следующей формулой: 10 где заданное число итераций . 11 Результаты: В вычислениях использовался следующий алгоритм: Задаём количество итераций полагаем тогда шаг сетки =01. Полученный ответ с точностью до...

Русский

2013-10-24

103.07 KB

44 чел.

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

учреждение высшего профессионального образования

Уфимский государственный авиационный технический университет

Лабораторная работа №5

по дисциплине «Численные методы»

На тему: «Решение задачи Дирихле для уравнения Пуассона

методом Чебышева»

Выполнил:

студент группы ПМ-335

Ямилев И.М.

Проверил:

Голичев И.И.

Уфа

2012

Отчёт по лабораторной работе № 5.

 

Задача:

Явным методом Чебышева требуется найти приближённое решение уравнения

                                                                                  (1)

в квадрате с краевыми условиями

                                                           ,                                                        (2)

где – граница квадрата .

Выбираем функцию, удовлетворяющую краевым условиям (2)

.

Вычислим

.

Возьмём по определению в качестве правой части уравнения (1)

,

тогда нам известно точное решение задачи (1), (2).

Теория:

От задачи (1), (2) перейдём к разностной. Вводим на плоскости прямоугольную сетку с шагом по направлению и по направлению . Получим , . Обозначим .

Обозначим через множество внутренних узлов сетки, то есть

,

а через – множество граничных узлов, то есть

.

Пусть далее

Рассмотрим конечномерное пространство функции , заданных на сетке . Здесь и будем обозначать . Обозначим

.

Тогда разностный оператор Лапласа записывается в виде

                                               .                                        (3)

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

Заменим исходную задачу разностной задачей. При этом будем считать, что , тогда . Разностная аппроксимация задачи (1), (2), принимает вид

                                               ,                                        (4)

или более подробно

                               ,                       (5)

.

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

                                                   .                                      (6)

В пространстве определим оператор

                                                .                                       (7)

Тогда уравнение (5) можно записать в операторной форме

                                                               ,                                               (8)

где – функция, заданная на сетке и . Сеточные функции и будем рассматривать как вектора – мерного пространства с координатами .

Наименьшее  и наибольшее собственные значения оператора равны

,

                 (9)

.

Разностную задачу (5) будем решать явным итерационным методом с чебышевским набором параметров, который выражается следующей формулой:

                                                   ,                                   (10)

где , -заданное число итераций,

             .         (11)


Результаты:

В вычислениях использовался следующий алгоритм:

  1.  Задаём количество итераций , полагаем , тогда шаг сетки =0,1.
  2.  По формулам (9), (11) вычисляем , .
  3.  Вычисляем и по формулам (11).
  4.  Полагая , последовательно применяя формулу (10), находим .
  5.  Пункт 4 повторяем, полагая
  6.  Итерационный процесс продолжаем до совпадения первых четырех знаков в последних итерациях по циклам.

Полученный ответ с точностью до четвертой цифры получен после 50 итераций:

Точка

Точное значение функции

Значение после 40 итераций

Значение после 50 итераций

(0.1, 0.1)

0.0081

0.0080999

0.0081

(0.1, 0.2)

0.0144

0.0144

0.0144

(0.1, 0.3)

0.0189

0.0189

0.0189

(0.1, 0.4)

0.0216

0.0216

0.0216

(0.1, 0.5)

0.0225

0.0225

0.0225

(0.1, 0.6)

0.0216

0.0216

0.0216

(0.1, 0.7)

0.0189

0.0189

0.0189

(0.1, 0.8)

0.0144

0.0144

0.0144

(0.1, 0.9)

0.0081

0.0080999

0.0081

(0.2, 0.1)

0.0144

0.0144

0.0144

(0.2, 0.2)

0.0256

0.0256

0.0256

(0.2, 0.3)

0.0336

0.033599

0.0336

(0.2, 0.4)

0.0384

0.038399

0.0384

(0.2, 0.5)

0.04

0.039999

0.04

(0.2, 0.6)

0.0384

0.038399

0.0384

(0.2, 0.7)

0.0336

0.033599

0.0336

(0.2, 0.8)

0.0256

0.0256

0.0256

(0.2, 0.9)

0.0144

0.0144

0.0144

(0.3, 0.1)

0.0189

0.0189

0.0189

(0.3, 0.2)

0.0336

0.033599

0.0336

(0.3, 0.3)

0.0441

0.044099

0.0441

(0.3, 0.4)

0.0504

0.050399

0.0504

(0.3, 0.5)

0.0525

0.052499

0.0525

(0.3, 0.6)

0.0504

0.050399

0.0504

(0.3, 0.7)

0.0441

0.044099

0.0441

(0.3, 0.8)

0.0336

0.033599

0.0336

(0.3, 0.9)

0.0189

0.0189

0.0189

(0.4, 0.1)

0.0216

0.0216

0.0216

(0.4, 0.2)

0.0384

0.038399

0.0384

(0.4, 0.3)

0.0504

0.050399

0.0504

(0.4, 0.4)

0.0576

0.057599

0.0576

(0.4, 0.5)

0.06

0.059999

0.06

(0.4, 0.6)

0.0576

0.057599

0.0576

(0.4, 0.7)

0.0504

0.050399

0.0504

(0.4, 0.8)

0.0384

0.038399

0.0384

(0.4, 0.9)

0.0216

0.0216

0.0216

(0.5, 0.1)

0.0225

0.0225

0.0225

(0.5, 0.2)

0.04

0.039999

0.04

(0.5, 0.3)

0.0525

0.052499

0.0525

(0.5, 0.4)

0.06

0.059999

0.06

(0.5, 0.5)

0.0625

0.062499

0.0625

(0.5, 0.6)

0.06

0.059999

0.06

(0.5, 0.7)

0.0525

0.052499

0.0525

(0.5, 0.8)

0.04

0.039999

0.04

(0.5, 0.9)

0.0225

0.0225

0.0225

(0.6, 0.1)

0.0216

0.0216

0.0216

(0.6, 0.2)

0.0384

0.038399

0.0384

(0.6, 0.3)

0.0504

0.050399

0.0504

(0.6, 0.4)

0.0576

0.057599

0.0576

(0.6, 0.5)

0.06

0.059999

0.06

(0.6, 0.6)

0.0576

0.057599

0.0576

(0.6, 0.7)

0.0504

0.050399

0.0504

(0.6, 0.8)

0.0384

0.038399

0.0384

(0.6, 0.9)

0.0216

0.0216

0.0216

(0.7, 0.1)

0.0189

0.0189

0.0189

(0.7, 0.2)

0.0336

0.033599

0.0336

(0.7, 0.3)

0.0441

0.044099

0.0441

(0.7, 0.4)

0.0504

0.050399

0.0504

(0.7, 0.5)

0.0525

0.052499

0.0525

(0.7, 0.6)

0.0504

0.050399

0.0504

(0.7, 0.7)

0.0441

0.044099

0.0441

(0.7, 0.8)

0.0336

0.033599

0.0336

(0.7, 0.9)

0.0189

0.0189

0.0189

(0.8, 0.1)

0.0144

0.0144

0.0144

(0.8, 0.2)

0.0256

0.0256

0.0256

(0.8, 0.3)

0.0336

0.033599

0.0336

(0.8, 0.4)

0.0384

0.038399

0.0384

(0.8, 0.5)

0.04

0.039999

0.04

(0.8, 0.6)

0.0384

0.038399

0.0384

(0.8, 0.7)

0.0336

0.033599

0.0336

(0.8, 0.8)

0.0256

0.0256

0.0256

(0.8, 0.9)

0.0144

0.0144

0.0144

(0.9, 0.1)

0.0081

0.0080999

0.0081

(0.9, 0.2)

0.0144

0.0144

0.0144

(0.9, 0.3)

0.0189

0.0189

0.0189

(0.9, 0.4)

0.0216

0.0216

0.0216

(0.9, 0.5)

0.0225

0.0225

0.0225

(0.9, 0.6)

0.0216

0.0216

0.0216

(0.9, 0.7)

0.0189

0.0189

0.0189

(0.9, 0.8)

0.0144

0.0144

0.0144

(0.9, 0.9)

0.0081

0.0080999

0.0081


 

А также другие работы, которые могут Вас заинтересовать

36842. ИЗУЧЕНИЕ И ПОВЕРКА МАНОМЕТРОВ 298 KB
  Давлением называется физическая величина характеризующая интенсивность нормальных распределенных сил с которыми одно тело действует на поверхность другого. Если силы распределены вдоль поверхности равномерно то давление на любую часть поверхности определяется следующим образом: P=G F 1 где F площадь поверхности G сумма приложенных сил. Давление в один ньютон на квадратный метр в системе СИ получила название Паскаль [Па]. = 1013 бар = 0101 МПа Классификация манометров По принципу действия манометры делятся на : Жидкостные...
36843. ПОСТРОЕНИЕ ДИАГРАММ В MS EXCEL 318.5 KB
  Освоить технологии построения диаграммы различных типов. Научиться работать с компонентами диаграммы и настраивать параметры диаграммы. Задания для выполнения и методические рекомендации: С помощью Microsoft Excel можно создавать сложные диаграммы для данных рабочего листа. Прежде чем начать построение диаграммы рассмотрим два важных определения.
36844. Основные определения и критерии классификации угроз 223.2 KB
  Потенциальные злоумышленники называются источниками угрозы. Нарушение безопасности это реализация угрозы. Естественные угрозы это угрозы вызванные воздействием на АС объективных физических процессов стихийных природных явлений не зависящих от человека. Искусственные делят на: непреднамеренные совершенные по незнанию и без злого умысла из любопытности или халатности преднамеренные Каналы проникновения в систему и их классификация: По способу: прямые косвенные По типу основного средства для реализации угрозы: человек...
36845. Подготовка грунтовой площадки к строительству 570.5 KB
  Свойства и технологические характеристики грунтов Любое здание или инженерное сооружение возводится на подстилающем слое грунта. От физикомеханических свойств подстилающего слоя грунта зависит величина осадочных деформаций и долговечность сооружения в целом. К скальным однородным грунтам относят массивы изверженных пород с кристаллической структурой которые характеризуются значительной плотностью и малой влагоемкостью. К скальным слоистым грунтам относят породы сложенные из песчаников доломитов и глинистых сланцев.
36846. КОМПЬЮТЕРНАЯ СИСТЕМА PROJECT EXPERT. АНАЛИЗ ФИНАНСОВЫХ РЕЗУЛЬТАТОВ ПРОЕКТА 64 KB
  Нижняя граница обусловлена тем что оборотных средств должно быть достаточно для погашения краткосрочных обязательств иначе компания окажется под угрозой банкротства. Превышение оборотных средств над краткосрочными обязательствами более чем в три раза также является нежелательным поскольку свидетельствует о нерациональной структуре активов. Показывает отношение наиболее ликвидной части оборотных средств денежных средств дебиторской задолженности краткосрочных финансовых вложений к краткосрочным обязательствам. Чистый оборотный капитал...
36847. Массивы и матрицы. Решение задач линейной алгебры 121.5 KB
  9000 Ввод элементов матрицы также осуществляется в квадратных скобках при этом элементы строки отделяются друг от друга пробелом или запятой а строки разделяются между собой точкой с запятой: nme=[x11 x12 . xmn;] Обратиться к элементу матрицы можно указав после имени матрицы в круглых скобках через запятую номер строки и номер столбца на пересечении которых элемент расположен: nmeиндекс1 индекс2 Листинг 3. Пример обращения к элементам матрицы =[1 2 3;4 5 6;7 8 9] = 1 2 3 4 5 6 7 8 9 12^22 33 ns = 3.
36848. Система автоматического регулирования температуры 488 KB
  Лабораторная работа Система автоматического регулирования температуры. Система автоматического регулирования температуры. Цель работы: Ознакомление с принципами построения системы автоматического регулирования и принципами работы такой системы. Экспериментальное получение переходных процессов системы автоматического регулирования.
36849. Логическая организация оперативной памяти 236.2 KB
  Определить объем основной памяти 2.Определить объем дополнительной памяти 3.Определить объем отображаемой памяти 4.
36850. КОНСОЛИДАЦИЯ ДАННЫХ В MS EXCEL 421 KB
  Создайте три однотипные таблицы по образцу на одном листе или на разных листах MS Excel рис. Проведите консолидацию 3х таблиц аттестации в одну с вычислением среднего балла по каждому предмету и разместите консолидированную таблицу на листе Консолидация для чего: перейдите на чистый лист в книге и установите маркер мыши в левый верхний угол будущей таблицы; на панели Данные выберите Консолидация; в окне Консолидация рис. 2 Диалоговое окно Консолидация перейдите в строку Ссылка затем выделите на листе Данные для консолидации...