神奇的分形藝術(shù): Mandelbrot集和Julia集

這段時(shí)間看了一個(gè)關(guān)于維度的視頻介紹,嘆于其驚艷的多維幾何體和分形的視覺(jué)動(dòng)畫(huà)效果。其實(shí)關(guān)于分形,已經(jīng)有很成熟的分形軟件和應(yīng)用場(chǎng)景,可以參看目前流行的分形軟件一覽,不過(guò)沒(méi)有及時(shí)更新,有些鏈接已經(jīng)進(jìn)不了,還有一個(gè) 分形藝術(shù)網(wǎng)。
關(guān)于分形視頻的制作,有不同的制作辦法,例如,這位仁兄的作品:他的思路是通過(guò)Shader渲染,然后通過(guò)ffmpeg命令把圖片生成mp4視頻;也可以通過(guò)專業(yè)的分形軟件來(lái)做也可以在分形軟件基礎(chǔ)上定制開(kāi)發(fā);本文采用純Python實(shí)現(xiàn),借助MoviePy這個(gè)包來(lái)生成視頻或者gif動(dòng)畫(huà),本文只是提供有興趣的朋友,了解這個(gè)領(lǐng)域,并可以動(dòng)手實(shí)踐體驗(yàn)。
本文討論范圍局限在 Mandelbrot集合和Julia集合。

Mandelbrot集合

如果英語(yǔ)還行,建議去維基百科讀英文 Mandelbrot集合

Mandelbrot集合是一個(gè)復(fù)數(shù)c的集合,c由 z_{n+1}=z_{n}^2 + c公式在z_0=0開(kāi)始迭代而得到。得到的值可以組成一個(gè)數(shù)列,依次為 c,c^2+c,(c^2+c)^2+c……。當(dāng)該數(shù)列發(fā)散到無(wú)窮時(shí),對(duì)應(yīng)的點(diǎn)就屬于Mandelbrot集合。Mandelbrot集合是分形中最經(jīng)典例子。

c=0 時(shí),顯然數(shù)列永遠(yuǎn)是0,并不發(fā)散,因此 c=0 不屬于Mandelbrot集合。

又如c=3i 時(shí),對(duì)應(yīng)的數(shù)列為3i, -9+3i, 63-51i, 1431-6477j …. ,數(shù)字越來(lái)越龐大,因此3i就屬于Mandelbrot集合。

在二維平面上,將所有不屬于Mandelbrot集合的點(diǎn)標(biāo)記為黑色,將所有屬于Mandelbrot集合的點(diǎn)按照其發(fā)散速度賦予不同的顏色,就可以得到Mandelbrot的經(jīng)典圖像:


在這里插入圖片描述

請(qǐng)注意: Mandelbrot集合是在z_0=0時(shí),不斷的迭代c得到。

在生成的Mandelbrot集合中,我們可以將圖像放大,選取某些區(qū)域進(jìn)行生成,就可以得到格式各樣造型迥異的分形圖案。在Mandelbrot集合中,有很多地方圖案比較奇特,如下圖中的3個(gè)位置。


在這里插入圖片描述

其中編號(hào)為2的地方被稱為“Elephant Valley”,因?yàn)榇颂幍膱D案與大象很像,直接運(yùn)行可以得到該區(qū)域的圖像:

 # Elephant Valley
 fractal.gen_mandelbrot(fractal.set_data(x_tuple=(.275, .28),
                                     y_tuple=(.006, .01)),
                        rgb=(.9, .6, .6)).save("mandelbrot_area.png")
大象山谷

編號(hào)為3的地方被稱為“Triple Spiral Valley”(三重螺旋),在mandelbrot_area.py修改一下坐標(biāo)位置為(ratio調(diào)整的是顏色):

# 三重螺旋
fractal.gen_mandelbrot(fractal.set_data(x_tuple=(-.09, -.086),
                                        y_tuple=(.654, .657)),
                       rgb=(.2, .6, .6)).save("mandelbrot_3.png")
   

就可以得到該處的圖案:


三重螺旋

最后編號(hào)為1的地方被稱為“Seahorse Valley”(海馬山谷),對(duì)應(yīng)的坐標(biāo)為:

 # Seahorse Valley(海馬山谷)
  fractal.gen_mandelbrot(fractal.set_data(x_tuple=(-.75, -.747),
                                          y_tuple=(.099, .102)),
                          rgb=(.1, .1, .3)).save("mandelbrot_sea.png")

圖像如下,確實(shí)和海馬有一點(diǎn)神似:


海馬山谷
Julia集合

如果英語(yǔ)還行,建議去維基百科讀英文 Julia集合

Julia集合和Mandelbrot集合差不多,但這次我們固定c,轉(zhuǎn)而計(jì)算發(fā)散的z的值。即c是固定的常數(shù)(可以任取),數(shù)列變成z,z^2+c,(z^2+c)^2 +c,…..。如果該數(shù)列發(fā)散,對(duì)應(yīng)的z集合就屬于Julia集合。我們看看不同的c值下Julia圖案的差別:

c = -0.835-0.2321i 時(shí):

在這里插入圖片描述

c = 0.285 - 0.01i 時(shí)::

在這里插入圖片描述

c = -0.8 + 0.156i,圖案又變得完全不同:

在這里插入圖片描述

c = -0.835 - 0.2321i,圖案又變得完全不同:

在這里插入圖片描述

c = 0.70176 - 0.3842i,圖案成為樹(shù)狀:

分形樹(shù)

生成Julia集合的動(dòng)畫(huà)

在Julia集合中,每次都對(duì)c的值進(jìn)行微小的改變,并將依次生成圖片制作為gif,就可以生成如下所示的動(dòng)畫(huà): 從這里下載觀看,
文件超過(guò)尺寸了。程序中提供了一個(gè)width參數(shù),可以修改它以生成更大尺寸,質(zhì)量更高的動(dòng)畫(huà)圖像。

代碼實(shí)現(xiàn)

代碼采用TensorFlow實(shí)現(xiàn), 比較簡(jiǎn)單,可以從這里下載代碼。希望你能點(diǎn)贊,反饋,謝謝!

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書(shū)系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • 令人目眩的萬(wàn)花筒,螺旋紋路的西蘭花,它們之間存在什么相似之處? 我們說(shuō)“一花一世界,一樹(shù)一菩提”,說(shuō)的是以小見(jiàn)大,...
    牛油果進(jìn)化論閱讀 4,150評(píng)論 0 16
  • 20170301 協(xié)方差矩陣的理解 協(xié)方差矩陣的幾何解釋(Very Good) [LDA 線性判別分析](1)(2...
    NoneLand閱讀 474評(píng)論 0 1
  • 早上好,對(duì)自己打個(gè)招呼 真勇敢,太棒了,對(duì)自己豎起大拇指 能量滿滿的遇見(jiàn)最好的自己 在頹廢的時(shí)候,在梨花帶雨痛哭流...
    輪滑公主閱讀 197評(píng)論 0 0
  • 上個(gè)月反反復(fù)復(fù)的看這幾本書(shū),看到激動(dòng)處會(huì)立即起來(lái)扔掉一樣糾結(jié)很久依然不舍得扔的物品,心中涌起想扔?xùn)|西的欲望。因?yàn)閻?ài)...
    榮榮手工閱讀 927評(píng)論 5 10
  • 常會(huì)杰就是這樣一個(gè)嫁給了愛(ài)情的人,她和丈夫付偉超是在大學(xué)時(shí)期自由戀愛(ài)的。大學(xué)畢業(yè)后,付偉超分配到了當(dāng)?shù)氐男W(xué)當(dāng)老師...
    bright_e596閱讀 386評(píng)論 0 0

友情鏈接更多精彩內(nèi)容