Matematické Fórum

Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.

Nástěnka
22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.

Nejste přihlášen(a). Přihlásit

#1 26. 01. 2021 16:05

gag
Příspěvky: 44
Reputace:   
 

Matlab - kinetika

AHoj, prosím, kde je problém ve vykreslení 3 grafů pro katalyzátory C,CaCO3 a Al2O3. Zobrazuje se mi pouze jeden graf a nemůžu najít chybu. Označení k značí rychlostní konstanty pro katalyzátor C(první řádek, Caco3(2.řádek) a Al2O3(3.řádek).


function Pd
t_span=[0 1000];
c0=[1 1 1 0 0 0 0 0 0];
[t,c]=ode45(@kinetics,t_span,c0);
plot(t,c)
xlabel('Time,[min]');
ylabel('Concentration,[mol/cm^3]');
legend('cct-CDT','ctt-CDT','ttt-CDT','cc-CDD', 'ct-CDD', 'tt-CDD', 'c-CDE', 't-CDE', 'CDA')

end
function dz=kinetics(t,c);
k=[0.105 0.097 0.065 0.071 0.060 0.056 0.065 0.250 0.061 0.07 0.0110 0.120;%%%C
   0.015 0.007 0.010 0.008 0.025 0.016 0.005 0.003 0.007 0.015 0.001 0.035;%%CaCO3
   0.090 0.078 0.070 0.058 0.125 0.045 0.069 0.007 0.073 0.080 0.006 0.039];%%Al2O3
dz=zeros(9,1);
r1=k(1)*c(1);
r2=k(2)*c(1);
r3=k(3)*c(2);
r4=k(4)*c(2);
r5=k(5)*c(3);
r6=k(1)*c(1)-k(6)*c(4);
r7=k(2)*c(1)+k(3)*c(2)-k(7)*c(5)-k(8)*c(5)-k(9)*c(5);
r8=k(4)*c(2)+k(5)*c(3)-k(10)*c(6);
r9=k(6)*c(4)+k(7)*c(5)-k(11)*c(7);
r10=k(9)*c(5)+k(10)*c(6)-k(12)*c(8);
r11=k(11)*c(7)+k(8)*c(5)+k(12)*c(8);
dz(1)=-r1-r2;
dz(2)=-r3-r4;
dz(3)=-r5;
dz(4)=r6;
dz(5)=r7;
dz(6)=r8;
dz(7)=r9;
dz(8)=r10;
dz(9)=r11;
end

Offline

 

Zápatí

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson