情人节,送你一朵小红花~

文章目录

    • 前言
    • 1. 含苞待放
    • 2. 灼灼其华

前言

  在github上看到一个有趣的代码,正值情人节,分享给大家^_^

1. 含苞待放

  3D模型的绘制需要网格点,关于网格点的作用,在这篇博客中已经介绍,这里不再赘述。

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure()
ax = Axes3D(fig)

x = np.linspace(0, 1, num=30)
t = np.linspace(0, 1, num=1200) * 20 * np.pi + 4 * np.pi
x, t = np.meshgrid(x, t)

p = 0.5 * np.pi * np.exp(-t / (8 * np.pi))
change = np.sin(15 * t) / 150
u = 1 - (1 - np.mod(3.3 * t, 2 * np.pi) / np.pi)**4 / 2 + change
y = 2 * (x**2 - x)**2 * np.sin(p)
r = u * (x * np.sin(p) + y * np.cos(p)) * 1.5
h = u * (x * np.cos(p) - y * np.sin(p))

# PiYG_r
# RdBu_r
surf = ax.plot_surface(r * np.cos(t), r * np.sin(t), h,
                       rstride=1, cstride=1, cmap=plt.cm.RdPu_r)
plt.savefig('img/img1.png')
plt.show()

情人节,送你一朵小红花~_第1张图片
情人节,送你一朵小红花~_第2张图片
情人节,送你一朵小红花~_第3张图片

2. 灼灼其华

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure()
ax = Axes3D(fig)
# plt.axis('off')

x = np.linspace(0, 1, num=30)
t = np.linspace(0, 1, num=1200) * 50 * np.pi - 4 * np.pi
x, t = np.meshgrid(x, t)

p = 0.5 * np.pi * np.exp(-t / (8 * np.pi))
change = np.sin(20 * t) / 50
u = 1 - (1 - np.mod(3.3 * t, 2 * np.pi) / np.pi)**4 / 2 + change
y = 2 * (x**2 - x)**2 * np.sin(p)
r = u * (x * np.sin(p) + y * np.cos(p)) * 1.5
h = u * (x * np.cos(p) - y * np.sin(p))

ax = ax.plot_surface(r * np.cos(t), r * np.sin(t), h,
                     rstride=1, cstride=1, cmap=plt.cm.RdPu_r)

plt.savefig('img/img2.png')
plt.show()

情人节,送你一朵小红花~_第4张图片

情人节,送你一朵小红花~_第5张图片
情人节,送你一朵小红花~_第6张图片

  有关mpl_toolkits.mplot3d的使用可以参考官方文档;

  更多的颜色搭配可参考matplotlibcolormap官方手册。

  3D花的绘制原理请参考我的这篇博客。

你可能感兴趣的:(闲了没事,python,算法,线性代数,可视化,matplotlib)