[转]矩阵、随机化与分形图形

转自 http://www.51xuewen.com/group/513/topic_18097.htm

or: http://www.matrix67.com/blog/archives/500

里头的代码是

1. list = {{0, 0}};
   
2. last = {{0}, {0}};
   
3. For[i = 0, i < 50000, i++, r = Random[];
   
4.    If[r < 0.85, last = {{0.83, 0.03}, {-0.03, 0.86}}.last + {{0}, {1.5}},
   
5.      If[r < 0.91, last = {{0.2, -0.25}, {0.21, 0.23}}.last + {{0}, {1.5}},
   
6.        If[r < 0.99, last = {{-0.15, 0.27}, {0.25, 0.26}}.last + {{0}, {0.45}},
   
7.          last = {{0, 0}, {0, 0.17}}.last + {{0}, {0}}
   
8.        ]
   
9.      ]
   
10.    ];
    
11.    list = Append[list, First[Transpose[last]]];
    
12. ]
    
13. ListPlot[list, PlotStyle -> PointSize[0.002]]

复制代码