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

愛無止境

我想根據(jù)一篇文獻(xiàn)里的關(guān)系式重現(xiàn)齒輪圓錐曲線，前面公式及命令如下：
fi1=[0:0.1:2*pi];
$ e0 a9 l$ J% i7 B2 ec=0.16;
- u  G, N8 q4 w! g! d8 HR=43.6;
4 k: B/ v6 B# Dfi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接換成0.75代入關(guān)系式了,。
1 {5 _# M4 ~) M# _delta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
delta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
8 S2 `% Z# [0 x$ j5 r" A
6 c" P9 K, B/ O  y( u) ?球面方程如下：
x1=R*(sin(delta1).*cos(fi1));
) Y# [+ b( Y' X; Q& x( N0 y9 Dy1=R*sin(delta1).*sin(fi1);
z1=R*cos(delta1);
% v' B3 l5 K  tr1=x1*i+y1*j+z1*k     %這里的i,j,k前的系數(shù)是節(jié)曲線上的點(diǎn)在三個(gè)坐標(biāo)上的坐標(biāo),。寫到程序里面的不是這樣的,，下面也是一樣。這樣寫出來大家理解直觀一些,。
x2=R*(sin(delta2).*cos(fi2));
  w$ Y* ^+ e! Oy2=R*sin(delta2).*sin(fi2);
z2=R*cos(delta2);
/ E& N" i' ]  h0 G7 C% }& mr2=x2*i+y2*j+z2*k

r1和r2的值驗(yàn)算出來的都是43.6,，沒有問題。但是我是不知道用什么樣的命令得到下面的球面圖形,。由于本人對(duì)matlab不很了解,，到這里用命令老是出錯(cuò)，各位英雄女俠幫幫忙,。感激不盡�,。�,！


; ?. u: Z) N% m6 y6 [5 H