Example: Function call: kb_err(0.65,matrix_1). Parameter = 0.65 and matrix_1:
x1 |
x2 |
x3 |
x4 |
scenario_benchmark |
scenario_probability |
1 7 2 0 |
4 5 8 3 |
8 4 1 4 |
3 6 0 9 |
2 11 6 10 |
0.2 0.2 0.3 0.3 |
MATLAB code (./PSG/Example/Func_value_kb_err.m) for function minimization:
%Define data:
H1 = [1 4 8 3;7 -5 4 -6;-2 8 -1 0;0 -3 -4 9];
c1 = [-2;11;6;10];
p1 = [0.2;0.2;0.3;0.3];
%Create the PSG matrix 'matrix_1'(structure containing header and matrix body) and pack it to structure 'iargstruc_arr':
iargstruc_arr(1) = tbpsg_matrix_pack('matrix_1', H1, [], c1, p1);
%Create problem statement:
problem_statement = sprintf('%s\n',...
'minimize',...
' kb_err(0.65,matrix_1)',...
'');
%Minimize kb_err function with parameter 0.65 :
[solution_str, outargstruc_arr] = tbpsg_run(problem_statement, iargstruc_arr);
%Extract value of function from solution report:
val = tbpsg_function_data(solution_str, outargstruc_arr);
%Extract values of output points:
output = tbpsg_point_data(solution_str, outargstruc_arr);
%Display function value:
disp(' ')
disp(sprintf('kb_err function value = %g', val(1)));
disp(' ')
disp(sprintf('Pseudo R2 = %g', val(2)));
disp(' ')
disp('Contributions (normed increments) = ');
disp(output{2})
Program output:
kb_err function value = 3.86358e-015
Pseudo R2 = 1
Contributions (normed increments) =
0.3602 0.2864 0.2234 0.1300