机械必威体育网址

标题: 求助：matlab里球面坐标系下的非圆齿轮渐开线命令 [打印本页]

---

作者: 爱无止境    时间: 2010-1-19 20:30
标题: 求助：matlab里球面坐标系下的非圆齿轮渐开线命令
我想根据一篇文献里的关系式重现齿轮圆锥曲线，前面公式及命令如下：
0 P0 Y' {( A+ h  G  Lfi1=[0:0.1:2*pi];
4 c# y( {: U4 M6 Cc=0.16;
+ {, v9 I8 W. |& Q, v* \+ o' `R=43.6;
) q1 r7 q( |, a# j4 Sfi2=0.75*fi1+(c/3)*0.75*cos(3*fi1+3*pi/2);      %0.75=z1/z2=3/4,我直接换成0.75代入关系式了。
delta1=atan(0.75*(1-c*sin(3*fi1+3*pi/2)));
' g7 |2 ~; V+ l0 Z; O: ~3 N5 ~delta2=acot(0.75*(1-c*sin(3*fi1+3*pi/2)));
: c/ F- t; j; N8 q2 X6 `3 @+ s
. Q* l* M5 H) r! ?- b' d& `  [球面方程如下：
/ i6 e2 {( r3 h9 a. _* |) }x1=R*(sin(delta1).*cos(fi1));
/ q% f$ h6 j8 S; D$ K) G* jy1=R*sin(delta1).*sin(fi1);
* J: |" m+ h7 G" A# Lz1=R*cos(delta1);
: f& s. X2 o- r) f) C. Fr1=x1*i+y1*j+z1*k     %这里的i,j,k前的系数是节曲线上的点在三个坐标上的坐标。写到程序里面的不是这样的，下面也是一样。这样写出来大家理解直观一些。
6 l; E. Z# z, H( m- K6 p" c4 q2 kx2=R*(sin(delta2).*cos(fi2));
y2=R*sin(delta2).*sin(fi2);
+ R2 b, x* G( H) Y9 \! L# [% Q1 kz2=R*cos(delta2);
3 z- ^! e& l0 |. L" q) K3 r5 v0 {8 rr2=x2*i+y2*j+z2*k
, o1 j3 f; U0 o# U) Z
; H5 a! C( u5 o; s6 |8 tr1和r2的值验算出来的都是43.6，没有问题。但是我是不知道用什么样的命令得到下面的球面图形。由于本人对matlab不很了解，到这里用命令老是出错，各位英雄女侠帮帮忙。感激不尽！！！
3 b  D; w5 T0 `* j% V [attach]163490[/attach]
& J" E! k* l% k2 W- Y. u0 u( J
/ y* h% ~  }0 j7 L

---

欢迎光临 机械必威体育网址 (//www.szfco.com/)

Powered by Discuz! X3.4