Главная » Просмотр файлов » Д.П. Костомаров, А.П. Фаворский - Вводные лекции по численным методам (pdf)

Д.П. Костомаров, А.П. Фаворский - Вводные лекции по численным методам (pdf) (1113729), страница 8

Файл №1113729 Д.П. Костомаров, А.П. Фаворский - Вводные лекции по численным методам (pdf) (Д.П. Костомаров, А.П. Фаворский - Вводные лекции по численным методам (pdf)) 8 страницаД.П. Костомаров, А.П. Фаворский - Вводные лекции по численным методам (pdf) (1113729) страница 82019-05-05СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Метод итераций (метод последовательных приближений).В этом параграфе мы познакомимся еще с одним численным методом решенияуравнений. Предположим, что уравнение можно записать в виде(19)x = ϕ ( x).Возьмем произвольную точку x0 из области определения функции ϕ ( x ) и будемстроить последовательность чисел { xn } , определенных с помощью рекуррентнойформулы(20)xn+1 = ϕ ( xn ) .Последовательность { xn } называется итерационной последовательностью. При ееизучении встают два вопроса:- 35 1.

Можно ли процесс построения последовательности { xn } продолжатьнеограниченно, т. е. будут ли эти числа принадлежать области определения функцииϕ ( x) ?2. Если итерационная последовательность (20) бесконечна, то как ведут себя еечлены при n → ∞ ?Ответ на оба вопроса дает следующая теорема.Теорема о сходимости итерационной последовательности.Пусть c - корень уравнения (19) и пусть функция ϕ ( x ) удовлетворяет на отрезке[c − δ , c + δ ] условию Липшица с константой L < 1y2 − y1 = ϕ ( x2 ) − ϕ ( x1 ) ≤ L x2 − x1 , L < 1 .(21)Тогда при любом выборе x0 на отрезке [ c − δ , c + δ ] существует бесконечнаяитерационная последовательность { xn } (20), сходящаяся к корню уравнения (19)x = c . Этот корень является единственным на отрезке [ c − δ , c + δ ] .Напомним известный факт из математического анализа: выполнение условияЛипшица (21) будет заведомо обеспечено, если предположить, что функция ϕ ( x )имеет на отрезке [ c − δ , c + δ ] непрерывную производную, ϕ ′ ( x ) модуль которойменьше единицы: ϕ ′ ( x ) ≤ m < 1 .

В этом случае согласно формуле конечныхприращений Лагранжа будем иметь(22)y2 − y1 = ϕ ′ (ξ )( x2 − x1 ) ≤ m x2 − x1 .Мы получили неравенство (21) с константой Липшица L = m .После этого замечания перейдем к доказательству теоремы. Число c являетсякорнем уравнения (19), так что c = ϕ ( c ) . Возьмем произвольную точку x0 на отрезке[c − δ , c + δ ] .

Она отстоит от точки c не больше, чем на δ : x0 − c ≤ δ .Вычислим x1 = ϕ ( x0 ) . При этом будем иметьx1 − c = ϕ ( x0 ) − ϕ ( c ) ≤ L x0 − c ≤ Lδ .(23)Неравенство (23) показывает, что точка x1 принадлежит отрезку [ c − δ , c + δ ] ирасположена ближе к точке c чем x0 .Продолжим построение итерационной последовательности. Вычислимx2 = ϕ ( x1 ) .

При этомx2 − c = ϕ ( x1 ) − ϕ ( c ) ≤ L x1 − c ≤ L2 x0 − c ≤ L2δ .Точка x2 тоже принадлежит отрезку [ c − δ , c + δ ] и расположена ближе к точке c чемx1 . На второй итерации мы опять приблизились к c .По индукции легко доказать, что все последующие итерации также существуюти удовлетворяют неравенствам- 36 Отсюда следует, чтоxn − c ≤ Ln x0 − c ≤ Lnδ .(24)lim ( xn − c ) = 0 , т.

е. lim xn = c .(25)n →∞n→∞Нам остается доказать, что корень x = c является единственным решениемуравнения (19) на отрезке [ c − δ , c + δ ] . Действительно, предположим, что существуетеще один корень x = c1(26)c1 = ϕ ( c1 ) , c − δ ≤ c1 ≤ c + δ .Примем c1 за нулевое приближение и будем строить итерационнуюпоследовательность (20). С учетом (26) получим xn = c1 , n = 0,1, 2,K . С другойстороны, по доказанному lim xn = c , т.

е. c1 = c . Никаких других решений, кроме x = c ,n→∞уравнение (19) на рассматриваемом отрезке не имеет.Доказанная теорема имеет простой смысл. Будем говорить, что функция ϕ ( x )осуществляет отображение точки x отрезка [ c − δ , c + δ ] на точку y = ϕ ( x ) .Рассмотрим пару точек x1 , x2 и их образы y1 , y2 .Условие Липшица (21) приводит ктому, что расстояние между образами оказывается меньше расстояния междуисходными точками, т. е. отображение y = ϕ ( x ) является сжимающим. Корень c неподвижная точка отображения: c = ϕ ( c ) . В результате каждый шаг в итерационномпроцессе, сжимая расстояние, приближает очередную итерацию к корню.Центральная идея метода итераций – сжимающие отображения – являетсявесьма общей.

Например, одно из доказательств теоремы существования иединственности решения задачи Коши для обыкновенных дифференциальныхуравнений основано на методе последовательных приближений в условиях, когдадействует принцип сжимающих отображений. Для многих сложных нелинейных задачпринцип сжимающих отображений оказывается основным методом исследования.Задача 2.Записать уравнение (16) в видеx = cos x ,и найти приближенное значение его корня методом итераций.В данном случае(27)ϕ ( x ) = cos x , ϕ ′ ( x ) = − sin x .На отрезке [ 0,1] , на котором расположен корень уравнения (27), для модуляпроизводной справедлива оценкаϕ ′ ( x ) ≤ sin1 < 0.85Она обеспечивает выполнение условия Липшица с константой Липшица L = 0.85 .Результаты вычислений по рекуррентной формулеxn+1 = cos xn- 37 даны в таблице 2.

За нулевое приближение выбрана средняя точка отрезка x0 = 0.5 .Для удобства анализа итерационной последовательности ее члены расположены подва в строке. В результате образовались столбцы членов с четными и нечетныминомерами. Сравнивая их между собой, мы видим, что четные члены меньше нечетных:итерационная последовательность «скачет» то вверх, то вниз. С возрастанием номерачетные члены последовательности возрастают, а нечетные убывают, приближаясьдруг к другу. Такое поведение последовательности означает, что корень уравнениялежит между ними.

Четные члены дают его значение с недостатком, нечетные – сизбытком. Это позволяет легко контролировать точность, достигнутую после любогочисла итераций: погрешность не превышает разности между последним нечетным ичетным членами. Мы остановили процесс вычислений на 19-ой итерации и можемнаписать для корня c двойное неравенствоx18 = 0.738912449332 < c < x19 = 0.739201444135или, отбрасывая лишние десятичные знаки,0.73891 < c < 0.73921 .Таким образом, члены итерационной последовательности x18 и x19 определяют c снедостатком и с избытком с погрешностью, которая не превышает разности x19 − x18 :ε < ∆19 = x19 − x18 < 0.0003 .Точность, которой мы достигли после 19 итераций, примерно соответствует точности12 шагов метода вилки.

Причина такого различия ясна. В обоих методах погрешностьубывает по закону геометрической прогрессии. Для метода вилки знаменательпрогрессии равен 1/2. Он не зависит от вида функции f ( x) . Для метода итерацийзнаменатель прогрессии равен константе Липшица. В рассматриваемом примереL = 0.85 . Поэтому скорость сходимости метода итераций медленнее скоростисходимости метода вилки.

Метод итераций имеет преимущество перед методом вилкив скорости сходимости только при L < 1/ 2 .Таблица 2n0123456789x2n0,5000000000000,6390124941650,6947780267880,7191654459420,7300810631380,7350063090150,7372357254420,7382462383320,7387045393570,738912449332x2 n+10,8775825618900,8026851006820,7681958312820,7523557594220,7451203413510,7418265226430,7403296518780,7396499627700,7393414522810,739201444136- 38 §3. Метод касательных (метод Ньютона).Метод касательных, связанный с именем Ньютона, является одним из наиболееэффективных численных методов решения уравнений.

Идея метода очень проста.Предположим, что функция f ( x) , имеющая корень c на отрезке [ a, b ] ,дифференцируема на этом отрезке и ее производная f ′( x) не обращается на нем вноль. Возьмем произвольную точку x0 ∈ [ a, b ] и запишем уравнение касательной кграфику функции f ( x) в этой точке(28)y = f ( x0 ) + f ′( x0 ) ( x − x0 ) .График функции f ( x) и ее касательной близки около точки касания, поэтомуестественно ожидать, что точка x1 пересечения касательной с осью x будетрасположена недалеко от корня c (см.

рис. 2). Для определения точки x1 имеемуравнениеf ( x0 ) + f ′( x0 ) ( x1 − x0 ) = 0 ,согласно которомуf ( x0 ).x1 = x0 −f ′ ( x0 )Повторим проделанную процедуру: напишем уравнение касательной к графикуфункции f ( x) в точке x1 и найдем для нее точку пересечения x2 с осью x (см. рис. 2):f ( x1 ).x2 = x1 −f ′ ( x1 )Продолжая этот процесс, получим последовательность { xn } , определенную спомощью рекуррентной формулыf ( xn ).(29)xn+1 = xn −f ′ ( xn )При ее исследовании, как и при исследовании последовательности (20) методаитераций, встают два вопроса:1.

Можно ли процесс вычисления чисел { xn } по рекуррентной формуле (29)продолжать неограниченно, т. е. будут ли эти числа принадлежать отрезку [ a, b ] ?2. Если процесс (29) бесконечен, то как ведет себя последовательностьn → ∞?{ xn }приПри анализе этих вопросов предположим, что корень x = c является внутреннейточкой отрезка [ a, b ] , а функция f ( x) дважды непрерывно дифференцируема наданном отрезке, причем ее производные удовлетворяют неравенствам(30)f ′( x) ≥ m > 0 , f ′′( x) ≤ M , x ∈ [ a, b ] .- 39 Следует обратить внимание на то, что в неравенствах (30) величина m дает оценкумодуля первой производной f ′( x) снизу, а величина M оценку модуля второйпроизводной f ′′( x) сверху.Теорема о сходимости метода касательных.Если функция f ( x) удовлетворяет сформулированным условиям, то найдется такоеδ : 0 < δ ≤ min( c − a , b − c ) , что при любом выборе начального приближения x0 наотрезкесуществуетбесконечнаяитерационная[ c − δ , c + δ ] ⊂ [ a, b ]последовательность (29) и эта последовательность сходится к корню c .В силу предположения о дифференцируемости функции f ( x) и неравенственулю ее производной, уравнение (1) эквивалентно на отрезке [ a, b ] уравнениюf ( x),(31)f ′( x )так что корень x = c исходного уравнения является одновременно корнем уравнения(31).

Исследуем возможность отыскания этого корня с помощью метода итераций.Вычислим и оценим производную функции ϕ ( x ) (31):x = ϕ ( x ) , где ϕ ( x ) = x −( f ′ ( x ) ) − f ( x ) f ′′ ( x ) = f ( x ) f ′′ ( x ) ,ϕ ′( x ) = 1 −( f ′( x ))( f ′( x ))222(32)M(33)f ( x) .m2Теперь воспользуемся непрерывностью функции f ( x) и ее равенством нулю в точкеx = c . Возьмем ε = m 2 / ( 2 M ) .

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

Список файлов лекций

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