怎样用matlab表示如下方程

justlikethisis

ny=nu=2;L=2;其中ny，nu分辨表示系统的输出和输入的最大时滞；L为系统的自由度；
请问如何用matlab表示下列多项式：
y(k)=y(k-1)+y(k-2)+u(k-1)+u(k-2)+y^2(k-1)+y^2(k-2)+u^2(k-1)+u^2(k-2)+y(k-1)u(k-1)+y(k-1)u(k-2)+y(k-2)u(k-1)+y(k-2)u(k-2)+y(k-1)y(k-2)+u(k-1)u(k-2);

geronimo123

这就是时间序列里的ARMA模型啊