用matlab 某系统的传递函数为 G(s)=5s+100/s+2s+16s+18s+18,输入为u=1+e-tcos(5t)，当系统的初始条件为零时，求出其时域响应。

在MATLAB中，计算系统对给定输入的时域响应通常需要将传递函数转换为微分方程形式，然后解这个微分方程。首先，我们需要将传递函数G(s)转换成状态空间模型（State-Space Representation），因为MATLAB的`lsim`或`ode45`函数需要这种形式的模型。

传递函数 G(s) = (5s + 100) / (s^2 + 2s + 16s + 18s + 18) 可以简化一下：

G(s) = (5s + 100) / (s*(s + 14)(s + 9))

接下来，我们需要找到该系统的零点、极点以及增益矩阵。然后我们可以设置初始条件x0 = [0; 0]，输入信号u(t) = 1 + e^-t * cos(5t)，并使用`ss`函数将其转换为状态空间模型 (`sys`)。

假设我们完成了上述步骤，然后可以使用下面的MATLAB代码来计算时域响应y(t)：

% 系统状态空间模型
sys = ss(A,B,C,D); % A, B, C, D 分别对应系统矩阵

% 输入信号
tspan = [0, T]; % 时间范围
u = @(t) 1 + exp(-t).*cos(5*t); % 输入函数

% 初始条件
x0 = zeros(numStates, 1);

% 解微分方程得到时域响应
[t, y] = ode45(@(t,x) ODEfun(t,x,A,B,C,D,u(t)), tspan, x0);

% ODEfun是一个内部函数，用于定义状态方程
function dydt = ODEfun(t, x, A, B, C, D, u)
    dxdt = A*x + B*u;
    dydt = C*x + D*u;
end

% y(t) 就是时域响应

请注意，由于没有给出具体的A, B, C, D矩阵，这部分需要根据实际的传递函数进行计算。运行以上代码后，变量`y`将包含系统在给定初始条件下的时域响应。