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');