Главная » Просмотр файлов » Учебник - Практические занятия по вычислительной математике - Аристова

Учебник - Практические занятия по вычислительной математике - Аристова (1238839), страница 15

Файл №1238839 Учебник - Практические занятия по вычислительной математике - Аристова (Учебник - Практические занятия по вычислительной математике - Аристова) 15 страницаУчебник - Практические занятия по вычислительной математике - Аристова (1238839) страница 152020-10-28СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Можем заметить лишь, что в достаточно малой окрестности решения скорость сходимости метода квадратичная....IV.8. Критерии сходимости итерацийВ методах, использующих локализацию корня, останов итераций призаданной точности  > 0 вычисления положения корня осуществить просто: когда длина очередного отрезка локализации| an – bn | < ,итерационный процесс останавливается.Для одноточечных линейно сходящихся методов (МПИ) итерационный процесс следует остановить при выполнении оценки|| x*  x n ||  .(8.1)Так как точное решение x* неизвестно, то это условие на практике вявном виде проверить невозможно. Здесь поступают следующим образом.Введем величину ||xn – xn – 1||, которая называется итерационной поправкой.

Иногда требуют просто, чтобы итерационная поправка не превосходила заданного значения точности:|| xn  xn 1 ||  .(8.2)В реальных больших задачах скорость сходимости линейного итерационного процесса может быть весьма медленной, поэтому критерий (8.2)не является удовлетворительным. Для обеспечения критерия сходимости(8.1) по итерационной поправке в предположении сходимости процессаиспользуем следующий прием:94IV.9. Задачи на доказательство|| xn  x* ||  || xn  x ||  || xn  xn 1  xn 1  xn  2  xn  2  xn 3  ...

 x ||  || xn  xn 1 ||  || xn 1  xn  2 ||  || xn  2  xn 3 ||  ...  || x n  x n 1 ||  (1  q  q 2  ...)  || x n  x n 1 ||(8.3)11 qТаким образом, мы видим, что для достижения заданной точности критерий (8.1) в терминах итерационной поправки нужно переформулироватьследующим образом:|| xn  xn 1 ||   (1  q) ,(8.4)что при q, близких к единице, может существенно увеличить количествотребуемых итераций для достижения заданной точности по сравнению скритерием (8.2).Величина q в оценке (8.4) должна уточняться на итерациях:q || xn  xn 1 || / || xn 1  xn 2 || .(8.5)В случае достаточной малости величины q в (8.4) можно пренебречьмножителем (1 – q) и перейти к оценке (8.2).Для метода Ньютона нахождения простого корня нелинейного уравнения выполнение условия (8.2) на итерационную погрешность являетсядостаточным в силу квадратичной сходимости метода.

Для кратных корней необходимо переходить к оценке (8.4).Можно использовать также критерий сходимости по самой функции.Считается, что итерационный процесс сошелся, если выполнено условие||F(xn)|| < ε.Для кратных корней такой критерий является неудовлетворительным.IV.9. Задачи на доказательствоIV.9.1. Доказать формулу (4.4).IV.9.2. Доказать, что метод простой итерации x(n+1) = φ(x(n)) сходится длялюбого начального приближения x0  R:а) φ(x) = cos (x);б) φ(x) = a sin2x + b cos2x +c, где |a – b| < 1;в) φ(x) = a exp(– bx2) + c, где b  0, 2a 2b  e .95IV.

ЧИСЛЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙIV.9.3. Написать формулу метода Ньютона для решения уравнения x2 = a,a > 0. Доказать, что если за начальное приближение принять произвольное x0 > 0, то все последующие приближения xk, k ≥ 1, больше, чем a .IV.9.4. Известно, что уравнение F(x) = 0 имеет решение на отрезкеα ≤ x ≤ β. Причем на этом отрезке Fʹ(x) > 0, Fʹʹ(x) > 0. Показать, что задачуможно решить численно методом Ньютона, положив x0 = β. Построитьпример, в котором при выборе в качестве начального приближения числаx0 = α сходимость не имеет места.IV.9.5. Известно, что уравнение F(x) = 0 имеет решение на отрезкеα ≤ x ≤ β. Как следует выбирать начальные приближения x0, чтобы былагарантирована сходимость метода Ньютона в следующих случаях: на отрезке   x   всюду F   x   0 , F   x   0 , на отрезке   x   всюду F   x   0 , F   x   0 , на отрезке   x   всюду F   x   0 , F   x   0 .IV.9.6.

Доказать, что уравнение x + 0,5 sin x + a = 0 имеет единственныйкорень при любом a. Найти его значение с точностью ε = 10- 3 дляa = ±1, ±2, ±3.IV.9.7. Какие условия теоремы о сходимости простых итераций для уравнения x = φ(x), φ(x) = x/2 – 4/5, x0 = 0, на отрезке [–1, 1], не выполнены?Будет ли это уравнение иметь решение на этом отрезке?IV.9.8. Пусть уравнение f(x) = 0 имеет на отрезке [a, b] корень z кратностиp, причем f (x) – дважды дифференцируемая функция.а) Показать, что при этих условиях метод Ньютона сходится линейносо скоростью геометрической прогрессии со знаменателем (p – 1)/p.б) Построить модификацию метода Ньютона, имеющую квадратичную скорость сходимости.Указание. Искать указанную модификацию среди методов видаf ( x(n) )( n 1)и найти значение параметра α, обеспечивающееx x( n)  f ( x ( n ) )квадратичную сходимость.IV.9.9. Для нахождения простого нуля z функции f ( x)  C 4 используетсяитерационный процессx( n 1)  0.5(u ( n 1)  v( n 1) ) , где96IV.10.

Задачи с решениямиf ( x( n) )g ( x( n) )f ( x)( n 1)(n),vx, g ( x) .f ( x)f ( x( n ) )g ( x ( n ) )Доказать, что если метод сходится, то скорость – кубическая.u ( n 1)  x ( n ) IV.9.10. Пусть уравнение f (x) = 0 имеет единственный корень на отрезке[a, b]. Для его численного решения используются метод деления отрезкапополам и другой похожий метод, в котором отрезок делится на три равные части. Пусть каждое вычисление функции выполняется за время О(1),а сравнение знаков выполняется мгновенно. Доказать, что метод деленияпополам сходится к корню быстрее.IV.10.

Задачи с решениямиIV.10.1. (Е.Н. Аристова) Методом простой итерации найти ширину функции y = x e –x , x > 0 на полувысоте с точностью ε = 10 –2.Решение. Найдем максимальное значение функции:y' = e –x – x e –x обращается в ноль при x = 1, при этом максимальное значение функции ymax = e –1 (рис. 4.5).Рис. 4.5 К задаче IV.10.1Необходимо найти два значения аргумента, при которых функция принимает значения, равные половине от максимального:x e –x = 1/2e.(10.1)Уравнение (10.1) может быть очевидным образом приведено к форме метода простой итерации двумя способами:1 спо со б .

xk + 1 = 0.5exp(xk – 1), в этом случае правая часть метода простой итерации f(x) = 0.5exp(xk – 1). Для сходимости метода необходимовыполнение условия |f (x)| < 1:97IV. ЧИСЛЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙf (x) = 0.5exp(xk – 1), и эта величина меньше единицы при выполненииусловия x < 1 + ln2 ≈ 1.7.Этот итерационный процесс пригоден для поиска левого корня уравнения (10.1).2 спо со б . xk + 1 = ln(2exk), f (x) = 1/x, производная |f (x)| < 1 при x > 1.Этот итерационный процесс пригоден для поиска правого корняуравнения (10.1).Каждый из корней должен быть найден с точностью ε/2.

С учетомкритерия сходимости (8.4) имеем для обоих итерационных процессов,стартующих от точки максимума функции:Правый кореньx0 = 1.x1 = 1.69314x2 = 2.21973x3 = 2.49053x4 = 2.60564x5 = 2.65083x6 = 2.66802x7 = 2.67448x8 = 2.67690x9 = 2.67780x10 = 2.67814Левый кореньx0 = 1.x1 = 0.5x2 = 0.30326x3 = 0.24910x4 = 0.23597x5 = 0.23289x6 = 0.23217Таким образом, искомая величина 1/2  x2  x1  2.678  0.232  2.446 .IV.10.2. Оценить число итераций метода Ньютона для нахождения положительного корня уравнения sinx + x 2 –1 = 0 с точностью ε = 10 –4.Решение. Для начала локализуем корень.

Графически можно показать,что функции sin x и 1 – x 2 пересекаются на отрезке [0, 1]. Для более точного определения отрезка локализации для оценки производных заменимв уравнении sin x ≈ x:x 2 + x –1 = 0, положительный корень x = 0.5(–1 + 5 ) ≈ 0.618. В качествеотрезка локализации возьмем интервал [π/6, π/4] а в качестве начальногоприближения величину x0 ≈ 0,7.98IV.10. Задачи с решениямиТогдаM1 F ( x) ( F ( x))1   inf /6 x/4 /6 x /4sup  inf2 x  cos x  /6 x/4M 2  sup F ( x) /6 x/4111 0.4063 . / 3  cos  / 4sup 2  sin x  2  sin  / 4  1.2929 Тогда/6 x/4показатель квадратичного убыванияq = 0.5 M1M2 = 0.2626.По теореме 2 справедлива оценка| xn  x* |   q 2n 1n| x0  x* |2  .Подставляя грубую оценку |x0 – x*| ≈ 0,5, будем иметь для оценки необходимого числа итераций 0.5  0.26262n1 2 , или  0.13132n 1 2 104 , n  3 .Начиная от выбранного начального приближения x0 ≈ 0.7, получим последующие приближения метода Ньютона:x1 ≈ 0.63801, x2 ≈ 0.63673, x3 ≈ 0.63673.IV.10.3.

Доказать, что если методx n 1  x n F ( xn )Fx ( x n )F ( x n ) F 2 ( x n )2 Fx ( x n )3сходится, то порядок сходимости – третий.Решение. Используем представление этого метода в качестве методапростой итерации xn+1 = f(xn) с функциейf ( x)  x F ( x) F ( x) F 2 ( x).Fx ( x) 2  F  ( x) 3xВычислим производные этой функции в ряд Тейлора в корне нелинейного уравнения F(x*) = 0:99IV. ЧИСЛЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙf ( x)  1 F 2  FF F 2F F 2 F   2F F 2 F  3F 2 F 2 F 22 F 4,илиf ( x) FF F 2F F 22 F 3F F 2 F   2 F F 2 F  3F 2 F 22 F 4 F  3F 23F 2 F 2 F2  2 F 3 2 F 42 F 4.Множитель F 2 перед скобкой обеспечивает второй порядок нуляпроизводной f ′(x) в корне, поэтому сходимость не ниже квадратичной.Вторую производную правой части МПИ в корне можно не вычислять – она обращается в нуль при x = x* в силу множителя F2 перед скобкой.

Таким образом, сходимость не хуже кубической. Можно показать,что третья производная в корне в нуль не обращается, т.е. если сходимость имеет место, то с кубической скоростью.IV.10.4. В середине XIX века Ферхюльст для описания динамики популяционной системы предложил измерять ежегодно численность особейuk, где k – номер года. Относительная численность uk+1 полагалась пропорциональной численности в k год, однако она начинает убывать, когдаживотных становится много (uk сравнимо с 1):uk 1  uk (1  uk ), u0  aПусть численность популяции нормирована на максимальную возможную численность. Так как uk  [0,1] , то   [0, 4] .

Функция в правой части называется логистическим отображением.Исследовать свойства логистического отображения при значенияхпараметра 0    1  6. Представить графически последовательностьитераций логистического отображения при разных значениях параметра.Решение. Рассмотрим подробнее свойства отображенияuk 1  uk (1 uk ), u0  a.Заметим, что если f (0) = f (1) = 0 и max f (u) = f (0,5) = /4, то при0 < λ <4 интервал Х = [0, 1] отображается в себя, u  X.Рассмотрим вначале случай 0 < λ <1.

На Х = [0, 1] существует только одна предельная (или неподвижная) точка x = 0. Любая последова-100IV.10. Задачи с решениямительность, { f k (u0 )}k 0 сходится к предельной точке рассматриваемогоотображения x = 0. Если рассматривается популяционная модель, то этоозначает, что рассматриваемая популяция не может выжить.Из теоремы о сжимающем отображении следует, что последова-тельность {uk }n 0 сходится к своей предельной точке, если |fu′| ≤ 1.

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

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

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

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