Главная » Просмотр файлов » Дьяченко В.Ф. Основные понятия вычислительной математики

Дьяченко В.Ф. Основные понятия вычислительной математики (1185904), страница 16

Файл №1185904 Дьяченко В.Ф. Основные понятия вычислительной математики (Дьяченко В.Ф. Основные понятия вычислительной математики.djvu) 16 страницаДьяченко В.Ф. Основные понятия вычислительной математики (1185904) страница 162020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

(6, Исследовать также рааностную схему вида ем о-Ь е . о е о "з м "з +, "а+ь "а-ье +ь "е, ы 2т 2йх 2ае 2. Рассмотреть итерационные способы решения системы урав- пений (219), (220), Оценить количество итераций. 3. Оценить число арифметических операций, необходимых дая решения аадачн на конечном отрезке времени н в ограниченной об- ласти, при использовании: а) явной схемы (214), б) неявнойсхемм (217), решаемой методом матричной прогонки, в) неявной схемы (217), применяя итерационные способы, г) метода перемеаных напраелеяий (226), (22?). В последних трех случаях привять т а, 4. Исключить из формул (226), (227) промежуточную величину и.

Сравнить нолучениое ревностное уравнение с неявной схемой (2(7). б. Показать, что для решения аадачи (213) могут быть использованы разностйые алгоритмы, определяемые формулами "а+1, 2"а, т+ "1-1, ю + х и и и + "а, тш 2"а. т + "а, ю-1 ! и юа — и и+1 и, и+1, «+1 а 1в+1 а Фл + а ж-1 и+1 иа,1и а,ив 1 — 2ва ю+ ва и "а, "а, в "ам, т 2"а, + "1-1, т т аа„ и+1 и+1 «..ив1 и+1 .ва, ю ва, ю "а, ю+1 а, ю+ Иа, ив-1 Провести сравнение их с неявной схемой (217) и методом переменных направлеввй (226), (227) (исключвв и).

6. Рассмотреть воаможности обобщения всех упоминаемых в атом параграфе методов на соответствующие трехмерные задачи, 6 (2. СТАЦИОНАРНЫЕ ЗАДАЧИ даУ даУ да + (229) и принимаюшую на границе Г этой области заданные значения (7(г = я. (230) Этот термин применяется для задач, описывающих стационарные, не меняющиеся во времени, состояния различных систем. Типичным представителем такого рода задач является следующая. Требуется найти функцию У (х, р), удовлетворяющую, в некоторой ограниченной области 6 плоскости х, у, уравнению Построение соответствующей раэностнои задачи не и область 6 асчетной сеткой, для простоты с равными шагами по л и у (рис. ).

(229) аменяем на этой сетке раэностным соотношением у у — В,ии которое имеет смысл для каждой внутреннеи расчетной й, т, для которой все четыре соседние точки й +. 1, ау~ Рис. 20. используемые в (231), расположены внутри области 6. Остальные расчетные точки я, т, при д ина лежащие 6, б р ничными их совокупность обозначим через ,о ъявимг а Ф пе еносом значения у. Значения и„, на у получим просто переносо у иэ ближайшей точки границы Г и 1„=а1г. (232) Разностная задача (231), (232) принадлежит к рассмоткак (231) в пределе при й — ~ 0 переходит в уравнение ( (232) в (230) (расстояние между у и Г порядка Ь). а в Дон~жом ущойчшюсть пестр Р т оенной азностной задачи, т.

е. совпадение порядко ядков решения и правых частей (232) при лю ом . Для б я. Д этого воспользуемся следующим приемом. 107 Пусть решение задачи и существует. Рассмотрим две вспомогательные функции о, и с, положив ох = ~и + а (хз + у') + р, (233) где а, р — пока произвольные постоянные. Обозначим левую часть (231) через Пи. Подставим (233) в (231). Получим Юсх = ~7+ 4а, Очевидно, для этого достаточно положить 1 а = — шах1(а,,„~. а,т (235) Определим теперь р так, чтобы пь ~„< О.

(236) В соответствии с (232), (233) это будет при 'р = — шах ! у / — а шах (х' + у'). (237) г с Допустим, что шах ох доствгается в некоторой внутренней точкей, ж. Тогда хотя бы в одной из соседних с ней точек значение ах меньше максимального, и, как нетрудно увидеть из (231), (Псх)„, „( О. Зто противоречит (234) и, следовательно, шах ох может достигаться только на границе у. Но здесь, в силу (236), сх отрицательна. Значит, она отрицательна всюду. Таким образом, ~ ин ~ + а (х' + у') и ~ + р ( О~ т. е. в соответствии с (235), (237), 1 шах/иа „~с шах!у~+ 4 шах~7и /шах(х'-(-у~). (238) ьм г цт ' с Полученное неравенство означает устойчивость задачи (231), (232). Заодно мы доказали существование и единственность решения этой задачи.

Действительно, поскольку всякое решение должно удовлетворять (238), то при д = 7" = О 108 так как Юи = 7, П (хз + у') = 4, Пр = О. Выберем а таким, чтобы всюду в области 6 выполнялось неравенство Вс+ >О. (234) возможно только тривиальное решение и = О. Следовательно, решение неоднородной системы линейных уравнений (23»), (232) существует и единственно.

Перейдем к вопросу о способах решения системы (231), (232). Исторически первыми были способы итерационные. Простейший из нкх следующий. Разрешим каждое уравнение (23») относительно значения и» в центральной точке расчетной ячейки; 1 в»," 4 (и»-ь +к»,ь„+и»,„,+и„,„„— Ь»Л» ) (239) и используем эту формулу для проведения итераций. Вычислительный процесс весьма прост — на каждой т-й итерации вычисляем среднее арифметическое значений и»акт»» в точках, окружающих данную центральну»о, го и получаем следующее приближение и» [м~п Исследуем сходимость этого процесса. Положим оэ все ик„= и,,„+Ь,, где и», — точное решение системы (23»), (232). Тогда, очевидно, ошибка Ькм будет определяться следующим итерационным процессом: Ь).

в = — (Ь)"+'ь + Ь»"-'ь + Ь»",' + Ь»",' — ), Ь'""'~, = О. (24О) Обоаначим максимум модуля Ь~»",~ через Ь и будем рассуждать так. Посколькуб»'~ есть среднее арифметическое четырех значений ЬДь вм то~Ь»~',1 ~ также не превосходит Ь, это справедливо для всех точек, на всех итерациях. Но для точек, соседних с граничными, можно сделать более точную оценку. А именно, если хотя бы одна из соседних для данной точки й, л» вЂ” граничная, где Ь = О, то в этой точке й, т ~бк ! 4 4 Ь' и, о+зй з Последнее неравенство справедливо для всего приграничного слоя точек, Перейдем ко второй итерации — влияние границы распространится еще на один слой точек, 109 В которых з оз 4 + 15- Зта оценка заведомо справедлива и для первого приграничного слоя точек.

Продолжая рассуждение, будем продвигаться с каждой итерацией в глубь области 6, получая для пройденных слоев точек оценку Наконец, на какой-то и-й итерации, и ИЬ, мы исчерпаем все расчетные точки. Это значит, что за и итераций ошибка 6 уменьшится, по крайней мере, в (1 — 4 ") раз. За следующие л итераций — еще во столько же раз, и .. Ю т. д. Мы доказали, что при т ->- оо ошибка 6(">-~0, т. е. итерационный процесс сходится. За я итераций ошиб»'->з 4я ка убывает в (1 — 4 ") раз, следовательно за Ряс.

21. одну итерацию, в сред- нем, в(1 — 4 ")и" или(так как и 1/Ь) в (1 — 4 ма)" 1 — Ь4»" раз. Принято скоросл>ь слодилос>ви характеризовать величиной относительного убывания ошибки за одну итерацию, т. е. отношением (6<"> — бои>)/бо> = к. В данном случае для этой величины мы получили оценку к Ь4»а (241) Фактически для многих случаев эта оценка оказывается слишком завышенной, итерационный процесс сходится быстрее.

Так, если область б является прямоугольником (рис. 21), то нетрудно получить более точную оценку. Для этого заметим, что формула (240), описывающая эволюцию ошибки 6<"> от итерации к итерации, может быть интерпретирована как формула разностной задачи бе+>> >б(» ыо уже знакомого нам класса (3 6). Все отличие состоит лишь в том, что т есть теперь номер итерации, а не номер временнбго слоя. Поэтому для исследования эволюции ошибки Ж"> мы можем применить спектральный приэнак. Не будем пока принимать в расчет граничных условий б(„= 0 и положим б(ч> фм> цкоиоо> к, «= осе (242) Как всегда, бом> окааывается равным Лб<">, причем в дан- ном случае легко получаем, что сов ~р + соо $ '2 >243) Поскольку экспоненты, фигурирующие в (242), выражаются черве эш >с„, соэ Ьр, зш ш>(>, соэ т>р, то интересующие нас функции являются комбинациями последних.

Чтобы удовлетворить левому граничному условию б„= О, эти функции должны содержать множитель э)й Ьр. Требуя выполнения условия на правой границе при Й = К, приходим к равенству эш Коэ = О. Это возможно только при Ку = ря, где р — целое. Таким обраэом, мы должны рассматривать только дискретный набор р, и даже конеч- ный ~рр — — р — "., р=1,2,..., К вЂ” 1, (244) так как при остальных р мы получим те же сеточные функ- цииэш Ьр„, а при р = О или р = К вЂ” нулевую функцию на сетке. т.

е. ~ Л ) ( 1. Такая оценка обеспечивала нам устойчивость эволюционных задач, но сейчас она недостаточна. Нас устроит только строгое неравенство (Л (~ 1, только это гарантирует сходимость: бе> -о- 0 при т -~- оо. Легко видеть, что (Л ~ = 1 соответствует функциям (242), которые не удовлетворяют граничным условиям б(„= 0 (они получаются при >р = о(> = 0 или >р = о)> = я), и потому оценка ! Л ( может быть уточнена. Оставим только те комбинации функций (242), которые обращаются в нуль на границах нашей прямоугольной области 6 (рис.

21), т. е. удовлетворяют условиям бо бл 0 ш 0 1 М бк, о = б», и = О. й = О, 1,..., К. По аналогичным соображениям, наши функции Ьз,,„ должны содержать множитель зш тф„причем ф~=д —, 7=4,2,...,М вЂ” $. (245) Итак, сеточные функции Ь», = з(пйф„з(птф ' (246) при любых ф„, фю определяемых равенствами (244), (245), удовлетворяют граничным условиям. Подставим (246) в (240) вместо Ьз("' . Получим после несложных выкладок („соз ф + соз ф Ь„"+ = ", ' з(п йф„з(п тф, т. е.

Ьм (246) являются собственными функциями итерационного оператора, а соответствующие им собственные значения выражаются формулой соз ф„+ соз ф Р Ю (247) совпадающей с (243). Запас собственных функций (246) болыпой (можно показать, что достаточно большой), и по величине Хр з (247) можно судить о действительной скорости сходимости итерационного процесса. Наибольшие значения ! Хр, з( достигаются при крайних значениях фр, фю т. е. при ! соз фр! = соз (я/К) и ! соз фз! = соз (я/М). Так как КЬ и Мй определяют размеры области С, т.

е. порядка единицы, то мы можем написать шах~Ар с<в (я/К) + соз (я/М) у 2 — соей — $ — —. 2 Следовательно, введенная выше характеристика скорости сходимости итераций х Ь', (248) что, конечно, лучше, чем (24$), полученное грубой оценкой. Мы уже отметили аналогию между итерационным процессом и эволюционной разностной задачей. Фактически П2 она простирается гораздо дальше. Всякую стационарную задачу можно рассматривать как частный случай эволюционной, нестационарной, где нас интересует лишь конечное, установившееся состояние, а не сам процесс установления.

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

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

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