stick slip code for matlab
Show older comments
Problem plotting stick slip graph with this code, where am I going wrong?
[t,y] = ode45(@Stick_Slip,[0:0.01:10],[0.01; 0]);
figure
plot( y(:,1),y(:,2) )
xlabel( 'Position' )
ylabel( 'Velocity' )
title('Stick-Slip');
function dy = Stick_Slip(t,y)
m = 1;
k = 100;
V = 0.5;
Mu=0.5;
c = 2;
f = c*y(2) + k*y(1);
if f <= Mu*m*9.81
dy(1,1) = V;
dy(2,1) = 0;
else
dy(1,1) = y(2);
dy(2,1) = (f - k*y(1) - c*y(2)) / m;
end
end
3 Comments
John D'Errico
on 4 Mar 2020
And your question is what? The mind reading toolbox is fuzzy today,so we need for you to tell us what is on your mind.
Milind Khanna
on 4 Mar 2020
venkatesu u
on 12 Jun 2023
Moved: DGM
on 12 Jun 2023
Show me the Stick_Slip function
Answers (0)
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!