数控铣床凸面球的编程可以通过宏程序来实现。以下是一个示例程序，用于加工半径为SR20、台展角为67°的凸球面，使用半径为R8的立铣刀。
```plaintext
% O0001
N10 M6 T 1
N20 G54 G90 G0 G43 H1 Z200
N30 M3 S2000
N40 X0 Y0
N50 Z2
N60 M8
N70 G1 Z0 F50
N80 1=0
N90 2=20
N100 WHILE[1LE67]DO 1
N110 3=8
N120 4=2*(1-COS[1])
N130 5=3+2*SIN[1]
N140 G1 X5 Y0 F200
N150 Z -4 F50
N160 G2 I-5 F200
N170 1=1+1
N180 END1
N190 G0 Z200
M9
N200 G49 G90 Z0
N210 M30
```
程序解释
初始化
`N10 M6 T 1`：设置主轴转速为1转/分钟，冷却液开启。
`N20 G54 G90 G0 G43 H1 Z200`：设置工件坐标系原点，取消刀具半径补偿，将Z轴坐标设置为200。
`N30 M3 S2000`：设置主轴转速为2000转/分钟。
`N40 X0 Y0`：将X和Y坐标设置为0。
`N50 Z2`：将Z坐标设置为2。
`N60 M8`：准备刀具。
`N70 G1 Z0 F50`：以速度F50开始加工，Z轴坐标设置为0。
`N80 1=0`：初始化角度变量。
`N90 2=20`：初始化球面半径变量。
循环加工
`N100 WHILE[1LE67]DO 1`：当角度小于等于67°时，执行循环。
`N110 3=8`：初始化铣刀半径变量。
`N120 4=2*（1-COS[1]）`：计算当前角度对应的铣刀半径。
`N130 5=3+2*SIN[1]`：计算当前角度对应的球面加工半径。
`N140 G1 X5 Y0 F200`：以速度F200沿X方向移动到计算出的X坐标。
`N150 Z -4 F50`：沿Z轴负方向移动到计算出的Z坐标，速度为F50。
`N160 G2 I-5 F200`：以速度F200沿I方向移动到计算出的I坐标。
`N170 1=1+1`：角度变量加1。
结束循环
`N180 END1`：结束循环。
`N190 G0 Z200`：回到初始Z坐标。
`M9`：取消主轴转速和冷却液。
`N200 G49 G90 Z0`：取消刀具半径补偿，将Z轴坐标设置为0。
`N210 M30`：取消主轴转速。
建议
在实际应用中，可能需要根据具体的加工条件和机床参数进行调整。
使用宏程序可以大大提高编程效率，减少重复代码。
在编写宏程序时，确保变量的定义和使用正确无误，以避免加工错误。