Главная » Просмотр файлов » Численные методы. Ионкин (2013)

Численные методы. Ионкин (2013) (1160444), страница 18

Файл №1160444 Численные методы. Ионкин (2013) (Численные методы. Ионкин (2013)) 18 страницаЧисленные методы. Ионкин (2013) (1160444) страница 182019-09-19СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 18)

Запишем попеременно-треугольный итерационный метод:( + 1 )( + 2 ) (+1) − ()+ () = , ∈ ℎ , ∈ Z+ , (0) — задано,где действие матриц 1 и 2 на сеточную функцию определяется равенствами(1 ) = − −1, − ,−1+,2ℎ1ℎ22(2 ) = − +1, − ,+1+.2ℎ1ℎ22110Глава 4. Разностные методы решения задач математической физикиЗаметим, воспользовавшись результатами, полученными в §8 главы 1, что этот итерационный метод сходится при > 4 , а каждую следующую итерацию данного метода, как и врассмотренных ранее методах, можно вычислить по явным формулам.В заключение отметим, что при ℎ1 = ℎ2 = ℎ количество итераций 0 , необходимое дляполучения решения с заданной точностью > 0 обратно пропорционально шагу ℎ:(︀)︀0 () = O ℎ−1 .Таким образом, попеременно-треугольный итерационный метод на порядок быстрее методов Якоби и Зейделя и тем самым является эффективным методом, широко применяемымв настоящее время при численном решении разностной задачи Дирихле.§9Основные понятия теории разностных схем: аппроксимация, устойчивость, сходимостьРассмотрим задачу:(()) = (), ∈ ,(1)где – линейный дифференциальный оператор, – произвольная, в общем случае многомерная, область.Произвольную линейную дифференциальную задачу, в том числе, содержащую краевыеи начальные условия, можно привести к виду (1).Мы рассматриваем алгоритмы решения только для корректно поставленных задач.

Задача (1) называется корректно поставленной, если1. Решение () рассматриваемой задачи существует и единственно.2. Решение () непрерывно зависит от входных данных, к которым относится праваячасть уравнения, краевые и начальные условия.Построим в области сетку ℎ с некоторым обобщенным шагом ℎ = max ℎ — максимальным среди всех используемых на сетке ℎ шагов ℎ . Заметим, что при ℎ → 0 числоузлов в данной сетке стремится к бесконечности.Поставим в соответствие непрерывным функциям () и () их дискретные аналогиℎ () и ℎ (), ∈ ℎ соответственно. Пусть также оператор ℎ является дискретным аналогом оператора .

Тогда поставим в соответствие исходной дифференциальной задаче еедискретный аналог:ℎ ℎ () = ℎ (), ∈ ℎ .(2)Замечание. Заметим, что дискретных аналогов каждого дифференциального уравнениявида (1) существует бесконечно много, поэтому при выборе такого аналога необходиморуководствоваться потребностями конкретной задачи и выбирать приближение оптимальным для данного случая образом.Определение. Дискретная задача вида (2) называется разностной схемой для уравнения (1).Замечание.

Заметим, что разностная схема является, фактически, системой линейных алгебраических уравнений, и при малых значениях шага ℎ является системой высокихпорядков.§9. Основные понятия теории разностных схем111При рассмотрении приближений заданных в непрерывной области задач их дискретными аналогами встает вопрос о том, каким образом измерять близость решений обеих задач.Существует два способа измерения близости решений уравнений (1) и (2).Введем два линейных пространства: 0 – пространство непрерывных функций, котороеудовлетворяет уравнению (1), и ℎ – пространство дискретных функций, удовлетворяющихуравнению (2).

Близость функций из пространств 0 и ℎ можно измерять:1. В норме пространства ℎ .Пусть ℎ – оператор проектирования пространства 0 на пространство ℎ . Тогдаблизость функций () и ℎ () будем измерять в норме ‖·‖ℎ пространства ℎ . Если‖ℎ (()) − ℎ ()‖ℎ → 0,то в этой норме ℎ () сходится к ().2. В норме пространства 0 с помощью восстановления дискретной функции ℎ () донепрерывной функции пространства 0 .Замечание. При рассмотрении норм в пространствах 0 и ℎ необходимо учесть, чтонормы в этих пространствах должны быть согласованными, то есть должен существовать пределlim ‖ℎ (())‖ℎ = ‖‖0 ,(3)ℎ→0где ‖·‖0 — норма в пространстве 0 . Как будет показано после доказательства теоремыФилиппова, согласованность норм гарантирует сходимость решения разностной схемыименно к решению исходного уравнения — если согласованность не выполнена, мы можемполучить сходящуюся разностную схему, но сходиться она будет не к решению исходногоуравнения.Пример.

Рассмотрим область = { : 0 6 6 1} и зададим на этой области равномернуюсетку с шагом ℎ и числом узлов, равным ∈ N:ℎ = { = ℎ, = 0, , ℎ = 1}.Введем нормы в пространствах 0 и ℎ . В пространстве 0 :‖‖0 = ‖‖ = max |()|,0661() ∈ 0 .В пространстве ℎ :∀ℎ = (0 , 1 , . . . , ) ∈ ℎ‖ℎ ‖ℎ = ‖ℎ ‖ = max | |.066Для введенных таким образом норм выполняется условие согласованности.Рассмотрим еще несколько примеров:1. Пусть 0 — пространство функций, интегрируемых с квадратом. Введем норму⎛ 1⎞ 12∫︁‖‖0 = ⎝ 2 ()⎠ = ‖‖2 , () ∈ 0 .0Тогда рассмотрим в дискретном пространстве ℎ следующую норму:(︃ )︃ 21∑︁2 ℎ‖ℎ ‖ℎ == ‖ℎ ‖2 (ℎ ) , ℎ ∈ ℎ .=0Для этих норм выполнено условие согласованности.112Глава 4.

Разностные методы решения задач математической физики2. Введем норму в дискретном пространстве ℎ(︃‖ℎ ‖ℎ =∑︁)︃ 212, ℎ ∈ ℎ .=0Докажем, что эта норма не согласована ни с одной нормой пространства 0 .Пусть () ∈ 0 : () ≡ 1. Тогда(︃ )︃ 21∑︁√‖ℎ (())‖ℎ == 2 + 1.1=0Следовательно, при ℎ → 0lim ‖ℎ (())‖ℎ = ∞.ℎ→0Очевидно, что ‖‖0 не может равняться бесконечности, так как норма должна бытьконкретным числом. Значит норма ‖·‖ℎ не согласована ни с одной нормой пространства 0 .Рассмотрим подробнее оператор проектирования ℎ пространства 0 на пространствоℎ . Ранее в этом параграфе мы предполагали, что этот оператор задан следующим образом:(ℎ ())( ) = ( ), ∈ ℎ ,() ∈ 0 .Однако оператор проектирования можно ввести бесконечным числом способов.

Например,оператор проектирования1(ℎ ())( ) =ℎ∫︁+0.5ℎ(),() ∈ 0 , −0.5ℎзадает среднее значение функции () в узле , = 1, . Значения оператора в граничныхточках области ℎ определяются следующим образом:1(ℎ ())(0 ) =0.5ℎ0.5ℎ∫︁(),0∫︁11(ℎ ())( ) =0.5ℎ().1−0.5ℎПерейдем к основным понятиям теории разностных схем.Определение. Сеточная функцияℎ () = ℎ () − ℎ () = ℎ () − ℎ (()),ℎ () ∈ ℎ , () ∈ 0 ,(4)называется погрешностью разностной схемы (2).Подставим выражение ℎ () = ℎ () + ℎ () в разностную схему (2) и получим разностную схему относительно погрешности ℎ ():ℎ ℎ () = ℎ (), ∈ ℎ ,(5)гдеℎ () = ℎ () − ℎ ℎ ().(6)§9.

Основные понятия теории разностных схем113Определение. Сеточная функция, задаваемая соотношением (6), называется погрешностью аппроксимации на решении исходной задачи (1).Замечание. Погрешность аппроксимации можно представить в виде суммы погрешности приближения оператора и погрешности приближения правой части.Определение. Разностная схема (2) аппроксимирует исходную задачу, еслиlim ‖ℎ ‖ℎ = 0.ℎ→0Определение. Разностная схема (2) имеет -ый порядок аппроксимации, если существуют положительные константы 1 и , не зависящие от шага ℎ, такие, что‖ℎ ‖ℎ 6 1 ℎдля достаточно малых ℎ.Определение. Разностная задача (2) называется корректно поставленной, если при достаточно малых ℎ выполнено:1.

При любых погрешностях аппроксимации ℎ и при любых правых частях ℎ решениезадачи (2) ℎ () существует и единственно.2. Существует константа 2 , не зависящая от шага ℎ, для которой выполнена априорная оценка:‖ℎ ‖ℎ 6 2 ‖ℎ ‖ℎ .Это оценка означает устойчивость в норме ‖·‖ℎ решения разностной схемы по правой части уравнения.Замечание 1. Свойства существования и единственности решения задачи определяютсуществование оператора −1ℎ .Замечание 2. Второе условие корректности постановки задачи означает равномернуюограниченность по ℎ оператора −1ℎ .Определение. Решение разностной задачи (2) сходится к решению исходной дифференциальной задачи (1), еслиlim ‖ℎ ‖ℎ = lim ‖ℎ − ℎ ‖ℎ = 0.ℎ→0ℎ→0Определение.

Разностная схема (2) имеет -ый порядок точности, если существуетконстанта 3 , не зависящая от шага ℎ, такая, что‖ℎ ‖ℎ 6 3 ℎ .Теорема 1. (Филиппова). Пусть исходная задача (1) и разностная схема (2) поставленыкорректно, и пусть разностная схема аппроксимирует исходную задачу. Тогда решениеразностной задачи сходится к решению исходной задачи, причем порядок точности разностной схемы совпадает с порядком аппроксимации.Доказательство. Рассмотрим задачу для погрешности разностной схемы (5). Так как поусловию разностная схема корректна, то выполнено условие‖ℎ ‖ℎ 6 2 ‖ℎ ‖ℎ ,114Глава 4. Разностные методы решения задач математической физикигде константа 2 не зависит от ℎ.

Заметим, что задача для погрешности ℎ ∈ ℎ тоже является корректно поставленной, значит, для погрешности выполнена аналогичная оценка:‖ℎ ‖ℎ 6 2 ‖ℎ ‖ℎ .Так как разностная схема аппроксимирует исходную задачу, тоlim ‖ℎ ‖ℎ = 0.ℎ→0Пусть разностная схема имеет -ый порядок аппроксимации:‖ℎ ‖ℎ 6 1 ℎ ,где константа 1 не зависит от ℎ. Тогда аналогичная оценка справедлива и для погрешности ℎ :‖ℎ ‖ℎ 6 3 ℎ ,где 3 = 2 1 > 0, не зависит от ℎ. Это и означает, что разностная схема (2) сходится с-ым порядком точности по ℎ.Замечание. При доказательстве теоремы условие согласованности норм (3) не использовалось. Это условие нужно для того, чтобы гарантировать единственность предельнойфункции.Покажем, что возможно отсутствие сходимости решения задачи (2) к исходной задаче (1), если нормы пространств 0 и ℎ не согласованы.

Пустьlim ‖ℎ − ℎ ‖ℎ = 0.ℎ→0Мы хотим выяснить, существует ли какая-либо другая функция () ∈ 0 , для которойlim ‖ℎ − ℎ ‖ℎ = 0.ℎ→0Покажем, что если норма ‖·‖0 пространства 0 согласована с нормой ‖·‖ℎ пространства ℎ ,то функция (), если существует, то тождественно совпадает с функцией (). Оценимнорму ‖ℎ − ℎ ‖ℎ :‖ℎ − ℎ ‖ℎ = ‖ℎ − ℎ − (ℎ − ℎ )‖ℎ 6 ‖ℎ − ℎ ‖ℎ + ‖ℎ − ℎ ‖ℎ .Так какlim ‖ℎ − ℎ ‖ℎ = 0,ℎ→0lim ‖ℎ − ℎ ‖ℎ = 0,ℎ→0тоlim ‖ℎ − ℎ ‖ℎ = 0.ℎ→0В случае согласованности норм отсюда следует, что‖() − ()‖0 = 0,а это значит, что функции () и () совпадают: () ≡ (), и решение разностнойсхемы (2) ℎ однозначно определяет решение исходной задачи (1).Этот факт гарантируется согласованностью норм пространств 0 и ℎ .

Характеристики

Тип файла
PDF-файл
Размер
1005,37 Kb
Тип материала
Высшее учебное заведение

Список файлов лекций

Свежие статьи
Популярно сейчас
А знаете ли Вы, что из года в год задания практически не меняются? Математика, преподаваемая в учебных заведениях, никак не менялась минимум 30 лет. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6381
Авторов
на СтудИзбе
308
Средний доход
с одного платного файла
Обучение Подробнее