Главная » Просмотр файлов » Самарский А.А. Гулин А.В. - Численные методы

Самарский А.А. Гулин А.В. - Численные методы (1078412), страница 47

Файл №1078412 Самарский А.А. Гулин А.В. - Численные методы (Самарский А.А. Гулин А.В. - Численные методы) 47 страницаСамарский А.А. Гулин А.В. - Численные методы (1078412) страница 472018-01-10СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Могут оказаться равными нулю и другие коэффициенты Ь„, но не все, так как по условию Ь,+Ь,+...+Ь =1 (см. (3) из 33). Пусть Ь,=Ь,=... ...=-Ь,,=О, Ь;ФО, 0<)~т, Тогда пз (20) получим а„д~ + а,ч'" '+ ... + е,„ ";4 + -т",а Отсюда видно, что прн больших о функция 1л(о) ведет себя как "о ь / Следовательно, для любого достаточно большого по модулю числа р (в том числе н для 1л, лежащих в левой полуплоскости) найдется корень о уравнения (13) с ~д~ )1, Доказано (см. [37] и указанную там литературу), что среди неявных линейных лгногошаговых л1етодов нет А-устойчивых методов, имеющих порядок точности выше второго.

Примером А-устойчивости метода второго порядка точности является симметричная схема (19). 254 имеет четвертый порядок точности и А(а)-устойчива при некотором а>0. 5. Чисто неявные разностные методы. В настоящее время при интегрировании жестких систем уравнений широко используется метод Гира (37), в основу которого положены чисто неявные многошаговые ревностные методы высокого порядка точности. Разностный метод л ~~; а„ул, = т1'(1„, у,) (22) называется чисто неявная Он является частным случаем не|ода (2), когда Ь,=Ь,=...=Ь =О, Ь,=1.

Для отысканвя у„получаем из (22) нелинейное уравнение »1 аоул — т~(1л, у„) = — '~~', алу, в, А=1 (23) которое можно решать тем или иным итерационным методом. Условия р-го порядка аппроксимации (9), (10) из 3 3 в случае метода (22) принимают вид л3 ~л Ш а = — ~~~ ~ам '~ ~йах= — 1, '~~ и'ил=О, 1=2, 3, ..., р. (24) Отсюда видно, что наинысший достижимый порядок аппроксимации чисто неявного гп-шагового метода равен т. Упомянутый выше метод Гира использует чисто неявные схемы наивысшего порядка аппроксимации. Система уравнений (24) для определения ж.

В связи с этим было введено еще несколько определений устойчивости, которые являются менее ограничительными, чем определе. ние А-устойчивости. Разностный метод называется А(а)-устойчивым, если область его устойчивости содержит угол (агд ( — р) (<а, р=тХ. гл 1 В частности, А~ — ~-устойчивость совпадает с А-устойчивостью.

(2 Доказано, что ни для какого а не существует явного А (а) -устойчивого линейного многошагового метода. Построены А(а)-устойчивые неявные методы третьего н четвертого порядка точности. К ннм относятся, в частности, чисто неявные многошаговые ревностные схемы, у которых правая часть 1(1, и) вычисляется только при 1= =1„, а производная и'(1) аппроксимируется в точке 1„по нескольким предыдущим точкам. Например, схема ул ул 1+ ул л ул-л+ ул-л ~(1 ) (21) 1хт козффициентов а„..., а,„метода наивысшего порядка имеет вид а,+2а,+...+та = — 1, а, + 2'а,+... + т-'а „= а,„ (25) а, + 2"а,+...+ т"а„=О. Эта система однозначно разрешима, так как ее определитель отличен от нуля.

При т=1 метод (23), (25) совпадает с неявным методом Эйлера. При т=2 и т=З получаем методы з ! — у, — 2у„, + — у. = т) (1,„ул), 2 2 (26) !! з ! — у.— Зу.-.+ — у.— — — у- =т! (!. у.) (27) где 9=т).. Ему соответствует характеристическое уравнение — — р)ра — 2а+ — =О.

3 ! (29) 2 Наы нужно найти множество точек 6 комплексной плоскости ц = ц, +!йь для которых оба корня дь, (!т) уравнения (29) нс пре- восходят по модулю единицу, Границей области 6 является множе- ство таких точек и, для которых ~д~ =1, Выразим нз уравнения (29) параметр р, через переменное д, т. е. запишем = — — 2д '+ — 3 2 (30) Отсюда видно, что если ~д~ =1, т. е. а=а-", то и = — — 2есч + — ен!.

3 ! ! 2 2 (31) При изменении ар!.умента гГ от 0 до 2п точка 9 описывает замкнутую кривую Г, симметричную относительно действительной оси имеющие, соответственно, второй и третий порядок точности, При и! =4 из (23), (25) получим схему (21).

Для практических расчетов используются аналогичные методы вплоть до десятого порядка точности. Важна отметить, что чисто неявные разностные методы обладают хорошими свойствами устойчивости, позволяющими использовать их для решения жестких систем уравнений.

Рассмотрим более подробно метод второго порядка (26) и найдем область его устойчивости. Для модельного уравнения (4) метод (26) принимает вид 3 ! — у, — 2у,, + — у„, = иу„, 2 2 (см. рис. 7). Для точек )А(д), расположенных снаружи от этой кривой, выполнено условие !с)~ <1, поэтому область устойчивости 6 метода (28) представляет собой внешность кривой Г. Точки, расположенные внутри Г, составляют область неустойчивости.

Обозначая х=соэсс,можно переписать (81) в виде !А= (1 — х)*-~-!У! — х'(2--х), откуда следует, что вся кривая Г расположена в правой полу- плоскости. Поэтому область устойчивости метода (26) целиком содержит левую полуплоскость и тем самым метод (26) является А-устойчивым. Исследуем аналогичным образом область устойчивости метода четвертого порядка (21).

/та Записывая характеристическое уравнение в виде р = — (28 — 484 '+Збд ' — 164 адни ') Рис. 7. Граница устойчивости метода (26) и полагая с)=е ', находим уравнение границы, разделяющей области устойчивости и неустойчивости: Р = — — (1 — х)'(Зх+ 1) +- 1 (бх' — 16х'+ 1бх — 8), г'! — к' 3 3 где х=созср (см. рис. 8). В отличие от предыдущего примера, имеются точки границы, расположенные в левой полуплоскости. Поэтому метод (21) не является А-устойчивым. Найдем теперь значение а, м,,,' С':,эск.,' .;;;, при котором метод (21) А (а)- Р~ ',, ',;.;,' Я'~~~:,.

устойчив. Для этого достаточно .'Ф : "".' найти угол се, который образует касательная 7 (см. рис. 8), профст; ,,;,". ходящая через точку (О, 0), с й' ' . Фо отрицательным направлением ОСИ )А,. Обозначим ра(х) = — — (1 — х)'(Зх+ 1), (32) Рис. 8. Граница устойчивости мото. да (2!) р,(х)= к (бха — 16х'+15х — 8). 3 Условие касания у с границей определяется следующим уравнением относительно параметра х: Р,(х) И, (к! Ра(х) Ра(х! 257 9 А. А. Самарский, А.

В. Гуана Из (32) получим ро(х) =8х(1 — х)', — 8к~ + ! 6хз — 4к~ — 8к + 5 р, (х)— г ! — х' (34) Уравнение (33) после подстановки в него выражений для р,, р„рм р, из (32), (34) и приведения подобных членов сводится к линейному уравнению х — 0,2=0. Прн х=0,2 из (32) получим 4' )Г24 699 р 3 ° 5~ 3 ° 5а так что 1ясг= — '= ))'6 Иг 4Г6 699 ра 5Г2 Поэтому метод (21) А(и)-устойчив, где а=агс16 16=68', В заключение отметим, что подробное изложение численных методов решения жестких систем дифференциальных уравнений содержится в книгах 126, 371.

ЧАСТЬ П! РАЗНОСТНЫЕ МЕТОДЫ РЕШЕНИЯ ЗАДАЧ МАТЕМАТИЧЕСКОЙ ФИЗИКИ Здесь излагаются приближенные методы решения краевых задач для уравнений с частными производными. В основе этих методов лежит сведение дифференциальной задачи к системе линейных алгебраических уравнений путем замены дифференциального оператора разностным. гллвл ВНОДНЫЕ ПОНЯТИЯ 9 1. Примеры разностных аппроксимаций Различные способы приближенной замены одномерных диффе- ренциальных уравнений разностными изучались в 5 4 ч.

1. Напом- ним примеры разностных аппроксимаций и введем необходимые обозначения. Будем рассматривать равномерную сетку с шагом А, т, е. множество точек со<,=(х< —— й, <=О, +-1, <-2,...). Пусть и(х) — достаточно гладкая функция, заданная на отрезке и,— и, и,.— и,, (х, „х,+,). Обозначим и<=и(х<), и,,= ", и-. = з и. = м< й< 2Ь Разностные отношения и„„и-., и; называются соответственно .<,< правой, левой и центральной разностнымн производными функции и(х) в точке х,. Каждое из этих разностных отношений аппроксими- рует и'(х) в точке хь т. е. при фиксированном х< и при А-<-0 (тем са- мым при 1-<-«и) пределом этих о~ношений является и'(х,), Проводя разложение по формуле Тейлора, получим и„х — и'(х,) = 0,5ди" (х,) + 0(/Р), и-,, — и' (х<) = — 0,5йи" (х<) + 0 (<<'), и — и' (х<) = 0 (1<').

Отсюда видно, что левая и правая разностные производные ап- проксимируют и'(х) с первым порядком по А, а центральная раз- 299 ностная производная — со вторым порядком. Нетрудно показать, что вторая разностная производная иих — и- с ис 2и,+ ис икх,с л л* аппрокснмирует и" (х,) со вторым порядком по й, причем справед- ливо разложение и-„,— и" (х;) = — иссс(хс) + 0(й'). Рассмотрим дифференциальное выражение 1.и = — с Й(х) — 1 сСх ~ ах 1 с переменным коэффициентом Й(х), Заменим выражение (1) раз- ностным отношением ис — ис л (2) где а=а(х) — функция, определенная на сетке ас,.

Найдем условия, которым должна удовлетворять функция а(х) для того, чтобы от- ношение (аи-„), с аппроксимировало (Йи')' в точке х, со вторым по- рядком по й. Подставляя в (2) разложения и,,с =и; + — и;+ — ис + 0(й ), л " л' 2 б л " лс и-„,. = и,'. — — и; + — й + 0 (/са), где и,'=сс'(х,), получим Л 2 б С другой стороны, 1.и= (Йи')'=Йсс"+Й'и', т, е. С ас„— ас,;, с ас„+ас сс 2 б Отсюда видно, что 1„и — 1лс=О(й'), если выполнены условия — '-=Й'(х;)+О(/сс), '"' ' =1с;+0(йа). (3) л Условия (3) называются достаточньсхси условияхси второго порядка аппроксимации.

Прн ссх выводе предполагалось, что функция и(х) имеес непрерывную четвертую производную и Й(х) — диффе- 260 ренцируемая функция, Нетрудно показать, что условиям (3) удо- влетворяют, например, следующие функции: ил = 0,5 (й(х;) + й (х;,)), а; =й(х; — 0,5й), а;=''г'/г(х;) й(хт,), Введем на плоскости (х„х,) прямоугольную сетку с шагом й, по направлению х, и с шагом й, по направлению х„т.

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

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

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

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