Cet article traite de l’identification des systèmes par des modèles d’ordre fractionnaire multivariables. Une méthode à erreur de sortie, basée sur une technique d’optimisation non linéaire est développée. Elle repose sur un modèle d’état fractionnaire multivariable des fonctions de sensibilité paramétrique, dont la simulation est basée sur l’opérateur d’intégration fractionnaire borné en fréquences. Différentes simulations illustrent les performances de la méthode.