太极图和八卦

fengaas

原文来自CTEX论坛：
http://bbs.ctex.org/viewthread.p ... hlight=%B0%CB%D8%D4

用asympteto所画的太极图及先天八卦和后天八卦。

1. texpreamble("\usepackage{CJK}\AtBeginDocument{\begin{CJK*}{GBK}{song}}
   
2.     \AtEndDocument{\clearpage\end{CJK*}}");  //中文支持
   
3. 
   
4. size(10cm);
   
5. 
   
6. real r1=0.15;
   
7. 
   
8. guide g1=arc((0,-0.5),0.5,-90,90)..arc((0,0.5),0.5,270,90)..arc((0,0),1,90,-90)..cycle;
   
9. 
   
10. filldraw(g1^^circle((0,0.5),r1));
    
11. draw(unitcircle);
    
12. filldraw(circle((0,-0.5),r1));
    
13. 
    
14. void gua(int nn, int pos=nn, real initangle=0, bool showname=false)
    
15. {
    
16.   int n,n1,n2,n3;
    
17.   real ang;
    
18.   n=nn%8;
    
19.   ang=initangle+45*(pos%8);
    
20.   n1=quotient(n,4);
    
21.   n2=quotient(n%4,2);
    
22.   n3=quotient(n%2,1);
    
23. // write(n1,n2,n3);
    
24. 
    
25.   pair p;
    
26.   path m0[];  // --- 阳
    
27.   path m1[];  // - - 阴
    
28.   string name[]={"乾","巽","离","艮","兑","坎","震","坤"};
    
29.   path g1[],g2[],g3[];
    
30.   real r,len=0.65,inter=0.15,interh=len/6,lw=6;  
    
31.   // len 长度，interh 阴的间隔长度，inter 三道之间的间隔, lw 线宽
    
32.   r=1+3*inter;
    
33.   p=(r*Cos(ang),r*Sin(ang));
    
34.   m0=(-len/2,0)--(len/2,0);
    
35.   m1=(-len/2,0)--(-interh/2,0)^^(interh/2,0)--(len/2,0);
    
36.   if(n1==0)
    
37.        g1=shift(p)*rotate(ang+90)*m0;
    
38.   else
    
39.        g1=shift(p)*rotate(ang+90)*m1;
    
40.   if(n2==0)
    
41.        g2=shift(p-(inter*Cos(ang),inter*Sin(ang)))*rotate(ang+90)*m0;
    
42.   else
    
43.        g2=shift(p-(inter*Cos(ang),inter*Sin(ang)))*rotate(ang+90)*m1;
    
44.   if(n3==0)
    
45.        g3=shift(p-(2inter*Cos(ang),2inter*Sin(ang)))*rotate(ang+90)*m0;
    
46.   else
    
47.        g3=shift(p-(2inter*Cos(ang),2inter*Sin(ang)))*rotate(ang+90)*m1;
    
48.   draw(g1,linewidth(lw)+linecap(0));
    
49.   draw(g2,linewidth(lw)+linecap(0));
    
50.   draw(g3,linewidth(lw)+linecap(0));
    
51.   if(showname)
    
52.    {
    
53.     label(scale(2)*rotate(ang-90)*name[n],p+(2inter*Cos(ang),2inter*Sin(ang)));
    
54.    }
    
55. }
    
56. 
    
57. bool showname;
    
58. real rot;
    
59. rot=-45;
    
60. showname=false;
    
61. //-------后天八卦--------
    
62. gua(0,0,rot,showname); //乾
    
63. gua(4,1,rot,showname); //兑
    
64. gua(7,2,rot,showname); //坤
    
65. gua(2,3,rot,showname); //离
    
66. gua(1,4,rot,showname); //巽
    
67. gua(6,5,rot,showname); //震
    
68. gua(3,6,rot,showname); //艮
    
69. gua(5,7,rot,showname); //坎
    
70. 
    
71. 
    
72. //-------先天八卦--------
    
73. rot=90;
    
74. showname=false;
    
75. //gua(0,0,rot,showname); //乾
    
76. //gua(1,1,rot,showname); //巽
    
77. //gua(2,2,rot,showname); //离
    
78. //gua(3,3,rot,showname); //艮
    
79. //gua(7,4,rot,showname); //坤
    
80. //gua(6,5,rot,showname); //震
    
81. //gua(5,6,rot,showname); //坎
    
82. //gua(4,7,rot,showname); //兑

复制代码

得到的效果图如下：
先天八卦

太极先天八卦1

太极先天八卦.png (4.89 KB, 下载次数: 0, 售价: 1 枚金币)

2009-9-10 13:19 上传

点击文件名下载附件
太极先天八卦2

后天八卦
太极后天八卦2.png (4.32 KB, 下载次数: 0, 售价: 1 枚金币)

2009-9-10 13:19 上传

点击文件名下载附件
太极后天八卦1

太极后天八卦.png (5.06 KB, 下载次数: 0, 售价: 1 枚金币)

2009-9-10 13:19 上传

点击文件名下载附件
太极后天八卦2

wayne

哈哈，我可是在八卦图中浸泡着长大的~~

无心人

楼上会武？

wayne

八卦与武何相关？

楼主的软件我很喜欢，用楼主的代码，我在此贴出一个大图来：

wayne

以前我用Mathematica画了一个阴阳鱼此消彼长的动画
刚好在硬盘里搜到了代码，运行出来了，也奉上

wcvery

谢谢！