- •Классификация методов решения слау
- •[Править]Прямые методы
- •[Править]Итерационные методы
- •Погрешность численного дифференцирования
- •Итерационные методы решения слау Итерационные методы решения линейных алгебраических систем Метод простой итерации или метод Якоби
- •Метод Гаусса – Зейделя
- •Метод Гаусса решения слау
- •[Править]Условие совместности
- •Метод Эйлера
- •[Править]Оценка погрешности
- •[Править]Значение метода Эйлера
- •Метод прогонки решения слАу
- •Метод Эйлера с пересчетом
- •Метод Гаусса-Зейделя
- •Одношаговые методы
- •Решение нелинейных уравнений
- •Многошаговые методы
- •Метод деления отрезка пополам для решения нелинейного уравнения
- •Методы Рунге-Кутта
- •Метод хорд для решения нелинейных уравнений Метод хорд
- •Задача Коши для разрешения обыкновыеннх дифф. Уравнений
- •Метод Ньютона Метод Ньютона (метод касательных)
- •Краевая задача для уравнений
- •2. Теоретическая справка
- •2.1. Пример краевой задачи
- •Метод простой итерации Метод простых итераций
- •Основные понятия метода сеток
- •Комплексные корни
- •Многошаговый метод Адамса
- •Метод простой итерации для решения системы
- •Задача Дирихле для уравнения Лапласа
- •Приближения функции. Аппроксимация.
- •Явные и неявные разностные схемы
- •Явные схемы
- •[Править]Неявные схемы
- •Полунеявные схемы
- •Интерполяция
- •Определения
- •[Править]Пример
- •Интерполяция методом ближайшего соседа
- •[Править]Интерполяция многочленами
- •[Править]Обратное интерполирование (вычисление X при заданном y)
- •[Править]Интерполяция функции нескольких переменных
- •[Править]Другие способы интерполяции
- •Аппроксимация, сходимость, устойчивость разностной схемы
- •1.2.7 Аппроксимация.
- •Линейная интерполяция
- •[Править]Геометрическая интерпретация
- •[Править]Применение
- •Уравнение Пуассона
- •Квадратичная интерполяция
- •Уравнение теплопроводности
Интерполяция методом ближайшего соседа
Простейшим способом интерполяции является интерполяция методом ближайшего соседа.
[Править]Интерполяция многочленами
На практике чаще всего применяют интерполяцию многочленами. Это связано прежде всего с тем, что многочлены легко вычислять, легко аналитически находить их производные и множество многочленов плотно в пространстве непрерывных функций (теорема Вейерштрасса).
Линейная интерполяция
Интерполяционная формула Ньютона
Метод конечных разностей
ИМН-1 и ИМН-2
Многочлен Лагранжа (интерполяционный многочлен)
По схеме Эйткена
Сплайн-функция
Кубический сплайн
[Править]Обратное интерполирование (вычисление X при заданном y)
Полином Лагранжа
Обратное интерполирование по формуле Ньютона
Обратное интерполирование по формуле Гаусса
[Править]Интерполяция функции нескольких переменных
Билинейная интерполяция
Бикубическая интерполяция
[Править]Другие способы интерполяции
Рациональная интерполяция
Тригонометрическая интерполяция
Аппроксимация, сходимость, устойчивость разностной схемы
Сходимость разностных схем.
Будем говорить, что разностная схема (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_________________________________