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

UptoLike

ВУЗ: 

ТГТУ | Тамбов

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

Рубрика: 

% вызов солвера от файл-функции plugreact,
% начальной и конечной точке реактора
% и вектора начальных условий
[L, Y] = ode45('plugreact', [0 10], Y0);
% вывод графика решен. исходных дифф. уравнений
plot (L, Y(:, 1), 'ko')
hold on
plot (L, Y(:, 2), 'k>')
plot (L, Y(:, 3), 'k')
plot (L, Y(:, 4), 'k--')
hold off
% оформление графика
title('Изменение концентраций веществ по длине реактора') xlabel('Длина')
ylabel('Концентрации')
legend('A', 'B', 'D', 'E')
grid on
% проверка материального баланса
Ca = Y0(1)-Y0(2)/2+Y(10, 2)/2-Y(10, 4)*2
Y10 = Y(10, 1)
Результаты при вызове программы solvdem:
>> solvdem
Ca =
1.9975
Y10 =
1.9975
Визуализация решения приведена на рис. 6.1.
Рис. 6.1 Решение системы дифференциальных уравнений
6.2 Решение уравнений

Страницы