Интерполяция ии аппроксимация функций Таблично-заданная функция:

Интерполяция ии аппроксимация функций

Таблично-заданная функция: Достоинство: Нет необходимости проводить дополнительные расчеты значения функции. Недостаток: Невозможность определения значений функции при значениях аргумента, отличных от табличных. XX YY xx 11 yy 11 xx 22 yy 22. . . xx ii yy ii …… …… xx nn yy nn

Графическая интерпретация методов: Интерполяция Аппроксимация xy y=f(x) x 1 x iy 1 y i Узлыy*=f(x)

1. Интерполяция функций А) А) Глобальная интерполяция. . Б) Б) Локальная интерполяция. .

Локальная интерполяция Необходимо определиться с видом интерполирующей функции. На практике чаще всего использую следующие виды локальной интерполяции: 1. 1. Линейная интерполяция 2. 2. Квадратичная интерполяция 3. 3. Сплайн-интерполяция

Глобальная интерполяция Необходимо, чтобы одна интерполирующая функция проходила через все узлы таблично заданной функции. Обычно в этом случае интерполирующую функцию ищут в виде полинома: y=cy=c 00 +c+c 11 xx ++ cc 22 xx 2 2 +…++…+ cc ii xxii ++ …+c…+c n-1 n-1 xxn-1 n-1 Через nn точек может проходить только один график функции полинома n-1 n-1 порядка. .

Интерполяция методом канонических полиномов В уравнение полиномиальной зависимости подставляются все значения точек xx — yy (( аргумент — функция ) ) табличной зависимости. Таким образом, получается система линейных уравнений nn -го порядка: nn nni ninn in ini iiii n ni i yxcxcxcxcc

Решаем полученную СЛАУ одним из известных методов, получаем значения параметров полиномиальной зависимости: cc 00 , , cc 11 , , cc 22 , …, ccii , …, ccn-1 n-1 В полученное уравнение полиномиальной функции с найденными параметрами подставляем значение аргумента xx цц и и находим значение функции yy цц. . Данный метод удобен в случаях, когда необходимо найти несколько значений функции при нескольких значениях аргумента.

Интерполяция методом полинома Лагранжа Данный метод удобен в случаях, когда необходимо найти только одно значение функции при одном значении аргумента. В этом методе полином n-1 степени представляется в виде: где ll 11 , , ll 22 , , …, …, llii , , …, llnn — — являются также полиномиальными зависимостями. . nniiцlylyy. . .

Данные полиномы рассчитываются по следующим формулам: Коэффициент llii будет равен 1, если xxцц == xxii , и будет равен 0 (нулю), если xx цц == xx 11 , , xx 22 и так далее. ni nц ii iц iц iц i n nц i iццц xx xx xx l xx xx xx xx l . . . .

Блок-схема интерполяции полиномом Лагранжа НАЧАЛО ввод n i=1, n ввод x(i), y(i) ввод xx i=1, n L=L*(xx-x(j))/(x(i)-x(j))Количество пар точек yy=yy+y(i)*L вывод yy КОНЕЦ дане т. Ввод табличной функции Ввод аргумента, при котором необходимо рассчитать значение функции yy=0 L=1 j=1, n i≠j

Достоинства метода интерполяции: Если исходная таблично заданная функция не содержит в себе погрешности, то метод интерполяции дает достаточно точные и достоверные результаты. Недостатки метода интерполяции: 1. В случае локальной интерполяции получается несколько интерполирующих функций, с которыми неудобно работать. 2. В случае глобальной интерполяции при достаточно большом количестве пар точек (узлов) получается достаточно сложная и громоздкая интерполирующая функция. 3. Если исходная таблично заданная функция содержит в себе погрешность, то интерполирующая функция также будет «копировать» эту погрешность и не даст необходимого достоверного результата.

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

Пример: Рассматривается линейная функция. Метод интерполяции в этом случае неприменим. Применяется метод аппроксимации. xy

Возникает вопрос: Каков критерий близости аппроксимирующей функции к исходной таблично заданной? МЕТОД НАИМЕНЬШИХ КВАДРАТОВ (( МНКМНК )). . Суть: Аппроксимирующая функция должна быть такой, чтобы сумма квадратов отклонений аппроксимирующей функции от исходной таблично заданной была минимальной.

Согласно МНК только одна функция будет лежать максимально близко к исходным узлам. xy xy 1 3 4 5 2 7 6 y=y*(x)min))(*( 1 2 n i ii n i ixyy. S

Замечание: МЕТОД НАИМЕНЬШИХ КВАДРАТОВ не не позволяет определить видвид аппроксимирующей функции. МНКМНК позволяет определить только параметры выбранного вида аппроксимирующей функции. Вид аппроксимирующей функции определяет пользователь исходя из теоретических знаний, логических предпосылок или других разумных соображений.

Аппроксимация линейной функцией В качестве аппроксимирующей выбирают линейную функцию: yy ** =c=c 00 +c+c 11 xx. . Здесь cc 00 ии cc 11 — параметры функции, которые определим с использованием МНК.

В формулу суммы квадратов отклонений подставим аппроксимирующую линейную функцию, получим: Находим минимум функции SS , приравняв частные производные функции S S по по параметрам аппроксимирующей функции нулю: min)())(( 1 2 10 1 2* 1 2 n i iin i ixccyxyy. S 0)(2 1 10 0 1 10 1 n i iii xccy c S xxccy c S

Преобразуем полученные выражения: В результате получаем СЛАУ, где в качестве неизвестных выступают параметры линейной функции cc 00 ии cc 11 : : 0)( 1 01 1 0 2 1 n i iiii cxcy xcxcyx n i i in i i yncxc yxxcxc

Решаем полученную СЛАУ методом Крамера: Таким образом, вычислив параметры cc 00 ии cc 11 , мы получаем линейную аппроксимирующую функцию: yy ** =c=c 00 +c+c 11 xx. . 2 1111 2 0 2 111 1 n i iin i i n i in i i xnx xyxyx c xnx yxnyx c

Блок-схема аппроксимации линейной функцией НАЧАЛО ввод n i=1, n ввод x(i), y(i) ввод xx i= 1 , n c 1 =(S 3 *n-S 2 *S 4 )/(S 1 *n-S 2 2 )) c 0 =(S 1 *S 4 -S 2 *S 3 )/(S 1 *n-S 2 2 ))Количество пар точек yy=c 0 +c 1 *xx вывод c 1 , c 0 , yy КОНЕЦВвод табличной функции Ввод аргумента, при котором необходимо рассчитать значение функции S 1 =0; S 2 =0 S 3 =0; S 4 =0 S 1 =S 1 +x i *x i S 2 =S 2 +x i S 3 =S 3 +x i *y i S 4 =S 4 +y i

Аппроксимация квадратичной функцией В качестве аппроксимирующей выбирают квадратичную функцию: yy ** =c=c 00 +c+c 11 xx ++ cc 22 xx 22. . Здесь cc 00 , , cc 11 ии cc 22 — параметры функции, которые определим с использованием МНК.

В формулу суммы квадратов отклонений подставим аппроксимирующую квадратичную функцию, получим: Для нахождения минимума функции SS приравняем частные производные функции S S по параметрам аппроксимирующей функции нулю: min)())(( 1 22 210 1 2* 1 2 n i iii n i ixcxccyxyy. S 0)(2 1 2 210 0 1 2 210 1 1 22 210 2 n i iiii xcxccy c S xxcxccy c S

Преобразуем полученные выражения следующим образом: В результате получаем систему линейных алгебраических уравнений, где в качестве неизвестных выступают параметры квадратичной функции cc 00 , , cc 11 ии cc 22 : : 0)( 0)( 1 01 2 2 1 0 2 1 3 2 1 2 0 3 1 4 2 2 n i iiiii cxcxcy xcxcxcyx n i in i in i in i i yncxcxc yxxcxcxc

Решаем полученную СЛАУ методом Крамера для систем 3 -го порядка и вычисляем параметры cc 00 , , cc 11 ии cc 22. . Таким образом, мы получаем квадратичную аппроксимирующую функцию: yy ** =c=c 00 +c+c 11 xx ++ cc 22 xx

Блок-схема аппроксимации квадратичной функцией НАЧАЛО ввод n i=1, n ввод x(i), y(i) ввод xx i= 1 , n a 11 =S 1 ; a 12 =S 2 ; a 13 =S 3 a 21 =S 2 ; a 22 =S 3 ; a 23 =S 4 a 31 =S 3 ; a 32 =S 4 ; a 33 =n b 1 =S 5 ; b 2 =S 6 ; b 3 =S 7 Количество пар точек yy=c 0 +c 1 *xx+с 2 * xx 2 вывод c 2 , c 1 , c 0 , yy КОНЕЦВвод табличной функции Ввод аргумента, при котором необходимо рассчитать значение функции S 1 =0; S 2 =0 S 3 =0; S 4 =0 S 5 =0; S 6 =0 S 7 =0 S 1 =S 1 +x i 4 S 2 =S 2 +x i 3 S 3 =S 3 +x i 2 S 4 =S 4 +x i S 5 =S 5 +x i 2 *y i S 6 =S 6 +x i *y i S 7 =S 7 +y i c 2 = вычисляются c 1 = по формулам c 0 = Крамера