用MATLB软件完成对压床六杆机构的运动学分析，编写程序代码，求出所有从动件的位移曲线，速度曲线，加速度曲线。

以下是一份MATLAB程序代码，用于完成对压床六杆机构的运动学分析，并求出从动件的位移曲线、速度曲线和加速度曲线：

clc;
clear all;
close all;

% 输入参数
a = 0.2; % AB杆长
b = 0.3; % BC杆长
l = 0.5; % CD杆长
h = 0.4; % AD杆长
omega = 2; % 转速
theta0 = 0; % 初始角度
t = linspace(0,5,500); % 时间数组

% 计算AB杆的角速度和角位移
thetaAB = omega*t + theta0;
omegaAB = omega*ones(1,length(t));

% 计算BC杆的角速度和角位移
thetaBC = acos((a*cos(thetaAB) - b)/l);
omegaBC = -(a*omegaAB.*sin(thetaAB))./(l*sin(thetaBC));

% 计算CD杆的角速度和角位移
thetaCD = acos((a*cos(thetaAB) + h)/l) + thetaBC;
omegaCD = -(a*omegaAB.*sin(thetaAB) + h*omegaAB.*cos(thetaAB))./(l*sin(thetaCD));

% 计算从动件的位移、速度和加速度
x = a*cos(thetaAB) + l*cos(thetaCD);
y = a*sin(thetaAB) + l*sin(thetaCD) - h;
v = sqrt((a*omegaAB.*sin(thetaAB) - l*omegaCD.*sin(thetaCD)).^2 + (a*omegaAB.*cos(thetaAB) + l*omegaCD.*cos(thetaCD)).^2);
a = sqrt((a*omegaAB.^2.*cos(thetaAB) - l*omegaCD.^2.*cos(thetaCD)).^2 + (a*omegaAB.^2.*sin(thetaAB) + l*omegaCD.^2.*sin(thetaCD)).^2);

% 绘制曲线图
figure;
subplot(3,1,1);
plot(t,x);
title('从动件的位移曲线');
xlabel('时间 (s)');
ylabel('位移 (m)');
grid on;

subplot(3,1,2);
plot(t,v);
title('从动件的速度曲线');
xlabel('时间 (s)');
ylabel('速度 (m/s)');
grid on;

subplot(3,1,3);
plot(t,a);
title('从动件的加速度曲线');
xlabel('时间 (s)');
ylabel('加速度 (m/s^2)');
grid on;

该程序首先输入了压床六杆机构的各个参数，包括AB杆长、BC杆长、CD杆长、AD杆长、转速和初始角度等。然后程序通过计算AB杆、BC杆和CD杆的角速度和角位移，最终求出从动件的位移曲线、速度曲线和加速度曲线，并用MATLAB绘图函数绘制出来。