机械工程郑佳文学号：1508520102 任务：根据杆长判断机构类型，模拟仿真曲柄滑块机构运动，并绘制滑块速度及加速度图像。

源代码：
clear all
clc
l1=4,l2=10,e=0,w1=10;
if(e==0) & (l1>l2)
% load handel
% sound(y,Fs)
disp('当e=0时,l1需小于l2。

不满足杆长条件，请重新输入l1,l2,e的值')
elseif((e~=0)&(l1+e>l2))
% load handel
% sound(y,Fs)
disp('当e~=0时,l1+e需小于l2。

不满足杆长条件，请重新输入l1,l2,e的值') else
end
fai=linspace(0,2*pi);
for i=1:10:length(fai)
JGT2(fai(i),l1,l2,e,w1)
end
%%%%%%%%%%%计算滑块速度并绘制图像
figure (2)
fai=0:0.01:2*pi;
v=-l1.*sin(fai).*w1-(l2^2-(l1.*sin(fai)-e). ^2).^(-0.5)*(l1^2).*(l1.*sin(fai)-e).*sin(f ai).*cos(fai).*w1;
plot(fai,v,'-b')
title('滑块角度-速度图像')
xlabel('角度\phi/°')
ylabel('速度/m/s')
figure (3)
A=- l1.*cos(fai) -
(l1^2.*cos(fai).^2)./(l2^2 - (e -
l1.*sin(fai)).^2).^(1/2) - (l1.*sin(fai).*(e - l1.*sin(fai)))./(l2^2 - (e -
l1.*sin(fai)).^2).^(1/2) ...
- (l1^2.*cos(fai).^2.*(e -
l1.*sin(fai)).^2)./(l2^2 - (e -
l1.*sin(fai)).^2).^(3/2);
plot(fai,A,'-b')
title('滑块角度-加速度图像')
xlabel('角度\phi/°') ylabel('速度/m/s')
结果图：
-2
2
468
10
12
14
滑块位移图像
角度φ/°
滑块位移/
m
0123
4567
滑块角度-速度图像
角度φ/°
速度/m /s
0123
4567
-6
-5-4-3-2-1012
3滑块角度-加速度图像
角度 /°
速度/m /s。