求助：matlab里球面坐標(biāo)系下的非圓齒輪漸開(kāi)線(xiàn)命令

愛(ài)無(wú)止境

我想根據(jù)一篇文獻(xiàn)里的關(guān)系式重現(xiàn)齒輪圓錐曲線(xiàn)，前面公式及命令如下：
' T, v2 b6 G& q# a  ^& Ofi1=[0:0.1:2*pi];
c=0.16;
R=43.6;
fi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接換成0.75代入關(guān)系式了。
& F" Y9 E3 D- _5 Ldelta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
delta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
2 o3 v, f: e3 x4 p( `
球面方程如下：
x1=R*(sin(delta1).*cos(fi1));
, B: i3 E/ o( k6 W) {5 }2 Zy1=R*sin(delta1).*sin(fi1);
0 W; C1 g( D6 v- W- L/ u: Uz1=R*cos(delta1);
0 E% \! r' n9 R$ L' \r1=x1*i+y1*j+z1*k     %這里的i,j,k前的系數(shù)是節(jié)曲線(xiàn)上的點(diǎn)在三個(gè)坐標(biāo)上的坐標(biāo)。寫(xiě)到程序里面的不是這樣的，下面也是一樣。這樣寫(xiě)出來(lái)大家理解直觀一些。
" _3 W  p2 \  D! D2 L) W! Ux2=R*(sin(delta2).*cos(fi2));
3 Q) Y2 _3 s* c+ `& Ey2=R*sin(delta2).*sin(fi2);
& b3 Q% @3 K$ x1 z# tz2=R*cos(delta2);
7 g  L0 ]2 m/ D) {r2=x2*i+y2*j+z2*k

r1和r2的值驗(yàn)算出來(lái)的都是43.6，沒(méi)有問(wèn)題。但是我是不知道用什么樣的命令得到下面的球面圖形。由于本人對(duì)matlab不很了解，到這里用命令老是出錯(cuò)，各位英雄女俠幫幫忙。感激不盡?。?！
4 c7 X: d6 ?: {: K# b