I want to implement a filter on discrete system in Matlab. The filter is A(z)= z^3 + a1 z^2 + a2 z +a3, where as I need to apply it on y(k)=(1/A(z))*u(k). Here 'k' is discrete variable and 'z' is time shift operator. Here u(k) is input vector of discrete values i.e. a vector.