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.

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