- 积分
- 0
- 注册时间
- 2007-5-23
- 仿真币
-
- 最后登录
- 1970-1-1
|
我在Maple 11写了一个动坐标的程序。但是最后总出错。请各位大哥帮个忙。文件也上传了。谢谢!
c:=t->[cos(t),sin(t),t/6]:
ll:=plots[spacecurve](c(t), t = 0 .. 4*Pi, thickness = 1, numpoints = 100, color = COLOR(RGB, 0, 0, 1)):
plots[spacecurve](c(t), t = 0 .. 4*Pi, thickness = 1, numpoints = 100, color = COLOR(RGB, 0, 0, 1),axes = BOXED):
curve
e11:=t->evalf(eval(diff(c(u),u)/sqrt(add(diff(c(u),u)*diff(c(u),u),i= 1..3)),u = t)):
d1:=u->diff(c(u),u):
d2:=u->diff(c(u),u$2):
cross22:=u->convert(LinearAlgebra[CrossProduct](Vector(d1(u)),Vector(d2(u))),list):
e33:=t->evalf(eval(cross22(u)/sqrt(add(cross22(u)*cross22(u),i=1..3)),u=t)):
e22:=t->convert(LinearAlgebra[CrossProduct](Vector(e33(t)),Vector(e11(t))),list):
ll1:=t->plots[display](plottools[arrow](c(t),c(t)+e11(t),0.03,0.08,0.1,cylindrical_arrow,color=COLOR(RGB,0,1,0))):
ll2:=t->plots[display](plottools[arrow](c(t),c(t)+e22(t),0.03,0.08,0.1,cylindrical_arrow,color=COLOR(RGB,0,0,1))):
NULL;
ll3:=t->plots[display](plottools[arrow](c(t),c(t)+e33(t),0.03,0.08,0.1,cylindrical_arrow,color=COLOR(RGB,1,0,0))):
point1:=t->plots[display](plottools[curve]([c(t)+e11(t),c(t)+e11(t)+e22(t),c(t)+e22(t)]),plottools[curve]([c(t)+e22(t),c(t)+e22(t)+e33(t),c(t)+e33(t)]),plottools[curve]([c(t)+e33(t),c(t)+e33(t)+e11(t),c(t)+e11(t)]),linestyle = dash,color=COLOR(RGB,1,0,0)):
dzb:=t->plots[pointplot3d](evalf(c(t)),symbol = Solidcircle,symbolsize = 30,color=COLOR(RGB,0,0,0)):
hsh:=t->plots[display]({ll,ll1(t),ll2(t),ll3(t),point1(t),dzb(t)}):
hsh(2*Pi);
plots[animate3d](hsh(t), [fw], fw = 0 .. 4*Pi); |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|