Force = [13 469.11 406.6 239.91 104.5 73.24 62.83 38.09 31.58 10.74 0.33 -124.69 -291.37 -364.27];
Velocity = [13 -263.46 -157.61 -79.39 -40.52 -21.78 -11.48 0 10.08 19.92 39.13 78.96 156.28 262.66]/1000;

p = polyfit(Velocity,Force,1)
figure (1)
hold on
box on
plot(Velocity, Force, 'or')
xlabel('Velocity mm/sec')
ylabel('Force N')



Part 2

m=
a=
b=
alpha=

p = polyfit(Velocity,Force,1)
c=-p(1)
ce=c*((a/b)*cos(alpha))^2
omega_n=ce/(2*1*m)
ke=(omega_n)^2*m
k=ke/((a/b)*cos(alpha))^2


Part 3


t = 0:0.01:T;
T=30

H=0.5;
y = H.*(1-(1+(omega_n.*t)).*exp(-omega_n.*t));

figure(1)
hold on
box on
plot(t,y,'y')
xlabel('t, [sec]')
ylabel('Step Input, H, [m^2/s^2]')
title('Response of the Second order System')