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

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

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

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

1), так как условия устойчивости прогонки выполнены. 276 Для исследования устойчивости разностной схемы (12) будем искать частные решения уравнения лн л л~!л лгг 1'-~-1 г ' г-г О г г аг г имеющие вид (10). Тогда получим д = !' ! + 4у гйп — ), у = —, ,!ге ~,-г а' следовательно, 1д)(! при любых гр, т, й. Таким образом, схема (12) абсолютно устойчива, т. е. устойчива прн любых шагах т и Ь. Абсолютная устойчивость является основным преимуи(еством неявных схем. Теперь уже не надо брать шаг т слишком малым, можно взять, например, т=й=10 '. Величина шагов сетки т, Ь определяется теперь необходимой точностью расчета, а не соображениями устойчивости.

Шеститочечной симметричной схемой называется разностная схема лгг л Уг Уг 1 л.г л л = — (У;„; + У„-„;) + %, (14) (! 5) 1=1, 2,..., )ч' — 1, и=О, 1,..., К вЂ” 1, у"„" =)гг(Г„„„), ун" =рг(!л.г), п=О, 1,, К вЂ” 1, У,'=ил(х;), Е=О, 1, ..., М. При а=0 получим отсюда явную схему, при а=! — чисто неявную схему и при а=0,5 — симметричную схему (14). Исследуем погрешность аппроксимации схемы (15) на решении исходной за.

дачи (1) — (3). Представим решение задачи (15) в виде уг"= =и(х,.1„)+ а",, где и(х„1„) — точное решение дифференциальной задачи (1) — (З). Тогда для погрешности получим систему для которой начальные и граничные условия задаются так же, как и в схеме (12). Эта схема использует шеститочечный шаблон, изображенный на рис. 11, в. Предлагаем читателю самостоятельно доказать, что эта схема имеет второй порядок аппроксимации как но Л, так и по т (если только гр",.

=1(хь 1.+0,5т)+0(т'+и')), она абсолютно устойчива и ее можно решать методом прогонки. Обобщением трех рассмотренных схем является однопараметрнческое семейство схем с весами. Зададим произвольный действительный параметр а и определим разностную схему л+г л = ау-„', + (! — а) у"-„, + Чг„ уравнений (!б) и, перегруппнровывая слагаемые, получим, что ~>; = (и" — и + ср;) + (и — 0,5) ти" + — и'~ + 0 (т' + Ь'). 12 Учитывая уравнение (1) и" — и= — !' и следствие из него и'~ — и"= = — !", окончательно можем записать, что ф";= ~(а — О 5) т+ — 1и" +~р' — )(хь 7„.л)— !2! — — !'"(хь г„,и) + 0(та+ и'), (1о) Из формулы (18) можно сделать следующие выводы.

Если и' а а = а, = — — — ср". = 1 (хь 7„, „) + — !" (хь 7„.к ) + 0 (тг + й), гч и — г" 1=1,2,...,М вЂ” 1, п=0,1,...,К вЂ” 1, г"" = гчы = О, и = О, 1,, К вЂ” 1, г'. = О, ! = О, 1, ..., !Ч. Сеточная функция ф,"., входящая в правую часть уравнения (16) и равная ф; = пи„-, 1 + (1 — о) и„-„, — и,л + ср;, (17) называется погрешностью аппроксимации схемы (15) на решении задачи (1) — (3). Получим первые члены разложения функции по степеням й и т. Будем разлагать все функции, входящие в выражение для ф";, по формуле Тейлора в точке (хв 1„+О,бт).

Учитывая разлогкепия и~ с =и(хь г„„п)+0(т'), и„-„,. =й(хс) + +,— и" (х,)+0(п'), где й'=д'и)дх', и=ди)д1, 1„+,и — — Г„+0,5т, получим л l ~ч ф = а (й (хь 1„,Д+ — и" (хь 7„п)) + (1 — о) ~и" (хь 1„) + 12 + — и'"(ХЬ 7„)) — и(Х„7„~И) + ~рб+ 0(тз) + 0(84).

Отсюда, проводя разложение в точке (х„г„.„п) и обозначая и= = и (х„7„.„, и), будем и меть Ь~ т /Р 2 12 2 12 — и+ р",+О( ')+ 0(й4) то схема (15) имеет второй порядок аппроксимации по т и четвер- тый — по й, Такая схема называется схемой повышенного порядки 278 аппроксимации. Если а=0,5, (р,".=1(хп «„еи) + 0(т~+ ««~), то схема (15) имеет второй порядок аппроксимации по т н по ««. Прн остальных значениях а и при ~р~ =1(х„«„„)+0(т+««') схема (15) имеет первый порядок аппроксимации по т н второй — по й. Опуская выкладки, отметим, что если искать решение уравнения (!5) с ~рр=— 0 в виде (!0), то получим иа ! — 471! — а) вп'— 2 у= ар 1+ 4та яз†2 н ) д~ -. 1 прн всех ~р, если л а '=- .— —— 2 4т (19) где т л л и к у = —, г'; = у~ + (! — а) ту;,, + т<р; .

Система (20) решается методом прогонки. Условия устойчивости прогонки (условия (47), (48) нз 9 4 ч. 1) при а~О сводятся к неравенству )1+2сг~) э2)а)7 и выполнены прн а) — 1/(47). Последнее неравенство следует из условия устойчивости (19) разностной схемы. 4. Уравнения с переменными коэффициентами н нелинейные уравнения.

Рассмотрим первую краевую задачу для уравнения теплопроводности с переменными коэффициентами р(х, «) — "= — 1««(х, «) — 1+~(х, «), 0(х(1, 0(«(Т,(21) д«дя1 ' дх/ и(х, 0) =и,(х), й(0, 1) =)«,(1), и(1, 1) =)«,(1), где р(х, «), й(х, «), )(х, «) — достаточно гладкие функции, удовлетворяющие условиям 0<с,(««(х, «) ~с„р(х, «) )с,>0. (22) чтя Отсюда видно, в частяости, что все схемы с а)0,5 абсолютно устойчивы. Схема повышенного порядка аппроксимации (а=а.) также абсолютно устойчива, что проверяется непосредственно. Прн ачьО разностная схема (15) является неявной схемой. Для нахождения решения у',"' по заданным у,".

требуется решать систему уравнений аууй~~~ — (1+ 2ау)у,""+ ауу,"",'= — г,", «= 1, 2, ..., М вЂ” 1, (20) ен км уе =Мт(«ан), уи =ух(«п~1)~ 'д ! ди! Дифференциальное выражение Еи = — (н(х, !) — ) прп каждом д~ фиксированном ! аппроксимируем в точке (хи !) так же, как н в стационарном случае (см. В !), разпостным отношением Л (!) у; = (а (хп !) у;)хи = = — ( а(хи„!) ' — а(хи !) ' ' ~, (23) где разностный коэффициент тенлопроводности а(хи !) должен удовлетворять условиям второго порядка аппроксимации а(хьп !)+а(хи !) =2й(хи !)+ 0(!!'), =й'(хи !) + 0(й').

и Наиболее употребительны следующие выражения для а(хь !): а(хь !) =0,5(й(хи !)+а(х; и !)), а(хп !) =й~х; — —,т!), й 2х (х~ „ !! х (хи !! а(хи !)= А(х; !)+А(хи О Здесь в качестве ! можно взять любое значение !ен(1„, („„), например 1=!.+О,бт.

Если в уравнении (24) !=!„+0,5т, о=0,5, то схема (24) имеет второй порядок аппроксимации по т и по й. Прн остальных значениях о и ! выполняется первый порядок аппроксимации по т н второй — по Й. При исследовании устойчивости разностных схем с переменными коэффициентами иногда применяется принцип замороженных коэффициентов, сводящий задачу к уравнению с постоянными коэффициентами. Рассмотрим явную схему, соответствующую уравнению (24! с п=О н )(хи !) =— О, т. е.

схему лх1 и р(хп !) ' ' =(а(хи !)у"-,),х. (25) Предположим, что коэффициенты р(хь !), а(х„!) — постоянные, р(хь!)=р=сопз(, а(хь!)=а=сопя!. Тогда уравнение (25) можно записать в ниде у" ,"' — у'," =ау- 280 Разностная схема с весами для задачи (2!) имеет вид ~» л р (хи !) =Л (!) (оу",.'1+ (1 — о) у",.) + ~ (хи !), 2= (, 2, ...,Л! — ), (24) У,",=Р1((х) У;"~=Рх(4и), У,'=их(х). или «+1 « = у-„,, та Р Из п. 2 известно, что последнее уравнение устойчиво при т'(О,бй', т.

е. прн (26) р з Принцип замороженных коэффициецтов утверждает, что схема (25) устойчива, если условие (26) выполнено при всех допустимых значениях а(х;,1), р(х„(), т. е. если при всех х, 1 выполнены неравенства та(ки 0 а« ( р(ки О к (27) Если известно, что О~с,(а(х„() (с„р(хь 1)=»с,)0, то неравен- ство (27~ будет выполнено при Строгое обоснование устойчивости схемы (25) будет дано в примере 2 нз ~ 4 гл. 2. Если параметр о)ОД то из принципа замороженных коэффициентов следует абсолютная устойчивость схемы (24).

Рассмотрим теперь первую краевую задачу для нелинейного уравнения теплопроводности — = — (й(и) — ) +1(и). (28) В случае нелинейных уравнений, когда заранее неизвестны пределы изменения функции й(и), избегают пользоваться явными схемами. Чисто неявная схема, линейная относительно у~+', 1=1, 2,... ..., Л' — 1, имеет вид «и «( ««««и «~1 « 1 ' ' = — „~и;„'" ' —; * ' ' ~+1Ьс), (29) где й; — й (у",.).

яа! где а,= 0,6(й(у";) +й(у,",)). Эта схема абсолютно устойчива, имеет первый порядок аппроксимации по т и второй- по й. Решение у,"+', 1=1,2,...,М вЂ” 1, находится методом прогонки. Заметим, что схему (29) можно записать в виде у"; = — ((йу„-) . + (йук)-,;) + ~(ут) 1 Часто используется нелинейная схема Л+1 Л ЛМ .Л11 — а (ул") '" ' — а (у л~ " ь а+1 а+1 ")" -'- '+~(") а(и~~-11 ! ), („а+11 а(ул") = 2 (зо) Для реализации этой схемы необходимо применить тот или иной итерационный метод. Например такой: (1+1) Л / (5+1) (111) (111) ($+1) ') М 1 У(а) Ул у(м) ул+1 ! 1 ' 1 из которой находятся промежуточные значения у("4, (=0,1,...

...,()). Затем на втором этапе используется симметричная шести- точечная схема для уравнения (28), в которой нелинейные коэффициенты а(у), 1(у) вычисляются при у= у,""', т. е. схема +1 ( ) (( ( л+и) лм) + ( ( л„и) л) ) +( ( л+!1) )=1,2, ..., Ф вЂ” 1, Уо" =Р1((л~ ). Уя' =))1((л+1).

282 Здесь э — номер итерации. Как видим, нелинейные коэффициенты берутся с предыдущей итерации, а в качестве начального приближения для у,"." выбирается у( Это начальное приближение тем лучше, чем меньше шаг т. Число итераций М задается из соображений точности. В задачах с гладкими коэффициентами прн я(и) ~с,)0 часто бывает достаточно провести две — трн итерации. Значения у('"" на новой итерации находятся из системы (31) методом прогонки. Прн М=1 итерационный метод (3!) совпадает с разностной схемой (29). Для приближенного решения нелинейного уравнения (28) применяются также схемы предиктор — корректор второго порядка точности, аналогичные методу Рунге — Кутта для обыкновенных дифференциальных уравнений (см. п.

2 $ 1 гл. 6 ч. П). Здесь переход со слоя п на слой и+1 осуществляется в два этапа. Приведем пример такой схемы. На первом этапе решается неявная линейная система уравнений +Н л ' =(а(у";) у-„"'и)ля +)(у,"), (=1, 2, ..., л( — 1, уа""М =)(1 (1,+ О ~бт) у~~"' = ра (1л+ О,бт), 1.

Разностные схемы для уравнения колебаний. Рассмотрим первую краевую задачу для уравнения колебаний — — 0». х». 1, 0(1(Т, дм дх» и(0, Г) =1х,(1), и(1, 1) =1х,(1), О~С~Т, и(х, 0)=и,(х), — (х, 0)=и,(х), 0 ='х(1. й 5. Трехслойные разностные схемы (1) (2) (3) Известно (см., например, (41)), что эта задача поставлена корректно, т. е. ее решение существует, опо единственно и непрерывно зависит от начальных н граничных данных.

Будем использовать ту же сетку ы„, что н в ф 4, т. е. а„= =.х „ ыь=(хс=й, 1=0, 1, ..., У, йМ= 1), ы,=(1„=пт, и =О, 1, ..., К, Кт =Т). Очевидно, минимальный шаблон, на котором можно аппроксимировать уравнение (1), это пятнточечный шаблон, изображенный па рнс ! 1, г. Таким образом, в отличие от схем для уравнения теплопроводности, в которых использовалось только два временных слоя (слон и и а+1), здесь требуется использовать три слоя: и — 1, п, и+1. Такие схемы называются трехслойными. Их применение предполагает, что при нахождении значений у',." на верхнем слое значения на предыдущих слоях у',-', уе 1=0,1,...,Ф хранятся в памяти ЭВМ. Видимо, простейшей разностной аппроксимацией уравнения (1) и граничных условий (2) является следующая система уравнений: тх »х 1=1,2,...,Л' — 1, п=!,2,...,К вЂ” 1, у,""=р,(1,), у""=р,(1»ы), п=О, 1,..., К вЂ” 1.

(5) Разностиое уравнение (4) имеет второй порядок погрешности аппроксимации по т н по и. Решение у,"." выражается явным образом через значения на предыдущих слоях: у»ы — 2у» у»-ь + у (у» 2у» + у» ) 1=1, 2,..., Ф вЂ” 1, (=та/й', п=1, 2,..., К вЂ” 1. (6) Для начала счета по формулам (6) должны быть заданы значения у,'., уе 1=0, 1„...,У вЂ” 1,У. Из первого начального условия (3) сразу получаем у,'=ио(хс), 1=1, 2, ..., Й вЂ” 1. (Т) 283 Простейшая замена второго из начальных условий (3) уравнением (ут — уа)/т=йО(х,) имеет лишь первый порядок аппроксимации по т.

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

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

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

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