Python绘制牛奶冻曲线(高木曲线)案例_F11 - 专业站长和开发者的学习网站

分享到

Python绘制牛奶冻曲线(高木曲线)案例

python 来源：互联网作者：佚名发布时间：2022-08-26 20:06:22 人浏览

摘要

牛奶冻曲线(blancmange curve)，因在1901年由高木贞治所研究，又称高木曲线。在单位区间内，牛奶冻函数定义为：分形曲线的轮廓会随着阶数的增多而填充细节，即对于下面的来说， N的

牛奶冻曲线(blancmange curve)，因在1901年由高木贞治所研究，又称高木曲线。

在单位区间内，牛奶冻函数定义为：

分形曲线的轮廓会随着阶数的增多而填充细节，即对于下面的来说， 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()

如图所示：

牛奶冻曲线是一种典型的分形曲线，即随着区间的不断缩小，其形状几乎不发生什么变化，例如更改自变量的范围，令

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

最终得到的牛奶冻曲线在观感上是没什么区别的。

接下来绘制一下，当区间发生变化时，牛奶冻曲线的变化过程

绘图代码为：

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()

您可能感兴趣的文章 :

原文链接 : https://blog.csdn.net/m0_37816922/article/details/126517261

Tag : python(762)曲线(4)

Python Django教程之实现新闻应用程序

Django是一个用Python编写的高级框架，它允许我们创建服务器端Web应用程序。在本文中，我们将了解如何使用Django创建新闻应用程序。我们将
书写Python代码的一种更优雅方式(推荐!)

一些比较熟悉pandas的读者朋友应该经常会使用query()、eval()、pipe()、assign()等pandas的常用方法，书写可读性很高的「链式」数据分析处理代码
Python灰度变换中伽马变换分析实现

1. 介绍伽马变换主要目的是对比度拉伸，将图像灰度较低的部分进行修正伽马变换针对的是对单个像素点的变换，也就是点对点的映射形
使用OpenCV实现迷宫解密的全过程

一、你能自己走出迷宫吗？如下图所示，可以看到是一张较为复杂的迷宫图，相信也有人尝试过自己一点一点的找出口，但我们肉眼来解谜
Python中的数据精度问题的介绍

一、python运算时精度问题 1.运行时精度问题在Python中（其他语言中也存在这个问题，这是计算机采用二进制导致的），有时候由于二进制和
Python随机值生成的常用方法

一、随机整数 1.包含上下限：[a, b] 1 2 3 4 import random #1、随机整数：包含上下限：[a, b] for i in range(10): print(random.randint(0,5),end= | ) 查看运行结
Python字典高级用法深入分析讲解

一、 collections 中 defaultdict 的使用 1.字典的键映射多个值将下面的列表转成字典 l = [(a,2),(b,3),(a,1),(b,4),(a,3),(a,1),(b,3)] 一个字典就是一个键对
Python浅析多态与鸭子类型使用实例

什么多态：同一事物有多种形态为何要有多态=》多态会带来什么样的特性，多态性多态性指的是可以在不考虑对象具体类型的情况下而直
Python字典高级用法深入分析介绍

一、 collections 中 defaultdict 的使用 1.字典的键映射多个值将下面的列表转成字典 l = [(a,2),(b,3),(a,1),(b,4),(a,3),(a,1),(b,3)] 一个字典就是一个键对
Python淘宝或京东等秒杀抢购脚本实现(秒杀脚本

我们的目标是秒杀淘宝或京东等的订单，这里面有几个关键点，首先需要登录淘宝或京东，其次你需要准备好订单，最后要在指定时间快速