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

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

Mandelbrot集合

如果英語還行伯铣,建議去維基百科讀英文 Mandelbrot集合

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

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

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

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


在這里插入圖片描述

請注意: Mandelbrot集合是在z_0=0時偶器,不斷的迭代c得到。

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


在這里插入圖片描述

其中編號為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")
大象山谷

編號為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")
   

就可以得到該處的圖案:


三重螺旋

最后編號為1的地方被稱為“Seahorse Valley”(海馬山谷),對應(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集合

如果英語還行轰传,建議去維基百科讀英文 Julia集合

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

c = -0.835-0.2321i 時:

在這里插入圖片描述

c = 0.285 - 0.01i 時::

在這里插入圖片描述

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

在這里插入圖片描述

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

在這里插入圖片描述

c = 0.70176 - 0.3842i,圖案成為樹狀:

分形樹

生成Julia集合的動畫

在Julia集合中实蓬,每次都對c的值進(jìn)行微小的改變茸俭,并將依次生成圖片制作為gif,就可以生成如下所示的動畫: 從這里下載觀看安皱,
文件超過尺寸了调鬓。程序中提供了一個width參數(shù),可以修改它以生成更大尺寸酌伊,質(zhì)量更高的動畫圖像袖迎。

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

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

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末托慨,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子暇榴,更是在濱河造成了極大的恐慌厚棵,老刑警劉巖,帶你破解...
    沈念sama閱讀 222,104評論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件蔼紧,死亡現(xiàn)場離奇詭異婆硬,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)奸例,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,816評論 3 399
  • 文/潘曉璐 我一進(jìn)店門彬犯,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人查吊,你說我怎么就攤上這事谐区。” “怎么了逻卖?”我有些...
    開封第一講書人閱讀 168,697評論 0 360
  • 文/不壞的土叔 我叫張陵宋列,是天一觀的道長。 經(jīng)常有香客問我评也,道長炼杖,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 59,836評論 1 298
  • 正文 為了忘掉前任盗迟,我火速辦了婚禮坤邪,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘诈乒。我一直安慰自己罩扇,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,851評論 6 397
  • 文/花漫 我一把揭開白布怕磨。 她就那樣靜靜地躺著喂饥,像睡著了一般。 火紅的嫁衣襯著肌膚如雪肠鲫。 梳的紋絲不亂的頭發(fā)上员帮,一...
    開封第一講書人閱讀 52,441評論 1 310
  • 那天,我揣著相機(jī)與錄音导饲,去河邊找鬼捞高。 笑死氯材,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的硝岗。 我是一名探鬼主播氢哮,決...
    沈念sama閱讀 40,992評論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼型檀!你這毒婦竟也來了冗尤?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,899評論 0 276
  • 序言:老撾萬榮一對情侶失蹤胀溺,失蹤者是張志新(化名)和其女友劉穎裂七,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體仓坞,經(jīng)...
    沈念sama閱讀 46,457評論 1 318
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡背零,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,529評論 3 341
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了无埃。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片徙瓶。...
    茶點(diǎn)故事閱讀 40,664評論 1 352
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖录语,靈堂內(nèi)的尸體忽然破棺而出倍啥,到底是詐尸還是另有隱情禾乘,我是刑警寧澤澎埠,帶...
    沈念sama閱讀 36,346評論 5 350
  • 正文 年R本政府宣布,位于F島的核電站始藕,受9級特大地震影響蒲稳,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜伍派,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 42,025評論 3 334
  • 文/蒙蒙 一江耀、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧诉植,春花似錦祥国、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,511評論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至灼擂,卻和暖如春壁查,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背剔应。 一陣腳步聲響...
    開封第一講書人閱讀 33,611評論 1 272
  • 我被黑心中介騙來泰國打工睡腿, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留语御,地道東北人。 一個月前我還...
    沈念sama閱讀 49,081評論 3 377
  • 正文 我出身青樓席怪,卻偏偏與公主長得像应闯,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子挂捻,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,675評論 2 359

推薦閱讀更多精彩內(nèi)容

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