Составители:
10
>> x=[170 180 190 200 210 220 230 240 250 260];
>> y=[2.382 2.321 2.265 2.211 2.161 2.113 2.067 2.024 1.982 1.943];
>> p=polyfit (x, y, 3)
После нажатия кнопки <Enter> получим ответ в следующем виде:
p =
-0.0000 0.0000 -0.0161 4.1757
Тогда функцией интерполяции будет следующий полином третьей
степени:
φ(x) = -0.0161 x + 4.1757,
или
φ(t) = -0.0161 t + 4.1757.
По этой формуле можно определять значение плотности пара в
зависимости от температуры.
В системе Matlab имеется функция вычисления математического
выражения при заданных значениях аргументов. Функция имеет вид:
polyval (p, x).
Здесь:
p ─ вычисляемая функция;
x ─ вектор аргументов функции.
Воспользуемся этой функцией для проверки достоверности результатов
интерполяции.
Введем функцию f=polyval (p,x) и нажмем кнопку <Enter>.
Откликом будет следующее решение:
f =
2.3817 2.3215 2.2648 2.2113 2.1607 2.1127 2.0672 2.0239
1.9825 1.9427
Сравнивая это решение с вектором у исходных данных, видим, что они
отличаются несущественно, а значит интерполяционный полином f(x) третьей
степени хорошо отображает исходную функцию.
Вычислим теперь абсолютную среднеквадратическую погрешность
аппроксимации по формуле:
2
1
.
n
i
i
n
ε
=
∆
=
∑
(1.1)
В нашем случае:
Страницы
- « первая
- ‹ предыдущая
- …
- 9
- 10
- 11
- 12
- 13
- …
- следующая ›
- последняя »
