Advertisement
Guest User

Untitled

a guest
Apr 27th, 2017
72
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.97 KB | None | 0 0
  1. % Dane
  2. C0 = [2000; 0; 0; 0; 0; 0];
  3. % Czas
  4. tR1_span = [0 u(2)];
  5. % Sterowanie
  6. Temp = 298;
  7. tx = 75;
  8. S = 0.05;
  9. u = [Temp, tx, S];
  10. % Reaktor 1
  11. [T1, C1] = ode45(@(t, C)R1(t, C, u), tR1_span, C0);
  12. function dC = R1(t, C, u)
  13. k1 = 0.0444*exp(-2500/u(1));
  14. k2 = 6889*exp(-1*5000/u(1));
  15. dC = [ -2*k1*(C(1)*C(1));
  16. k1*(C(1)*C(1)) - k2*C(2);
  17. k2*C(2);
  18. 0;
  19. 0;
  20. 0; ];
  21. end
  22. % Czas
  23. tR2_span = [u(2) 180];
  24. % Sterowanie
  25. Temp = 298;
  26. tx = 75;
  27. S = 0.05;
  28. u = [Temp, tx, S];
  29. % Reaktor 2
  30. [T2, C2] = ode45(@(t, C)R2(t, C, u), tR2_span, C0_2);
  31. function dC = R2(t, C, u)
  32. dC = [ 0;
  33. -1*0.02*C(2)-0.05*C(2)-2*4*(10^(-5))*(C(2)*C(2));
  34. 0;
  35. 0.02*C(2);
  36. 0.05*C(2);
  37. 4*(10^(-5))*(C(2)*C(2)); ];
  38. end
  39.  
  40. % Dane
  41. V1 = 0.1;
  42. % Sterowanie
  43. S = 0.05;
  44. u = [Temp, tx, S];
  45. C0B = 600;
  46. % Mieszanie
  47. V2 = V1 + u(3);
  48. C0_2 = M(C1(end,:), u, V1, C0B);
  49. function Cn = M(C, u, V1, C0B)
  50. Cn = [ (V1*C(1))/(V1+u(3));
  51. (V1*C(2) + u(3)*C0B)/(V1+u(3));
  52. (V1*C(3))/(V1+u(3));
  53. 0;
  54. 0;
  55. 0;
  56. ];
  57. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement