MeteoInfoLab脚本示例：绘制山魈脸

MeteoInfo

用到了很多复数运算。

1. Z=+0.04+.020211004j
   
2. col=9
   
3. m=1200
   
4. cx=-1
   
5. cy=-0.5
   
6. l=2
   
7. x=linspace(cx-l,cx+l,m)
   
8. y=linspace(cy-l,cy+l,m)
   
9. [X,Y]=meshgrid(x,y)
   
10. C=X+1j*Y
    
11. for k in range(col):
    
12.     Z=conj(Z)**(Z)+C
    
13.     W=exp(-abs(Z))
    
14. 
    
15. axes(aspect=1.5, axis=False)
    
16. imshow(rot90(W))
    
17. axis([100, 1200, 200, 1000])

Z=+0.04+.020211004j<br /> col=9<br /> m=1200<br /> cx=-1<br /> cy=-0.5<br /> l=2<br /> x=linspace(cx-l,cx+l,m)<br /> y=linspace(cy-l,cy+l,m)<br /> [X,Y]=meshgrid(x,y)<br /> C=X+1j*Y<br /> for k in range(col):<br />     Z=conj(Z)**(Z)+C <br />     W=exp(-abs(Z))<br /> <br /> axes(aspect=1.5, axis=False)<br /> imshow(rot90(W))<br /> axis([100, 1200, 200, 1000])