分形图基本图形以及源程序.docVIP

分形图基本图形以及源程序 第一部分 本人新手，如有错误请指正。程序完成于2011/6/17晚间到2011/6/18。 很多变量名称采用的是同学的姓名拼音，为的是告诉大家这些都是可以随意命名的变量或函数名，一般大写字母开头的是系统定义的变量不可以随意更改。 (*雪花*) 源程序 lovelyduwangen[zhengguojie_List]:=Block[{weihuayan={},i,wuxiaonan=Length[zhengguojie],gengping=60Degree,sa=Sin[gengping],ca=Cos[gengping],c,d,e,T={{ca,-sa},{sa,ca}}}, For[i=1,i wuxiaonan,i++,c=zhengguojie[[i]]*2/3+zhengguojie[[i+1]]/3; e=zhengguojie[[i]]/3+zhengguojie[[i+1]]*2/3; d=c+T.(e-c); weihuayan=Join[weihuayan,{zhengguojie[[i]],c,d,e,zhengguojie[[i+1]]}]]; weihuayan] dongquanfa={{0,0},{1/2,Sqrt[3]/2},{1,0},{0,0}}; Show[Graphics[Line[Nest[lovelyduwangen,dongquanfa,0]],AspectRatio?Sqrt[3]/2]] Show[Graphics[Line[Nest[lovelyduwangen,dongquanfa,5]],AspectRatio?Sqrt[3]/2]] 基本生成元 分形图 (*上三角下三角1*) lovelyduwangen[zhengguojie_List]:=Block[{weihuayan={},i, wuxiaonan =Length[zhengguojie],gengping=60Degree,sa=Sin[gengping],ca=Cos[gengping],c,d,g,f,e,T={{ca,-sa},{sa,ca}},S={{ca,sa},{-sa,ca}}}, For[i=1,i wuxiaonan,i++,c=zhengguojie[[i]]*3/4+zhengguojie[[i+1]]/4; e=zhengguojie[[i]]/2+zhengguojie[[i+1]]/2; f=zhengguojie[[i]]/4+zhengguojie[[i+1]]*3/4; d=c+T.(e-c); g=e+S.(f-e); weihuayan=Join[weihuayan,{zhengguojie[[i]],c,d,e,g,f,zhengguojie[[i+1]]}]]; weihuayan] dongquanfa={{0,0},{1,0}}; Show[Graphics[Line[Nest[lovelyduwangen,dongquanfa,1]],AspectRatio?Sqrt[3]/2]] Show[Graphics[Line[Nest[lovelyduwangen,dongquanfa,4]],AspectRatio?Sqrt[3]/2]] 基本生成元 分形图 (*上三角下三角2*) lovelyduwangen[zhengguojie_List]:=Block[{weihuayan={},i, wuxiaonan =Length[zhengguojie],gengping=60Degree,sa=Sin[gengping],ca=Cos[gengping],c,d,g,f,e,T={{ca,-sa},{sa,ca}},S}, For[i=1,i wuxiaonan,i++,c=zhengguojie[[i]]*3/4+zhengguojie[[i+1]]/4; e=zhengguojie[[i]]/2+zhengguojie[[i+1]]/2; f=zhengguojie[[i]]/4+zhengguojie[[i+1]]*3/4; d=c+T.(e-c); S=Transpose[T]; g=e+S.(f-e); weihuayan=Join[weihuayan,{zhengguojie[[i]],c,d,e,g,f,zhengguojie[[i+1]]}]]; weihuay