| 查看: 775 | 回复: 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 |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
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
200
