Guest User

pha_contd

a guest
Sep 9th, 2019
39
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. y0=[8.73317E-06 0.00059068 1.11715E-06 6.70578E-05 7.24048E-05 2.77752E-08 7.65432E-06 8.77345E-05];
  2. [t,y]=ode15s(@pha_cont,[0,15],y0);
  3.  
  4.  
  5. %plot glucose and biomass
  6. % figure;
  7. % subplot(1,2,1)
  8. % plot(t,[y(:,1),y(:,5)]);
  9. % legend('R1 Glucose','R2 Glucose');
  10. % xlabel('Time (hrs)');
  11. % ylabel('Concentration (g/L)');
  12.  
  13. %plot nh3
  14. % subplot(1,2,2)
  15. % plot(t,[y(:,2),y(:,6)]);
  16. % legend('R1 NH3','R2 NH3');
  17. % xlabel('Time (hrs)');
  18. % ylabel('Concentration (g/L)');
  19.  
  20. %plot Biomass and PHA
  21. figure
  22. plot(t,[y(:,4),y(:,3),y(:,8),y(:,7)]);
  23. legend('R1 PHA','R1 Biomass','R2 PHA','R2 Biomass');
  24. xlabel('Time (hrs)');
  25. ylabel('Concentration (g/L)');
  26.  
  27. %plot PHA fraction
  28. % figure
  29. % frac=[y(:,4)./(y(:,3)+y(:,4)),y(:,8)./(y(:,7)+y(:,8))];
  30. % plot(t,frac)
  31. % legend('R1','R2');
  32. % xlabel('Time (hrs)');
  33. % ylabel('PHA DCW Fraction');
  34.  
  35. % Ss_in = 206;
  36. % Snh_in_1 = 3;
  37. % r1_glu_end = y(end,1)./Ss_in
  38. % r1_nh3_end = y(end,2)./Snh_in_1
  39. % r2_glu_end = y(end,5)./Ss_in
RAW Paste Data