绘制三维水母

MeteoInfo

MeteoInfoLab脚本程序绘制三维水母的例子，用到了airy函数。

1. from mipylib.numeric.special import airy
   
2. 
   
3. [x,y,z] = sphere(80)
   
4. c = lambda x : cos(x * pi)
   
5. r = c(arange(-4, 4.1, .1))
   
6. d = lambda x, y, z : surf(x, y, z, alpha=.5, facecolor='interp',
   
7.     edgecolor=None, cmap='MPL_spring_r')
   
8. 
   
9. figure(facecolor='k', newfig=False)
   
10. axes3d(aspect='equal', axis=False, clip_plane=False)
    
11. lighting(position=[1,1,1,0], ambient=1, mat_specular=1)
    
12. 
    
13. d(x, y, z*cos(z)*r)
    
14. d(x*2, y*2, 2*abs(z+.4)-.5)
    
15. r = arange(1, 3.1, .1)
    
16. for o in r[:]:
    
17.     for n in r:
    
18.         d(x/4/n+n*c(o), c(o)*y/4/n+n*c(o+.5), z/9/n+airy(4*n-9)[0]/4-.7)

from mipylib.numeric.special import airy<br /> <br /> [x,y,z] = sphere(80)<br /> c = lambda x : cos(x * pi)<br /> r = c(arange(-4, 4.1, .1))<br /> d = lambda x, y, z : surf(x, y, z, alpha=.5, facecolor='interp',<br />     edgecolor=None, cmap='MPL_spring_r')<br /> <br /> figure(facecolor='k', newfig=False)<br /> axes3d(aspect='equal', axis=False, clip_plane=False)<br /> lighting(position=[1,1,1,0], ambient=1, mat_specular=1)<br /> <br /> d(x, y, z*cos(z)*r)<br /> d(x*2, y*2, 2*abs(z+.4)-.5)<br /> r = arange(1, 3.1, .1)<br /> for o in r[:]:<br />     for n in r:<br />         d(x/4/n+n*c(o), c(o)*y/4/n+n*c(o+.5), z/9/n+airy(4*n-9)[0]/4-.7)

MgLily

赞！！喜欢水母