Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- % Dane
- C0 = [2000; 0; 0; 0; 0; 0];
- % Czas
- tR1_span = [0 u(2)];
- % Sterowanie
- Temp = 298;
- tx = 75;
- S = 0.05;
- u = [Temp, tx, S];
- % Reaktor 1
- [T1, C1] = ode45(@(t, C)R1(t, C, u), tR1_span, C0);
- function dC = R1(t, C, u)
- k1 = 0.0444*exp(-2500/u(1));
- k2 = 6889*exp(-1*5000/u(1));
- dC = [ -2*k1*(C(1)*C(1));
- k1*(C(1)*C(1)) - k2*C(2);
- k2*C(2);
- 0;
- 0;
- 0; ];
- end
- % Czas
- tR2_span = [u(2) 180];
- % Sterowanie
- Temp = 298;
- tx = 75;
- S = 0.05;
- u = [Temp, tx, S];
- % Reaktor 2
- [T2, C2] = ode45(@(t, C)R2(t, C, u), tR2_span, C0_2);
- function dC = R2(t, C, u)
- dC = [ 0;
- -1*0.02*C(2)-0.05*C(2)-2*4*(10^(-5))*(C(2)*C(2));
- 0;
- 0.02*C(2);
- 0.05*C(2);
- 4*(10^(-5))*(C(2)*C(2)); ];
- end
- % Dane
- V1 = 0.1;
- % Sterowanie
- S = 0.05;
- u = [Temp, tx, S];
- C0B = 600;
- % Mieszanie
- V2 = V1 + u(3);
- C0_2 = M(C1(end,:), u, V1, C0B);
- function Cn = M(C, u, V1, C0B)
- Cn = [ (V1*C(1))/(V1+u(3));
- (V1*C(2) + u(3)*C0B)/(V1+u(3));
- (V1*C(3))/(V1+u(3));
- 0;
- 0;
- 0;
- ];
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement