小珠子在环内的凹槽中滚动仿真图

%% 小珠子在环内的凹槽中滚动仿真图

clc;clear;figure('numbertitle','off','name','小珠子在环内的凹槽中滚动仿真图');

global K R 

K=4;w=1.3;R=0.2;r=1;

t1=[0,80];H1=[0;0.3];

[t,H] = ode45(@dH,t1,H1);

z=r-cos(H(:,1));

x=sin(H(:,1)).*cos(w*t);

y=sin(H(:,1)).*sin(w*t);

plot3(x,y,z,'b','linewidth',1.2);

line([-1,1],[0,0],[0,0],'linestyle','-','LineWidth',1.1);

line([0,0],[-1,1],[0,0],'linestyle','-','LineWidth',1.1);

line([0,0],[0,0],[0,2],'linestyle','-','LineWidth',1.1);

axis([-1 1 -1 1 0 2]);legend('小珠子在环内的凹槽中滚动位置');

xlabel('x');ylabel('y');zlabel('z');title('小珠子在环内的凹槽中滚动仿真图');