Глобальная и кусочно-полиномиальная интерполяция - раздел Математика, ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА Пусть Функция F(X) Интерполируется На Отрезке [A, B]. Метод Реш...
Пусть функция f(x) интерполируется на отрезке [a, b]. Метод решения этой задачи с помощью единого многочлена 
для всего отрезка называют глобальной полиномиальной интерполяцией. В вычислительной практике такой поход применяется редко в силу различных причин. Одна из причин v необходимо задать стратегию выбора узлов при интерполяции функции f многочленами все возрастающей степени n.
Теорема Фабера. Какова бы ни была стратегия выбора узлов интерполяции, найдется непрерывная на отрезке [a, b] функция f, для которой 
при 
. Таким образом, теорема Фабера отрицает существование единой для всех непрерывных функций стратегии выбора узлов. Проиллюстрируем сказанное примером.
Предположим, что выбираем равноотстоящие узлы, то есть 
, i = 0, 1,...n, где 
. Покажем, что для функции Рунге такая стратегия является неудачной.
ПРИМЕР 1. Глобальная интерполяция функции Рунге (рис. 10.1).
% Интерполяция функции Рунге
% Введём функцию Рунге
f = inline('1./(1+25*x.^2)');
% Вычислим таблицу значений
x = linspace(-1, 1, 10);
y = f(x);
% Проинтерполируем функцию Рунге многочленами Лагранжа
p = polyfit(x, y, 10);
xx = linspace(-1, 1, 100);
yy = polyval(p, xx);
axes('NextPlot', 'Add');
% Покажем, что глобальная аппроксимация плохо работает для функции Рунге
plot(x, y);
plot(xx, yy, 'Color', 'r');
Рис. 10.1 - Глобальная интерполяция функции Рунге
% С увеличением узлов сетки, ситуация только ухудщается
% Вычислим таблицу значений. 20 узлов сетки
x = linspace(-1, 1, 20);
y = f(x);
% Проинтерполируем функцию Рунге многочленами Лагранжа
p = polyfit(x, y, 20);
xx = linspace(-1, 1, 100);
yy = polyval(p, xx);
figure
axes('NextPlot', 'Add');
% Покажем, что глобальная аппроксимация плохо работает для функции Рунге
plot(x, y);
plot(xx, yy, 'Color', 'r');
Рис. 10.2 - интерполяция функции Рунге многочленами Лагранжа
На практике чаще используют кусочно-полиномиальную интерполяцию: исходный отрезок разбивается на части и на каждом отрезке малой длины исходная функция заменяется многочленом невысокой степени. Система Mathcad предоставляет возможность аппроксимации двумя важными классами функций: кусочно-линейной и сплайнами.
При кусочно-линейной интерполяции узловые точки соединяются отрезками прямых, то есть через каждые две точки 
, 
проводится полином первой степени.
ПРИМЕР 2. Кусочно-линейная интерполяция функции Рунге.
% Кусочно-линейная интерполяция функции Рунге
% Введём функцию Рунге
f = inline('1./(1+25*x.^2)');
% Вычислим таблицу значений
x = linspace(-1, 1, 10);
y = f(x);
% Начертим график кусочно-линейной аппроксимации (рис. 10.3)
plot(x, y);
Рис. 10.3 - график кусочно-линейной аппроксимации
Как видно из приведенного примера этот способ приближения также имеет недостаток: в точках «стыка» двух соседних многочленов производная, как правило, имеет разрыв.
Если исходная функция была гладкой и требуется, чтобы и аппроксимирующая функция была гладкой, то кусочно-полиномиальная интерполяция неприемлема. В этом случае применяют сплайны v специальным образом построенные гладкие кусочно-многочленные функции.
Все темы данного раздела:
Таганрог 2007
УДК 681.326.3 – 181.48(075.8)
Синютин С.А. Вычислительная математика. Учебное пособие по дисциплине «Вычислительная математика» – Таганрог: Изд-во ТТИ ЮФУ, 2007. – 123
Ввод матриц
Вы можете вводить матрицы в MATLAB несколькими способами:
• вводить полный список элементов;
• загружать матрицы из внешних файлов;
• генерировать матрицы, используя встр
Индексы
Элемент в строке i и столбце j матрицы А обозначается A(i,j). Например, А(4,2) - это число в четвертой строке и втором столбце. Таким образом, можно вычислить сумму элементов в четвертом столбце м
Функции
MATLAB предоставляет большое количество элементарных математических функций, таких как abs, sqrt, exp, sin. Вычисление квадратного корня или логарифма отрицательного числа не является ошибко
Редактор командной строки
Различные стрелки и управляющие клавиши на вашей клавиатуре позволяют вам вызывать, редактировать и многократно использовать команды, набранные ранее. Например, предположим, что вы допустили ошибку
Создание графика
Функция plot имеет различные формы, связанные с входными параметрами, например plot(y) создает кусочно-линейный график зависимости элементов у от их индексов. Если вы задаете два вект
Управление осями
Функция axis имеет несколько возможностей для настройки масштаба, ориентации и коэффициента сжатия.
Обычно MATLAB находит максимальное и минимальное значение и выбирает соответствую
Печать графики
Опция Printв меню Fileи команда print печатают графику MATLAB. Меню Printвызывает диалоговое окно, которое позволяет выбирать общие стандар
Команда help
Команда help – это самый основной способ определения синтаксиса и поведения отдельных функций. Информация отображается прямо в командном окне. Например,
help magic
Сценарии и функции
MATLAB - это мощный язык программирования, также как и интерактивная вычислительная среда. Файлы, которые содержат код на языке MATLAB, называются М-файлами. Вы создаете М-файлы, используя текстово
Сценарии
Когда вы вызываете сценарий, MATLAB просто вызывает команды, содержащиеся в файле. Сценарии могут оперировать существующими данными в рабочем пространстве или они могут сами создавать эти данные. Х
Функции
Функции - это М-файлы, которые могут иметь входные и выходные возвращать. Имя М-файла и функции должно быть одним и тем же. Функции работают с переменными в пределах их собственного рабочего простр
Теория погрешностей и машинная арифметика.
Пусть - точное значение, - приближенное значение некоторой величины
Локализация корней
ПРИМЕР 1. Локализация корней (рис. 4.1).
% Локализовать корни уравнения f(x)=0, где f(x)= x^3 - cos(x) + 1
% Введём функцию f(x)
f = inline('x.^3 - cos(x)
Метод бисекции
Пусть[a,b] v отрезок локализации. Предположим, что функция f(x) непрерывна на [a,b] и на концах принимает значения разных знаков .
Метод Ньютона (метод касательных)
Расчетная формула метода Ньютона имеет вид:
. Геометрически метод Ньютона означает, что следующее приближение к кор
Метод простой итерации (метод последовательных повторений)
Для применения метода простой итерации следует исходное уравнение преобразовать к виду, удобному для итерации
Обусловленность задачи нахождения корня.
Пусть v корень, подлежащий определению. Будем считать, что входными данными для задачи вычисления корня являются значения функции
Интервал неопределенности корня.
Если функция непрерывна, то найдется такая малая окрестность корня
Применение метода Ньютона для нахождения кратного корня.
Метод Ньютона для случая кратного корня обладает лишь линейной скоростью сходимости. Чтобы сохранить квадратичную сходимость его модифицируют следующим образом:
Методика решения алгебраического уравнения
Мы остановимся здесь подробнее на методике решения алгебраического уравнения, т.е. уравнения вида: , левую часть которого будем обозначать т
Нормы векторов и матриц
Обозначим через - точное решение системы, а через - приближенное решение системы. Д
Обусловленность задачи
Так же как и другие задачи, задача вычисления решения системы может быть как хорошо обусловленной, так и плохо обусловленной.
Теорема об оценке погрешности решения по погрешностям входн
Метод Гаусса
Рассмотрим метод Гаусса (схему единственного деления) решения системы уравнений. Прямой ход состоит из m-1 шагов исключения.
1. Шаг. Исключим неизвестное
LU разложение матрицы
Представим матрицу A в виде произведения нижней треугольной матрицы L и верхней треугольной U.
Введем в рассмотрение матрицы
Метод Гаусса с выбором главного элемента по столбцу
На k-ом шаге прямого хода в качестве ведущего элемента выбирают максимальный по модулю коэффициент при неизвестной
Метод Холецкого
Если матрица системы является симметричной и положительно определенной, то для решения системы применяют метод Холецкого (метод квадратных корней). В основе метода лежит алгоритм специального LU
Метод прогонки
Если матрица системы является разреженной, то есть содержит большое число нулевых элементов, то применяют еще одну модификацию метода Гаусса - метод прогонки. Рассмотрим систему уравнений с трехдиа
Метод Якоби
Самый простой способ приведения системы к виду удобному для итерации состоит в следующем: из первого уравнения системы выразим неизвестное x1, из второго уравнения системы выраз
Метод Зейделя
Метод можно рассматривать как модификацию метода Якоби. Основная идея состоит в том, что при вычислении очередного (n+1)-го приближения к неизвестному xi при i >1 используют
Метод простой итерации.
Если A - симметричная и положительно определенная матрица, то систему уравнений часто приводят к эквивалентному виду: x = x -
Постановка задачи приближения функции по методу наименьших квадратов
Пусть функция y=f(x) задана таблицей своих значений: , i=0,1,-n. Требуется найти многочлен фиксированной степени m
Определение параметров эмпирической зависимости
Часто из физических соображений следует, что зависимость между величинами хорошо описывается моделью вида
Многочлены Бернштейна
Предположим, что функция задана в отрезке [0,1] в точках , при некотором фиксирован
Постановка задачи интерполяции функций
Пусть функция y = f(x) задана таблицей своих значений:
, i=0,1,...n. Требуется найти многочлен степени n
Оценка погрешности интерполяции
Если функция n+1 раз на отрезке [a,b] , содержащем узлы интерполяции , i=0,1,...n, то для погрешности интерполяц
Интерполяция сплайнами
Пусть отрезок [a,b] разбит точками на n частичных отрезков . Сплайном степени m называется функция
Первая производная. Двухточечные методы
Для двухточечных методов при вычислении производных используется значение функции в двух точках. Приращение аргумента задается тремя способами, откладывая Δx = h вправо, влево и в обе стороны
Вычисление производных второго порядка
Вторая производная вычисляется как первая производная от первой производной. Для следующей пятиточечной схемы (рис. 11.2):
Вычисление производных третьего порядка
Производные третьего порядка вычисляются как первая производная от производной второго порядка. Для рассмотренной пятиточечной схемы расчетная формула имеет вид:
Численное интегрирование
Определенным интегралом функции f(x), взятом в интервале от a до b, называется предел, к которому стремится интегральная сумма при стремл
Постановка задачи Коши для дифференциального уравнения первого порядка
Пусть дано дифференциальное уравнение первого порядка . Требуется найти функцию , у
Численное решение задачи Коши методом Эйлера
Численное решение задачи Коши состоит в построении таблицы приближенных значений в точках
Оценка погрешности метода Эйлера
Локальной погрешностью метода называется величина . Найдем величину локальной погрешности метода Эйлера:
Модификации метода Эйлера
Метод Эйлера обладает медленной сходимостью, поэтому чаще применяют методы более высокого порядка точности. Второй порядок точности по имеет
Решение систем дифференциальных уравнений методом Эйлера.
Пусть требуется решить нормальную систему дифференциальных уравнений:
с начальными условиями:
Численное интегрирование систем дифференциальных уравнений
Системой дифференциальных уравнений называется система вида
где x - независимый аргумент;
yi
Аппроксимация функции по Фурье
Пусть функция задана в интервале . В этом случае (при наличии у нее соответствующих
Преобразование Фурье
Так называется действие, с помощью которого по заданной в интервале функции строитс
Быстрое преобразование Фурье.
В предыдущем пункте было описано дискретное преобразование Фурье - сопоставление набору значений функции набора коэффициентов
Уравнение Лапласа (эллиптическое уравнение)
Это уравнение описывает некоторое установившееся стационарное состояние в пространстве (x,y) в некоторой области G (16.1):
Уравнение теплопроводности (параболическое уравнение)
Типичная задача, описывающая теплоперенос вдоль одномерного стержня единичной длинны при заданном начальном распределении температуры от
Новости и инфо для студентов