请写出一份matlab代码，要求如下。做出凸轮机构型线的优化仿真设计

以下是一个凸轮机构型线的优化仿真设计的MATLAB代码示例：

% 凸轮机构型线优化仿真设计

% 定义凸轮轴线转动的角度范围
theta = linspace(0, 2*pi, 360);

% 定义凸轮参数
R = 20; % 凸轮半径
r = 10; % 从动件半径
e = 5; % 凸轮偏心距

% 定义从动件的位置和初始角度
x = zeros(size(theta));
y = r * ones(size(theta));
phi = zeros(size(theta));

% 计算凸轮型线
for i = 1:length(theta)
    x(i) = (R + e)*cos(theta(i)) - e*cos((R+r)/r*theta(i));
    y(i) = (R + e)*sin(theta(i)) - e*sin((R+r)/r*theta(i));
    phi(i) = atan2(y(i), x(i));
end

% 画出凸轮型线
figure;
plot(x, y);
axis equal;
title('凸轮机构型线');

% 计算从动件的运动轨迹
dx = r*cos(phi);
dy = r*sin(phi);
x2 = x + dx;
y2 = y + dy;

% 画出从动件的运动轨迹
figure;
plot(x2, y2);
axis equal;
title('从动件的运动轨迹');

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

% 画出从动件的速度和加速度曲线
figure;
subplot(2,1,1);
plot(theta(2:end), v);
title('从动件的速度');
subplot(2,1,2);
plot(theta(2:end-1), a);
title('从动件的加速度');

该程序可以绘制凸轮机构的型线，并计算从动件的运动轨迹、速度和加速度曲线。您可以根据需要进行修改和优化。