stick slip code for matlab

Milind Khanna

2020 3 3

0 个回答

5 次查看（30 天）

0 个投票

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 个评论
显示 1更早的评论隐藏 1更早的评论

Milind Khanna 2020-3-4

I am actually having problem in plotting the graph showing stick slip but I am not getting the desired output with the above code. Could you tell where I am doing the mistake?

Thank you.