第12章-主题画典例分析.docVIP

第十二章 主题数学动画典例分析 图12-01~02 圆柱面渐伸展开为平面* 图 1 圆柱面渐伸展开为平面 restart:with plots : N: 24:M: 4:R: N:w: 1:r: 5: a: 0:b: 2*Pi:c: -r:d: r: u: i- 2*Pi*i/N:v: j- c+2*d*j/M: f: i,j - r*cos u i : g: i,j - -r*sin u i : h: i,j - v j : F: k- [MESH [seq [seq evalf [f i,j ,g i,j ,h i,j ] ,j 0..M ],i 0..R-w*k ],COLOR RGB,1,0,1 ]: PLOT3D ANIMATE seq F k ,k 0..R ,AXES BOX :; display ANIMATE seq F k ,k 0..R ,insequence true,axes frame,labels [x,y,z] :;W1: %: Q: i- [polygon [[r*cos u i ,r*sin u i ,c],[r*cos u i +r*sin u i * u i ,r*sin u i -r*cos u i * u i ,c],[r*cos u i +r*sin u i * u i ,r*sin u i -r*cos u i * u i ,d],[r*cos u i ,r*sin u i ,d]],color blue ]:Q i : display ANIMATE seq Q i ,i 0..N ,insequence true,axes frame : PLOT3D ANIMATE seq Q i ,i 0..N ,AXES BOX :WW: %: f: i,j - r*cos u i : g: i,j - -r*sin u i : h: i,j - v j : FF: PLOT3D MESH [seq [seq evalf [f i,j ,g i,j ,h i,j ] ,j 0..M ],i 0..R ],STYLE LINE ,COLOR RGB,0,1,0 : display W1,WW,FF,view [-25..8,-32..16,c..d],orientation [20,72],scaling constrained,axesfont [TIMES,ROMAN,5],scaling constrained ; 本主题动画，如果采用新版animate函数，可以使用如下实现方案。 restart:with plots :r: 5:c: -r:d: r: cylinderplot r,u 0..2*Pi,v -c..c,style line,grid [13,3] :B: %: FF: proc t display polygonplot3d [[r*cos t ,r*sin t ,c],[r*cos t +r*sin t * t ,r*sin t -r*cos t * t ,c],[r*cos t +r*sin t * t ,r*sin t -r*cos t * t ,d],[r*cos t ,r*sin t ,d]] ,cylinderplot r,u t..2*Pi,v -c..c,color red,grid [13,3],thickness 0,style PATCHNOGRID :end: animate FF,[t],t 0..2*Pi,frames 9,background B,scaling constrained,orientation [40,70],axes box ; 图 2 新方案动画的一帧 图12-03~05 圆柱面的截口曲线* 图 3 截平面位置变动时圆柱面的截口曲线 restart:with plots :with plottools : a: -12:b: 12:N: 48:M: 36:r: 1:R: 18: t: i- a+ b-a /N*i: s: j- 2*Pi*j/M: x: i,j - r*cos s j : y: i,j - r*sin s j : z: i,j - t i : plot3d [x i,j ,y i,j ,z i,j ],i 0..N,j 0..M,style wireframe,color green :ZM: %: k: n- tan -0.0001+Pi*n/2/R : PM: n- polygonplot3d [[-1,-1,-k n ],[1,-1,-k n ],[1,1,k n ],[-1,1,k n ]]