| 查看: 739 | 回复: 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 |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
filecode=后面第一个是大写字母
已经有4人回复
-
时间戳又变了8-15
已经有20人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
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
