用MATLAB编写基圆为70，近休止角为120度，推程和回程角为105度，远休止角为30度，推程为40mm的用极坐标表示的凸轮廓线代码，

以下是MATLAB代码：

% 参数定义
r0 = 70; % 基圆半径
delta = deg2rad(120); % 近休止角，单位为弧度
beta1 = deg2rad(105); % 推程角1，单位为弧度
beta2 = deg2rad(-105); % 推程角2，单位为弧度
gamma = deg2rad(30); % 远休止角，单位为弧度
h = 40; % 推程高度

% 极坐标变量定义
t = linspace(0, 2*pi, 1000); % 角度变量
r = zeros(size(t)); % 半径变量

% 计算凸轮廓线
for i = 1:length(t)
    if t(i) <= delta || t(i) >= 2*pi-delta % 近休止角
        r(i) = r0;
    elseif t(i) >= delta+beta1 && t(i) <= 2*pi-delta-beta2 % 推程
        r(i) = r0 + h/2*(1-cos((t(i)-delta-beta1)*(pi/(beta2-beta1))));
    elseif t(i) >= 2*pi-delta-beta2+gamma && t(i) <= 2*pi-delta-gamma % 回程
        r(i) = r0 + h/2*(1-cos((t(i)-2*pi+delta+beta2-gamma)*(pi/(beta2-beta1))));
    else % 远休止角
        r(i) = r0;
    end
end

% 极坐标图像绘制
polarplot(t, r)

运行后可得到凸轮廓线的极坐标图像。