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

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

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

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

Поэтому численно решать уравнение (15) проще, чем исходное уравнение (6). 3 а м е ч а н и е 2. Уравнение Вольтерра (7) формально сводится к уравнению Фредгольма (6), но ядро прн этом имеет особенность (обычно разрыв) на диагонали х=$. Поэтому для уравнения Вольтерра следует выбирать обобщенную формулу трапеций и проводить уточнение способом Рунге. Многомерные задачи допускают, в принципе, применение описанного 'метода; надо только в (10) и других формулах под х„подразумевать узлы многомерной кубатурной формулы х„. Однако получить удовлетворительную точность' при умеренном объеме расчетов удается лишь для достаточно гладких К(х, «) и Г(х), когда можно использовать кубатурные формулы высокого порядка точности (напрнмер, произведение одномерных формул Гаусса с небольшим числом узлов л по каждой переменной).

В более сложных случаях развивают специальные методы; многие из них используют симметрию задачи и слабую зависимость' решении от части переменных. 3. Метод последовательных приближений. Это простейший приближенный лдетод. Запишем для неоднородного уравнения Фредгольма (6) итерационный процесс: ис (х) =О, и,чд(х) =1'(х)+й')К(х, е) и„(е) Щ.

(16) й 5=Х (рис. 102). Затем в качестнс (10) выбирают обобщенную формулу трапеций (4.7), причем в интервалах, прилдыкающнх к особой липни, используют соответствующие односторонние пределы функций. Вели вне особых линий все функции непрерывны вместе с достаточным числом своих производных, то при сгущении специальной сетки ьюжно уточнять решение способом Рунге.

Полезна предварительно так преобразовать исходное уравнение, чтобы гладкость решения повысилась. Например, если ядро непрерывно, а )(х) разрывна, то и(х) тоже разрывна. Полагая г(х) =-и (х)— — 7'(х), получим вместо (6) уравнение когягктно постлвленпые злдхяи 5 и Нетрудно показать, что при ограниченном ядре и достаточно малом значении !Х ! этот процесс сходится к решению уравнения (6). Д о к а з а т е л ь с т в о.

Обозначим погрешность п-й итерации через г„(х) =и, (х) — и(х). Вычитая (6) из (16), получим га, (х) =-) ~ К (х, ~) га Я) Л$. а (17) Отсюда следует неравенство !!г„,!!с=-:!)',(Ь вЂ” и) !!К(х, Е)!!с!!г,(х)!!с. (16) Тем самым, если выполнено условие Ч=,'Х (Ь вЂ” и) !!К(х, Е)!!с(1, (19) г„,(х) =1~ К(х, с) га (а)~Ц. а (20) Выкладки, полностью аналогичные доказательству сходимости метода Пикара (гл. ЧП1, э 1, и, 3), приводят к оценке !!га (х)!!с-== ', (! л ! (Ь вЂ” и)!!К (х, Е)!!с)"!!га (х)!!с (21) При п- оо правая часть этого неравенства стремится к нулю прн любых значениях )., что доказывает паше утверждение. Замечание 2.

Оценку (18) можно переписать вследующем виде: !!г„(х) !!с ==- д" !!г, (х)!!с д" (22) Отсюда видно, что метод последовательных приближений для уравнения Фредгольыа эквивалентен разложению в ряд по степеням параметра Х. Это можно строго показать, выражая и„(х) через и,(х) =1(х) при помощи рекуррентного соотношения (16). то итерации (16) сходятся равномерно по х, причем сходимость линейная. При достаточно малом !).! условие (19) выполняется. В практических вычислениях квадратуры, возникающие в этом методе, редко удается выразить через элементарные функции.

Поэтому обычно ограничиваются нахождением первых приближений. Замечание 1. Для уравнения Вольтерра (7) метод последовательных приближений сходится равномерно по х при любых значениях Х. Действительно, в этом случае вместо (17) справедливо соотношение 460 ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ [ГЛ. Х!У Пример. Рассмотрим уравнение и(х) — Л~ е-!"!ми(е) ие=х. о (23) Применяя процесс (16), получим и„(х) =О, иь(х) =х, и,(х) =х+Лг-", и, (х) = х + Ле + — !Ре -, и, (х) = х + ~) + — А- + 4 Л~) е 1 ь 1 .„1 н т. д. В этом случае удается найти точное решение Ль Ла и (х) = х+(Л+ — + — + ...1е-" = х+ — е-".

(24) 2 4 '''1 2 — Л Нетрудно заметить, что последовательные приближения здесь сходятся только при ~Л) ( 2. 4. Замена ядра вырожденным. Ядро уравнения Фредгольма называется вырожденным, если оно является суммой конечного числа членов вида К(х, Ц= ~Ч~~ Ал(х) В„($) л= ! (25) (для уравнения Вольтерра ядро не может быть вырожденным, иначе оно тождественно равнялось бы нулю). Решение уравнения с вырожденным ядром находится за конечное число действий.

В самом деле, подставляя ядро (25) в неоднородное уравнение (6), представим решение в виде суммы конечного числа членов: и (х) =Г(х)+Л ~Ч" „а„А (х), л= ! ь сьл = ~ В„ ($) и ($) г(е. а (26а) (266) Г ь 1 ь ~, '~6„— Л) В„(З) А,л(~) ьча,л=~ В„(З)1($) ьф, 1 =.и =!тг. (27) ~л .=. ! а а Решая эту систему и подставляя найденные значения ал в (26а), найдем искомое решение. Для однородного уравнения Фредгольма (8) надо положить во всех формулах Г(х)=0. Тогда система (27) становится однородной и представляет собой задачу на нахождение собственных Подставляя (26а) в (26б), получил! линейную систему для нахож- дения коэффициентов сь„: 461 4 и кОРРектнО постлвленные злдлчи значений матрицы У-го порядка.

Отсюда видно, что вырожденное ядро (25) имеет ровно 1Е' собственных значений Х1. Произвольное ядро нередко удается хорошо аппроксимировать вырожденным ядром. Например, разложим К(х, $) в ряд Фурье по некоторой полной ортонормированной системе функций В„($); коэффициенты этого разложения будут функциями от х: ОЗ ь К(х, 5)= У„А„(х)В„($), А„(х)=~В„"Я)К(х, $)|ф. (28) л =- 1 В качестве (25) можно взять отрезок разложения (28). Тогда формулы (25) — (27) позволяют найти приближенное решение. Оценки точности таких приближений мы не рассматриваем, поскольку они громоздки и неудобны в практических вычислениях.

Замечание. Пусть в неоднородном уравнении (6) свырожденным ядром (25) правая часть 7(х) =„йО и такова, что выполняется ') ) (с) В„$) 11Ч = О, 1 ~ 11 ( 11(. (29) О Тогда при Х, равном одному из собственных значений ядра Ц, система (27) имеет нетривиальное решение, причем не единственное. Тем самым, в данном случае существует решение и(х) уравнения (6). Г1р имер. Рассмотренное в п. 3 уравнение (23) имеет вырожденное ядро К(х, $)=е-"е-». У него должно быть ровно одно собственное значение; Определим его. Полагая в (27) й| = 1 и Г(х) =О, легко получим со 1 — Х,( е '»|(а|с»1=0, о ! откуда Л,=2. Заметим, что точное решение (24) неоднородного уравнения (23) при Х=-Х1 не существует.

5. Метод Галеркина (который для интегральных уравнений Обычно называют методом л|ол|ентов). Будем искать решение в виде разложения по полной системе функций фл(х): и(х) =) (х)+Х ';1, 'с»,фл(х); (30) Л=| поскольку от и (х) не надо специально требовать удовлетворения каким-либо краевым условиям, то от системы фл(х) ничего, кроме полноты, требовать не надо. Подставляя разложение (30) в неоднородное уравнение Фредгольма (6) и требуя ортогональности невязки ко всем функциям 46а интеГРАльные уРАВнения |гл. х|ч «р» (х), 1 =- л == Г|Г, получим линейную алгебраическую систему уравнений для нахождения сь»! 'У'„а»|ь» —— -Г«„, 1~л«="й(, »=1 ь ьь Я»=) «р (х) «р»(х) г(х — ).~ ~К(х, $) «р (х) «р»($) г(ХЩ, (31) а аа ьь Ь =') )К (х, 2) «(! (Х)Г'д) »ГХ й.

а В случае задачи на собственные значения (8) надо полагать в (30) и (3!) Г'(х) =О, Г«(етод применим и к нелинейному уравнению (1), но тогда он приводит к нелинейной алгебраической системе. Основной трудностью, препятствующей применению метода моментов, является сложность вычисления двукратных интегралов, входящих в (31). Поэтому обычно ограничиваются малым числом членов суммы (30).

3 а м е ч а н и е. Если система «р» (х) ортогональна, то метод моментов эквивалентен приближенной замене ядра на специальное вырожденное ядро: К (х, $) = ~', «р» (х) Ч'» (3), (32) Ч'» а = $ К (х, 2) р» (х) Гх. 5 2. Некорректные задачи 1. Регуляризация. Если в интегральном уравнении (1) правая часть г" (х, и (х)) не зависит от решения, т. е, и (х) входит только под знак интеграла, то задача оказывается некорректно поставленной, Классическими примерами некорректных задач являются уравнение Фредгольма первого рода: ь ) К(х, $) и (~)г(а=~(х), с(х(|Г, а (33) и уравнение Вольтерра первого рода: ')К(х, В)и($)«1$=((х), с~х~ь(. (34) а В отличие от уравнений второго рода, ядро уравнения Фредгольма (33) задано на прямоугольнике (с~х(ь(, а"=1(«1 а уран НЕКОРРЕКТНЫЕ ЗАДАЧИ ненни Вольтерра (34) — на трапеции [с==х-.-ег, а($ (х]*), при- чем функции гь(2) и 1(х) определены на разных отрезках и при- надлежат разным классам функций У и г".

Покажем, что задача (33) неустойчива по правой части и, тем самым, некорректна. Для этого рассмотрим высокочастотное возмущение с конечной амплитудой би Д) = ехр (1гое), оэ ,=~ 1. Еьгу соответствует возмущение правой части ь ь 61 (х) =~К(х, 4) би(е) с$= ~К(х, Е) ег"ЕЕ$. а а Интегрируя по частям, получим ь 67(х) = —,. ег"З К(х, с)1 — —. ~ ' ег"е.гЦ=О( — ). (35) а Это означает, что для достаточно больших частот величина з 6(~1 = — 0 (1!го) оказывается сколь угодно малой. Следовательно, су- ществуют такие сколь угодно малые возмущения правой части 61 (х), которым соответствуют большие возмущения решения би Я), т.

е. задача (33) неустойчива. Для уравнения Вольтерра (34) справедливы те же рассужде- ния. Напомним, что в главе П1 мы уже сталкивались с некор- ректностью задачи численного дифференцирования функции 1(х); эта задача сводится к решению уравнения к ~ и (с) г( $ = г (х), (36) а т. е. является частным случаем уравнения Вольтерра первого рода, с ядром К(х, с) = 1 (при 3 ~х), Кроме того, задачи (33), (34) имеют решение не при любых непрерывных правых частях 1" (х). Так, задача (36) имеет реше- ние только для днфференцируемых 1(х).

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

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

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

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