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

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

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

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

е. множество точек тех — ((хт х ) ! хт — (йт ха — /йа' ( ) 0 ,' '2,...) (см. рнс. 9), н обозначим и;,; — 2ип -)- и; А,х,.н йт 1 и;;, — 2ип + и,.;, «,х„н Ьа Рис. 9. Сетка ыл и патпточеч- ный шаблон Из предыдущих рассуждений следует, что разностное выраже- ние Ь,ин =и-„х, + и„-, н (5) аппроксимирует дифференциальное выражение (4) со вторым порядком, т. е. А„ио — Еи (х,', ха) = О (й,') + О (/~,*). Более того, для функций и(х„х,), обладающих непрерывными шестыми производными, справедливо разложение йт даи (х', х~) йа д'и (х'„ х~а) дх~ ! 2 дхаа Разностное выражение (5) называется пятпточечным разностныи оператором Лапласа, так как оно содержит значения функции и(х„ х,) в пяти точках сетки, а именно в точках (х'„ х.,), (х,' ',х,'), тхт (х,', х, ) (см, рис.

9). Указанное множество точек называется шаблоном разностного оператора. Возможны разностные аппроксимации оператора Лапласа и на шаблонах, содержащих большее число точек. 261 Заметим, что если положить а,=й(х,), то получим только первый порядок аппроксимации. В 2 2 будет рассмотрен регулярный метод получения разностных аппроксимаций вида (2). В качестве следующего примера рассмотрим разностную аппроксимацию оператора Лапласа да и д'и 1и = —,,-+ —, (4) дхт дх„' 2 2. Построение разностных схем интегро-интерполяционным методом 1.

Построение разностной схемы. Задачи математической физики формулируются в виде основного дифференциального уравнения и дополнительных (начальных, граничных) условий, обеспечивающих существование и единственность решения. Типичными примерами являются: задача Дирихле для уравнения Лапласа, задача Коши и смешанные задачи для уравнений параболического и гиперболического типов. Под разносгной схемой понимают совокупность разностных уравнений, аппроксимирующих основное уравнение и дополнительные условия исходной дифференциальной задачи.

Существуют различные способы построения разностных схем. В этом параграфе будет рассмотрен один из способов, носящий название интееро-интерполяиионного метода (или метода баланса) построения разностных схем. В качестве примера рассмотрим применение интегро-интерполяционного метода к построению разностной схемы следующей краевой задачи для обыкновенного дифференциального уравнения второго порядка: (й(х)и'(х))' — д(х)и(х)+1(х) =О, 0<х<1, (1) (2) — И(0) и'(0) + би(0) = р„и(1) = рн где н(х), д(х), 1(х) — заданные достаточно гладкие функции, удовлетворяющие условиям й(х) '= А.)0, д(х) = О, и (1= О, р,„р,,— заданные числа При сформулированных предположениях существует и единственно решение и(х) задачи (1), (2).

Будем считать, что это решение является достаточно гладким. Уравнение (1) можно трактовать как уравнение установившегося распределения температуры и(х) в стержне длины 1, на одном конце которого (х= 1) поддерживается заданная температура рь а на другом (х=О) происходит теплообмен с окружающей средой по закону Ньютона (см. [41)). При этом й(х) — коэффициент теплопроводности, — н(х) и'(х) — тепловой поток, коэффициенты а(х), 1(х) характеризуют плотность тепловых источников, Для построения разностной схемы введем прежде всего на отрезке [О, 1) равномерную сетку с шагом и, т.

е. множество точек в =(х,=й, 1=0, 1,..., )У, ЬИ=1). Обозначим х,. а=х,-4-0,5й, ш(х) =й(х)и'(х), ш,,а=в(х; „,) и проинтегрируем уравнение (1) на отрезке [х; „„х,„„,), Тогда получим уравнение .~и'/з Й+~, пчхи — и, и — ~ г)(х)и(х)йх+ ~ 1(х)йх=О, (3) к, .ч у которое представляет собой уравнение баланса тепла на отрезке 262 [х, и„х<+„Д. Далее, заменим интеграл кои д(х) и(х) ах К! у Х,+у его приближенным значением и; ~ !7(х)дх и введем обозначения сч кои йии ! с а!= — ! а(х)ах, ср;= — ! 7(х)г(х.

! (4) х;-н сч В результате вместо уравнения (3) получим уравнение ь (5) Выразим теперь ю, „, через значения функции и(х) в точках сетки. Для этого проинтегрируем соотношение и'(х) =!е(х)/й(х) на отрезке [х, ь х,). Тогда получим к! Х! Г в(х) Г дх и; — и;!,= ) Йх из ,) а « '* ,) ь (.! ' сп 1 Обозначая (5) и — и. 4-1 получим ю, и, ~ а! =а;и-„и и;„.„,=а,+,икь Подставляя эти выражения в уравнение (5), получим разностное уравнение, содержащее значения искомой функции в точках х„ х„,: ! — (а;„и„х — а;и- ) — !(;и!+<р;=О хх или в сокращенной записи (аи„-), ! — Аи, + ~р! = О. (7) Это уравнение по своему построению является разностным аналогом основного дифференциального уравнения (1).

Записывая уравнение (7) во всех точках сетки, в которых оно определено, т. е. при(=1, 2,..., й! — 1, получим систему из !У вЂ” 1 линейных алгебраических уравнений относительно У+ 1 неизвестных и„и„..., и„. Два недостающих уравнения получаются путем аппроксимации граничных условий (2). Одним из этих уравнений является условие и = заз =!ко, а второе может быть получено интегро-пнтерноляцнонным методом. Для этого проинтегрируем основное уравнение (!) на отрез ке (О, хч,), где хи=0,56: км км ши — юо — ') а (х) и (х) к(х + ) ! (х) дх = О. (8) о о Полагая, как и ранее, в, „,=а~и;о получим при 1=1, что шоо†,и; „, Выражение для ш, следует из граничного условия при х= =0: ои,= — !к,+рио. Наконец, полагая к,, ку ) д(х)и(х)бх=и ~ д(х)дх, получим из тождества (8) разностное уравнение ки к~ а,и-„, — !)ио + р, — и, ~ д (х) бх+ ~ 1(х) г(х = О.

(9) Обозначая ку к/ б = ~ д(х)дх, сро= — ( г'(х)йх, (10) 2э4 перепишем уравнение (9) в виде — а,и„,+ (и+0,5йг(„) и,= р„+0,бар, Из этой записи видно, что полученное уравнение является разностным аналогом граничного условия — й(0) и'(0) + ри(0) = иь В дальнейшем решение разностной задачи в отличие от решения дифференциальной задачи будем обозначать буквой у, так что у,=у(х,), х,~ы,, Объединяя все уравнения (7), (9), получаем следующую разностную схему для задачи (1), (2): (ау-,)кх — Ау + чч= О, 1= 1, 2, ..., л! — 1, — а~у.,о+ (Р+О 5ЙЫо)уо=!к~+О 56~?о, Ук= Ро. При анализе разностной схемы (10)„как впрочем и любой другой разностной схемы, возникают следующие вопросы: а) существование и единственность решения системы линейных алгебраических уравнений (10); б) каким методом надо отыскивать это решение; в) какое отношение имеет система разностных уравнений (10) к исходной задаче (1), (2), иначе говоря, переходит ли разностное уравнение (10) в уравнение (1), если шаг сетки й стремится к нулю? Это вопрос об аппроксимации дифференциальной задачи (1), (2) разностной схемой (10); г) сходится ли решение у(х) разностной задачи к решению и(х) дифференциальной задачи при Й- О? На первые два вопроса можно ответить немедленно.

Разностная задача (10) является типичным примером задачи, которая решает- ся методом прогонки, изложенным в и. 7 3 4 ч. 1. Систему уравнений (10) можно записать в виде АУ,,— СУ,+ВУ,,= — Рн 1=1, 2,, )Ч вЂ” 1, Уо = н1У1 + Р ! Ум = нэух+, Рм где А;=аь В,=а; „С,=а,+а,+,+й иь Р,=й2ср„н«=0, ! й (Н1-(- О зрю) ! + йа, ~ (й + О 5йио) Из условий а,)0, ))~0, а,)0 следует, что С;= А,+В,>0, т. е.

выполнены условия устойчивости прогонки. Поэтому разностная задача (10) однозначно разрешима и ее можно решать методом прогонки по формулам (44) — (46) из 3 4 ч. 1. Вопросы аппроксимации и сходнмости обсуждаются в следующем параграфе. в 3. Исследование аппроксимации и сходимости 1. Аппроксимация дифференциального уравнения. В ~ 2 рассма- тривалась краевая задача (й(х)и'(х))' — у(х)и(х)+1(х)=0, 0<х<1, (1) — й(0)и'(0)+()и(0) =)хо и(1) =)хм (2) й(х) )с,>0, р~О, для которой ннтегро-интерполяцнонным методом была построена разностная схема (ау„-), ~ — Ау; + <р; = О, 1 = 1, 2, ..., Ф вЂ” 1, а1у«,а+ах=)хм ум =)хм (3) (4) где «~ — 1 а~= ! Г дх 1=1,2, ..., А(, й,~ й (х) «и« (6) 265 «ли «1>и — д(х)дх, гр;= — ( 1(х)с(х, 1=1,2...,, )ч' — 1, (6) й й «; «~-у )) =6+ 0,6М., р1=!1+ О,бйр., о,м О.Ы с(,= — ( а(х)((х, <рв= ') )(х)с(х.

0,5й,) 0,5й о о Обозначим через (.и(х) левую часть уравнения (1) и через х«у, — левую часть уравнения (3), т. е. Еи(х)=(й(х)и')' — у(х)и(х)+)(х), (,йу;=(ау„-)«,; — Ау;+!рь выполняется соотношение (ав„-)„и — (й (х) в' (х))' )„= 0 (й'). Если кроме того, докажем, что с(,=д(х,)+0(й'), <р,=~(х;)+0(й'), (9) то тем самым будет установлено, что оператор Е, аппроксимирует Е со вторым порядком по й, т.

е. Е„п,— Еп(х<) =0(й'), 1=1, 2,..., й! — 1. (1О) Итак, доказательство второго порядка аппроксимации сводится к проверке условий (8), (9) для коэффициентов (б), (б). Проверим сначала выполнение условий (8). Обозначая р(х)=й '(х), получим к! 1 1 С !Р— — 1 р(х)дх=р! и+ — р! „+ 0(Ч), а; й у ° тн — 1 следовательно, Р; у а,=й! и —— !2 Р, и Аналогично +0(й!)=й, ', +О(/!).

Ьа Р! 12 Р! а„, = йи и — — —, + 0 (й ). й' Р~ з 12 Р~ О~сюда получим ""*+" =й(х!)+О(й), 2 ' и + 0 (й') = й' (к,) + 0 (1 '), т. е. условия (В) выполнены, Условия (9) выполнены в силу того, 266 Пусть п(х) — достаточно гладкая функция н п(х,) — ее значение в точке х, сетки ам=(х,=!й, 1=О, 1,..., Л', ЕМ=1).

(7) Говорят, что разностный оператор Е, аппраксимирует дифференциальный оператор Е в точке х=х„если разность Е,п,— Ев(х;) стремится к нулю при й — «О. В этом случае говорят также, что разностное уравнение (3) аппроксимирует дифференциальное уравнение (1). Чтобы установить наличие аппроксимации, достаточно разложить по формуле Тейлора в точке х=х, значения в;„=в(х,-«й), входящие в разностное выражение Е,п«Большая часть этой работы проделана в 5 1, где показано, что при условиях й '=й'(х;)+0(й'), '" ' =й(х!)+0(й') (8) 2 что замена интегралов (б) значениями дь !', соответствует приближенному вычислению этих интегралов по формуле прямоугольников с узлам в середине отрезка интегрирования.

2. Аппроксимация граничного условия. Исследуем погрешность аппроксимации разностного граничного условия (4). Обозначим 1еи(0) = — а,п„,+рп,. Если п(х) — произвольная достаточно гладкая функция, то очевидно 1ьп(0) = — й(0) и'(0)+ра(0) + 0(й), т. е. имеет место аппроксимация первого порядка по й. Однако если п=и(х) — решение задачи (1), (2), то разностнае граничное условие (4) имеет второй порядок аппроксимации, т.

е. — а,и,, + ~3и, — ц, = — й (0) и' (О) + ри (0) — р, + 0 (62). Докажем последнее утверждение, Используя разложения и„д —— (и,— и,)16=и'(х„,) + 0(ЕР), х,и — — 0,5й, а,=й„,+0(й-'), получим а,и„, = юи + 0 (й""1 = ц, + 0,56ш, + 0 (6') = =й(0)и'(О)+ 0,56(йи')' (0)+ 0(й'). Отсюда имеем 1„и (О) = — й (0) и' (0) — 0,56 (йи') ' (0) + ри, — и, + 0 (6') = = 1 — й (0) и' (0) + (Зи (0) — р,)+ 0,56 ! — (йи') ' (О) + Н,и,— АД + 0 (6') Учитывая граничное условие (2), получаем 1,и(0) =0 56( — (йи') '(0) +д,и,— ср,!+0(й'). Выражение, стоящее в квадратных скобках, преобразуем, учитывая уравнение (1), к виду — (йи') ' (0) + сК„и,— аь = — (йи') ' (0) + д (0) и (О) — ) (0) + + (г1,— а(0) ) и,— Ц(0) — ч4) = (д,— д(0) ) и,— (1(0) — ф,) . Из соотношений к,, Н,= ( д(х)йх=д(0)+0(6), п,=)(0)+0(6) о получаем 1ьи (0) = — о,и,, + (Ти, — рт = 0 (6'), что и требовалось доказать.

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

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

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

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