Главная » Просмотр файлов » Н.Н. Калиткин - Численные методы

Н.Н. Калиткин - Численные методы (1133437), страница 68

Файл №1133437 Н.Н. Калиткин - Численные методы (Н.Н. Калиткин - Численные методы) 68 страницаН.Н. Калиткин - Численные методы (1133437) страница 682019-05-12СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

53: ! - ° Л л1 !3 2!2'1 2т (Уп+Уп-2) (Зс+ Ап/ (Ул-1+Упп2)+(41 Лп) Ул. (29) УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ [гл. 1х в этих узлах нетрудно написать разностные уравнения (17)1 и лог Уо = К ((т), УА1 = Рг (Ге) которые являются точными (их невязка равна нулю). Более сложен случай второй краевой задачи для того же уравнения (дальше будем рассматривать только левое условие): ти и-т и„(0, ()=р,((), и,(а,()=р,(().

(30) Можно аппроксимировать производную односторонней разностью: Л (У1 Уо) = Р1(гто1). (31) Однако невязка этого разностного уравнения равна 1 . Ь оро = (й„) о — — „(й, — й,) = — — и„= 0 (й), (32) т. е. имеет меньший порядок малости, чем невязка (25) в регулярных узлах. Это приводит к понижению общей точности расчета. Рассмотрим способы написания разностного краевого условия нормальной точности 0 (Ло). Сделаем это на примере явной схемы (18). С п о с о б ф и к т и в н ы х т о ч е к очень нагляден. Введем вне отрезка 0 =хк-а фиктивную точку х,=х,— й и будем считать исходное уравнение справедливым прн х 1~х.

Тогда разностное Возможны случаи, когда часть коэффициентов схемы типа (28) определяют из условия наивысшего порядка малости невязки, а часть коэффициентов выбирают из других соображений. Метод неопределенных коэффициентов (как и метод разностной аппроксимации) применйм к уравнениям с непрерывными и достаточное число раз дифференцируемыми коэффициентами и решениями.

Из-за сравнительной громоздкости он применяется реже двух ранее описанных методов. К р а е в ы е у с л о в и я. Остановимся на записи разностной схемы в нерегулярных узлах (на границе илн вблизи нее). В этих узлах для записи разностных уравнений необходимо привлекать краевые условия. Например, в разностных схемах (16) и (18) для уравнения теплопроводности и, = йи„„ нерегулярными являются граничные узлы и =О, и =уо'. Лля первой краевой И-1 и ы задачи и(0, Т)=111(~), и(а, ()=ро(() зот лппгоксимлция уравнение (18) можно написать при и =0: 1 /г —,, (Уо — Уо) =9 (У-т 3Уо+Уд. Заменим в левом краевом условии (30) производную симметричной разностью: 1 2а (Ут Исключая из последних двух уравнений фиктивную точку, получим разностный аналог краевого условия: 1 Ь Ь(У1 УО) Р1 (~ ) + 2 (Уо Уо)' (33) Заметим, что это уравнение содержит только одно значение с нового слоя у„ т. е.

оно явное. Метод уменьшения невязки менее нагляден, но более универсален. Выразим и(х„т) при помощи формулы Тейлора: и (х„1) = и (х„1)+Й и, (х„()+- Ь'и„„+... На основании краевого условия (30) положим и,(х„1)=р,(1), а из уравнения теплопроводности (15а) найдем и, -ийй. Под- ставляя эти величины в формулу Тейлора, получим и(х1, 1) = и(хо, 1)+Ьрг Е+,— и1+". А и(х) — )(х) =О, хеп б, (34а) Ри(х) — р(х)=0, хенГ.

(34б) Введем в области О+Г сетку с шагом и, состоящую из множества внутренних (регулярных) узлов ю„и множества граничных (нерегулярных) узлов у„. Заменим задачу (34) в регулярных узлах разностным аналогом уравнения.(34а): А ьуь (х) — ~рь (х) = О, х ев ы„, (35а) Заменяя здесь и~ — (у,— у,)!т, снова придем к краевому условию (33). В последнем способе можно учесть большее число членов ряда Тейлора и получить краевые условия не только нормальной, но и повышенной точности. 5. Аппроксимация и ее порядок. Пусть имеется область 6 переменных х = (х„..., х ) с границей Г и поставлена корректная задача для некоторого уравнения с граничными условиями: УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ )ГЛ. 1Х (37) (ь ) 112 Ци (х)Ц =(((р,(х) и„'(х)+р, (х) и'(х)) с)х~ а р„(х) ) О, р, (х) ~ О.

(38в) Употребляются и другие нормы. Напомним (см. главу !), что из локальной близости функций следует их среднеквадратичная близость; поэтому Ц Ц, называют более сильной, чем Ц Цс, Нетрудно проверить, что энергетическая норма (Ззв) сильней Ц Ц . Выбор той или иной нормы в конкретной задаче определяется двумя соображениями. Желательно, чтобы разностное решение у а в нерегулярных узлах — разностным аналогом краевых усло- вий (34б): )с»у» (х) )(» (х) =О х е— = 7» (35б) (индексом й отмечены величины, определенные только на сетке; мы будем опускать его там, где это не вызовет недоразумений). Близость разностной схемы (35) к исходной задаче (34) будем определять по величине невязки: )р» (х) = (Аи — )) — (А»и — гр„), х ен ю», У» (х) =(Рль — р) — Я»и — )(») х ен у».

Определение. Разностная схема (35) аппраксимирует за- дачу (34), если Ц 1Р Це»-ьО, ЦУ Цх„-ьО пРи )т-+О; (38) аппроксимация имеет р-й порядок, если )) РЦ»=О()гл), Цу Ц,„=О()ьл) при )г — О. Обсудим вопрос о выборе норм в этом определении. Функции и (х), ((х), р(х) определены обычно на отрезке а(х =(т или во всех точках некоторой области пространства большего числа измерений, Для них можно ввести такие нормы, как чебышевская (локальная): Ци(х) Цс= шах !и(х)), (38а) а(к(Ь или гильбертова (среднеквадратичная): гь ч1/2 ))и (х)))с, = р (х) и'(х) с(х~, р (х) ) О (38б) (выражения написаны для одномерного случая).

Часто исполь. зуют связанные с оператором А энергетические нормы, напоми- нающие формулы для полной энергии колебательной системы, например: % г1 АППРОКСИМАЦИЯ ) Аг 'цг !!у!1,= шах !у„!, !!у!1г, =( ~', р„у,"-/г„~ 0(а~% го= 1 (39) В выборе разностных аналогов норм существует некоторый произвол.

Например, сумму в (39) козино брать по л=б, 1, ..., У вЂ” 1, что соответствует выбору другой квадратурнсй формулы для интеграла (386). Этим пользуются, определяя сеточные нормы так, чгобы облегчить доказательство сходнмости. Как надо понимать Ь -ь Ор Лля равномерной сетки это не требует пояснений. На неравномерных сетках рассматривают совокупность шагов Ьа как некоторую сеточную функцию и вводят какую-либо норму шага, например: ум — ~ 'пг Ь=!!Ь„!1,=снах Ь„или Ь=!!Ьз!!ь =~ ~„"Ь,',/ а —. о Эту норму считают «величиной шага» в определениях аппроксимации, порядка аппроксимации н т.

д. Если невязку оценивают в 11 11„то аппроксимацию называют локалььог". Для уравнений с достаточно гладкими решениями наличие локальной аппроксимации и ее порядок легко проверяются; в таких задачах нередко ограничиваются установлением локальной аппроксимации. Однако наиболее сильные результаты по сходимости разностных схем связаны с использованием более слабых норм для невязки (но сильных норм для решения).

Замечание 1. Факт наличия или отсутствия аппроксимации и порядок аппроксимации зависят не только от операторов А и А„"), но также от классов, к которым принадлежат и(х), )(х), и от выбора норм. Чем сильнее норма или чем шире классы функций, тем ниже, вообще говоря, порядок аппроксимации *) Операторы краевых условий для краткости обычно будем опускать. было близко к точному решению в возможно более сильной норме; например, в задачах на разрушение конструкций малость деформаций в 11 11с, не гарантирует сохранения конструкции, а малость в 11 1!с — гарантирует. С другой стороны, чем слабее 11 11„, тем легче построить сходящуюся в этой норме разностную схему и исследовать ее.

Заметим, что функции и(х), )(х), р(х) принадлежат, вообще говоря, разным классам. Например, если и(х) есть четырежды дифференцируемая функция и А = (г(г)с(хг), то ) (х) является дважды дифференцируемой функцией. Поэтому каждую из этих функций можно оценивать в своей норме: 11 11в, 11.!1), 11 11„.

ФуНКцИИ уь, Чга, уа ОПрЕдЕЛЕНЫ ТОЛЬКО На СЕТКЕ, ПОЭтОМу дпя них надо ввести сеточные нормы 11 11 „, 11.1!ть, 11 1!х„. Их вводят так, чтобы при й — 0 они переходили в выбранные 11 11„, 11 1!), 11 11„. За разностные аналоги чебышевской и гильбертовон норм можно принять соответственно зш (ГЛ. (Х УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ (последнее видно по замечанию 2 к п.

3, если оценивать невязку в (( ((с). Замечание 2. Как правило, решение и(х) исходной задачи (34) неизвестно, так что использовать его для получения невязки затруднительно. В этом случае берут достаточно широкий класс У функций и(х), которому и (х) заведомо принадлежит (обычно это класс функций, непрерывных вместе с достаточным числом своих производных). Если на всех функциях класса У имеется аппроксимация порядка р: ((Ап(х) — Аьп(х)+(ра(х) — )(х)ЦА=О(йл) при й- О, то, очевидно, аппроксимация на решении и (х) имеет порядок ие ниже р. В подобных случаях аппроксимация на решении и (х) может иметь порядок выше р. В замечании 3 к п.

3 мы видели, что для уравнения из=йи„„явная разностная схема (18) при т = = из((бя) имеет в классе сколь угодно гладких функций п(х) аппроксимацию 0(6я), а на решений — более высокого порядка (четвертого, как нетрудно проверить). Сл у ч а и м н о г и х и е р е м е н ц ы х имеет некоторые особенности.

Определение аппроксимации остается в основном прежним; надо только требовать стремления к нулю шагов по всем переменным. Порядок аппроксимации может быть разный по разным переменным. Например, для двух переменных соотношение ((зр()е =0(тл+йт) при т-ьО, Ь-ьО (40) означает р-й порядок по времени и д-й по пространству. Это хорошо видно на примере схемы (18) с иевязкой (25). Аппроксимация вида (40), погрешность которой стремится к нулю при любом законе стремления шагов к нулю, называется безусловной или абсолютной.

Если же погрешность аппроксимации стремится к нулю при одних законах убывания шагов и не стремится при других, то аппроксимацию называют условной. Например, если ((фЦь= 0,тя+йт+ — „—,) при т-ьО, Ь-ьО, (41) то аппроксимация условная: кроме т -ь О, Ь-ь О, надо дополнительно требовать, чтобы (т'/Ь') -~- О. Если аппроксимация условная, то разностный оператор АА при разных законах изменения т(Ь) может аппроксимировать разные дифференциальные операторы. Например, можно проверить, что 1 1 ! "'(ьу= (Уа — Ул)+ (Узы Ул-т) (Ул.ьт — аул+ Уя-т) (4х) зп э з1 нстоичивость при т=сл аппраксимирует сператор д д А= — + —, дт дх' а при т=сла — оператор д д 1 У А = — + — — — —. дт+дх с дхз' Поэтому, если нет специальных соображений, лучше пользоваться разностными схемами с беаусловной аппроксимацией.

5 3. Устойчивость 1. Неустойчивость. Для некоторых разностных уравнений малые ошибки, допущенные на каком-либо этапе вычисления решения, при дальнейших выкладках сильно возрастают и делают невозможным получение сколько-нибудь пригодного результата. При численном дифференцировании и суммировании рядов Фурье мы встречались с некорректными задачами, где бесконечно малая ошибка входных данных может привести к большой ошибке решения. Теперь рассмотрим пример неустойчивой разностной схемы решения задачи Коши для уравнения и'(х) =аи(х). Выберем следующую схему: а 1 — а Т (Ух:т Уа) + Ь (Уз — Ул-з) = с'Уа. (43) При а=1 она переходит в схему ломаных (8.15), устойчивость которой была доказана. 1зассмотрим случай а~ 1.

Для простоты отбросим ошибку аппроксимации и исследуем только рост ошибки начальных данных. Тогда ошибка бу„будет удовлетворять тому же уравнению (43), ибо оно линейное однородное. Удобно исследовать рост ошибки специального вида бу„=г". Подставляя ее в (43), получим аг'+ (1 — 2а — аЬ) г — (1 — а) = О. Если Ь ~1, то корни этого квадратного уравнения равны г, = ==- 1+ О (Ь), ге =- 1 — (1!а) + О (Ь) . Тогда при с ( з/з будет ~ г, ~ ) 1, т. е. ошибка такого вида. возрастает за шаг в несколько раз.

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

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

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

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