| 查看: 645 | 回复: 3 | |||
[求助]
VB 分形算法实现科赫格子
|
|
自己写的，加大写A的地方也就是角度的变化不太清楚如何实现？请教高手啊！新虫，金币不多，谅解。 Const pi = 3.14159 Private Sub Form_Click() ScaleTop = 50 ScaleLeft = 0 ScaleWidth = 100 ScaleHeight = -50 Call koch(20, 20, 80, 20, 1) End Sub Sub koch(ax As Single, ay As Single, bx As Single, by As Single, c As Single) Dim cx As Single Dim cy As Single Dim dx As Single Dim dy As Single Dim ex As Single Dim ey As Single Dim fx As Single Dim fy As Single Dim l As Single Dim alpha As Single If (bx - ax) * (bx - ax) + (by - ay) * (by - ay) < c Then Line (ax, ay)-(bx, by) Else cx = ax + (bx - ax) / 3 cy = ay + (by - ay) / 3 ex = bx - (bx - ax) / 3 ey = by - (by - ay) / 3 l = Sqr((ex - cx) * (ex - cx) + (ey - cy) * (ey - cy)) AA”If (ex - cx) > 0 Then AA”alpha = pi / 2 AA”Else AA”alpha = 0 / pi AA” End If dy = cy + Sin(alpha) * l dx = cx + Cos(alpha) * l fx = ex + Cos(alpha) * l fy = ey + Sin(alpha) * l Call koch(ax, ay, cx, cy, c) Call koch(cx, cy, ex, ey, c) Call koch(ex, ey, bx, by, c) Call koch(cx, cy, dx, dy, c) Call koch(dx, dy, fx, fy, c) Call koch(fx, fy, ex, ey, c) End If End Sub  0C99Z%H8}531TL`$N)FUAA6.jpg |
回复此楼» 猜你喜欢
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
-
本人考085602 化学工程 专硕
已经有8人回复
2楼2013-01-02 10:32:43
|
Const pi = 3.14159 Private Sub Form_Click() ScaleTop = 50 ScaleLeft = 0 ScaleWidth = 100 ScaleHeight = -50 Call koch(5, 0, 20, 20, 80, 20, 1) End Sub Sub koch(n As Integer, s As Integer, ax As Single, ay As Single, bx As Single, by As Single, c As Single) Dim cx As Single Dim cy As Single Dim dx As Single Dim dy As Single Dim ex As Single Dim ey As Single Dim fx As Single Dim fy As Single Dim l As Single Dim alpha As Single If n = 2 Then Line (ax, ay)-(bx, by) Else cx = ax + (bx - ax) / 3 cy = ay + (by - ay) / 3 ex = bx - (bx - ax) / 3 ey = by - (by - ay) / 3 l = Sqr((ex - cx) * (ex - cx) + (ey - cy) * (ey - cy)) If s = 0 Then alpha = pi / 2 ElseIf s = 1 Then alpha = 2 * pi ElseIf s = -1 Then alpha = pi End If dy = cy + Sin(alpha) * l dx = cx + Cos(alpha) * l fx = ex + Cos(alpha) * l fy = ey + Sin(alpha) * l Call koch(n - 1, 0, ax, ay, cx, cy, c) Call koch(n - 1, 1, cx, cy, dx, dy, c) Call koch(n - 1, 0, dx, dy, fx, fy, c) Call koch(n - 1, -1, fx, fy, ex, ey, c) Call koch(n - 1, 0, ex, ey, bx, by, c) Call koch(n - 1, 2, cx, cy, ex, ey, c) End If End Sub 迭代次数为一时正确，但是大于1就有问题了 |
3楼2013-01-04 13:24:16
4楼2013-01-04 13:24:51
20