import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
defheart_3d(x,y,z):return(x**2+(9/4)*y**2+z**2-1)**3-x**2*z**3-(9/80)*y**2*z**3defplot_implicit(fn,bbox=(-1.5,1.5)):
xmin,xmax,ymin,ymax,zmin,zmax = bbox*3
fig = plt.figure('HEART')
ax = fig.add_subplot(111,projection ='3d')
A = np.linspace(xmin,xmax,80)
B = np.linspace(xmin,xmax,30)
A1,A2 = np.meshgrid(A,A)for z in B:
X,Y = A1,A2
Z = fn(X,Y,z)
cest = ax.contour(X,Y,z+Z,[z],zdir='z',colors=('r',))for y in B:
X,Z = A1,A2
Y = fn(X,y,Z)
cest = ax.contour(X,Y+y,Z,[y],zdir ='y',colors =('red',))for x in B :
Y,Z=A1,A2
X = fn(x,Y,Z)
cest = ax.contour(X+x,Y,Z,[x],zdir ='x',colors =('red',))
ax.set_zlim3d(zmin,zmax)
ax.set_xlim3d(xmin,xmax)
ax.set_ylim3d(ymin,ymax)
plt.show()if __name__=='__main__':
plot_implicit(heart_3d)
print('\n'.join([''.join([('ForChange'[(x-y)%8]if((x*0.05)**2+(y*0.1)**2-1)**3-(x*0.05)**2*(y*0.1)**3<=0else' ')for x inrange(-30,30)])for y inrange(15,-15,-1)]))
效果
love-you
代码
import matplotlib.pyplot as plt
import seaborn
import numpy
l = numpy.arange(0,4,0.01)
L =1.0/ l
theta = numpy.arange(-4,4,0.01)
o =3.0* numpy.cos(theta)
O =3.0* numpy.sin(theta)
v = numpy.arange(-4,4,0.01)
V = numpy.abs(-2.0* v)
e = numpy.arange(-3,3,0.01)
E =-1.0* numpy.abs(numpy.sin(e))
y = numpy.arange(-10,10,0.01)
Y = numpy.log2(numpy.abs(y))
u = numpy.arange(-4,4,0.01)
U =2.0* u **2
points =[]for heartY in numpy.linspace(-100,100,500):for heartX in numpy.linspace(-100,100,500):if((heartX*0.03)**2+(heartY*0.03)**2-1)**3-(heartX*0.03)**2*(heartY*0.03)**3<=0:
points.append({"x": heartX,"y": heartY})
heart_x =list(map(lambda point: point["x"], points))
heart_y =list(map(lambda point: point["y"], points))
fig = plt.figure(figsize=(13,7))
ax_L = fig.add_subplot(2,4,1)
ax_O = fig.add_subplot(2,4,2)
ax_V = fig.add_subplot(2,4,3)
ax_E = fig.add_subplot(2,4,4)
ax_Y = fig.add_subplot(2,4,5)
ax_O_2 = fig.add_subplot(2,4,6)
ax_U = fig.add_subplot(2,4,7)
ax_heart = fig.add_subplot(2,4,8)
plt.plot(colos='tomato')
ax_L.plot(l, L)
ax_O.plot(o, O)
ax_V.plot(v, V)
ax_E.plot(E, e)
ax_Y.plot(y, Y)
ax_Y.axis([-10.0,10.0,-10.0,5.0])
ax_O_2.plot(o, O)
ax_U.plot(u, U)
ax_heart.scatter(heart_x, heart_y, s=10, alpha=0.5)
seaborn.set_style('whitegrid')
plt.show()