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

愛無止境

我想根據(jù)一篇文獻(xiàn)里的關(guān)系式重現(xiàn)齒輪圓錐曲線,，前面公式及命令如下：
fi1=[0:0.1:2*pi];
c=0.16;
- z) C  [3 L* |5 l0 c; r& ER=43.6;
4 |, O4 X  f/ }& A/ r5 tfi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接換成0.75代入關(guān)系式了,。
7 a) ?' M  B7 Y  C( {delta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
delta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
) m0 k" j/ c' g  X, M. J4 n0 {
3 f1 {. w; B% D# ?球面方程如下：
; a4 _1 I& @/ Nx1=R*(sin(delta1).*cos(fi1));
$ S6 Z8 V8 u- ?7 g) ^/ ]' N. D! xy1=R*sin(delta1).*sin(fi1);
z1=R*cos(delta1);
2 x4 D/ G* x( t& i; A4 \/ Fr1=x1*i+y1*j+z1*k     %這里的i,j,k前的系數(shù)是節(jié)曲線上的點(diǎn)在三個(gè)坐標(biāo)上的坐標(biāo),。寫到程序里面的不是這樣的,，下面也是一樣,。這樣寫出來大家理解直觀一些,。
x2=R*(sin(delta2).*cos(fi2));
y2=R*sin(delta2).*sin(fi2);
; B  @) J: }: M& C4 N$ o. a% Rz2=R*cos(delta2);
; B2 ~  A) P& Er2=x2*i+y2*j+z2*k
- O" l2 O% |3 h2 K* [
2 l2 e/ t) c" i* Z* Or1和r2的值驗(yàn)算出來的都是43.6，沒有問題,。但是我是不知道用什么樣的命令得到下面的球面圖形,。由于本人對matlab不很了解，到這里用命令老是出錯(cuò),，各位英雄女俠幫幫忙,。感激不盡！�,�,！


" x) k; ~( \* N- i