Les systèmes bilinéaires jouent un rôle central dans la représentation des systèmes non-linéaires en tant qu'approximant des systèmes non-linéaires analytiques généraux. Une classe particulière de ces systèmes bilinéaires, dits homogènes, permet de représenter la fonctionnelle entrée-sortie comme une série de Volterra à horizon infini mais de degré fini. L'algorithme d'identification héréditaire jusqu'alors limité aux systèmes linéaires est ici étendu à l'identification de ces systèmes en utilisant une structure canonique de base pour l'approximant de degré deux (quadratique). Une forme NARX multidimensionnelle (Nonlinear AutoRegressive eXogenous input) est exploitée ici pour autoriser l'identification du système par la méthode héréditaire. On montre, sur un exemple académique quadratique, la supériorité de cette approche comparée aux algorithmes classiques.