SIC^2: Simulation and Integration of Control for Canals
The LP method
This method does not calculate U control action variables to be applied to structures, but calculates the Lp-norm of a chosen Z measured variable.
For the LP method we read 2 parameters in the format (F8.2,1X,F8.2):
* P N
PS= 2. 1.
P: power coeff of the LP norm
N: flag for different options (0: **1/P, 1: /n, 2:*DTZ, otherwise no modification)
C-----Calculation of the power p, for the norm LP:
C-----|Z|lp = (Sum((Z)**p))**(1/p) (for p>0 integer)
C-----|Z|inf = Max(abs(Z)) (for p=0)
IF(IP.NE.0) THEN
PARA(3)=PARA(3)+ABS(Z(I))**IP
ELSE
PARA(3)=MAX(PARA(3),ABS(Z(I)))
ENDIFThe flag N is used at the end of the calculation:
IF(INT(PARA(2)).EQ.0) THEN
C------Classical norm LP
IF(IP.NE.0) PARA(3)=PARA(3)**(1./IP)
ELSEIF(INT(PARA(2)).EQ.1) THEN
C------We divide by the number of samples
IF(PARA(4).NE.0) PARA(3)=PARA(3)/PARA(4)
ELSEIF(INT(PARA(2)).EQ.2) THEN
C------Standard used for PID optimization for example
PARA(3)=PARA(3)*DTZ
ENDIF