1.4.3 分形_ 雪花曲线_ Mathematica 程序.docVIP

1.4.3 Koch雪花曲线的MATHEMATIC程序设计 （王积社） 1.4.3.1 试验问题： Koch雪花曲线的MATHEMATIC程序设计. 1.4.3.2 试验目的： ① 掌握MATHEMATIC迭代程序设计方法； ② 用MATHEMATIC作出Koch雪花曲线； ③ 通过试验，认识分形的妙处，掌握分形程序设计基本方法。 1.4.3.3 问题分析 1.4.3.3.1 算法分析： 1o 雪花曲线的生成元为： 可见，从几何作图方法来看，生成的方法为：把线段AB三等分，得到分点C、E，再把向量CE绕C逆时针旋转60o得到D点，这样图形就生成了。 2o 为了便于计算和处理，我们把图形放在坐标系中。设线段AB的两端点的坐标为A(x1,y1)、B(x2,y2),则Koch雪花曲线的生成元的关系可用向量表示为： 可见，要想确定雪花曲线的生成元，关键是确定C、D、E点的坐标，为此，应该分析图中的向量关系. 3o 分析图中的向量关系： ① 先分析向量： 这样即得： ② 再分析向量： 是由旋转60o而得，于是令矩阵:,则有: 而又 ， 所以 从而 于是我们可用来表示. 4o 以上过程可以说是从向量出发，求出了向量、，相当于从点集｛A、B｝出发求得了点集｛A、C、D、E、B｝. 5o 对点集｛A、C｝、{C、D}、{D、E}、{E、B}分别重复上面的过程（即迭代），…，如此下去（有限次） 6o 将所上面所生成的相应的点对连线，即可做出Koch雪花曲线. 1.4.3.3.2 Mathematica的画线功能： 要作出雪花曲线，需要将点连线，即需要Mathematica的画线功能，经查可知Mathematica的“画线”功能函数为： Line[{{x1,y1},{x2,y2},…}](P224)—— 依次连接相邻两点的线段，作如下试验，输入： a={{0,0},{1,0},{1,1},{0,1},{0,0}}; Line[a] 如果连线成功的话，应该是一个矩形，然而运行结果如下： 未得所想，原因是Line[a]只是一种“意连”，但不真正画出线条（或者说没有显示出线条），然而我们想到Show[]可以显示图形，于是试输入： Show[Line[a]] 希望显示出图形，可是运行结果却为： 可见有错，错误提示的意思是“Line”不是图形元素，于是需要将“Line”转换为图形元素，这需要Graphics[]来完成，总之于是应输入以下命令： Show[Graphics[Line[a]]] 这样可得到： 请大家注意这种做法。 1.4.3.3.3 算法大意： 1) 取点a(0,0)、b(0,1)，以其坐标做集合ab=｛a、 b｝ 记： L=|ab|=2 （ab的长度，即元素个数） 2) 求点c、e 3) 求点d 4) 令：tmp={a,c,d,e,b}。 5) 对点集{｛a、c｝、{c、d}、{d、e}、{e、b}}中的每对点分别重复上面的过程(迭代). 6) 依次连接上面所生成的相邻的点，即可划出雪花曲线。 将1)～4)设计成一个从点集ab出发，构作点集tmp的模块(定义为函数“koch”)，当“ab”中点较多时，应循环实现。 1.4.3.3.4 算法框图： 由上面的分析，可得以下的算法框图 1.4.3.4 程序设计 koch[ab_List]:= Block[{tmp={},i,L=Length[ab], jd=60Degree, sa=Sin[jd],ca=Cos[jd],c,d,e,T={{ca,-sa},{sa,ca}}}, For[i=1,iL,i++, c=ab[[i]]*2/3+ab[[i+1]]/3; e=ab[[i]]/3+ab[[i+1]]*2/3; d=c+T.(e-c); tmp=Join[tmp,{ab[[i]],c,d,e,ab[[i+1]]}] ]; Return[tmp] ] pt={{0,0},{1,0}}; Show[Graphics[Line[Nest[koch,pt,5]],AspectRatio-Sqrt[3]/6]] ⑸ 运行结果： 　　迭代１次： 迭代２次： 迭代３次： 迭代４次： 迭代５次： 迭代６次： 1.4.3.6 结果分析 上面只是做出了Koch雪花曲线的一支，做法是对点A(0,0)、B(0,1)迭代生成的，而整个Koch雪花曲线应该是对正三角形： 的三边生成，请自己考虑写出做整个Koch雪花曲线的Math