www.gusucode.com > simbio 案例源码程序 matlab代码 > simbio/CalculateSensitivitiesExample.m

    %% Calculate Sensitivities
% 
%% Load the Example Model
% 

% Copyright 2015 The MathWorks, Inc.

sbioloadproject gprotein_norules m1
%% Get the configuration set
%
csObj = getconfigset(m1);
%% Select Ga species
% 
Ga = sbioselect(m1,'Type','species','Where','Name','==','Ga');

%% Set the OutPuts property
% 
csObj.SensitivityAnalysisOptions.Outputs = Ga;
%% Select parameters
%
pif = sbioselect(m1,'Type','parameter');
%% Set Inputs property
% 
csObj.SensitivityAnalysisOptions.Inputs =  pif;
%% Turn on SensitivityAnalysis
% 
csObj.SolverOptions.SensitivityAnalysis = true;
%% Set Normalization
%
csObj.SensitivityAnalysisOptions.Normalization = 'Full';
%% Perform sensitivity analysis via simulation
%
simDataObj = sbiosimulate(m1);
%% Extract sensitivity data
%
[T, R, snames, ifacs] = getsensmatrix(simDataObj);
%% Reshape R
% 
R2 = squeeze(R);
%% Plot data
% 
figure;
plot(T,R2);
title('Normalized Sensitivity of Ga With Respect To Various Parameters');
xlabel('Time (seconds)');
ylabel('Normalized Sensitivity of Ga');
leg = legend(ifacs, 'Location', 'NorthEastOutside');
set(leg, 'Interpreter', 'none');