Example: Calculation of L1 Distance between points x = (3,4) and y = (2,3).

 

Function call: polynom_abs(matrix_1).

 

matrix_1:

 

 

at point_1:

 

 

Note. For this example y1=2, y2=3 are set in matrix_1.

 

MATLAB code (.\Aorda\PSG\MATLAB\Examples\Func_value_L1_Distance.m) for function calculation:

 

%Define data:

H1 = [1 1;2 3];

a = [3;4];

 

%=======================================================================

%Calculation of function with PSG Subroutine 'functionvalue'

%Calculate polynom_abs function at point 'a':

val1 = functionvalue('polynom_abs', [], H1, [], [], a);

%Display function value:

disp(sprintf('functionvalue = %g', val1));

%=======================================================================

%Calculation of function with PSG Subroutine 'tbpsg_function_value'

%Create the PSG matrix 'matrix_1'(structure containing header and matrix body) and pack it to structure 'toolboxstruc_arr':

toolboxstruc_arr(1) = tbpsg_matrix_pack('matrix_1', H1, [],  [], []);

%Create the PSG point 'point_1'(structure containing header and matrix body) and pack it to structure 'toolboxstruc_arr':

toolboxstruc_arr(2) = tbpsg_point_pack('point_1', a, []);

%Calculate polynom_abs function at point 'point_1':

val2 = tbpsg_function_value('polynom_abs(matrix_1)','point_1',toolboxstruc_arr);

%Display function value:

disp(sprintf('tbpsg_function_value = %g', val2));

%=======================================================================

%Creation of User Subroutine for calculating function value

%User subroutine is similar to the 'functionvalue'

%Create problem statement:

problem_statement = sprintf('%s\n',...

'calculate',...

'Point: point_1',...

'  polynom_abs(matrix_1)',...

' ');

%Create Subroutine:

tbpsg_create_user_subroutine(problem_statement,toolboxstruc_arr,'.\','calculate_polynom_abs');

%Uncomment section to call the user subroutine 'calculate_polynom_abs':

%{

%Calculation of function with User Subroutine 'calculate_polynom_abs':

[solution_str,outargstruc_arr] = calculate_polynom_abs(H1,matrix_1_vars,a,point_1_vars);

%Extract function value from the solution report:

val = tbpsg_function_data(solution_str, outargstruc_arr);

%}

 

 

Program output:

 

functionvalue = 2

tbpsg_function_value = 2