gas1388
عضو جدید
سلام
میتونید توی اررور های این برنامه کمکم کنید؟
function f=ex6_18(c,tau)
k1=0.12;k2=0.046;k3=0.02;k4=0.034;k5=0.04;cc0=2;
v0=10;v=300;
tau=v./v0;
eq1='c(1)+(k1+k2)*tau.*c(1)=cc0';
eq2='c(2)-(k1.*c(1)-k3.*c(2)+k4*tau.*c(3))=0';
eq3='c(3)-(k2.*c(1)+k3.*c(2)-k4.*c(3)-k5.*c(3)).*tau=0';
eq4='c(5)-k5*tau.*c(3)=0';
[c(1);c(2);c(3);c(4)]=solve(eq1,eq2,eq3,eq4,'c(1),c(2),c(3),c(4)')
for i=1:4
f(i)=c(i).*v0;
end
plot(tau,f,1),'+r',tau,f,2),'*g',tau,f,3),'-m',tau,f,4),'.b')
title('F-Tau(CSTR) \n [N.CHAVOSHY & A.PARSE]')
legend ('Fcoal','Fpreasfalin','Fasfalin','Foil')
xlabel ('tau(min)')
میتونید توی اررور های این برنامه کمکم کنید؟
function f=ex6_18(c,tau)
k1=0.12;k2=0.046;k3=0.02;k4=0.034;k5=0.04;cc0=2;
v0=10;v=300;
tau=v./v0;
eq1='c(1)+(k1+k2)*tau.*c(1)=cc0';
eq2='c(2)-(k1.*c(1)-k3.*c(2)+k4*tau.*c(3))=0';
eq3='c(3)-(k2.*c(1)+k3.*c(2)-k4.*c(3)-k5.*c(3)).*tau=0';
eq4='c(5)-k5*tau.*c(3)=0';
[c(1);c(2);c(3);c(4)]=solve(eq1,eq2,eq3,eq4,'c(1),c(2),c(3),c(4)')
for i=1:4
f(i)=c(i).*v0;
end
plot(tau,f,1),'+r',tau,f,2),'*g',tau,f,3),'-m',tau,f,4),'.b')
title('F-Tau(CSTR) \n [N.CHAVOSHY & A.PARSE]')
legend ('Fcoal','Fpreasfalin','Fasfalin','Foil')
xlabel ('tau(min)')