For an autoregressive model, the current output is dependent on past outputs and the present input. For the given equation, the present output depends on past output as well as past input.
As a workaround, the VARMA model can be used to model dynamical systems as described in the question. The link here shows functionalities like arma2ar which might be helpful to realise such models.
Hope this helps!