Python繪制牛奶凍曲線(xiàn)(高木曲線(xiàn))案例
前言:
牛奶凍曲線(xiàn)(blancmange curve),因在1901年由高木貞治所研究,又稱(chēng)高木曲線(xiàn)。
在單位區(qū)間內(nèi),牛奶凍函數(shù)定義為:
分形曲線(xiàn)的輪廓會(huì)隨著階數(shù)的增多而填充細(xì)節(jié),即對(duì)于下面的來(lái)說(shuō), N的變化會(huì)增添曲線(xiàn)的自相似特性
import numpy as np import matplotlib.pyplot as plt s = lambda x : np.min([x-np.floor(x), np.ceil(x)-x],0) x = np.arange(1000).reshape(-1,1)/1000 N = np.arange(30).reshape(1,-1) #2^N已經(jīng)很大了,精度足夠 b = np.sum(s(2**N*x)/2**N,1) plt.plot(b) plt.show()
如圖所示:
牛奶凍曲線(xiàn)是一種典型的分形曲線(xiàn),即隨著區(qū)間的不斷縮小,其形狀幾乎不發(fā)生什么變化,例如更改自變量的范圍,令
x = np.arange(0.25,0.5,1e-3).reshape(-1,1)
最終得到的牛奶凍曲線(xiàn)在觀(guān)感上是沒(méi)什么區(qū)別的。
接下來(lái)繪制一下,當(dāng)區(qū)間發(fā)生變化時(shí),牛奶凍曲線(xiàn)的變化過(guò)程
繪圖代碼為:
from aniDraw import * # 三角波函數(shù) s = lambda x : min(np.ceil(x)-x, x-np.floor(x)) s = lambda x : np.min([x-np.floor(x), np.ceil(x)-x],0) x = np.arange(1000).reshape(-1,1)/1000 N = np.arange(30).reshape(1,-1) #2^N已經(jīng)很大了,精度足夠 b = np.sum(s(2**N*x)/2**N,1) fig = plt.figure(figsize=(12,8)) ax = fig.add_subplot() # n為坐標(biāo)軸參數(shù) def bcFunc(n): st = 1/3 - (1/3)**n ed = 1/3 + (2/3)**n x = np.linspace(st,ed,1000).reshape(-1,1) b = np.sum(s(2**N*x)/2**N,1) return (x,b) line, = ax.plot([],[],lw=1) def animate(n): x,y = bcFunc(n) line.set_data(x,y) plt.xlim(x[0],x[-1]) plt.ylim(np.min(y),np.max(y)) return line, Ns = np.arange(1,10,0.1) ani = animation.FuncAnimation(fig, animate, Ns, interval=125, blit=False) plt.show()
到此這篇關(guān)于Python繪制牛奶凍曲線(xiàn)(高木曲線(xiàn))案例的文章就介紹到這了,更多相關(guān)Python 牛奶凍曲線(xiàn)內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
Keras搭建孿生神經(jīng)網(wǎng)絡(luò)Siamese?network比較圖片相似性
這篇文章主要為大家介紹了Keras搭建孿生神經(jīng)網(wǎng)絡(luò)Siamese?network比較圖片相似性,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進(jìn)步,早日升職加薪2022-05-05Python+Tkinter創(chuàng)建一個(gè)簡(jiǎn)單的鬧鐘程序
這篇文章主要為大家詳細(xì)介紹了如何使用 Python 的 Tkinter 庫(kù)創(chuàng)建一個(gè)簡(jiǎn)單的鬧鐘程序,它可以在指定的時(shí)間播放一個(gè)聲音來(lái)提醒你,感興趣的可以學(xué)習(xí)一下2023-04-04Python實(shí)現(xiàn)E-Mail收集插件實(shí)例教程
這篇文章主要給大家介紹了關(guān)于Python實(shí)現(xiàn)E-Mail收集插件的相關(guān)資料,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面來(lái)一起看看吧2019-02-02Linux下使用python調(diào)用top命令獲得CPU利用率
這篇文章主要介紹了Linux下使用python調(diào)用top命令獲得CPU利用率,本文直接給出實(shí)現(xiàn)代碼,需要的朋友可以參考下2015-03-03python實(shí)現(xiàn)爬取千萬(wàn)淘寶商品的方法
這篇文章主要介紹了python實(shí)現(xiàn)爬取千萬(wàn)淘寶商品的方法,涉及Python頁(yè)面抓取的相關(guān)技巧,需要的朋友可以參考下2015-06-06Django 開(kāi)發(fā)環(huán)境配置過(guò)程詳解
這篇文章主要介紹了Django 開(kāi)發(fā)環(huán)境配置過(guò)程詳解,文中通過(guò)示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-07-07Python 實(shí)時(shí)獲取任務(wù)請(qǐng)求對(duì)應(yīng)的Nginx日志的方法
本文給大家分享Python 實(shí)時(shí)獲取任務(wù)請(qǐng)求對(duì)應(yīng)的Nginx日志的方法,本文通過(guò)實(shí)例代碼給大家介紹的非常詳細(xì),需要的朋友參考下吧2021-07-07