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

Н.И. Ионкин - Численные методы (1163659), страница 16

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

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

Разностные методы решения задач математической физики{︃0+1 = 0+1= 0,0 = 0,(8)+1 ∈ ,(9) ∈ ℎ ,где — погрешность аппроксимации на решении:+1+1+ +1+1− −1 − 2+1++ ( , +1 ).= ( , ) = −ℎ2Задача.Доказать, что(10)(11)(︀)︀ = O + ℎ2 .Для оценки погрешности воспользуемся нормой ‖·‖ в пространствесеточных функций на слое, которую мы ввели в предыдущем параграфе.Теорема.

Пусть функция (, ) имеет достаточную гладкость (четырераза дифференцируема по и два раза по ). Тогда чисто неявная разностнаясхема сходится к решению исходной задачи в норме ‖·‖ с первым порядкомточности по и вторым порядком точности по ℎ.Доказательство. Пусть 0 ∈ ℎ — узел, на котором достигается максимумпогрешности на ( + 1)-м слое:⃒ +1 ⃒⃒⃒ ⃦⃦⃒⃒ = max ⃒ +1 ⃒ = ⃦ +1 ⃦ .0066Для доказательства теоремы воспользуемся, фактически, принципом максимума. Запишем уравнение (7) относительно узла 0 :(︀)︀(1 + 2) +1= 0 + +1+ +1+ 0 ,00 −10 +1=> 0.ℎ2Оценим левую и правую части равенства по модулю с учетом того, что(1 + 2) > 0:⃒⃒ ⃒ ⃒⃒ ⃒⃒)︀⃒ ⃒(︀⃒⃒ 6 ⃒ ⃒ + ⃒ +1 ⃒ + ⃒ +1 ⃒ + ⃒ ⃒ .(1 + 2) ⃒+1−1+100000Перейдем в правой части неравенства от модулей слагаемых к нормам соответствующих функций.

При таком переходе правая часть неравенства можеттолько увеличиться:⃒⃒ ⃦ ⃦⃦⃦⃦ ⃦⃒ 6 ⃦ ⃦ + 2 ⃦ +1 ⃦ + ⃦ ⃦ .(1 + 2) ⃒+10§3. Чисто неявная разностная схема123⃒⃒⃒ =Так как по предположению ⃒+10имеет вид⃦⃦⃦ ⃦(1 + 2) ⃦ +1 ⃦ 6 ⃦ ⃦⃦ +1 ⃦⃦⃦ , то полученное неравенство⃦⃦⃦ ⃦+ 2 ⃦ +1 ⃦ + ⃦ ⃦ .Отсюда следует, что⃦ +1 ⃦⃦ ⃦⃦ ⃦⃦⃦ 6 ⃦ ⃦ + ⃦ ⃦ .Из этого неравенства вытекает:∑︁⃦ +1 ⃦⃦ ⃦⃦ ⃦⃦⃦ 6 ⃦ 0 ⃦ + ⃦ ⃦ .=0Учитывая, что начальная погрешность равна нулю, получаем оценку∑︁⃦ +1 ⃦⃦ ⃦⃦⃦ 6⃦ ⃦ .=0Из (11) следует, что⃦ ⃦(︀)︀⃦ ⃦ 6 + ℎ2 ,где > 0 — константа, не зависящая от и ℎ, и∑︁ = +1 6 .=0Таким образом получим окончательную оценку:(︀)︀‖ +1 ‖ 6 1 + ℎ2 ,(12)где 1 = > 0 — константа, не зависящая от и ℎ.

Устремив и ℎк нулю, получим:⃦⃦lim ⃦ +1 − +1 ⃦ = 0. →0ℎ→0Равенство предела разности нулю означает, что решение разностной схемысходится к решению исходной задачи.Наличие оценки (12) означает, что схема имеет первый порядок точностипо и второй — по ℎ. Чисто неявная схема является абсолютно сходящейсяразностной схемой, так как оценка (12) была получена без всяких ограничений на и ℎ.124Глава IV . Разностные методы решения задач математической физикиЗамечание.вияЕсли в разностной задаче (4) – (6) взять нулевые краевые усло+10+1 = = 0,то для можно вывести оценку, аналогичную полученной выше:∑︁⃦ ⃦⃦ ⃦⃦ +1 ⃦⃦ ⃦ .⃦ 6 ⃦0 ⃦ + ⃦=0Эта оценка означает, что решение разностной схемы устойчиво по начальному условию и по правой части уравнения.§4Симметричная разностная схема. Задача на собственные значения.

Сходимость, устойчивость внорме 2 (ℎ )Рассмотрим уравнение теплопроводности с краевыми условиями первого рода:(, ) 2 (, )=+ (, ), (, ) ∈ = {(, ) : ∈ (0, 1), ∈ (0, ]}, (1)2{︃(0, ) = 1 () ∈ [0, ],(2)(1, ) = 2 (),(, 0) = 0 (), ∈ [0, 1].(3)Воспользуемся сетками ℎ и ℎ , введенными в первом параграфе даннойглавы на множествах и соответственно.Введем вторую разностную производную для дискретной функции =( , ), определенной на множестве ℎ :,= − 2 + −1+1.2ℎЭта производная является дискретным аналогом второй производной по функции (, ).Поставим в соответствие уравнению (1) его дискретный аналог в виде+1,+ ,+1 − =+ ( , + 1 ),22(︀)︀где ( , + 1 ) = , + 2 ∈ ℎ .2(4)§4. Симметричная разностная схема.

Задача на собственные значения. Сходимость, устойчивостьв норме 2 ( )Определение.125Слой + 1 = +22называется полуцелым слоем.Добавим краевые и начальное условия:{︃0+1 = 1 (+1 )+1 ∈ ,+1= 2 (+1 ),0 = 0 ( ), ∈ ℎ .(5)(6)В рассматриваемой разностной схеме использован шеститочечный шаблонвида−1+1−1+1+1Заметим, что данная схема, с точки зрения нахождения численного решения, похожа на ту, которую мы рассматривали в предыдущем параграфе, вчастности, матрица системы, соответствующей этой схеме, является трехдиагональной со строгим диагональным преобладанием. Это значит, что решениеразностной схемы (4) – (6) такой задачи существует, единственно и находитсяс помощью метода прогонки.Введем погрешность решения разностной схемы: = ( , ) = − ,где = ( , ), ( , ) ∈ ℎ .Выразив из этого выражения и подставив его в уравнение (4), получимзадачу относительно :+1,+ ,+1 − =+ , ( , ), ( , +1 ) ∈ ℎ ,2{︃0+1 = 0+1 ∈ ,+1= 0,0 = 0,где ∈ ℎ ,(7)(8)(9)— погрешность аппроксимации на решении исходной задачи (1)–(3): = ( , ) = −+1+1− , + ,++ ( , + 1 ),22( , ), ( , +1 ) ∈ ℎ .(10)126Глава IV .

Разностные методы решения задач математической физикиЗадача.Доказать, что(11)(︀)︀ = O 2 + ℎ2 .Переходим к изучению вопросов сходимости и устойчивости разностнойзадачи (4)–(6).Рассмотрим вещественное пространство −1 сеточных функций , заданных на одномерной сетке ℎ , содержащей ( − 1) узел и обращающихсяв нуль на границе (0 = = 0).Значение функции ∈ −1 в -м узле сетки, = 1, ( − 1), обозначимчерез . Заметим, чтоdim −1 = − 1.Введем скалярное произведение и норму в пространстве −1 :(, ) =−1∑︁ ℎ, ‖‖2 (ℎ ) =(︃ −1∑︁=1)︃ 212 ℎ,, ∈ −1 ..(12)=1Заметим, что если взять значения сеточной функции , рассматриваемойна сетке ℎ , принадлежащие одному слою, пусть -ому, то эти значения образуют функцию , принадлежащую пространству −1 .

Тогда, если будетверна оценка⃦ +1 ⃦(︀)︀⃦⃦6 2 + ℎ2 ,2 ( )ℎгде константa не зависит от и ℎ, то это будет означать сходимость рассматриваемой разностной схемы к решению исходной задачи в норме 2 (ℎ )со вторым порядком точности по и ℎ.Наряду с вещественным пространством −1 будем рассматривать гильбертово пространство 2 — линейное пространство функций, интегрируемыхс квадратом на интервале (0, 1):∫︁1 2 () < ∞.0Введем скалярное произведение и норму в пространстве 2 :∫︁1⎛∫︁1 ()(), ‖ ‖2 = ⎝(, ) =00⎞ 21 2 ()⎠ , (), () ∈ 2 .§4. Симметричная разностная схема. Задача на собственные значения.

Сходимость, устойчивостьв норме 2 ( )127Задача на собственные значенияРассмотрим задачу на собственные значения (задачу Штурма-Лиувилля) дляфункции () ∈ 2 , обладающей достаточной гладкостью:⎧ 2⎨ + () = 0, ∈ (0, 1),(13)2⎩(0) = (1) = 0,причем () ̸≡ 0.Решениями данной задачи являются собственные значения и собственные функции (): = 2 2 , ∈ N,0 < 1 < 2 < .

. . < < . . . , = ̸= 0. () = sin(),Одним из свойств собственных функций задачи Штурма-Лиувилля являетсятот факт, что эти√функции образуют ортогональный базис пространства 2 .Положим = 2 и получим:√ () = 2 sin().Тогда функции { ()}∞=1 образуют ортонормированный базис в пространстве 2 :( , ) = .Значит, произвольную функцию ()∈2 можно разложить∑︀∞по базису { ()}∞:()=(),гдекоэффициенты = (, )=1 =1называются коэффициентами Фурье.Тогда справедливо равенство Парсеваля:‖ ‖22 =∞∑︁2 .=1Рассмотрим теперь разностный аналог задачи Штурма-Лиувилля для сеточной функции ∈ −1 :{︃, + ( ) = 0, ∈ ℎ , = 1, ( − 1),(14)0 = = 0,128Глава IV .

Разностные методы решения задач математической физикипричем () ̸≡ 0. Будем искать собственные функции в виде( ) = sin( ), ∈ R, = 1, ( − 1).Распишем уравнение (14) подробнее:+1 − 2 + −1+ = 0ℎ2и перенесем слагаемые, содержащие , в правую часть:(︀)︀+1 + −1 = 2 − ℎ2 , = 1, ( − 1).Очевидно, что+1 +−1 = ( +ℎ)+( −ℎ) = sin ( +ℎ)+sin ( −ℎ) = 2 sin( ) cos(ℎ).Следовательно,2 sin( ) cos(ℎ) = (2 − ℎ2 ) sin .sin( ) ̸= 0, так как собственные функции не могут быть нулевыми, значит(︂ )︂ℎ22 ℎ= 1 − cos ℎ = 2 sin.22Отсюда следует, что(︂ )︂42 ℎ = 2 sin.ℎ2Для того, чтобы найти , воспользуемся краевым условием для : = sin = 0,откуда следует, что = , ∈ N. Тогда собственные значения равны)︂(︂42 ℎ, = 1, ( − 1), = 2 sinℎ2а соответствующие им собственные функции имеют вид = sin( ), ∈ 1, ( − 1).Системафункций ( ), = 1, ( − 1) ортогональна,√√ а если положить = 2, то совокупность сеточных функций ( ) = 2 sin( ) образует ортонормированный (в смысле скалярного произведения (12)) базис пространства −1 .

Следовательно, любая сеточная функция ( ),−1 = 1, ( − 1), однозначно разложима по базису { }, то есть1 ( ) =−1∑︁=1 ( ),§4. Симметричная разностная схема. Задача на собственные значения. Сходимость, устойчивостьв норме 2 ( )129где = (, ), = 1, ( − 1) — коэффициенты Фурье. Имеет место равенство Парсеваля:−1∑︁‖ ‖22 (ℎ ) =2 .(15)=1Воспользуемся рассмотренной задачей Штурма-Лиувилля для доказательства следующей теоремы.Пусть функция (, ), являющаяся решением задачи для уравнения теплопроводности (1) – (3), имеет достаточную гладкость. Тогдасимметричная разностная схема (4) – (6) сходится к решению исходной задачи со вторым порядком по и вторым порядком по ℎ в 2 (ℎ )-нормепространства сеточных функций.Теорема.Обратимся к рассмотрению задачи (7) – (9) для погрешности решения разностной схемы :Доказательство.+1+ ,,+1 − =+ ,2{︃0+1 = 0+1= 0,0 = 0,где (︀ — погрешность)︀ = O 2 + ℎ2 :( , ), ( , +1 ) ∈ ℎ ,(17)+1 ∈ ,(18) ∈ ℎ ,аппроксимациина(16)решениизадачи(1) – (3),+1+1− , + ,= ( , ) = −++ ( , + 1 ), ( , ), ( , +1 ) ∈ ℎ .22(19)Будем искать погрешность в виде=−1∑︁ ( ) ( ), ∈ ℎ ,(20)=1где ( ), = 1, ( − 1) — дискретные функции только аргумента , а , =1, ( − 1) — собственные функции задачи, зависящие только от ∈ ℎ :, + ( ) = 0, = 1, ( − 1),0 = = 0.(21)130Глава IV .

Разностные методы решения задач математической физикиЗадача (21) была рассмотрена выше.Функции имеют вид√ ( ) = 2 sin( ),, = 1, ( − 1)и образуют ортонормированный базис в −1 . Эти функции соответствуютсобственным значениям , равным(︂)︂42 ℎ, = 1, ( − 1), = 2 sinℎ2 −1Так как функции { }=1образуют ортонормированный базис пространства −1 , то любой элемент пространства −1 можно разложить по этим функциям, следовательно, представление (20) корректно.Разложим по базисным функциям погрешность аппроксимации на решении:−1∑︁ () ( ) ( ),(22) ==1где ) — дискретные функции только аргумента .Подставим выражения (20) и (22) в уравнение (7): () (−1∑︁ −1( (+1 ) − ( ))1 ∑︁ ( ) =( (+1 ) + ( )) ( ), +2=1−1∑︁+=1 () ( ) ( ).=1Принимая во внимание уравнение (21), получаем−1 (︂∑︁=1)︂−1∑︁ (+1 ) − ( ) + ( (+1 ) + ( )) ( ) = () ( ) ( ).2=1−1Так как { }=1 — система линейно независимых функций, то полученноеравенство выполняется тогда и только тогда, когда коэффициенты при соответствующих функциях ( ), = 1, ( − 1) равны: (+1 ) − ( ) + ( (+1 ) + ( )) = () ( ),2 = 1, ( − 1).Разрешим это уравнение относительно ( + 1)-го слоя, домножив обе частина ̸= 0 и сгруппировав слагаемые с (+1 ) и ( ):(1 + 0.5 ) (+1 ) = (1 − 0.5 ) ( ) + () ( ).§4.

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

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

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

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