Team:DTU-Denmark/plot mutant2
From 2013.igem.org
(Difference between revisions)
(Created page with "{{:Team:DTU-Denmark/Templates/StartPage|}} ==Timer series plots of Mutant2== <pre> code </pre> {{:Team:DTU-Denmark/Templates/EndPage}}") |
(→Timer series plots of Mutant2) |
||
Line 4: | Line 4: | ||
<pre> | <pre> | ||
- | + | clear all | |
+ | % create model | ||
+ | model=SBmodel('mutant_2.txt') | ||
+ | SBPDmakeMEXmodel(model,'modelmex'); | ||
+ | |||
+ | %calculate steady state | ||
+ | SBsteadystate(model) | ||
+ | |||
+ | conc = [5.872e-2 5.872e-1 5.872 29.358]; | ||
+ | enz = [2.115e-4 1.0575e-03 2.1115e-3 2.115e-1]; | ||
+ | count = 0; | ||
+ | for i=1:4 | ||
+ | count = 0; | ||
+ | for j=1:4 | ||
+ | count = count+1; | ||
+ | % make parametervector | ||
+ | parametervector = makeparamvecSBPD('modelmex',{'NO20','E'},[conc(i); enz(j)]); | ||
+ | output=modelmex([0:0.001:10],[],parametervector); | ||
+ | % save results in separate vectors | ||
+ | % check in which order they are (double click on output in workspace) | ||
+ | NO(:,count) = output.statevalues(:,1); | ||
+ | N2O(:,count) = output.statevalues(:,2); | ||
+ | end | ||
+ | t = output.time'; | ||
+ | % plot concentrations over time: | ||
+ | figure | ||
+ | plot(t,N2O(:,1),t,N2O(:,2),'--',t,N2O(:,3),':',t,N2O(:,4),'-.','LineWidth',2) | ||
+ | if i == 1 | ||
+ | title('1 mg/l nitrite (58.72 uM)') | ||
+ | elseif i == 2 | ||
+ | title('10 mg/l nitrite (587.2 uM)') | ||
+ | elseif i == 3 | ||
+ | title('100 mg/l nitrite (5.872 mM)') | ||
+ | elseif i == 4 | ||
+ | title('500 mg/l nitrite (29.358 mM)') | ||
+ | end | ||
+ | |||
+ | xlabel('Time in [min]','FontSize',16) | ||
+ | ylabel('N2O concentration in [mM]','FontSize',16) | ||
+ | hleg1 = legend('100 enzymes per cell','500 enzymes per cell','1000 enzymes per cell','100 000 enzymes per cell'); | ||
+ | set(hleg1,'Location','SouthEast','FontSize',12) | ||
+ | end | ||
</pre> | </pre> | ||
{{:Team:DTU-Denmark/Templates/EndPage}} | {{:Team:DTU-Denmark/Templates/EndPage}} |
Latest revision as of 18:44, 4 October 2013
Timer series plots of Mutant2
clear all % create model model=SBmodel('mutant_2.txt') SBPDmakeMEXmodel(model,'modelmex'); %calculate steady state SBsteadystate(model) conc = [5.872e-2 5.872e-1 5.872 29.358]; enz = [2.115e-4 1.0575e-03 2.1115e-3 2.115e-1]; count = 0; for i=1:4 count = 0; for j=1:4 count = count+1; % make parametervector parametervector = makeparamvecSBPD('modelmex',{'NO20','E'},[conc(i); enz(j)]); output=modelmex([0:0.001:10],[],parametervector); % save results in separate vectors % check in which order they are (double click on output in workspace) NO(:,count) = output.statevalues(:,1); N2O(:,count) = output.statevalues(:,2); end t = output.time'; % plot concentrations over time: figure plot(t,N2O(:,1),t,N2O(:,2),'--',t,N2O(:,3),':',t,N2O(:,4),'-.','LineWidth',2) if i == 1 title('1 mg/l nitrite (58.72 uM)') elseif i == 2 title('10 mg/l nitrite (587.2 uM)') elseif i == 3 title('100 mg/l nitrite (5.872 mM)') elseif i == 4 title('500 mg/l nitrite (29.358 mM)') end xlabel('Time in [min]','FontSize',16) ylabel('N2O concentration in [mM]','FontSize',16) hleg1 = legend('100 enzymes per cell','500 enzymes per cell','1000 enzymes per cell','100 000 enzymes per cell'); set(hleg1,'Location','SouthEast','FontSize',12) end