netgraph_movie.m

％The Metabolic Networks Toolbox contains functions to create, ％modify, display, and simulate bioche

% [M,ttt]=netgraph_movie(network,t,s_t,data_type,n_frames);
%
% s_t must be #variables x #timepoints !!!
%
% data_type {'concentrations','reactions','both','compute reactions'}
% default: concentrations
%
% usage: movie(M);

function [M,ttt] = netgraph_movie(network,t,s_t,data_type,n_frames);

if ~exist('data_type','var'), data_type='concentrations'; end 
if ~exist('n_frames','var'),  n_frames = 10; end

ttt = max(t)*(0:1/(n_frames-1):1);

s_t = interp1(t,s_t',ttt,'cubic')';

for j = 1:length(ttt),
  metvalues(:,j) = zeros(length(network.metabolites),1);
  actvalues(:,j) = zeros(length(network.actions),1);
  switch data_type,
    case 'concentrations',
      metvalues(:,j) = s_t(:,j);
    case 'reactions',
      actvalues(:,j) =  s_t(:,j);
    case 'both',
      metvalues(:,j) = s_t(1:length(network.metabolites),j);
      actvalues(:,j) = s_t(length(network.metabolites)+1:end,j);
    case 'compute reactions',
      metvalues(:,j) = s_t(:,j);
      actvalues(:,j) = network_velocities(metvalues(:,j),network,network.kinetics);    
  end
  end

jj=j;

for j=1:jj,
  netgraph_concentrations(network,metvalues(:,j),actvalues(:,j));
  noticks;
  M(j)=getframe;
end