clear;clc;n=750;l=0.975;R=0.0381;h=0.2;omiga=n.*pi/30;tmax=2.*pi/omiga;t=0:0.001:tmax; %计算曲柄转一圈的总t值alpha1=atan(h+R.*sin(omiga.*t)./sqrt(l.*l-(h+R.*sin(omiga.*t)+pi;alpha1p=-(R.*omiga.*cos(omiga.*t)./(l.*cos(alpha1);vb=-R.*omiga.*sin(omiga.*t)+R.*omiga.*cos(omiga.*t).*tan(alpha1);ab=-R.*omiga.2.*cos(omiga.*t)-(R.*omiga.*cos(omiga.*t).2./(l.*(cos(alpha1).3) -R.*omiga.2.*sin(omiga.*t).*tan(alpha1);subplot(1,2,1);plot(t,vb);title(曲柄滑块机构的滑块v-t图);x
