Simulink : How to simulate response using the transfer function matrix with the time-domain input data?

Hello, I have a LTI system with 9 inputs and 2 outputs. I want to simulate the output of the system from the 'experimentally obtained' input data in Simulink. (Each of P1, P2 and P3 are in [3, time_length]. )

And the system's dynamic equation is given in state-space form, which is

                                  x_dot = Ax + Bu;

                                      y = Cx + Du;

Because of my experimental set-up, I cannot use 'State-Space' block to simulate the response. Instead, I utilized 'MATLAB Function' block and implement the part which solves the equation in it.

Since I have all the (A,B,C,D) matrices, I can fully formulate the (2 x 9) transfer function matrix from it. (using 'ss2tf' function.)

And there gonna be a clear expression which relates the output from the input, in the s-domain.

(i=1,2; r=1,2,...,9)

But I don't know how to solve this equation with the time-domain input data, in real-time sense. Looking forward to any help. Thanks. (I can't use 'lsim' function neither, because input data is given in time-by-time, which means the whole input history is not available, and output data should also be obtained time-by-time.)

[Birdman]

Send me your simulink model(export to R2017a) and all other necessary workspace variables in a zip file.

[jongmin]

Ok, I just attached the files. Thanks for your help.