Question: Solve ode with time and input vector from CSV file

I solve ode for simply DC motor system:

J:= 0.01;b:=0.1;K:= 0.01;R:=1;L:=1;
eq1:=J*diff(theta(t),t,t)+b*diff(theta(t),t)=K*i(t):
eq2:=L*diff(i(t),t)+R*i(t)=V-K*diff(theta(t),t):

ICs:= theta(0)=0, D(theta)(0) = 0, D(theta)(0) = 0, i(0) = 0:
sol:=dsolve({eq1,eq2,ICs},numeric,parameters=[V],output=listprocedure):
sol(parameters=[10]):
ode_x1:=sol[2];
ode_x1:=rhs(ode_x1);
ode_x2:=sol[3];
ode_x2:=rhs(ode_x2);

p1:=plot([ode_x1(t),ode_x2(t)],t=0..50,gridlines=true):
plots[display](array([p1]));

I have time wektor and input wektor V in csv file. How to change this code to simulate 

ode solution ode_x1 and ode_x2 with data from file? Now I set V as 

sol(parameters=[10])

and t in 

plot([ode_x1(t),ode_x2(t)],t=0..50,gridlines=true)

Best