帮我用matlab设计一平底直动从动件盘形凸轮机构。凸轮的基圆半径𝑟0 =40mm从动件升程ℎ =10mm推程运动角120度远休止角30度 回程运动角120度近休止角90度从动件在推程和回程均作余弦加速度运动。绘制出完整凸轮曲线并给出代码

首先，根据给出的凸轮参数，可以计算出凸轮的压力角、迎角、侧隙等。这里假设压力角为20度，迎角为5度，侧隙为0.01mm。

然后，根据从动件的升程、推程、远休止角、回程、近休止角等参数，可以计算出从动件在推程和回程上的余弦加速度运动规律。

接下来，可以通过计算得到凸轮上各点的坐标，并绘制出完整的凸轮曲线。

最后，给出MATLAB代码如下：

% 凸轮参数
r0 = 40; % 基圆半径
phi = 20; % 压力角
alpha = 5; % 迎角
e = 0.01; % 侧隙

% 从动件参数
h = 10; % 升程
theta1 = 120; % 推程运动角
theta2 = 30; % 远休止角
theta3 = 120; % 回程运动角
theta4 = 90; % 近休止角

% 计算从动件在推程和回程上的余弦加速度运动规律
t1 = linspace(0, pi/2, 100);
t2 = linspace(0, pi/3, 100);
t3 = linspace(0, pi/2, 100);
t4 = linspace(0, pi/2.5, 100);
a1 = h/2*(1-cos(t1));
v1 = h/2*sin(t1);
s1 = h/2*(t1-sin(t1));
a2 = 0;
v2 = h/2;
s2 = h/2*pi/3;
a3 = -h/2*(1-cos(t3));
v3 = h/2*sin(t3);
s3 = h/2*(pi/2+t3-sin(t3));
a4 = 0;
v4 = 0;
s4 = h/2*(pi/2.5+pi/2-sin(pi/2.5));

% 计算凸轮上各点的坐标
theta = linspace(0, 2*pi, 360);
r = r0 + h + e;
xc = r*cos(theta);
yc = r*sin(theta);
x = (r-h-e)*cos(theta+theta1/180*pi).*cos(phi/180*pi) + h/2*cos(theta1/180*pi).*cos(phi/180*pi-theta+alpha/180*pi);
y = (r-h-e)*sin(theta+theta1/180*pi) + h/2*sin(theta1/180*pi).*cos(phi/180*pi-theta+alpha/180*pi);

% 绘制凸轮曲线
figure;
plot(xc, yc, 'k', 'LineWidth', 1.5);
hold on;
plot(x, y, 'r', 'LineWidth', 1.5);
axis equal;
grid on;
xlabel('x');
ylabel('y');
title('平底直动从动件盘形凸轮机构');