Главная » Просмотр файлов » Тихонов, Свешников - Дифференциальные уравнения

Тихонов, Свешников - Дифференциальные уравнения (947323), страница 31

Файл №947323 Тихонов, Свешников - Дифференциальные уравнения (Тихонов, Свешников - Дифференциальные уравнения) 31 страницаТихонов, Свешников - Дифференциальные уравнения (947323) страница 312013-09-15СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Ошибкой решения схемы (6.46) по-прежнему назовем выражение б 'л'у = = и'у — 'л'у. Величина б'л'<р, как и <л»<р» оценивается в норме 16 <л<,з))< = и<эх(16 <л'<1 16 <л<у)<з). (6.48) Определение. Схема (6.16) называется устойчивой по начальным условиям и правой части уров н ения (или просто устойчиво й), если существует такая постоянная Ьс, что при Ь ( Ьс для нормы о<иибки решения выполняется неравенство )а <юу1 С С16 <л< !), (6.46) где постоянная С не зависит от Ь. Данное определение, очевидно, означает непрерывную зависимость решения разностной схемы (6.16) от входных данных задачи — начальных условий и правой части уравнеш<я: малое изменение входных данных приводит к малому изменению решения, Легко убедиться, что схема (6.31) при а) 1/2 устойчива по начальным данным.

Устойчивость по начальным данным можно проверить положив равной нулю б <"'( — компоненту в выражении для б'л'<р (так и было в нашем примере). Действительно, как было показано выше, бу( равномерно ограяичеио. Что касается (л< бу(л)„то, как видно из (6.41), при о ) 1/2 имеет место неравенство Решкние нАКАльнОЙ злдлчн 1)<л1«Ф и поэтому бу< тоже будет равномерно ограниченным и (Я) даже бесконечно малым с уменьшением )<. Поэтому из (6.44), (6.45) и (6.46) следует , «6 <л> (6.5О) где е=шахцео!, )е!!)Ф (6.5() что и доказывает высказанное утверждение, Перейдем теперь к доказательству устойчивости схемы Эйлера (6.22).

Соответствующая возмущенная задача имеет вид <лг жг У< З< — ) г ( (Л> — ~ л<Л)г (6.52) Ьл(л)у = А (л) ул = ул + Ел- Будем предполагать, что решение '"'у этой задачи принадлежит области П, в которой выполнены неравенства (6.4), Для ошибки решения бачу разностной задачи (6.53), очевидно, получим 6 (л)у б(л) б(л) з«- / <л) л < до (>о ) ~я< — и у<-)) + ) (х<-и )у<-г)= 6 Л-ж 6(л)у, =е,. Здесь ( 6.53) (6.54) Из (6.53) в силу оценок (6.4) получим 16 Ф>у<1 л 16 <л>у <! + ЬК16 <му< 11 + Ыб <л>)!1 = = (4+ йй) 16 <">у,,)+ й))6 <">7() ~...

((4 ) <>й)Р16<л>у ! ) (( ! яй;)о((6<л>)() ) я!)6<л>)!1 (655) Применяя этот процесс последовательных оценок, после <-го шага найдем 16(")У<1«(4+ ЬК)'1, !+ Ь 1(Ф+ йК)< '+ ... + (1~!6<л>У1« «((+ ЬК)<16(")у1, ++((+ йК)*'~6<">У1. (6.56) В силу (6.48) и определения чебышевской нормы отсюда получим 16(")4«((+)<К)"(~6(")у~,+ — '.

!!6(")У~). (6.5У) И* чнолкннык мктоды <гл. е Так как Х вЂ” хв=пЬ, то в силу неравенства (1+ Ьй)<~ "')<л ." ~~ е окончательно будем иметь к <к-з~) Пб<~)УИ~~ек<х )(Дб<~)ун + ~ Пб<~)УИ) < Снб< ) И (6.58) прн любом Ь, что и доказывает устойчивость разностной схемы (6.53). Введенные выше понятия порядка аппроксимапии и устойчивости разностной схемы позволяют доказать следующую теорему, играющую фундаментальную роль при исследовании сходимости раеностных схем. Теорема 6.1.Если разностная схема (6Л6) устойчива и аппраксимирует задачу (6Л4) с порядком Ь, то решения оэу при Ь - О сходятся к реи<ению у(х) дифференциальной задачи, причем имеет место оценка <ул'У вЂ” (У)л)! «СЬ", (6.59) зде постоянная С не зависит от Ь.

Доказательство. Обозначим разность значений сеточных фУнкций 'мУ и (У)л чеРез б 'л'У: б <л) <л) (6.60) В силу определения порядка аппроксимации Ыу)л = "~у+ 6 ~>ц, (6.60 причем «б <л)<р!)) «С)Ь'. Из устойчивости разностной схемы (6Л6) в силу (6.49) получим ))б <л<У)) «С<)б <л<<р<) «СС<йл (6.62) Сохранив для произведения постоянных СС) обозначение С, мы получим соотношение (6.59). Теорема доказана.

Выше было показано, что схема Эйлера устойчива и аппраксимирует задачу (6Л) с первым порядком. Иэ доказанной теоремы следует, что семейство решений "'у, полученных по схеме Эйлера при Ь- О, сходится к точному решению зада ш (6Л) также с первым порядком. Этот факт был установлен в п. 1 с помощью прямых оценок.

3. Метод Рунге — Кутта. Как мы уже отмечали, схема Эйлера имеет лишь первый порядок сходимости, и для получения высокой точности по этой схеме приходится вести вычисления с очень мелким шагом, что привод)гг к значительному росту ввемени счета. Поэтому естественно искать схемы, обладающие более высокими порядками сходимости. Одним иэ классов таких схем являются схемы метода Рунге — Кутта. В основе этого метода лежат следующие соображения. Рассмотрим тождество гишкпик н<дчальнои з<(дачи у'(х) =)(х, у(х)). Сопдаспо формуле Лагранжа о конечных приращениях ка отрезке (х<-д, х») существует такая точка х*, что у(х<) — у(х» д) = (х» — х» д)у'(х*) = Ц(х*, у(х<')).

(6.63) Однако значение х~ нам ле известно. В схеме Эйлера положено х* = х<-д, что и дает всего лишь первый порядок точности. Основная идея метода Рунге — Кутта заключается во введении в разностную схему ряда дополнительных параметров, уточняющих приближевпое определение значения хз. Будем по-прежпему.

рассматривать начальную задачу (6.1)д — ", =Их, у), у(*;) =у.. (6 1) Условия гладкости функции 1(х, у) будут сформулированы ниже. Сопоставим задаче (6.1) разпостиую схему (ю„(ю„ Х = — (66) (ю д "д где К,=Кх, »,'"'у, д), Кд д(х — д + а»Ь (д)у — » + с(»ЬК») (6.65) К» д(х< (+а, »Ь< 'му, »+а, »ЬК, (), а рд, ..., р», ад, ..., а<» — некоторые параметры, выбором которых можно обеспечить требуемый порядок аппроксимации схемы (6,64) иа решевии задачи (6.1). В частном случае при 1 1, р=1 схема (6.64) переходит в схему Эйлера (6.22), ддмеющую первый порядок аппроксимации.

Покажем, что при 1=2 моншо так выбрать параметры рд, рд и а», что схема (6.64) будет иметь второй порядок аппроксимации. При атом будем предполагать, что функция ((х, у) имеет в О непрерывные частные производные до второго порядка по обоим аргументам, что достаточно для существования непрерывной третьей производной от решения.

Тогда, подставляя решение задачи (6.1) в схему (6.64) при 1= 2 и разлагая решение в строку Тейлора, получим ьд р .У (х;) — У (х, д) = ЬУ' (хд д) + —, У" (хд,) + б У (х~) = =Ь(рД(хд д, у(хд д)) + + рд( (х( д+ а»Ь, у (хд-д) + а»Ь( (хд д, у(х;,)))). (6.66) численные методы >гл з Воспользовавшись разложением ) (х« ~+ а>Ь, у(х«>) +а«Ь~(х«и у(х«>))) = 7(х«м У(х«>))+ 1Ь а (')+а Ь7( ) ду (*)+О(Ь ) (6.67) и очевидным соотношением у" ( —.) =т(-)+ — ( ) у'(; —.) = — ( ° )+1( ) — (.), (666) перепишем (6.66) в виде у'(х; >) — (р, + рз) ~(х«п у(х«,))+ + Ьу (х«->) ( 2 — а,р.,) = 0 (Ьз).

(6.69) Отсюда следует, что при рг + рз ' а'р' г (6.70) схема (6.64) при с = 2 имеет второй порядок аппроксимации. Прн атом а> может быть выбрано произвольно. Наиболее широко применяются схемы (6.64) с сс« =1 и а>=1/2, Заметим, что путем выбора сс«повысить порядок аппроксимации схемы (6.64) прп 7 = 2 нельзя. При а>=1 из (6.70) получим р«=рз=1/2 и схема (6.64) принимает внд «л>„«ь >„ — — (7(х« „«юу«,) + +У(х >+Ь, «>У, >+Ь7(х«м «ь>У«г))) = О. (6 717 Геометрическая интерпретация этой формулы ясна из рис.

19. Сначала по методу Эйлера находим точку у«='"'у«,+Ь|(х, ь '"'у« >), затем находим среднее значение угла наклона касау у «х>> телыюй к интегральным кривым «ю ° « нашагеЬ:«йа= ~ ( У«->+У«) > у! " и по нему уточнлем значение «с~ '"'у,. Аналогичные рассмотрения « легко провести и для случая а = $/2. Подобные схемы обычно >У-/ х; носят название «предиктор-корректора. Мы рассмотрели схемы со вторым порядком акпроксимации. Аналогичные рассмотрения могут быть проведены и для схем более высокого порядка П > 2).

При этом приходится нала- РЕШЕНИЕ НАЧАЛЬНОЙ ЗАДАЧИ <67 гать иа функцию >(х, у) все более высокие требования гладкости. Наиболее широкие практические применения находят схемы четвертого порядка. В большинстве библиотек стандартных программ для ЭВМ используется следующая схема: <л>„ <л>„ 7. (>> ) с (К~+2К~+2К~+К ) (л) «ю О, (6 72) «о где Кг = ) (х<< м (й)у(-<), К. =-)(х-.

+ —,, < >у< —, + —,К.).. й й й яг= 7((х<-1+ 2 1 у<-<+ 2 ч) Кй = )(х< — >+)<, (")у< — + йКЕ). (6.76) <л) (й) "<-> г (й) — ~(е< м «< г) (Ы е (6. 74) (й)у где функция <7(х, у) представляет собой линейную комбинацию функций )(х, у) от промежуточных значений аргумента. Поэтому оценки для функции С(х, у) н ее производных лед<о выразить через соответствующие оценки функции 7(х, у) и ее производных, используемые при доказательстве устойчивости метода Эйлера, откуда н следует сделанное утверждение. В частности, имеет ь<есто следующая Теорема 6.2.

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

Тип файла
DJVU-файл
Размер
2,41 Mb
Тип материала
Учебное заведение
Неизвестно

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

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