Board logo

标题: [数值计算] 求批处理版Mandelbrot集合 [打印本页]
作者: email10t    时间: 2015-7-14 07:09     标题: 求批处理版Mandelbrot集合

求批处理版Mandelbrot集合
作者: aa77dd@163.com    时间: 2015-7-14 15:15

本帖最后由 aa77dd@163.com 于 2015-7-14 15:42 编辑

近似算法, 以下是迭代上限 MaxRepeats 设为 30 次得到的图像
  1.                                                                       # #

  2.                                                                     #

  3.                                                                     #

  4.                                                                     #

  5.                                                                    ###

  6.                                                                ## #####

  7.                                                                 ########

  8.                                                                 #########

  9.                                                                 #########

  10.                                                                 #########

  11.                                                                 ########

  12.                                                   #       #      #######

  13.                                                     #     ## ############## ###      #

  14.                                                     ###   ####################    #

  15.                                                     #### ##################### #  # ##

  16.                                                    ##################################

  17.                                                     #################################

  18.                                                      ################################

  19.                                                 #  # ################################

  20.                                                #  ###################################

  21.                                                 ######################################

  22.                                                 #####################################

  23.                           #                     #####################################

  24.                                   #          ###########################################

  25.                                    #            ########################################

  26.                                  ###          ########################################

  27.                              ## #######       #########################################

  28.                              ############     #########################################

  29.                             #############     ##########################################

  30.                            ################# ##########################################

  31.                            ################# ##########################################

  32.                           ##############################################################

  33.                        #   ###########################################################

  34.                       ################################################################

  35.                       ###############################################################

  36.         ##        ##################################################################

  37. ##################################################################################

  38.         ##        ##################################################################

  39.                       ###############################################################

  40.                       ################################################################

  41.                        #   ###########################################################

  42.                           ##############################################################

  43.                            ################# ##########################################

  44.                            ################# ##########################################

  45.                             #############     ##########################################

  46.                              ############     #########################################

  47.                              ## #######       #########################################

  48.                                  ###          ########################################

  49.                                    #            ########################################

  50.                                   #          ###########################################

  51.                           #                     #####################################

  52.                                                 #####################################

  53.                                                 ######################################

  54.                                                #  ###################################

  55.                                                 #  # ################################

  56.                                                      ################################

  57.                                                     #################################

  58.                                                    ##################################

  59.                                                     #### ##################### #  # ##

  60.                                                     ###   ####################    #

  61.                                                     #     ## ############## ###      #

  62.                                                   #       #      #######

  63.                                                                 ########

  64.                                                                 #########

  65.                                                                 #########

  66.                                                                 #########

  67.                                                                 ########

  68.                                                                ## #####

  69.                                                                    ###

  70.                                                                     #

  71.                                                                     #

  72.                                                                     #
复制代码
如果将 MaxRepeats 的值设小一些, 耗时会减短, 但结果的近似程度相应会更差
  1. @echo off & setlocal enabledelayedexpansion & color f0 & chcp 437 & title Mandelbrot set

  2. rem 将控制台字体设置为 点阵 8X8 查看效果比较好

  3. 

  4. set "DOT=#"

  5. set /a "MaxRepeats = 30"

  6. set /a "ratio=36, EscapeSQ=ratio*ratio*2*2+1, wid=ratio*3, Cols=wid, lines=ratio*2+1, iMax=wid*(lines-1)"

  7. set /a "XC = Cols*2/3, YC = lines/2"

  8. set /a "xl=ratio*-2, xh=ratio-1, yl=-ratio+1, yh=ratio"

  9. mode !Cols!, !lines!

  10. 

  11. for /l %%i in (1 1 !iMax!) do set "scr= !scr!"

  12. 

  13. for /l %%a in (!xl! 1 !xh!) do for /l %%b in (!yl! 1 !yh!) do call :Mandelbrot %%a %%b

  14. cls &(<nul set /p=!scr:~0,-1!)

  15. >nul pause

  16. exit

  17. 

  18. 

  19. :Mandelbrot a b  rem c = a + bi

  20. title %1,%2

  21. set /a "repeats=za=zb=0" & rem z = za + zbi

  22. 

  23. :loop

  24.     REM z = z^2 + c

  25.     set /a "ta = za, za = (za*za - zb*zb) / ratio + %1,  zb = 2 * ta * zb / ratio + %2"

  26.     set /a "repeats += 1"

  27.     set /a "SQ = za*za + zb*zb, t = SQ - EscapeSQ & repeats - MaxRepeats"

  28. if !t! lss 0 goto :loop

  29. 

  30. if !repeats! geq %MaxRepeats% (

  31.     set /a "x= XC + %1, y= YC - %2"

  32.     set /a "ind=x+y*wid+1, lenL=ind-1, lenR=iMax-ind"

  33.     for /f "tokens=1,2,3" %%a in ("!lenL! !ind! !lenR!") do (set scr=!scr:~0,%%a!%DOT%!scr:~%%b,%%c!)

  34.     cls &(<nul set /p=!scr:~0,-1!)

  35. )

  36. exit /b
复制代码
作者: email10t    时间: 2015-8-9 07:40

http://pan.baidu.com/s/1c02jn2o



欢迎光临 批处理之家 (http://bbs.bathome.net/) Powered by Discuz! 7.2