电算的源程序（MATLAB）1.滑块的位移源程序x1=0:0.05:4*pi; %x1--原动件的角度变量b=150; %b--原动件的长度c=200; %c--连杆的长度k=b/c;l=150*cos(x1)+sqrt(200^2-(150*sin(x1).^2)) %l--连杆的位移plot(x1,l);title('滑块的位移图像')xlabel('\it角度','FontSize',8)ylabel('位移大小','FontSize',8)2.滑块的速度源程序x1=0:0.05:4*pi; %x1--原动件的角度变量b=150; %b--原动件的长度c=200; %c--连杆的长度k=b/c;w=(2*pi*650)/60; %w--原动件的角速度x2=asin(-k*sin(x1)); %x2--连杆的角度v=b*w*sin(x1-x2)./cos(x2);plot(x1,v);title('滑块的速度图像')xlabel('\it角度','FontSize',8)ylabel('速度大小','FontSize',8)3.滑块的加速度源程序：x1=0:0.05:4*pi; %x1--原动件的角度变量b=150; %b--原动件的长度c=200; %c--连杆的长度k=b/c;w=(2*pi*650)/60; %w--原动件的角速度x2=asin(-k*sin(x1)); %x2--连杆的角度a=b*(w^2)*(cos(x1-x2)./cos(x2)+k*cos(x1).^2./cos(x2).^3); plot(x1,a);title('滑块的加速度图像')xlabel('\it角度','FontSize',8)ylabel('加速度大小','FontSize',8)4.从动件滚子的位移源程序：e=5; %凸轮的偏心距r=35; %凸轮的基圆半径u=0:0.001:360;s0=sqrt(r^2-e^2);s1=(4*(1-cos(pi*u/50))).*(0<=u & u<50); %凸轮推程运动阶段s2=8.*(u>=50&u<55); %远休止阶段s3=4*(1+cos(pi*(u-55)/50)).*(u>=55&u<105); %凸轮回程阶段s4=0.*(u>=105&u<=360); %近休止阶段s=s1+s2+s3+s4;plot(u,s);axis([0,360,0,15]);title('滚子的运动线图');xlabel('角度');ylabel('位移');5.凸轮的轮廓线源程序：e=5; %凸轮的偏心距r=35; %凸轮的基圆半径r1=5; %滚子的半径s0=sqrt(r^2-e^2); %区间分段x1=50;x2=55;x3=105;x4=360;w=pi*650/60; %凸轮的角速度u=0:0.05:360;s1=(4*(1-cos(pi*u/50))).*(0<=u & u<x1); %凸轮推程运动阶段v1=pi*8.*sin(pi*u/50).*(0<=u & u<x1);dx1=((v1-e).*sin(u*pi/180)+(s0+s1).*cos(u*pi/180)).*(0<=u & u<x1);dy1=((v1-e).*cos(u*pi/180)-(s0+s1).*sin(u*pi/180)).*(0<=u & u<x1);X1=((s0+s1).*sin(u*pi/180)+e*cos(u*pi/180)).*(0<=u &u<x1);Y1=((s0+s1).*cos(u*pi/180)-e*sin(u*pi/180)).*(0<=u &u<x1);XS1=X1-r1*dy1./sqrt(dx1.^2+dy1.^2);YS1=Y1+r1*dx1./sqrt(dx1.^2+dy1.^2);s2=8; %远休止阶段v2=0;dx2=((v2-e).*sin(u*pi/180)+(s0+s2).*cos(u*pi/180)).*(x1<= u & u<x2);dy2=((v2-e).*cos(u*pi/180)-(s0+s2).*sin(u*pi/180)).*(x1<= u & u<x2);X2=((s0+s2)*sin(u*pi/180)+e.*cos(u*pi/180)).*(x1<=u &u<x2);Y2=((s0+s2)*cos(u*pi/180)-e.*sin(u*pi/180)).*(x1<=u &u<x2);XS2=X2-r1*dy2./sqrt(dx2.^2+dy2.^2);YS2=Y2+r1*dx2./sqrt(dx2.^2+dy2.^2);s3=4*(1+cos(pi*(u-x2)/50)).*(x2<=u&u<x3); %凸轮回程阶段v3=-pi*8.*sin(pi*(u-x2)/50).*(x2<=u&u<x3);dx3=((v3-e).*sin(u*pi/180)+(s0+s3).*cos(u*pi/180)).*(x2<= u&u<x3);dy3=((v3-e).*cos(u*pi/180)-(s0+s3).*sin(u*pi/180)).*(x2<= u&u<x3);X3=((s0+s3).*sin(u*pi/180)+e.*cos(u*pi/180)).*(x2<=u&u<x3 );Y3=((s0+s3).*cos(u*pi/180)-e.*sin(u*pi/180)).*(x2<=u&u<x3 );XS3=X3-r1*dy3./sqrt(dx3.^2+dy3.^2);YS3=Y3+r1*dx3./sqrt(dx3.^2+dy3.^2);%近休止阶段v4=0;s4=0;dx4=((v4-e).*sin(u*pi/180)+(s0+s4).*cos(u*pi/180)).*(x3<= u&u<x4);dy4=((v4-e).*cos(u*pi/180)-(s0+s4).*sin(u*pi/180)).*(x3<= u&u<x4);X4=(s0*sin(u*pi/180)+e.*cos(u*pi/180)).*(x3<=u & u<=x4); Y4=(s0*cos(u*pi/180)-e.*sin(u*pi/180)).*(x3<=u & u<=x4);XS4=X4-r1*dy4./sqrt(dx4.^2+dy4.^2);YS4=Y4+r1*dx4./sqrt(dx4.^2+dy4.^2);X=X1+X2+X3+X4; %合并分段函数Y=Y1+Y2+Y3+Y4;plot(X,Y,'--',XS1,YS1,'r',XS2,YS2,'r',XS3,YS3,'r',XS4,YS4 ,'r');title('凸轮理论轮廓线和实际轮廓线的比较');xlabel('X');ylabel('Y');6.凸轮的理论轮廓线的部分数据：X 5 8.017 11.072 14.222 17.492 20.862 24.26927.608 30.748 33.549 35.879 37.797 39.25940.067 40.230 39.812 38.922 37.698 36.29034.834 33.439 32.167 30.842 29.282 27.50025.508 23.323 20.959 18.437 15.774 12.99010.108 7.149 4.136 1.091 -1.962 -5.000 -8.000-10.939 -13.795 -16.546 -19.171 -21.651 -23.965-26.097 -28.030 -29.750 -31.244 -32.500 -33.509-34.262 -34.755 -34.983 -34.945 -34.641 -34.073-33.247 -32.167 -30.842 -29.282 -27.500 -25.508-23.323 -20.959 -18.437 -15.774 -12.990 -10.108-7.149 -4.136 -1.091 1.962 5.000Y34.641 34.268 33.999 33.759 33.439 32.90832.035 30.711 28.866 26.478 23.579 20.36216.892 13.166 9.322 5.491 1.786 -1.721 -5.000-8.067 -10.973 -13.795 -16.546 -19.171 -21.651-23.965 -26.097 -28.030 -29.750 -31.244 -32.500-33.509 -34.262 -34.755 -34.983 -34.945 -34.641-34.073 -33.247 -32.167 -30.842 -29.282 -27.500-25.508 -23.323 -20.959 -18.437 -15.774 -12.990-10.108 -7.149 -4.136 -1.091 1.962 5.000 8.00010.939 13.795 16.546 19.171 21.651 23.96526.097 28.030 29.750 31.244 32.500 33.50934.262 34.755 34.983 34.945 34.6417.凸轮实际轮廓线的数据(略)。