MeteoInfoLab脚本示例：光线下的花朵

MeteoInfo

参考MATLAB代码用MeteoInfoLab程序绘制花朵，并加上光照效果。

1. # Ported from MATLAT code:
   
2. # https://www.mathworks.com/matlab ... sts/5/entries/12418
   
3. 
   
4. m=arange(0,1.01,.01)
   
5. m=np.expand_dims(m, axis=0)
   
6. T=arange(-4,25.01,.01)
   
7. x=1-(5./4*(1-mod(3.6*T, 2.))**2-1/4.)**2/2.
   
8. P=exp(-T/5.)/2
   
9. s=sin(P*pi)
   
10. c=cos(P*pi)
    
11. y=2*m.T**5*(1.6*m.T-1)**2*s
    
12. S=x*(m.T*s+y*c)
    
13. X=S*sin(T*pi)
    
14. Y=S*cos(T*pi)
    
15. Z=m.T*c+y*s
    
16. 
    
17. figure(bgcolor='k', newfig=False)
    
18. ax = axes3d(aspect='equal', clip_plane=False, axis=False)
    
19. grid(False)
    
20. lighting(mat_specular=1)
    
21. surf(X, Y, x*Z, X**2+Y**2+Z**4, facecolor='interp',
    
22.     edgecolor=None, cmap='MPL_autumn_r')
    
23. antialias(True)

# Ported from MATLAT code: <br /> # <a href=\"https://www.mathworks.com/matlabcentral/communitycontests/contests/5/entries/12418\" target=\"_blank\">https://www.mathworks.com/matlab ... sts/5/entries/12418</a><br /> <br /> m=arange(0,1.01,.01)<br /> m=np.expand_dims(m, axis=0)<br /> T=arange(-4,25.01,.01)<br /> x=1-(5./4*(1-mod(3.6*T, 2.))**2-1/4.)**2/2.<br /> P=exp(-T/5.)/2<br /> s=sin(P*pi)<br /> c=cos(P*pi)<br /> y=2*m.T**5*(1.6*m.T-1)**2*s<br /> S=x*(m.T*s+y*c)<br /> X=S*sin(T*pi)<br /> Y=S*cos(T*pi)<br /> Z=m.T*c+y*s<br /> <br /> figure(bgcolor='k', newfig=False)<br /> ax = axes3d(aspect='equal', clip_plane=False, axis=False)<br /> grid(False)<br /> lighting(mat_specular=1)<br /> surf(X, Y, x*Z, X**2+Y**2+Z**4, facecolor='interp', <br />     edgecolor=None, cmap='MPL_autumn_r')<br /> antialias(True)

随缘

{:eb502:}{:eb502:}

little卓忘密码

牛哇

yuxili

little卓忘密码 发表于 2023-11-26 15:03
牛哇

6 six,这种是怎么弄出来的