AF长度图像，程序如下： plot(x1,p(:,1))
xlabel('\\theta1/rad'); ylabel('L1''/mm'); 输出图像如下：
角速度：
plot(x1,q(1,:),'--',x1,q(2,:),':',x1,q(3,:),'*') xlabel('W1 rad/s');
ylabel('W2、W3、W4 rad/S');
11
title('角速度');
F点速度： plot(x1,q(4,:)) xlabel('W1 rad/s'); ylabel('V mm/s'); title('F点速度');
角加速度：
plot(x1,f(1,:),'--',x1,f(2,:),':',x1,f(3,:),'*') xlabel('W1 rad/s');
12
ylabel('\\alpha2、\\alpha3、\\alpha4 rad/S^2'); title('角加速度');
F点的加速度：plot(x1,f(4,:)) xlabel('W1 rad/s'); ylabel('a mm/s^2'); title('F点的加速度');
13
机械原理大作业2
4-23 在图示的正弦机构中,已知lAB =100 mm,h1=120 mm,h2 =80 mm,W1 =10 rad/s(常
数)，滑块2和构件3的重量分别为 G2 =40 N和 G3 =100 N,质心S2 和S3 的位置如图所示，加于构件3上的生产阻力Fr=400 N，构件1的重力和惯性力略去不计。试用解析法求机构在Φ1=60°、150°、220°位置时各运动副反力和需加于构件1上的平衡力偶Mb 。
分别对三个构件进行受力分析如图：
14
构件3受力图
构件2受力图
构件1受力图
（1）滑块2：VS2 =LAB W1 ① as2 = LAB W12 ② 构件3：S=LAB sinΦ1 ③
15
