function yp=qcar(t,z) M=1000; m=100; Cs=1000; Ks=36000; Kt=360000; R=150; yp=zeros(4,1); yp(1)=z(2); yp(2)=-1/M*[Ks(z(1)-z(3))+Cs(z(2)-z(4))]; yp(3)=z(4); yp(4)=-1/m*[Ks(z(1)-z(3))+Cs(z(2)-z(4))-Kt(z(3)-R)];And in another .m filez0=[0,0,0,0];[t,z]=ode113(qcar,[0:0.2:10],f0);plot(t,z);MATLAB: Getting “Not enough Input arguments” for this code. Can some one help
ode
Related Solutions
This code might do what you are looking for
ix = sscanf( T(:,end), '%1d' );Y = cat( 1, Y, rand ); % add one element to get 10
R = Y(ix+1);inspect the result
>> R(1:6)'ans = 11.1665 11.1665 11.1665 11.5740 11.5740 11.5740where
T=['cluster2''cluster2''cluster2''cluster7''cluster7''etc ....'];See odetovectorfield :https://www.mathworks.com/help/symbolic/odetovectorfield.html. It also have example to convert the symbolic ODE to the system of first order ODE which can be used to write a ss-model.
If you want to do it by hand, then you first need to write the two ODEs in this form
and then convert it into the state-space form
where 
.
.Following shows the code using ode2vectorfield:
syms ms mu ks cs binrt kt z0(t) z1(t) z2(t)Dz1 = diff(z1);DDz1 = diff(Dz1);Dz2 = diff(z2);DDz2 = diff(Dz2);eq1 = ms*DDz2 == ks*(z1-z2) + cs*(Dz1-Dz2) + binrt*(DDz1-DDz2);eq2 = mu*DDz1 == ks*(z1-z2) + cs*(Dz1-Dz2) + binrt*(DDz1-DDz2) + kt*(z0-z1);[V,S] = odeToVectorField(eq1, eq2)Result
>> VV = Y[2] (binrt*kt*z0(t) + kt*ms*z0(t) - binrt*kt*Y[1] + cs*ms*Y[2] - cs*ms*Y[4] + ks*ms*Y[1] - ks*ms*Y[3] - kt*ms*Y[1])/(binrt*mu - binrt*ms + ms*mu) Y[4] (binrt*kt*z0(t) - binrt*kt*Y[1] + cs*mu*Y[2] - cs*mu*Y[4] + ks*mu*Y[1] - ks*mu*Y[3])/(binrt*mu - binrt*ms + ms*mu)>> SS = z1 Dz1 z2 Dz2You can also use this outout to write a state-space model.
Best Answer