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

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

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