Расчет и оптимизация процессов и аппаратов химических и пищевых производств в среде MatLab. Дворецкий Д.С - 65 стр.

UptoLike

% вывод графика табличной функции маркерами
plot (x, y, 'ko')
% построение графиков полиномов
t = [50:5:110];
p1 = polyval(P1, t);
p2 = polyval(P2, t);
p5 = polyval(P5, t);
hold on
plot (t, p1, 'k:', t, p2, 'k-', t, p5, 'k-.')
legend('табличные данные', 'n=1', 'n=2', 'n=5', 0)
title('Приближение табличной функции полиномами степени n')
xlabel('t, град.С')
ylabel('P, кг/см2')
Результаты приближения методом наименьших квадратов представлены на рис. 8.2.
Очевидно, что приближение полиномом 5-го порядка наилучшее из рассмотренных, поэтому можно
воспользоваться полученным полиномом для расчета.
function P = vap(temp)
% функция зависимости давления паров воды
% от температуры (приближена полиномом 5-ой степени)
% использование P = vap(t),
% где P - давление, Па;
% temp - температура, град.С
p=[0.00000000007089 -0.00000001398190 0.00000284748046 ...
-0.00017754514603 0.00743024544080 -0.09255158370997];
P1 = polyval(p, temp);
P = 98100*P1;
Расчет вальцевой сушилки производится в программе dryer. Данная программа использует файл-
функцию vap, которая должна быть записана в ту же папку, что и программа dryer. Результаты расчёта
выводятся в командное окно.
Рис. 8.2 Результаты приближения методом наименьших квадратов
% расчёт вальцевой сушилки
% исходные данные