how to feed a signal through a transfer function

I'm sure this is probably really simple but I can't seem to find it. I have a z transform transfer function, given from tf, and I have an input signal. What i need is to pass the signal though the system and get the output. This is a small excerpt from my very long code that is the part I am stuck on.

if true
 % code
 %the signal
 N = 250;
 n = 0:1:N-1;
 x1 = 3*cos(20*pi*n/N);
 %the transfer function
 num = [0.3881 0.3881];
 den = [1 -0.4452 0.2700 -0.0486];
 fs = 1/250;
 hz = tf(num, den, fs)
end

Any help or direction would be greatly appreciated.

Answer 1

[Azzi Abdelmalek]

You can use lsim function

N = 250;
n = 0:1:N-1;
x1 = 3*cos(20*pi*n/N);
num = [0.3881 0.3881];
den = [1 -0.4452 0.2700 -0.0486];
fs = 1/250;
hz = tf(num, den, fs)
[y,t]=lsim(hz,x1)
stem(t,y)

[Gavin Davis]

Thank you

Answer 2

[Dominik]

I have a 2x2 transfer function, 'tf', and want to simulate this system with 2 inputs. every Inputis an array. How can I plo the response of my system? Do i can use lsim too? I get a lot of errors...

thanks a lot for your comments!