用MATLB编写程序代码，求出压床机构中所有从动件的位移曲线，速度曲线，加速度曲线。

这里提供一种可能的MATLAB程序代码，假设压床机构为四杆机构，其中的连杆长度和转动角度已知，求解从动件的位移曲线、速度曲线和加速度曲线。

% 压床机构四杆机构的连杆长度和转动角度
L1 = 0.2;
L2 = 0.3;
L3 = 0.4;
L4 = 0.5;
theta1 = 0:0.01:2*pi;

% 计算从动件的位移曲线
x = L1*cos(theta1) + L2*cos(theta1+acos((L1^2+L2^2-L3^2)/(2*L1*L2))) + L4*cos(theta1-acos((L2^2+L4^2-L1^2)/(2*L2*L4)));
y = L1*sin(theta1) + L2*sin(theta1+acos((L1^2+L2^2-L3^2)/(2*L1*L2))) + L4*sin(theta1-acos((L2^2+L4^2-L1^2)/(2*L2*L4)));

% 计算从动件的速度曲线和加速度曲线
dx = diff(x);
dy = diff(y);
v = sqrt(dx.^2 + dy.^2);
a = diff(v);

% 绘制从动件的位移曲线、速度曲线和加速度曲线
figure
subplot(3,1,1)
plot(theta1, x, 'b')
hold on
plot(theta1, y, 'r')
xlabel('转动角度')
ylabel('位移')
legend('x', 'y')
subplot(3,1,2)
plot(theta1(1:end-1), v, 'g')
xlabel('转动角度')
ylabel('速度')
subplot(3,1,3)
plot(theta1(1:end-2), a, 'm')
xlabel('转动角度')
ylabel('加速度')

其中，第一部分计算从动件的位移曲线，第二部分计算从动件的速度曲线和加速度曲线，第三部分绘制曲线图像。需要注意的是，由于计算速度和加速度需要对位移曲线进行差分，因此在绘制速度和加速度曲线时需要将位移曲线的最后一个点和倒数第二个点去掉，以使曲线长度一致。