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

UptoLike

ВУЗ: 

СПбГЛТУ | Санкт-Петербург

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

Рубрика: 

12
Заменим p = x, ρ
=
y.
>> x=[0.392 0.412 0.431 0.450 0.470 0.490 0.539 0.588 0.637 0.686];
>> y=[1.973 2.074 2.177 2.279 2.382 2.485 2.744 3.007 3.271 3.537];
>> diff(y, 1)
ans =
0.1010 0.1030 0.1020 0.1030 0.1030 0.2590 0.2630 0.2640
0.2660
>> dif f(y, 2)
ans =
0.0020 -0.0010 0.0010 -0.0000 0.1560 0.0040 0.0010 0.0020
>> diff(y, 3)
ans =
-0.0030 0.0020 -0.0010 0.1560 -0.1520 -0.0030 0.0010
Табличные разности третьего порядка близки между собой, поэтому
интерполяционный полином будет иметь вид:
φ(x) = ax
3
+ bx
2
+ cx + d,
или
φ(t) = at
3
+ bt
2
+ ct+ d.
Определим теперь значения коэффициентов a, b, c, d, воспользовавшись
функцией polyfit.
Так как векторы аргумента x и функции y уже введены, то достаточно
ввести функцию polyfit и нажать клавишу <Enter>:
>> p=polyfit(x,y,3)
p =
0.9531 -1.1096 5.6690 -0.1370

Страницы