This function solve MIMO-Diophantine equation.
This equation should be in this form
I = E(j)*A_tild + z^(-j)*F(j)
This equation is based on chapter 6 of Model Predictive Control Book by
E.F Comacho and C. Bordons.
[G,Gp,F,E] = mimodioph(A,B,N)
where A and B are the polynomials are calculated from system transfer function as following example.
please contact me for any question or suggestion.