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

UptoLike

ВУЗ: 

ТГТУ | Тамбов

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

Рубрика: 

% программа ammiak расчета летучести аммиака
% в зависимости от давления
% с использованием оператора for
R = 8.31514; % универсальная газовая постоянная
T = 473; % температура, К
% табличные данные
P = [0.001 1 10 20 30 40 50 60 80 100 200 300 ...
400 500 600 800];
Z = [0 0.998 0.981 0.961 0.942 0.922 0.902 0.882...
0.841 0.801 0.551 0.462 0.495 0.557 0.621 0.755];
% Расчет
for i=1:length(P)
I = quad('ammiak1', 1.0e-10, P(i), [], [], Z(i));
f(i) = P(i)*exp(I); % летучесть аммиака
v(i) = (Z(i)*R*T)/P(i); % объем газа
end
% графический вывод результатов
subplot(3,1,1)
plot ( P, Z, 'k-')
xlabel('P')
ylabel('Z')
grid on
subplot(3,1,2)
plot ( P, v, 'k-')
xlabel('P')
ylabel('v')
grid on
subplot(3,1,3)
plot ( P, f, 'k-')
xlabel('P')
ylabel('f')
grid on
Результаты расчета на отдельных графиках представлены на рис. 6.3.
Реализованный алгоритм основан на квадратурной формуле Симпсона с автоматическим подбором
шага интегрирования для достижения требуемой относительной погрешности. Интегрирование гладких
функций лучше производить при помощи quad8, основанной на наиболее точных квадратурных форму-
лах Ньютона-Котеса с автоматическим подбором шага.
При вычислении в MatLab 5.2, 5.3 интегралов от функций, имеющих интегрируемую особенность,
например x1 , в нуле выводится сообщение об ошибке. В версии 6.0 данный алгоритм улучшен для
функций с особенностями, и, кроме того, появился quadl для нахождения определенного интеграла по
квадратурным формулам Гаусса-Лобатто.

Страницы