登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

徐小湛的博客

高等数学课件、函数图形动画、数学软件应用

 
 
 

日志

 
 

圆柱面  

2009-07-13 20:58:57|  分类: 高等数学图形动画 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

 

所有图形见数学图形目录

 

圆柱面

圆柱面 - Calculus - 高等数学

圆柱面 - Calculus - 高等数学

with(plots):
a:=2:b:=2:
qumian:=plot3d([a*cos(theta),b*sin(theta),t],t=-1..4,theta=0..2*Pi,lightmodel=light1):
x_axis:=plot3d([u,0,0],u=-3..3,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-3..3,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-2.5..6,v=0..0.01,thickness=3):
xyz:=display(x_axis,y_axis,z_axis):
display(qumian,xyz,tickmarks=[4,4,4],axes=none,grid=[39,63],scaling=constrained);
 
 

圆柱面 - Calculus - 高等数学

with(plots): with(plottools):
x:=t->1:y:=t->0:z:=t->t:
u:=t->sqrt(x(t)^2+y(t)^2):
alpha:=t->piecewise(x(t)>0,arctan(y(t)/x(t)),arctan(y(t)/x(t))+Pi):
a:=-1.3:b:=1.3:
Muxian:=spacecurve([x(t),y(t),z(t)],t=a-0.1..b+0.1,thickness=3,color=blue):
qumian:=plot3d([u(t)*cos(theta+alpha(t)),u(t)*sin(theta+alpha(t)),z(t)],t=a..b,theta=0..2*Pi,style=patchnogrid):
K:=20 for i from 1 to K do ti:=i*2*Pi/K:
muxian[i]:=rotate(Muxian,ti,[[0,0,0],[0,0,1]]) od:
muxian:=display(seq(muxian[i],i=1..K)):
x_axis:=plot3d([u,0,0],u=-1.5..1.5,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-1.5..1.5,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-2..2,v=0..0.01,thickness=3):
xyz:=display(x_axis,y_axis,z_axis):
display(qumian,muxian,text,xyz,tickmarks=[4,4,4],axes=none,lightmodel=light1,
scaling=unconstrained);
 

圆柱面 - Calculus - 高等数学

with(plots): with(plottools):
x:=t->1:y:=t->0:z:=t->t:
u:=t->sqrt(x(t)^2+y(t)^2):
a:=-1.3:b:=1.3:R:=1:
Muxian:=spacecurve([x(t),y(t),z(t)],t=a..b,thickness=3,color=blue):
yuan1:=spacecurve([u(a)*cos(t),u(a)*sin(t),z(a)],t=0..2*Pi,thickness=3,color=blue):
yuan2:=spacecurve([u(b)*cos(t),u(b)*sin(t),z(b)],t=0..2*Pi,thickness=3,color=blue):
K:=20:for i from 1 to K do ti:=i*2*Pi/K:
muxian[i]:=rotate(Muxian,ti,[[0,0,0],[0,0,1]]) od:
muxian:=display(seq(muxian[i],i=1..K)):
x_axis:=plot3d([u,0,0],u=-1.5..1.5,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-1.5..1.5,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-2..2,v=0..0.01,thickness=3):
xyz:=display(x_axis,y_axis,z_axis):
display(yuan1,yuan2,muxian,orientation=[40,70],tickmarks=[4,4,4],
axes=none,lightmodel=light1,scaling=unconstrained);

 

圆柱面 - Calculus - 高等数学

with(plots):with(plottools):
x:=t->1:y:=t->0:z:=t->t:
u:=t->sqrt(x(t)^2+y(t)^2):
a:=-1.3:b:=1.3:
Muxian:=spacecurve([x(t),y(t),z(t)],t=a-0.5..b+0.5,thickness=3,color=red):
K:=70:for i from 1 to K do ti:=i*4*Pi/K:
qumian[i]:=plot3d([u(t)*cos(theta),u(t)*sin(theta),z(t)],t=a..b,theta=0..ti,scaling=constrained,style=patch):
muxian[i]:=rotate(Muxian,ti,[[0,0,0],[0,0,1]]) od:
qumian:=display(seq(qumian[i],i=1..K),insequence=true):
muxian:=display(seq(muxian[i],i=1..K),insequence=true):
x_axis:=plot3d([u,0,0],u=-1.5..1.5,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-1.5..1.5,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-2..2,v=0..0.01,thickness=3):
xyz:=display(x_axis,y_axis,z_axis):
display(qumian,muxian,xyz,orientation=[40,70],tickmarks=[4,4,4],
axes=none,lightmodel=light1,scaling=unconstrained);
 
圆柱面 - Calculus - 高等数学
 
with(plots):
x:=t->1:y:=t->0:z:=t->t:
u:=t->sqrt(x(t)^2+y(t)^2):a:=-2:b:=2:
qumian:=plot3d([u(t)*cos(theta),u(t)*sin(theta),z(t)],t=a..b,theta=0..2*Pi,scaling=constrained,style=patch):
quxian:=spacecurve([x(t),y(t),z(t)],t=a-1..b+1,thickness=4,color=red):
x_axis:=plot3d([u,0,0],u=-2.2..2.2,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-2.2..2.2,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-.2..3,v=0..0.01,thickness=3):
xyz:=display(x_axis,y_axis,z_axis):
display(qumian,quxian,xyz,orientation=[40,70],tickmarks=[4,4,4],axes=none,grid=[54,76],
lightmodel=light2,scaling=unconstrained);
 
圆柱面
 
圆柱面 - Calculus - 高等数学
圆柱面 - Calculus - 高等数学
 
with(plots):
a:=2:b:=2:
qumian:=plot3d([a*cos(theta),t,b*sin(theta)],t=-1..4,theta=0..2*Pi,style=patch,lightmodel=light1):
x_axis:=plot3d([u,0,0],u=-3..3,v=0..0.01,thickness=3):
y_axis:=plot3d([0,u,0],u=-3..4,v=0..0.01,thickness=3):
z_axis:=plot3d([0,0,u],u=-2.5..3,v=0..0.01,thickness=3):
xyz:=display(x_axis,y_axis,z_axis):
display(qumian,xyz,orientation=[35,70],tickmarks=[4,4,4],title="circular cylinder",
axes=none,grid=[39,63],scaling=constrained);
 
 
  评论这张
 
阅读(3527)| 评论(0)

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018