Моделирование систем. Елизаров И.А - 70 стр.

UptoLike

ВУЗ: 

ТГТУ | Тамбов

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

Рубрика: 

68
0
0 .5
1
1 .5
2
2 .5
3
0
0 .1
0 .2
0 .3
0 .4
0 .5
0 .6
0 .7
0 .8
0 .9
1
С , к м о л ь /м
3
τ
, с
1
2
5
4
3
Рис. 10. Графические результаты моделирования реакции омыления
Программа на Matlab, решающая систему уравнений (5.25) представ-
лена ниже.
dt=0.001; % Шаг интегрирования по времени
time=[0:dt:3]; % Массив времени
n=length(time); % Размерность массива времени
% Создание массивов концентраций
Ca=zeros(1,n);
Cb=zeros(1,n);
Cc=zeros(1,n);
Cd=zeros(1,n);
Ce=zeros(1,n);
% Начальные условия
Ca(1)=1;
Cb(1)=0.5;
Cc(1)=0;
Cd(1)=0;
Ce(1)=0;
% Константы скоростей реакции
k1=2.24; k2=0.33;
% Интегрирование дифференциальных уравнений методом Эйлера
for i=2:n
Ca(i)=Ca(i-1)+dt*(-k1*Ca(i-1)*Cb(i-1));
Ce(i)=Ce(i-1)+dt*(k2*Cc(i-1)*Cb(i-1));
Cb(i)=Cb(1)+(Ca(i-1)-Ca(1))-(Ce(i-1)-Ce(1));
Cc(i)=Cc(1)-(Ca(i-1)-Ca(1))-(Ce(i-1)-Ce(1));

Страницы