Un système adaptatif de commande de processus industriels de type H comprend une unité de commande (21) qui effectue la régulation d'un processus industriel (20) (navire, laminoir, etc.), dont l'état doit être mesuré (y). Le signaux de sortie de l'unité de commande et du processus sont acheminés dans un estimateur paramètrique (34) qui estime les polynomes opérationels (fonctions de transfert) ou fonctions opératrices de retard) A, B et C du processus. Un organe (35) de conception de l'unité de commande traite les signaux de sortie de l'estimateur (34) afin de produire des coefficients pour l'unité de commande (21), qui établit un raport polynomial Cn/Cd. L'organe (35) utilise un algorithme de type H en ligne se fondant sur la réduction au minimum d'une forme "carrée de somme" de la fonction (Pc.(t)+Fc.e(t))2, ce qui a pour résultat des écuations qui sont plus faciles à résoudre que les équations provenant de la forme habituelle "somme de carrés" de la fonction. Un polynome L=PcnFcdB-FcnPcdA est calculé en (72), les zéros instables L- sont calculés en (73), les équations matricielles FAPcdlambda+LG=PcnCF*zn; FBFcdlambda-LH=FcnCF*zn sont calculées en (74), la valeur propre extrême est extraite en (75), le vecteur propre correspondant est extrait en (76) et Cn et Cd sont calculés sous la forme GFcd et HPcd en (77) pour être transmis à l'unité de commande (21). A la place de cette technique par vecteur propre/valeur propre, on peut utiliser une technique modifiée (de type F-itération), dans laquelle les équations FAPcdlambda+LG=PcnC; FBFcdlambda-LH=FcnC sont résolues, F*, qui est l'adjoint de F, est calculé, les polynomes se trouvant du côté droit de ces équations sont multipliés par cet adjoint et la procédure est soumise à une itération. On peut faire varier les paramètres de Fc et Pc.