through non-classical damping system's modal summation analysis, i get
Zi'(t) - Di*Zi(t) = (Si)T * F0 * r(t)
Zi'(t) is first differential of function Zi(t), (Si)T is transpose matrix of Si n=7, i=1,2,....,2n Zi(t) is the function i want to get, Di is 14x14 diagonal matrix, Si is 14x14 matrix , F0 is 14x1 matrix and r(t) is a known function Di and Si are complex matrixes. In this case, how can i get Zi(t) through MATLAB? Let me know if you need more informations. Sorry for my poor english, Thank you.