main (1160440), страница 8

Файл №1160440 main (Численные методы. Ионкин (2013) (github clone)) 8 страницаmain (1160440) страница 82019-09-19СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

QR-алгоритм применим к произвольной матрице .Замечание 3. QR-алгоритм является очень затратным по необходимому количествуопераций и объему памяти, используемому для хранения промежуточных матриц.§12Предварительное преобразование матрицы к ВПТФ. Неухудшение ВПТФ при QR-алгоритмеЛемма 1. Пусть = , где имеет ВТФ, а имеет ВПТФ. Тогда имеет ВПТФ.Доказательство. Выпишем элемент матрицы по определению произведения матриц: =∑︁ , , = 1, .=1Учтем, что = 0 при < и = 0 при > + 1: =∑︁= =+1∑︁ , , = 1, .=При > + 1 получим, что = 0. Таким образом, имеет ВПТФ и лемма доказана.Аналогичным образом доказывается следующая лемма (ее непосредственное доказательство предоставляется читателю).Лемма 2.

Пусть = , где — матрица с ВПТФ, а — матрица с ВТФ. Тогда —матрица с ВПТФ.Рассмотрим применение QR-алгоритма для матрицы . Приведем матрицу к верхнейпочти треугольной матрице 0 . Запишем QR-разложение матрицы 0 :0 = 0 0 .Поскольку 0 и 0−1 — матрицы, имеющие ВТФ, то матрица 0 , определяемая выражением0 = 0 0−1 ,в силу леммы 2 имеет ВПТФ.

Матрица 1 = 0 0 в силу леммы 1 также имеет ВПТФ.Таким образом, леммы 1 и 2 гарантируют на каждом шаге QR-алгоритма неухудшениеВПТФ матрицы , ∈ Z+ .Глава 2Интерполирование и приближениефункций§1Постановка задачи интерполированияРассмотрим некоторый технологический процесс, характеризуемый множеством параметров. Разместим в среде протекания процесса конечное число датчиков, позволяющих получать точные значения параметров процесса в ограниченном числе точек среды. Для получения исчерпывающей информации о протекании процесса необходимо уметь оцениватьзначения параметров процесса в точках, в которых нет возможности их измерить.Под интерполированием (точное определение будет дано ниже) понимается процесспоиска промежуточных значений величины по имеющемуся дискретному набору известных значений. В вычислительной математике интерполирование обычно рассматриваетсяв рамках задачи вычисления промежуточных значений функций, например, при вычислении значений специальных функций, являющихся решениями дифференциальных уравнений специального вида (функции Бесселя, Ханкеля и другие).

Как правило, значенияфункций такого рода задаются таблицами, шаг которых может оказаться слишком большим для конкретной задачи. В таком случае используют интерполирование для получениязначений функции с заданной точностью.Интерполирование функций используется при исследовании сходимости разностных методов решения дифференциальных задач. При исследовании сходимости необходимо уметьсравнивать сеточные и непрерывные функции. Эту задачу можно решить двумя методами. Первый метод состоит в проецировании непрерывной функции на сетку и последующемсравнении сеточных функций.

Второй способ состоит в восстановлении непрерывной функции по сеточной с помощью интерполирования и последующем сравнении непрерывныхфункций.Постановка задачи. Рассмотрим вещественную функцию (), ∈ [, ] ⊂ Rи произвольным образом заданное разбиение области определения этой функции: 6 0 < 1 < 2 < . . . < 6 .Точки { }=0 называются узловыми точками функции (). В этих точках задано значение функции: ( ) = , = 0, .Задача интерполирования состоит в нахождении значений функции () на всем отрезке[a,b] по ее значениям в узловых точках.48Глава 2.

Интерполирование и приближение функцийЗамечание. Далее будем считать термины «интерполирование функции» и «приближение функции» синонимами.Заметим, что в постановке задачи интерполирования не указан конкретный метод построения приближенных значений функции (). В силу этого задача допускает скольугодно много решений. В этой главе рассматривается задача приближения заданной функции вещественными полиномами: () = 0 + 1 + 2 2 + . . . + , ∈ R,∑︁2 ̸= 0.=0Определение.

Вещественный полином -ой степени () называется интерполяционным полиномом для функции (), построенным по узлам { }=0 , если его значения вузловых точках совпадают со значениями функции в этих точках: = 0, . ( ) = ,(1)Утверждение. Для любой функции () существует единственный интерполяционныйполином степени , построенный по ( + 1)-му узлу.Доказательство.

Распишем систему (1) покоординатно:⎧⎪0 + 1 0 + 2 20 + . . . + 0 = 0⎪⎪⎪⎨ + + 2 + . . . + = 01 12 1 11⎪...⎪⎪⎪⎩ + + 2 + . . . + = 01 2 .(2)Получим систему линейных уравнений с ( + 1)-им уравнением относительно коэффициентов полинома () с матрицей⎛⎞1 0 20 . . . 0⎜1 1 2 . . . ⎟11⎟⎜ = ⎜. ...

. ... ⎟ ...⎝. .. . ⎠.21 . . . Определитель матрицы — это определитель Вандермонда ( + 1)-ого порядка:|| =∏︁( − ).06<6Поскольку все узлы различны, матрица невырождена: || ≠ 0.Из невырожденности матрицы следует существование и единственность решения системы (2). Таким образом, для любой функции () существует интерполяционный полином (), и его коэффициенты однозначно определяются по значениям функции в ( + 1)ой узловой точке.Замечание. Помимо интерполирования иногда решают задачу экстраполирования — прогнозирования поведения функции за пределами отрезка. Задача экстраполирования имеетбольшую погрешность, чем задача интерполирования.§2.

Интерполяционная формула Лагранжа§249Интерполяционная формула ЛагранжаРассмотрим вещественную функцию ∈ [, ] ⊂ R, (),заданную в узловых точках произвольного разбиения отрезка [, ]: 6 0 < 1 < 2 < . . . < 6 , = 0, . ( ) = ,Определение. Интерполяционный полином для функции (), заданный формулой () =∑︁ = 0, , () ( ),(1)=0где () — полином степени , называется интерполяционным полиномом в форме Лагранжа.Из определения интерполяционного полинома следует, что ( ) = ( ) = , = 0, .Из этих равенств следуют условия, = 0, .

( ) = ,(2)Будем искать полиномы () с учетом этих условий.Рассмотрим полином ( + 1)-ой степени вида() =∏︁( − ).=0Вынесем за скобку множитель ( − ):⎛⎞∏︁⎜⎟() = ( − ) ⎝ ( − )⎠ ,=0̸=продифференцируем по :⎛⎜ ′ () = ( − ) ⎝∏︁⎞′⎟ ⎜( − )⎠ + ⎝=0̸=и подставим в полученное выражение = :⎛⎞⎜∏︁⎟ ′ ( ) = ⎝ ( − )⎠ ,=0̸=⎛∏︁⎞⎟( − )⎠=0̸= = 0, .50Глава 2. Интерполирование и приближение функцийИскомые полиномы () можно представить следующим образом: () =(),( − ) ′ ( ) = 0, .(3)Заметим, что условия (2) для полиномов () выполнены.

Учитывая равенства (1) и (3),запишем интерполяционный полином в форме Лагранжа: () =∑︁=0() ( ).( − ) ′ ( )Оценим точность приближения функции () интерполяционным полиномом в формеЛагранжа.Определение. Пусть () — интерполяционный полином в форме Лагранжа для функции (). Тогда функция () = () − ()(4)называется погрешностью интерполирования функции () интерполяционным полиномом ().Пусть существует ( + 1)-ая производная функции () на отрезке [, ]. Тогда () = (+1) ()(),( + 1)!где ∈ [, ].Обычно оценку погрешности аппроксимации (5) записывают в виде⃒⃒+1⃒ (+1) ⃒()⃒.| ()| 6|()| , где +1 = sup ⃒( + 1)!∈[,](5)(6)Замечание 1. Вывод формул (5) и (6) в данном курсе не рассматривается, его можнонайти в [1].Замечание 2. Если исходная функция является полиномом степени, не превышающей, то интерполяционный полином в форме Лагранжа приближает ее точно.1Замечание 3.

Наличие в оценке погрешности (6) быстро убывающего множителя (+1)!вовсе не гарантирует сходимость интерполяционного полинома в форме Лагранжа к заданной функции при увеличении числа узлов в разбиении. Более того, начальное разбиениеможет быть выбрано так, что мы вовсе не получим сходимости. Поэтому на практике лучше разбивать область определения функции на меньшие отрезки, на каждом изкоторых приближать функцию полиномом невысокой степени, и потом «сшивать» полученные приближения в одну функцию, определенную уже на всем отрезке.§3Разделенные разностиРассмотрим вещественную функцию (), ∈ [, ] ⊂ R,заданную в узловых точках произвольного разбиения отрезка [, ]: 6 0 < 1 < 2 < .

. . < 6 , ( ) = , = 0, .§3. Разделенные разности51Определение. Разделенной разностью первого порядка, построенной по несовпадающимузлам и , называется отношение ( , ) = ( ) − ( ), − 0 6 , 6 .(1)Обычно мы будем рассматривать разделенные разности, составленные по соседним узлам.Например, (1 ) − (0 ) (2 ) − (1 ) (0 , 1 ) =, (1 , 2 ) =.1 − 02 − 1Замечание. Отношение (1) является дискретным аналогом первой производной.Определение. Разделенной разностью второго порядка, построенной по несовпадающимузлам 0 , 1 , 2 , называется отношение (0 , 1 , 2 ) = (1 , 2 ) − (0 , 1 ).2 − 0(2)Замечание.

Отношение (2) является дискретным аналогом второй производной.Определение. Пусть даны ( , . . . , + ) и (+1 , . . . , ++1 ) — разделенные разности-ого порядка по соответствующим узлам, где 0 6 , 6 . Тогда разделенной разностью ( + 1)-ого порядка, построенной по несовпадающим узлам , +1 , . . . , ++1 ,называется отношение ( , +1 , . . . , ++1 ) = (+1 , +2 , .

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

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

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

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