Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Chislennye_metody 2.docx
Скачиваний:
15
Добавлен:
09.09.2019
Размер:
1.03 Mб
Скачать

Интерполяция методом ближайшего соседа

Простейшим способом интерполяции является интерполяция методом ближайшего соседа.

[Править]Интерполяция многочленами

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

  • Линейная интерполяция

  • Интерполяционная формула Ньютона

  • Метод конечных разностей

  • ИМН-1 и ИМН-2

  • Многочлен Лагранжа (интерполяционный многочлен)

  • По схеме Эйткена

  • Сплайн-функция

  • Кубический сплайн

[Править]Обратное интерполирование (вычисление X при заданном y)

  • Полином Лагранжа

  • Обратное интерполирование по формуле Ньютона

  • Обратное интерполирование по формуле Гаусса

[Править]Интерполяция функции нескольких переменных

  • Билинейная интерполяция

  • Бикубическая интерполяция

[Править]Другие способы интерполяции

  • Рациональная интерполяция

  • Тригонометрическая интерполяция

  1. Аппроксимация, сходимость, устойчивость разностной схемы

Сходимость разностных схем.

Будем говорить, что разностная схема (S) сходится (к решению  задачи (E) – (C)), если

||  Pt||t 0 при  0.

(4)

Сходимость разностной схемы означает, что при достаточно малом  значения сеточного (приближенного) решения  и точного решения  мало отличаются. Соотношение (4) на практике оказывается мало полезным, поскольку на основании его нельзя судить о том насколько малым мы должны выбрать шаг , чтобы в узлах сетки точное и приближенное решения отличались друг от друга не более, чем на  (заранее заданную точность). Если удается доказать, что при достаточно малых  > 0

||  Pt||t  Ck,

(5)

где C — не зависящая от  константа, то говорят, что схема (S) сходится с порядком k (или является схемой k-го порядка (сходимости)). Оценка (5), если в ней известна (для конкретной задачи (E) – (C)) константа C, позволяет по заранее выбранной точности  a priori выбрать шаг так, чтобы приближенное решение аппроксимировало решение данной (дифференциальной) задачи Коши с точностью :

||  P||   ;

достаточно взять   k/C.

1.2.7 Аппроксимация.

Явная схема Эйлера обладает двумя важными свойствами, из которых, как будет показано ниже, последует ее сходимость с первым порядком. Во-первых, при достаточно малых 

||Ft(Pt)||t  Ca,

(6)

где Ca — константа, не зависящая от , а  — как обычно, точное решение задачи (E) – (C). В этом случае говорят, что схема (S) имеет первый порядок аппроксимации на решении. (Если в правой части неравенства (6) стоит Cak, то, соответственно, говорят о схеме k-го порядка аппроксимации (на решении).) Другими словами, неравенство(6) эквивалентно тому, что ||Ft(Pt)||t = O() (в случае схемы k-го порядка аппроксимации — O(k)). Тот факт, что разностная схема обладает аппроксимацией на решении, означает, грубо говоря, что при подстановке точного решения дифференциальной задачи в разностную схему мы получаем невязку соответствующего порядка малости по. (Было бы идеально, если бы после такой подстановки мы получали в левой части нуль, однако в общем случае конструктивно такие схемы выписать нельзя.)

Часто вместо свойства аппроксимации на решении рассматривают формально более жесткое требование, которое называют свойством аппроксимации (в зарубежнойлитературе — согласованностью); именно, говорят, что схема (S) является схемой k-го порядка аппроксимации на функции x, если при достаточно малых 

||Ft(Ptx)  PtF(x)||t  Cak.

Обычно требуется, чтобы схема обладала свойством аппроксимации на множестве функций из некоторого класса гладкости. Очевидно, если решения дифференциального уравнения (E) m раз непрерывно дифференцируемы, а разностная схема обладает k-м порядком аппроксимации (согласованности) на m раз непрерывно дифференцируемых функциях, то она обладает k-м порядком аппроксимации на решении.

Второе важное свойство, которым обладает явная схема Эйлера, называется устойчивостью и определяется так: если z  St и, кроме того, ||z||t и  достаточно малы, то уравнение

Ft(y) = z

(7)

однозначно разрешимо и существует такая не зависящая от  и ||z|| константа Cs, что

||y  || = ||F1(z) ||  Cs||z||.

(8)

Устойчивость разностной схемы означает, что малые возмущения z в начальных данных и правой части разностной схемы приводят к равномерно малому по  изменению решения (напомним, что  — решение невозмущенной системы, а F1(z) — возмущенной). Поскольку  = F1(z),неравенство (8), переписанное в виде ||F1(z)F1(0)||  Cs||z|| , означает, в частности, непрерывность обратного к разностному оператору оператора F1в нуле.

Устойчивость — очень важное в приложениях свойство разностных схем. При практической реализации на ЭВМ разностных методов возникают, в частности, проблемы, связанные с невозможностью представления точных чисел в компьютере. В результате мы решаем не разностную схему (S), а несколько отличающееся от (S) уравнение. Все такие возмущения в разностной схеме, грубо говоря, можно "перенести в правую часть" и, таким образом, считать, что в ЭВМ ищется решение не разностной схемы(S), но решение возмущенного уравнения (7). Свойство устойчивости разностной схемы гарантирует близость при достаточно малых  между точным (теоретическим) решением  разностной схемы и его практической реализацией F1(z)(где z — суммарный вектор возмущений). Источником возмущений служит не только невозможность точного представления данных в ЭВМ, но и неточность определения физических параметров модели, погрешность измерений и т.п.

БИЛЕТ 20_________________________________

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]