main (1160446), страница 14

Файл №1160446 main (Численные методы. Ионкин (методичка) (2015) (LaTeX source)) 14 страницаmain (1160446) страница 142019-09-19СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

также [11]) известно, что в такойпостановке существует единственное решение (, ), которое непрерывно зависит от правойчасти уравнения (, ), начального условия 0 () и краевых условий (2).Чтобы решить эту задачу численно, поставим ей в соответствие разностную схему, тоесть дискретный аналог рассматриваемого уравнения и дополнительных условий.

Такимобразом мы сведем непрерывную задачу к конечной системе линейных уравнений, которыеуже можно решать с использованием вычислительных машин.Сначала введем в рассматриваемой области равномерную по переменным и сетку.Сеткой в заданной области называется совокупность конечного числаточек, принадлежащих данной области. Эти точки называются узлами сетки.Определение.В частности, равномерная сетка размера ( − 1) × , , ∈ N в рассматриваемойобласти вводится так:{︁}︁{︀}︀ℎ = = ℎ, = 1, ( − 1) , = = , = 1, ,1> 0, => 0.Величину ℎ назовем шагом по переменной , величину — шагом по времени.Тогда множество точек ℎ = × ℎ ⊂ ℎ=задает равномерную сетку с шагом ℎ по переменной и шагом по времени в области .Эта сетка изображена на рисунке.Tℎ1Аналогичным образом введем равномерную сетку размера ( + 1) × ( + 1) на замыкании области с теми же размерами шагов ℎ и по переменной и по переменной соответственно.

Эту сетку задает множество точек ℎ = × ℎ ⊂ = {(, ) : ∈ [0, 1], ∈ [0, ]},где{︀}︀{︀}︀ ℎ = = ℎ, = 0, , = = , = 0, .В дальнейшем везде, где мы рассматриваем уравнение теплопроводности, будем использовать введенные сетки, если не указано иное.82В общем случае сетки могут иметь более сложную структуру, например,использовать переменный шаг, который зависит от расположения конкретной пары узлов, или для многомерной области иметь более сложную структуру расположения узлов относительно друг друга (в рассматриваемом примере равномерная сетка являетсяпрямоугольной). В последнее время часто используются сетки, автоматически подстраивающиеся под решение конкретной задачи.Замечание.Совокупность всех узлов в фиксированный момент времени называется слоем.

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

Введем обозначения для ( , ) ∈ ℎ : = ( , ), = ( , ).Обозначим численное решение задачи через( , ) = ,( , ) ∈ ℎ .Здесь ( , ) является сеточной функцией, заданной на сетке ℎ .Поставим в соответствие производным функции (, ) их дискретные аналоги дляфункции ( , ): +1 − ( , )≈ , − 2 + −1 2 ( , )+1≈.2ℎ2В результате получаем дискретный аналог уравнения (1): − 2 + +1+1 − = −1+ ( , ),ℎ2( , ) ∈ ℎ .(4)§26.

Явная разностная схема. Погрешность, сходимость, устойчивость83Запишем дискретные аналоги краевых условий первого рода (2) и начального условия (3):{︃0+1 = 1 (+1 )(5)+1 ∈ ,+1= 2 (+1 ),0 = 0 ( ), ∈ ℎ .(6)Дискретным аналогом задачи (1) – (3), или ее разностной схемой, называется система линейных уравнений (4) – (6).Определение.+1В первой краевой задаче численные значения решения 0+1 и равнызначениям функций 1 () и 2 () соответственно при = +1 (хотя это и не обязательно).

В случае краевых условий иного типа, аппроксимация краевых условий должна бытьсогласована по порядку погрешности с порядком аппроксимации уравнения. Определениеаппроксимации и порядка погрешности аппроксимации будет дано ниже.Замечание 1.Заметим, что в уравнении (4) значения функции (, ) не обязательно брать именно в узлах рассматриваемой сетки, можно использовать значения этойфункции с некоторой «поправкой». Что именно имеется в виду под «поправкой», будетрассмотрено далее, а также будет показано, что выбор значений функции (, ) дляразностной схемы, использующих такую «поправку», позволит получить более высокийпорядок погрешности аппроксимации, а стало быть и более точное решение исходногоуравнения.Замечание 2.Замечание 3. Качество и скорость решения численной задачи (4) – (6) во многом зависит от выбора числа узлов сетки ℎ : чем меньше узлов в сетке, тем меньше уравненийсодержится в системе, тем проще и быстрее ее решать, но и приближение решенияисходной задачи в этом случае будет более грубым.При изучении разностных схем возникают следующие вопросы:1.

Погрешность аппроксимации на решении (невязка).Каждой задаче может быть сопоставлено бесконечное число разностных схем, оценка погрешности аппроксимации позволяет их сравнивать. Разностная схема должнааппроксимировать исходную дифференциальную задачу. Если же аппроксимация отсутствует, то не будет сходимости решения численной задачи к решению исходнойзадачи, и рассмотрение такой разностной схемы не имеет смысла.2. Существование и единственность решения разностной задачи.Построенная разностная задача должна быть корректной, то есть должно существовать единственное решение. В ряде случаев доказательство существования и единственности решения является нетривиальной задачей.3.

Алгоритм нахождения разностного решения.В разностных схемах матрица системы линейных уравнений как правило содержитбольшое число нулей. Для таких систем существуют более эффективные алгоритмырешения, чем универсальный метод Гаусса, например, для систем с трехдиагональнойматрицей разумно использовать метод прогонки.4. Сходимость разностной схемы.Необходимо изучить условия, при которых решение данной разностной схемы сходится к точному решению исходной задачи с наперед заданной точностью.845.

Устойчивость разностной схемы.Устойчивость в данном контексте является чисто внутренним свойством разностныхсхем: разностная схема называется устойчивой в норме ‖·‖, если выполнена априорнаяоценка‖‖ 6 ‖ ‖,где > 0 — константа, не зависящая от шагов сетки.Для построения разностной схемы, обладающей хорошими свойствами, необходимо изучитьвесь круг перечисленных проблем.Вопросы сходимости и устойчивости разностной схемы являются ключевыми, однако обычно достаточно рассмотреть только один из этих двух вопросов: в концекурса будет доказано, что из устойчивости разностной схемы следует ее сходимость крешению исходной задачи при условии, что разностная схема аппроксимирует исходнуюзадачу.Замечание.Совокупность узлов, которые участвуют в записи разностной схемы,называют шаблоном.Определение.Вернемся к изучению явной разностной схемы (4) – (6).В рассматриваемой разностной схеме использован четырехточечный шаблон, схематично изображенный на рисунке.+1−1+1Для построенной разностной схемы решение на ( + 1)-м слое находится явно, поэтомуи рассматриваемая разностная схема называется явной:+1 = + (− 2 + +1) + , = 1, ( − 1),ℎ2 −1{︃0+1 = 1 (+1 )+1 ∈ ,+1= 2 (+1 ),0 = 0 ( ), = 0, .Представленные явные формулы позволяют утверждать, что решение разностной схемы(4) –(6) существует и единственно, значит, мы получили положительный ответ на вопрос 2.Перейдем к исследованию оставшихся вопросов.

Как мы уже упоминали в главе «Интерполирование и приближение функций», существует два подхода к измерению близоститочного решения задачи (1) – (3) (непрерывной функции) и численного решения задачи (4) –(6) (сеточной функции):1. Спроектировать непрерывную функцию (, ) на дискретное пространство и измерять близость функций (, ) и в норме дискретного пространства.2. С помощью интерполирования восполнить функцию до непрерывной и сравниватьрассматриваемые функции в пространстве непрерывных функций.§26. Явная разностная схема.

Погрешность, сходимость, устойчивость85В этом курсе будем пользоваться первым подходом. Под проекцией функции (, ) непрерывных аргумнтов (, ) будем понимать сеточную функцию = ( , ), определенную вузлах сетки ℎ . Обзоначим через точное значение решения (, ) дифференциальнойзадачи (1) – (3) в узле ( , ).Определение.Сеточная функция вида(7) = ( , ) = − , ( , ) ∈ ℎназывается погрешностью решения разностной схемы (4) – (6).Выразим = + и подставим это выражение в разностную схему.

Получим системууравнений для , аналогичную разностной схеме, но с нулевыми краевыми условиями инулевой начальной функцией: − 2 + −1+1 − +1=+ ,ℎ2+10+1 = = 0,0 = 0,Здесь =(8)( , ) ∈ ℎ ,(9)+1 ∈ , ∈ ℎ .−1 − 2 + +1 +1− −+ ,ℎ2(10)(11)Определение. Сеточная функция, задаваемая равенством (11) называется погрешностью аппроксимации разностной схемы (4) – (6) на решении исходной задачи.(︀)︀2Задача. Доказать, что = O + ℎ .Здесь и далее ( , ) ∈ ℎ , = 0, , = 0, .

Далее всюду при использованииформулы Тейлора будем предполагать, что разлагаемая функция обладает нужной гладкостью, то есть имеет непрерывные производные до соответствующего по ходу разложенияпорядка.Разложим ( , +1 ) в узле ( , ) по формуле Тейлора:Решение.( , +1 ) = +1= ( , ) + ′ ( , ) + ( 2 ).Разложим (+1 , ) в узле ( , ) по формуле Тейлора:11(+1 , ) = +1 = ( , ) + ′ ( , )ℎ + ′′ ( , )ℎ2 + ′′′( , )ℎ3 + (ℎ4 ).26 Разложим (−1 , ) в узле ( , ) по формуле Тейлора :11(−1 , ) = −1 = ( , ) − ′ ( , )ℎ + ′′ ( , )ℎ2 − ′′′( , )ℎ3 + (ℎ4 ).26 Полученные разложения подставим(︀)︀в формулу (11) и после приведения подобных слагаемых получим оценку = O + ℎ2 .Введем норму в пространстве сеточных функций на -м слое, = 0, :‖ ‖ = max | |.066Мы рассматриваем решение разностной задачи по слоям, поэтому нет необходимости вводить норму как максимум модуля для всех слоев.86Пусть решение (, ) задачи (1) – (3) обладает достаточной гладкостью (четыре раза дифференцируема по и два раза по ).

Тогда для сходимости решения разностной схемы (4) – (6) к решению исходной задачи (1) – (3) в норме ‖·‖ необходимо идостаточно, чтобы выполнялось условие:Теорема.=6 0.5.ℎ2При этом условии, выполняется оценка:⃦⃦ +1(︀)︀⃦− +1 ⃦ 6 1 + ℎ2 , = 0, 1, . . .где 1 > 0 — константа, не зависящая от и ℎ.Докажем, что выполнения условий теоремы достаточно для сходимостиразностной схемы к решению исходной задачи.Запишем выражение для +1 в видеДоказательство.(︀ )︀+1 = (1 − 2) + −1+ +1+ и оценим левую часть равенства по модулю с учетом условия 1 − 2 > 0.

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

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

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

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

Численные методы
.gitignore
README.md
biblio.bib
circle.eps
main.aux
main.bbl
main.blg
main.log
main.out
main.synctex.gz
main.tex
main.toc
msu.eps
utf8gost705u.bst
Свежие статьи
Популярно сейчас
Почему делать на заказ в разы дороже, чем купить готовую учебную работу на СтудИзбе? Наши учебные работы продаются каждый год, тогда как большинство заказов выполняются с нуля. Найдите подходящий учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6363
Авторов
на СтудИзбе
310
Средний доход
с одного платного файла
Обучение Подробнее