机械原理作业-程序及设计.docx

一、凸轮程序设计（matlab）a.位移程序% t表示转角，s表示位移t=0:0.01:pi/4; %升程阶段加速s=80*(t/pi*2).^2;hold onplot(t/pi*180,s);t=pi/4:0.01:pi/2; %升程阶段减速s=40-80*4/pi.^2*(pi/2-t).^2;hold onplot(t/pi*180,s); t= pi/2:0.01:19*pi/18; %远休止阶段s=40;hold onplot(t/pi*180,s);t=19*pi/18:0.01:4*pi/3; %回程阶段s=40*(1-35*[(18*t-19*pi)/(5*pi)].^4+84*[(18*t-19*pi)/(5*pi)].^5-70*[(18*t-19*pi)/(5*pi)].^6+20*[(18*t-19*pi)/(5*pi)].^7);hold onplot(t/pi*180,s);t=4*pi/3:0.01:2*pi; %近休止阶段s=0;hold onplot(t/pi*180,s);grid onhold offb．速度程序（ω=1）% t表示转角，令ω1=1t=0:0.01:pi/4; %升程加速阶段v=160/(pi/2).^2*t;hold onplot(t/pi*180,v/12.03211);t=pi/4:0.01:pi/2; %升程减速阶段v=160/(pi/2).^2*(pi/2-t);hold onplot(t/pi*180,v/12.03211);t= pi/2:0.01:19*pi/18; %远休止阶段v=0;hold onplot(t/pi*180,v/12.03211);t=19*pi/18:0.01:4*pi/3; %回程阶段v=-40*18/5/pi*(140*[(18*t-19*pi)/(5*pi)].^3-420*[(18*t-19*pi)/(5*pi)].^4+420*[(18*t-19*pi)/(5*pi)].^5-140*[(18*t-19*pi)/(5*pi)].^6);hold onplot(t/pi*180,v/12.03211);t=4*pi/3:0.01:2*pi; %近休止阶段v=0;hold onplot(t/pi*180,v/2/12.03211);grid onhold offc．加速度程序% t表示转角，令ω1=1t=0:0.001:pi/4; %升程加速阶段a=160/(pi/2).^2;hold onplot(t/pi*180,a/(46.53378));t=pi/4:0.001:pi/2; %升程减速阶段a=-160/(pi/2).^2;hold onplot(t/pi*180,a/(46.53378));t= pi/2:0.001:19*pi/18; %远休止阶段a=0;hold onplot(t/pi*180,a/(46.53378));t=19*pi/18:0.001:4*pi/3; %回程阶段a=-40/pi/(5*pi/18).^2*(420*[(18*t-19*pi)/(5*pi)].^2-1680*[(18*t-19*pi)/(5*pi)].^3+2100*[(18*t-19*pi)/(5*pi)].^4-840*[(18*t-19*pi)/(5*pi)].^5);hold onplot(t/pi*180,a/(46.53378));t=4*pi/3:0.001:2*pi; %近休止阶段a=0;hold onplot(t/pi*180,a/(46.53378));grid onhold offd．ds/dφ-s程序% t表示转角，x(横坐标)表示类速度ds/dφ，y（纵坐标）表示位移st=0:0.001:pi/4; %升程加速阶段x=640/pi.^2*t;y=2*40*(2*t/pi).^2;hold onplot(x/10y/4,-r);t=pi/4:0.001:pi/2; %升程减速阶段x=640/pi.^2*(pi/2-t);y=40-320*(pi/2-t).^2/pi/pi;hold on plot(x/10y/4,-r);t=19*pi/18:0.001:4*pi/3; %回程阶段x=80/pi*(-140*[(18*t-19*pi)/(5*pi)].^3+420*[(18*t-19*pi)/(5*pi)].^4-420*[(18*t-19*pi)/(5*pi)].^5+140*[(18*t-19*pi)/(5