混沌之蝴蝶效应.doc

下载文档

2
0
约3.02千字
约 3页
2018-03-29 发布于河南
举报
版权申诉
保障服务

混沌之蝴蝶效应.doc

1、本文档共3页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

混沌之蝴蝶效应

混沌之蝴蝶效应 Public Function fx(x, y, z) fx = -10 * x + 10 * y End Function Public Function fy(x, y, z) fy = 28 * x - y - x * z End Function Public Function fz(x, y, z) fz = x * y - 8 * z / 3 End Function Private Sub option1_click() Picture1.Cls Picture1.Scale (-26, 51)-(26, -1) Picture1.Line (-25, 50)-(25, 0), vbRed, B h = 0.01 x = 1 y = 1 z = 1 For n = 0 To 5000 X1 = h * fx(x, y, z) Y1 = h * fy(x, y, z) Z1 = h * fz(x, y, z) X2 = h * fx(x + X1 / 2, y + Y1 / 2, z + Z1 / 2) Y2 = h * fy(x + X1 / 2, y + Y1 / 2, z + Z1 / 2) Z2 = h * fz(x + X1 / 2, y + Y1 / 2, z + Z1 / 2) X3 = h * fx(x + X2 / 2, y + Y2 / 2, z + Z2 / 2) Y3 = h * fy(x + X2 / 2, y + Y2 / 2, z + Z2 / 2) Z3 = h * fz(x + X2 / 2, y + Y2 / 2, z + Z2 / 2) X4 = h * fx(x + X3, y + Y3, z + Z3) Y4 = h * fy(x + X3, y + Y3, z + Z3) Z4 = h * fz(x + X3, y + Y3, z + Z3) xa = x ya = y za = z x = x + (X1 + 2 * X2 + 2 * X3 + X4) / 6 y = y + (Y1 + 2 * Y2 + 2 * Y3 + Y4) / 6 z = z + (Z1 + 2 * Z2 + 2 * Z3 + Z4) / 6 If n 2000 Then Picture1.Line (x, z)-(xa, za), vbYellow Picture1.PSet (x, z), vbRed End If Next End Sub Private Sub Option2_Click() Picture1.Cls Picture1.Scale (-26, 31)-(26, -31) Picture1.Line (-25, 30)-(25, -30), vbRed, B h = 0.01 x = 1 y = 1 z = 1 For n = 0 To 5000 X1 = h * fx(x, y, z) Y1 = h * fy(x, y, z) Z1 = h * fz(x, y, z) X2 = h * fx(x + X1 / 2, y + Y1 / 2, z + Z1 / 2) Y2 = h * fy(x + X1 / 2, y + Y1 / 2, z + Z1 / 2) Z2 = h * fz(x + X1 / 2, y + Y1 / 2, z + Z1 / 2) X3 = h * fx(x + X2 / 2, y + Y2 / 2, z + Z2 / 2) Y3 = h * fy(x + X2 / 2, y + Y2 / 2, z + Z2 / 2) Z3 = h * fz(x + X2 / 2, y + Y2 / 2, z + Z2 / 2) X4 = h * fx(x + X3, y + Y3, z + Z3) Y4 = h * fy(x + X3, y + Y3, z + Z3) Z4 = h * fz(x + X3, y + Y3, z + Z3) xa = x ya = y za = z x = x + (X1 + 2 * X2 + 2 * X3 + X4) / 6 y = y + (Y1 + 2 * Y2 + 2 * Y3 + Y4) / 6 z = z + (Z1 + 2 * Z2 + 2 * Z3 + Z4) / 6 If n 2000 Then Picture1.