Автоматизированные системы управления технологическими процессами. Втюрин В.А - 12 стр.

UptoLike

Составители: 

11
i
= y
i
f(x) = y
i
- f
i;
n = 10
Тогда вычислительные процедуры в системе Matlab будут иметь вид:
z = у - f;
w = z . * z ;
Здесь точка перед знаком умножения означает что операция умножения
выполняется поэлементно, т.е. с каждым элементом массива чисел (вектора
строки).
R = sum (w);
Е = sqrt(R)/10
После нажатия клавиши <Enter> получим ответ:
E =
1.0186e-004.
Полученная ошибка очень мала, что подтверждает хорошее совпадение.
Система Matlab позволяет обоснованно выбрать степень полинома при
полиномиальной интерполяции путем вычисления табличных разностей. Для
этой цели служит функция diff. Эта функция имеет вид:
diff (v, n).
Здесь:
v — вектор функции у(х);
n — порядок конечных разностей.
Полиномиальная аппроксимация имеет смысл лишь тогда, когда n -ая ко-
нечная разность функции у(х) при постоянном шаге изменения аргумента х
является постоянной. При этом значение n является степенью полинома. Если
это условие не выполняется, то многочлен степени п либо не может быть
функцией интерполяции, либо является основным источником погрешностей.
Найдем полиномиальную зависимость ρ (плотность перегретого пара) от
давления p при постоянной температуре t. Результаты приведены в табл. 1. 2.
Таблица 1. 2.
Зависимость плотности перегретого пара от давления при постоянной температуре
Значения переменных
p, МПа 0.392 0.412 0.431 0.450 0.470 0.490 0.539
0.588
0.637
0.686
ρ, кг/м
3
1.973 2.074 2.177 2.279 2.382 2.485 2.744
3.007
3.271
3.537