La méthode LP

Cette méthode ne calcule pas des commandes U à appliquer à des ouvrages, mais calcule la norme Lp d’une mesure Z choisie.

editregv5_lp

Pour la méthode LP on lit 2 paramètres au format (F8.2,1X,F8.2): * P N PS= 2. 1.

P: coeff puissance de la norme LP N: flag pour différentes options (0: **1/P, 1: /n, 2:*DTZ, sinon idem)

C-----Calcul de l'exposant P, pour la norme LP :
C-----|Z|lp  = (Somme((Z)**p))**(1/p) (pour p entier positif)
C-----|Z|inf = Max(abs(Z)) (pour p=0)

IF(IP.NE.0) THEN
     PARA(3)=PARA(3)+ABS(Z(I))**IP
ELSE
     PARA(3)=MAX(PARA(3),ABS(Z(I)))
ENDIF

Le flag N est utilisé en fin de calcul:

IF(INT(PARA(2)).EQ.0) THEN
C------Norme classique LP
     IF(IP.NE.0) PARA(3)=PARA(3)**(1./IP)
ELSEIF(INT(PARA(2)).EQ.1) THEN
C------On divise par le nombre d'échantillons
     IF(PARA(4).NE.0) PARA(3)=PARA(3)/PARA(4)
ELSEIF(INT(PARA(2)).EQ.2) THEN
C------Norme utilisée pour l'optimisation des PID par exemple
     PARA(3)=PARA(3)*DTZ
ENDIF