欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

Python繪制牛奶凍曲線(高木曲線)案例

 更新時間:2022年08月26日 17:00:21   作者:微小冷  
這篇文章主要介紹了Python繪制牛奶凍曲線(高木曲線)案例,文章圍繞主題展開詳細的內容介紹,具有一定的參考價值,需要的小伙伴可以參考一下

前言:

牛奶凍曲線(blancmange curve),因在1901年由高木貞治所研究,又稱高木曲線。

在單位區(qū)間內,牛奶凍函數定義為:

分形曲線的輪廓會隨著階數的增多而填充細節(jié),即對于下面的來說, 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已經很大了,精度足夠
b = np.sum(s(2**N*x)/2**N,1)
plt.plot(b)
plt.show()

如圖所示:

牛奶凍曲線是一種典型的分形曲線,即隨著區(qū)間的不斷縮小,其形狀幾乎不發(fā)生什么變化,例如更改自變量的范圍,令

x = np.arange(0.25,0.5,1e-3).reshape(-1,1)

最終得到的牛奶凍曲線在觀感上是沒什么區(qū)別的。

接下來繪制一下,當區(qū)間發(fā)生變化時,牛奶凍曲線的變化過程

繪圖代碼為:

from aniDraw import *

# 三角波函數
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已經很大了,精度足夠
b = np.sum(s(2**N*x)/2**N,1)
fig = plt.figure(figsize=(12,8))
ax = fig.add_subplot()
# n為坐標軸參數
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()

到此這篇關于Python繪制牛奶凍曲線(高木曲線)案例的文章就介紹到這了,更多相關Python 牛奶凍曲線內容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關文章希望大家以后多多支持腳本之家!

相關文章

  • Keras搭建孿生神經網絡Siamese?network比較圖片相似性

    Keras搭建孿生神經網絡Siamese?network比較圖片相似性

    這篇文章主要為大家介紹了Keras搭建孿生神經網絡Siamese?network比較圖片相似性,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進步,早日升職加薪
    2022-05-05
  • Python+Tkinter創(chuàng)建一個簡單的鬧鐘程序

    Python+Tkinter創(chuàng)建一個簡單的鬧鐘程序

    這篇文章主要為大家詳細介紹了如何使用 Python 的 Tkinter 庫創(chuàng)建一個簡單的鬧鐘程序,它可以在指定的時間播放一個聲音來提醒你,感興趣的可以學習一下
    2023-04-04
  • Langchain集成管理prompt功能詳解

    Langchain集成管理prompt功能詳解

    這篇文章主要為大家介紹了Langchain集成管理prompt功能示例詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進步,早日升職加薪
    2023-03-03
  • Python實現E-Mail收集插件實例教程

    Python實現E-Mail收集插件實例教程

    這篇文章主要給大家介紹了關于Python實現E-Mail收集插件的相關資料,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面來一起看看吧
    2019-02-02
  • Linux下使用python調用top命令獲得CPU利用率

    Linux下使用python調用top命令獲得CPU利用率

    這篇文章主要介紹了Linux下使用python調用top命令獲得CPU利用率,本文直接給出實現代碼,需要的朋友可以參考下
    2015-03-03
  • python實現爬取千萬淘寶商品的方法

    python實現爬取千萬淘寶商品的方法

    這篇文章主要介紹了python實現爬取千萬淘寶商品的方法,涉及Python頁面抓取的相關技巧,需要的朋友可以參考下
    2015-06-06
  • pydantic進階用法示例詳解

    pydantic進階用法示例詳解

    這篇文章主要為大家介紹了pydantic進階用法示例詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進步,早日升職加薪
    2023-03-03
  • 用python解壓分析jar包實例

    用python解壓分析jar包實例

    今天小編就為大家分享一篇用python解壓分析jar包實例,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧
    2020-01-01
  • Django 開發(fā)環(huán)境配置過程詳解

    Django 開發(fā)環(huán)境配置過程詳解

    這篇文章主要介紹了Django 開發(fā)環(huán)境配置過程詳解,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下
    2019-07-07
  • Python 實時獲取任務請求對應的Nginx日志的方法

    Python 實時獲取任務請求對應的Nginx日志的方法

    本文給大家分享Python 實時獲取任務請求對應的Nginx日志的方法,本文通過實例代碼給大家介紹的非常詳細,需要的朋友參考下吧
    2021-07-07

最新評論