Hi Ishit,
As per my understanding, you would like to store the calculated values “time”, “displacement”, “velocity” into a table. And you would like to change the equation as the “time-displacement” plot does not conform to the behaviour of the spring-mass-damper system.
Here is an example showing how you can implement the spring-mass-damper system.
% Calculate the derivatives
dxdt = v(i);
dvdt = -(k/m) * x(i) - (c/m) * v(i);
% Update the state variables using the fourth-order Runge-Kutta method
k1x = dxdt;
k1v = dvdt;
k2x = v(i) + 0.5 * dt * k1v;
k2v = -(k/m) * (x(i) + 0.5 * dt * k1x) - (c/m) * k2x;
k3x = v(i) + 0.5 * dt * k2v;
k3v = -(k/m) * (x(i) + 0.5 * dt * k2x) - (c/m) * k3x;
k4x = v(i) + dt * k3v;
k4v = -(k/m) * (x(i) + dt * k3x) - (c/m) * k4x;
% Update the state variables
t(i+1) = t(i) + dt;
x(i+1) = x(i) + (dt/6) * (k1x + 2*k2x + 2*k3x + k4x);
v(i+1) = v(i) + (dt/6) * (k1v + 2*k2v + 2*k3v + k4v);
You can store the calculated values by using the “array2table” function. Refer to below code for the same:
arrr = [t x v]
output_table = array2table(arrr);
Please refer to the following documentation to know more about “array2table” function:
I hope this resolves the issue you were facing.
