3 часть (1081356), страница 44
Текст из файла (страница 44)
16.126*. Определить порядок аппроксимации разностной схемы задачи 16.117. В задачах 16.127 — 1б.129 для заданных задач Коши и соответствующих им разностных схем выполнить следующие задания: 1) Выписать разностный оператор Хь и правую часть 7"ь, возникающие при записи втой схемы в виде (4).
2) Нарисовать шаблон для разностного оператора Ьь. 3) Показать, что разностная схема аппроксимируег соответствующую краевую задачу с первым относительно Ь порядком на решении й(х, 1), имеющем ограниченные вторые производные. 4) Исследовать разностную схему на устойчивость в зависимости от значений параметра Л: т = Лб. Гл. 16. Уравнения в частных производных 318 16.127.и', — и' = 1'(х, 1), — со < х < +со, 0 < 1 < Т, и(х, 0) = = р(х); 1 1 (ит, и-1-1 ит, и) я (игп, и ит-1, и) т т, п~ т т=О,х1,...,хМ, и о=!р, и=0,1,...,(Т/т) — 1.
16.128. и~~ + и'„. = у (х, 1), — со < х < +со, 0 < 1 < Т, и(х, 0) = = Ю(х); 1 1 (ит,пЧ-1 ит,п) + ~ (ит,п ит-!,п) — !т,1и т й т=О, х1, ..., хМ, и о=!р, и=0,1,...,(Т(т) — 1. 16.129. и', + и' = у(х, 1), †< х < +со, 0 < 1 < Т,и(т,, 0) = = !р(х); ит,п.!.! ит,п ит+!,п игпоп) 1т,п т= О, х1, ..., хМ, и, о = рт, и = О, 1, ..., (Тут) — 1. 16.130. Исследовать устойчивость разностной схемы из примера 1. 16.131.
Исследовать устойчивость 11азностной схемы из задачи 16.112. 16.132. Исследовать устойчивость разностной схемы из задачи 16.116. 16.133. Исследовать устойчивость разностной схемы 1 1 2т (ит,п1-! ит,п — 1) З(ит — 1,п 2иттп + ит1-1,п)1 т=0,1,...,М вЂ” 1, и=0,1,...,Л вЂ” 1, ипЬО ='Рт ит,1 = Фт, ис,п = им,п = О. 2. Численное решение краевых задач методом сеток. Все разностные схемы, применяемые при решении краевых задач математической физики, делятся на два больших класса — явных и неявных схем. Под слоем разностной схемы понимается совокупность точек сетки З, лежащих на некоторой горизонтальной (или вертикальной) прямой.
з 3. Приближенные методы решения уравнений 319 Ясли значения сеточной функции и пеы заданные на (и+ 1)-м слое, выражаются в явном виде через значения этой же функции на слоях с меньшими номерами, то такая схема называется леной. В противном случае схема называется неявной. Например, разностиая схема, полученная в примере 1, является явной, так как может быть записана следующим образом (см. выраже- ние (7)): атт ип1,п~1 = тит,п + — „. (итэцп — 2ип|,п+ ит-цп) + утин и = О, 1, ..., я — 1; гп = 1, ..., т — 1, ипь о = рт, тп = О, 1, ..., т, иоп — — Я(1п), та=1,2,...,в, и„п —— фт(1п), гп т 1, 2, ..., ж Решение получающейся системы линейных алгебраических уравнений трудностей не представляет и осуществляется последовательно, переходом от слоя к слою.
Сложнее обстоит дело с неявными схемами (как в задаче 16.112,например). Для решения соответствующих им систем уравнений удобно применять гзетод прогонки, т.е. модифицированный метод исключения для решения системы линейных уравнений. Суть этого метода разберем на основе неявной разностной схемы для уравнения теплопроводности (см. задачу 16.112). Введем обозначения атт 2атт ат,п+1 = бт,п+1 = 1+ —, дпп,+1~ = ип, 1+ т)т,п+Г. (20) 62 ьт Разностную схему из задачи 16.112 перепишем в виде ио,п = Ф~,п, ат, пт1 ит ц и+ г + бт, пипп пью + ат, пЧ т итз ц пьг — дт, и+ ы (21) иь птуухп, п=О,...,я — 1, ит,о = уэт, Система уравнений (21) при каждом фиксированном значении и совпа- дает с системой апи,, + Ьпип+ с„и.ьг —— дп, п = 1, ..., Ж вЂ” 1, (22) интф, для которой справедлива следующая 320 Гл.
16. Уравнения в частных производных Лемма. Система уравнений (22), коэуьдьициентьь которой удовлетворяют неравенстоам )Ь„) > 1+ )а„!+ (с„), разрешима при любьж правых часпьлх, и длл ее решенил справедлива оценка )и( < шах()эв(, (дь(, ..., (д„ь(, ф). Для системы (21) условия леммы выполнены. Систему уравнения (21) удобно решать с помощью метода, который называется методоль прогонки. Будем искать решение системы уравнений (21) при каждом фиксированном и в виде иьи-! и!.! = ььььл п!.1ит,лв1 + Ни ив! ь т = ьь ... ь 2. (23) Исключая и,„! иэ! из системы (21), получим а,п иь! ит, п.гь— ипь-Ы, М-1 + аль, ив1ьгт, пь1 + бт, и-ь-! + дьл, и-!.1 аьл, и-ь-1Нт, и-1-1 аини~ЫЯ,„~Ы + б„п„.гь ' СООТНОШЕНИЕ (24) СВЯЗЫВаЕт ЗиаЧЕНИЯ фУНКЦИй и,и иэь, и ~1 иЫ, ПО- этому можно записать и,„,и.„! = Сд +ь,пь.ьи .~ь,пч.ь + Н.„, „,, (25) Сравнивая соотношения (24) и (25), имеем ат,пеь Ят-11, и-ь-1 = ььт, иэ1 ьгль, и ! 1 + бт, и ! ! (26) д„, иь! — О п.ььН т-ь-ь, и-ь-! ат, л-гьЬиЬт, иьь + бпь, и-1-1 Соотношения (26) определяют значения всех прогоночных коэффициентов Я,„, Н,„,и.
С помошью этих соотношений сетка проходится вверх по т от значения 1 до значения Й вЂ” 1 при фиксированном значении и. ПРи этом опРсделЯютсн все значеннЯ Сдт „, Нилл на сспье (пРЯмаЯ прогонка). Определив все прогоночные коэффициенты Я „, Ни, „, проходят сетку вниз от значения й до 2, последовательно определяя значения ип,,и из уравнения (23) (обратнал прогонка). Граничное условие прн т = 0 определяет начальные значения Яь „.ы, Н! иеь, а граничное условие при т = Й в обшем случае определяет первое значение ия.
Метод прогонки обладает тем свойством, что ошибки округления, получаемые на каждом шаге, не нарастают. Это свойство служит основанием ее широкого применения. 16.134**. Определить начальные значения прогоночных коэффициентов и значение иь для решения параболического уравнения З 3. Приояиженные методы решения уравнений 321 в случае граничных условий ь ди 1 / ди с21 — + р1и~ = Фн (аг — + Аи = Фг. В задачах 16.135-16.147, используя разобранные выше разностные схемы, найти приближенное численное решение: дги дги 16.135.
2 —, — 2 — — 1хг = О, Р = ((х, 1)(0 < х < 2, 0 < 1 < дхг <Ц,и~, о — — х, — =зшх,и~ =е — 1,и~ =4соз1. с=о дги дги 16.136. — — — — 8ху = О, Р = ((х, у) )О < х < 1, 0 < у < Ц, дхг дуг ди и(0, х) = х, и(0, у) = у, — (х, О) = О, и(1, у) = 1 — у. у 16.137. и... — 4и, = ху, Р = ((х, у))0 < х < 2, 0 < у < Ц, .2 и) о= х'+1, и,'.( -о — +6, и~*=о — з1ву+1 и~*=2 — бсоьу. О, Р = ((х у)(0 < х < 1 0 < 16.138.
и „вЂ” 4脄— ху = <у<Ц,и~ =хг, и=о ду /ди = 8х — 2, ( — +уи~ = уг, (,дх з ди — +уги~ =5уз+4у 16.139. и — 4и,„= е*", Р = ди Ц,и~ =4х+5, — =1, и=о ' ду и=о ((х, у))0 < х < 2, 0 < у < с ди — + зшуи = 5у, дх ' /,о — + узи = сову.
д2,и — 2 —,г — х = О, Р = ((х, у))0 < х < 2, 0 < у < Ц, д г и~ = сову, и~ = з1пу. ~т=е ' ~х.=г 16.140.— ди ду ! -о Задачи 16.140-16.143 решить, используя явную или неявную схемы: 322 Гл. 16. Уравнения в частных производных ди ди 16.141. — — 2 — + 1 = О, Р = ((х, у))0 ( х ( 1, 0 ( у ( Ц, х у /ди и1х, 0) =хч,и~ =еЯ", ~ — +и =О.
/ =о 16142. 4ил — 2и +хй = О, Р = 11х, С))0 < х < 2, 0 < 1 < Ц, ~4.=, ди д2и 16.143. — = — + 8х214, Р = ((х, 1))0 < х < 1, 0 < 1 < Ц, ' д1 дхз и~, о — — з1пнх, ~ — — 0,5и) = О, ~ — +0,5и = О. Решить следующие задачи Дирихле: 16.144. 2и + и„, = х~ + 2ху, Р = 11х, у))0 ( х < 2, 0 ( у < < Ц, ~ =1+х, Ч =3 — у', ~ =1+5у ( -,=6 дои д2и 16.145. — + —. = — х~ — у, Р = Цх, у) ~0 ( х ( 1, 0 ( У < дх дуз <Ц,и~ =2у~,и~,=у +4у,и~ =О,и~,=х2+2т,+2.
ди ди 16.146. — +4 — з — — вшху, Р = Их, У)~0 < х < У+1, О < У < дхз ду~ < Ц, и! = х, и(,= ху+ 1, и(,= 5 — х, и( „= 5У. 16.147. 2и +4лл„„+ух = О, Р = 1(х, у)/О < х < 2, 0 < у < Ц, и!.=о= 5У и~ =о= лх и! =л=5 х' и!х=2= ~/ лс 2 Глава 17 МЕТОДЫ ОПТИМИЗАЦИИ В 1. Численные методы минимизации функций одной переменной 1. Основные панлтил. Прямые методы минимизации. Пусть на множестве (/ с К определена функпил /(х). Под минимизацией функции у(х) на множестве У будем понимать решения следук/щей задачи; найти хотя бы одну точку минимума х" и минимум У'* = /'(х*) этой функции на множестве (/.
Задача нахождения точки максимума и максимального значения функции )(х) сводится к задаче минимизации залюной /(х) на — ) (х), поэтому ниже будут рассматриваться талыш задачи на минимизацшо. Напомним, что число х' е (/ называгтсл точкой абсолютного (глобального),минимума или просто точкой .минимума функции /'(х) на множестве (/', если у (х*) < / (х) длл всех х е (/. Значение /" = шш ) (/с) (/ называется абсолютным (глобальным) минимумом или просто .минимума,м )(/г) на г/.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
