Главная » Просмотр файлов » Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu)

Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu) (1160088), страница 98

Файл №1160088 Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu) (Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы) 98 страницаН.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков - Численные методы (djvu) (1160088) страница 982019-09-19СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Таким образом, окончательно имеем (( "((< Ьа((+((о ((< < гпах ~ шал )рз(, шах )рз), ((и ))) -~-Т шах ((97 !). Так как зто неравенства справедливо при любом и, О < и < Ж, то это и означает устойчивость разностной сжмы в сеточной норме пространства С. Теорема доказана. Глава 10. Методы решеевя уравнений е честных прокшодных (Этот вывод можно было бы сделать, непосредственно используя соотношение (9) и пе вводя в рассмотрение функции 9' н ех.] Отметим, что постоянная в (8) получилась в данном случае зависящей аг Т. Если т)1г = ь =- гопя1, то условен а < 1/2 является необходимым и достаточным условием устойчивости. тфз 1~2 Задача 1. Пусть 1ип — =-еш Доказать, что схема (3), (5) нее, † устойчива. Указание Рагсмотреть частные решгни» гггпй и,"„= Л" ып д.

Доказатееьство устойчивости разностной схемы (3), (5) было получено при соотношении па шаги сетки т < 1Р/2. Разиостные схемы, которые обладшог 1ттойчнвостьк) лишь прн определенных стютвошсниях па шаги сетки, называются условно ргглгючвеилш. Соответственно, если схема устойчива при лкгбьгх игггю1пеннях между шагами сетки, то такая схеьш называется безусловно угшоячиеой. Покаягем, что схеыа (4), (5) отвоснтсв к квассу бшусловно устойчивых схем.

С и рвведгп1ва Теорема 2. При любит А, и т длл решеигш зппп ги (4), (5) имсшл место оценке (8). Докиэетельсшео. По аналогии с доказательством предыдущей поремы преобразуем (4) к аиду Из вгвх значений и'™, по модулю равных ))йь~)), всзьмем то, у которого индекс т принимшн наименьшее значение. Если ш = О или же л~. = 54, то (9) выполнено.

Пусть теперь щ отлично от О и М. ПосколькУ )йе, ~) > )и +~ ) (по опР<щеле1гиго гл) и )гг,"„тг) > )и,„~~г), то (2йы~( > )и" Ь1г) + (и,"„~~г(; позтому вйл(2йьг — и"'~~1 — и,"',~~г) = е18пи'„'+г. Тогда )(ив ~ )( = )и"+ ) < )и'„'+ + р( — и ~1 + 2и,„+ — и" ьг|)( = = (и,„+ тр ~ ) < ()и~)) з-тфр + )); при лгобых 5 н т для схемы (4), (5) получена оценка (9). Заверпгение доказательства теоремы совпадыт с доказательсшом предыдущей теоремы. Таким образом, схема (4), (5) является безусловно устойчивой.

5 5. Рвзностные схемы д»я параболическою уравнения 533 (т — 1,л+1) (т,лт1) (т+1,п-г1) (т,л+1) (гл — !,л) (т,п) (ш+1,л) Рис. 10.5.1 (т,л) Рис. 10.5.2 Пусть в — функция, определенная на с»ое и принимающая в т-ы узле значение и . Обсеиачнм Лет = (емы — 2вв 4 ем — ~)Ф ЛЧ(жг) = п(в+ 5, 1) — 2в(к, 1) -~- и(х — Л, 1). Введем более общий, чем (3), (5] и (4), (5), ° ид схемы и" "~ — в" =пЛины 4-(1 — с)Ли +д,"„, гп = 1,..., И вЂ” 1. (11) Между явной (3), (5) и неявной (4), (5) схемамн имеется, таким образом, принципиальное отличие.

Явной схеме соответствуют явные формуЛы для вычисления функции на слав по известным значениям на прщььь дущих слоях. Однако эзв схема являсття ус»овио устойчивой. Это приводит к тому, что при малом шше Ь мы вынуждены выбирать слишком мелкий шаг по времени (г < Ьз/2), чтобы обеспечить устайчивгюгь. Это, в свою очередь, приводит к Значнтельнолзу уселнчениго затрат времони счета на ЭВМ и не может быть оправдано требованю~ыи точности, егин по временной переменной 1 решеане достаточно гладкое. С другой стороны, при использовании поленов стел!ы можно значительна увш!иппь шш по времени г, однако нри переходе от слоя к с»ою требуется каждый раз решать систему уравнений.

Впрочем, н одпоь~ерном случае это не представляет проблемы. В иктпости, испо»ьзуя мотод прогонки, можно получить опт~ при известном вд за 0(М) операций, т.е. количество арифметических операций при переходе ш слюя к слою по порядку будет теы же, что и в с»учае явной схемм. Это позволшг сделать вывод о том, по испозп,»осанне неявных схем в одномерном случае часто и»ляпни более предпочтиъюзьныы, так как ведет к уменьшению затрат времени счвп» па ЭВМ. Перейдем теперь к исследованию устойчивости разностных схем (3), (5) и (4), (б) в других нормах, в частности в сеточных аналогах норм Ьз и И'зг по слою.

Твк как исследование свойств схем (3), (5) и (4), (5) проводится по оцной и той жс мелодике, то имеет смыс» объединить ати две схемы. Как и ранее, для пшлядносги сопоставим развостной схеме ее шаблон. В нашем случае для схем (3), (5) и (4), (5) шаблоны имеюг внд, изображенный соответственно па рис. 10.5.1 и 10.5.2. Глава 10. Метали реюеввя ураваенип в часттнж пропзводвых 534 Постоянная а в (11] иазываегся весам и обычно берется в пределах 0 < гг < 1. В частности, при а = 0 выражение (11) переходит в [3), а прн т = 1 получаем (4).

Оп-1,п) (т,п) (т+1.п) Раанастную схему (11), (5) называют гтемой с еепьии. Она имеет ш~, ститочечный шаблон при о е (О, 1) (рис. 10.5.3). Схт1а (11), (5) явлнРслнастную схему (4), (5), чтобы отличить вггда (11) с 0 < о < 1, аачыюаот тзтпа (т-1,пе1) [т.п+1) (па+1,л+1) ется явной лишь прн а = О. ее ат других неявных схелг пеленой схемой. Пусть тг Ьа и [ = — *" — аЛи„, — (1 — о)Ливг л( е! Как и выше, назовем величину г = Ьг, [и)ь — ьг~', (12) и(:с,[+т) — и(т, [) да| тг д'гг! + з[ т д [[г,гт,(г) 24 дгз [[.ха 0 Лгг! = Ли( х-Л вЂ”,! +О(т ) = [ю1т-',з)) [,м-,) 2 дг [,,1ч-',) д гг! Лад'и~ т ди! дхг[[,ы-,") 12 дх~!р му д[ [~тир Отсюда т" =ба~„[и] — [с ~ =)~х,[+-)— (е,) ' 2 — )г,"„-~- т(о — 0,5)Л вЂ” [ + О(т + Л ).

ди! г г д[ [ю1тт) Таким образом, если [та = у(х, 1+ т/2), то г = О(бг+ т) при о т[ 0,5 и г' = О(бг + тг) при а =. 0,5. где [и)„— значения решения в узлах сетки Щ,, пазреилюетью огтроксимации разностной схемы (11) уравнения (1). Заьгнгим, что граничные и начальные условия для функции иь выполняются точно. Используя рюложение и в ряд Тейлора в точке (х, 1+у)'2) = (т)г, от+ т)'2), имеем 535 З 5 Рнзвгст„ого ~немы Ллн пнГгабохи*юского уГгеннения Будем исследовать устойчивость по начальным данным, т. е. будем оценивать чувствительнось репюння к возмущению начальных денных. Положим рь(1) ив е О, угг' = =О. Обознагилг ,глг г Л гуз ()о"((ь, „=- ~~, 1г(о,",,) ~ ,н=г Назонам разностную схому устой теса по напольным даппим е норме бз г„ сели супгествует постоянная сг, вс занисящая от шагов гтткгг Л и г, ти- кая, .гш для решения и" задачи (11), (5) с рг = рг = угг, = О сорансдлива оценка пшх )(он((с, „.

сг((гг")(ь,, о««ь (15) В теории разностных схем установилась традиция, когда не делают различия между лштрицей и порождаемым ею ливгйныы оператором. Обозначим через Л оператор (матриоу), конгрый функции е со значениими ое = О, нг,..., ом,, ом = О е узлах О,, М ставит н сгюгвежгтвио функцию со значениями в тех жс узлах равнымн О, Лог,..., Лом г, О. В снучно рг = рг = грь = О уравнение (11) можно записать н виде онь' = Егг" м-г о т 7гй(гггб) — ~ гйн1гг г=г Из равенств . з хй1г 4 нгпз— 2 "ь = яйтЬ хЬпб Лып — = -игжп —, Х Х следует, что чисва — ыь являлося сггркчненными значениями оператора Л.

Из (11) получаем (Š— огЛ)ои = (Е+ т(1 — о)Л)оа, т.о. Я = (Е— отЛ) '(Е + т(1 — о)Л). Поэтолгу (т хй(тЛ) 1 ч 1 — г(1 — о)оь тй(гпб) о =Е у сьн!в сь нш ь=г Х / ~ 1+гонг, ь=г Х Таким образом, собственные числа Ль матрицы Я имеют вид 1 — г(1 — о) ь Ль = 1+ель Матрица Я называется матроцен или оператором перехода от слоя к глого. В ггбщем шгучае Я может зависеть гп и. Пусть (Л„), г 1,..., М вЂ” 1, собственные числа мнгрицы Я.

Матрица Е симметрична и поэтому ()Яз = пзвх(Л,(. Функция ие может быть представлена в виде дискретной суммы Фурье Глава 10. ЛГепшы решени» уравнений в частных производных 536 Выясним, когда будет выполнено условие )Лг) < 1. Заменяя Ль езо выражением через 2ь, получаем условие — 1<- 1 — т [1 — и) иг. < 1. 1 й топь Так как пг., т > О, то эти нгравет2гзна эквивалентны аштн2лпениям — 1 — тпиь < 1+ топ! — тмь < 1+ тпит.

Правое неравенство выполнено всегда. а левоо — при т(1 — 2о)пь < 2. При гг Р 1/2 посчеднее неравенство будет выполняться при зпабоы т > О, а при и < 1/2 — если Х262 2 62 т< < <— 2(1 — 2гт] (1 — 2п) шзх ыа 2(1 — 2п) ! (14) 'Таким сбразом, нами получены достаточные условия устойчивости гхс мы (11), (5) по пачалыпаг данным. А именно, если ьа ш О и рм !вт ш О, то при и > 1/2 рази<к:тная схема (11), (Б) безусловно устойчива; при и < 1/2 схема угтойчива, если шаги й и т связаны гогггношепием (14), т.о, схема (11), (5) в этт2ьг случае ус~к~ало устойчива.

предположим, что условие )ль) < ! парт шветс», т. е. по )лз2 1) > 1 Р 6, гдг х(Л! — !)2пй д > О и не зависит от й и т. Положим 2з~'„ = ип „ . !авда Х п(3! — Ппй в,"„= (1-~-8)вин, и ))и")) = (! -~-6)"))пв)), гдв 2 > 4 > О. Л В этом <шучве швх ))и")) -т гю при т з О, т.е. 2хема неустойчива. кт Вногдв используют несколько отличнее опрелевеапе устойчивости по пвчаш пыла давным. Говорит, пю рззнпстнав схема угтюйчивв по начзльгшм данныи, если собствениыв числа оператора перехода лежат в круге рздпугз 1 ь ст.

Покажем, что и рас2матрнввсмом примере это определение согпвсуетсв с (13). Действительно, пусть Лв — гобсмшшые 2пслв лппрвцы о. Тогда и =.Попе = лт-~ Е" '21гп л,",ст вш — и ))гР)) < (! + ст)" ((пе)). В шюм случае при пт. = 2юпв! и т -з О с=2 имогм ()ш'() < е'~О2те)), пт < т. При исследовании разностньш схем длк более сложных задач, например при других типах граничных условий, более общем операторе в правой части уравнения (1) и т.пе доказательство устойчивости с использованием принципа максимума илн метода Фурье вызывает определенные затруднения, а иногда исследование устойчивости этими приемами является гц2осто невозможным. В этом случае исследование устойчивости разностных схем обычно проводится методом энергетических оценок.

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

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

Список файлов книги

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