MATLAB: Please tell how

Question

num=[303 4848 0];

den=[1 16 164];

ea=tf(num,den)

subplot(1,2,2)

step(ea)

express it with

Answer 1

Try this

num=[303 4848 0];
den=[1 16 164];
ea=tf(num,den);
subplot(1,2,1)
step(ea)
[t, y] = ode45(@odeFun, [0 1], [303; 0]);
subplot(1,2,2)
plot(t, y(:,1), '-')
function dydt = odeFun(t, y)
    % transfer function is equivalent to following ODE
    % y'' = -16y'-164y+303u''+4848u'
    u = 1; % step input
    du = 0;
    ddu = 0;
    dydt = zeros(2, 1);
    dydt(1) = y(2);
    dydt(2) = -16*y(2)-164*y(1)+303*ddu+4848*du+134*u;
end