Using ode45 to solve matrix/vector form of state space differential equation

I currently have two state space equations in matrix/vector form

The first is Ydot4=[a 0 0 0;0 0 0 0; 0 0 b 0;0 0 0 0 c]*[y1;y2;y3;y4]

The second is Ydot2= [d 0 0 0;0 e 0 0;0 0 f 0;0 0 0 g]*[y1;y2;y3;y4]+[h 0 0 0]*V0(t)

I would like to know how to solve these equations using ODE45 over a timespan of 0 to 1 seconds with intervals of 1ms. And V0(t)=5cos(40t)

Ydot4=V2doubledot

Y4=V2dot

Y3=V2

Y2=V1dot

Ydot2=V1doubledot

Y1=V1