Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- t=0:1:30; %s
- V0=1.5; %m3
- CA0=2.5; %mol/m3
- CB0=0; %mol/m3
- k=0.5
- function modelo6=modelo6(C) % C=[V CA CB CC CD]
- F1=3; %m3/s
- A=1; %m2
- F2=0.5*(C(1)/A); %m3/s
- CA1=1.5; %mol/m3
- CB1=1; %mol/m3
- modelo6(1)=F1-F2
- modelo6(2)=(F1*CA1-F2*C(2)-C(2)*(F1-F2)-k*C(2)*C(3)*C(1))/C(1)
- modelo6(3)=(F1*CA1-F2*C(3)-C(3)*(F1-F2)-k*C(2)*C(3)*C(1))/C(1)
- modelo6(4)=(-F2*C(4)-C(4)*(F1-F2)+k*C(2)*C(3)*C(1))/C(1)
- modelo6(5)=(-F2*C(5)-C(5)*(F1-F2)+k*C(2)*C(3)*C(1))/C(1)
- endfunction
- res=lsode("modelo5",[V0,CA0,CB0,CC,CD],t)
- figure(1)
- plot(t,res(:,2),"linewidth",1.5)
- hold on
- plot(t,res(:,3),"linewidth",1.5)
- plot(t,res(:,4),"linewidth",1.5)
- plot(t,res(:,5),"linewidth",1.5)
- xlabel("t, s","fontsize",12)
- ylabel("C_A, C_B, C_C, C_D, mol/m^3","fontsize",12)
- leyenda=legend("C_A","C_B","C_C","C_D");
- set(leyenda,"fontsize",12)
- figure(2)
- plot(t,res(:,1),"linewidth",1.5)
- xlabel("t, s","fontsize",12)
- ylabel("V, m^3","fontsize",12)
- leyenda2=legend("V");
- set(leyenda2,"fontsize",12)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement