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
! 2.11.2020 (L) Vykreslete si svůj první matematický výraz přes MathJax!
! 04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
17.01.2016 (Jel.) Rok 2016 s novými a novějšími krystaly od kolegy Pavla!
17.01.2016 (Jel.) Nabídka knih z oborů matematiky, fyziky, chemie
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: 39
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