如何在Matplotlib中的3D颤振器中旋转箭头?

时间:2019-09-07 09:10:22

标签: python matplotlib plot

我正在尝试使用python和matplotlib复制以下图。 Plot I want to replicate 但是,我能够生产的最好的产品如下 My attempt

这里的主要问题是不是平面箭头,即使我对总体图的质量不满意也是如此。我一直在寻找在3D图中使用2D颤动的解决方案,但是我没有找到有关如何执行此操作的任何有用信息。您是否知道其他实现平面箭头的方法?总的来说,任何有关如何改善该地块质量的建议都会受到欢迎。


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

params = {
   'font.family' : 'serif',
   'mathtext.fontset': 'stix',
   'axes.labelsize': 13,
   'legend.fontsize': 8,
   'xtick.labelsize': 13,
   'ytick.labelsize': 13,
   'text.usetex': True,
   'figure.figsize': [10, 5]
   }


plt.rcParams.update(params)

plt.close('all')


x_ax = np.linspace(-10, 10, 24)
y_ax = np.linspace(-10, 10, 24)

x, y = np.meshgrid(x_ax, y_ax, indexing='ij')

r = np.sqrt(x**2 + y**2)

j_x = -y/r*(- np.exp(-np.abs(r)) + np.exp(-np.abs(r)/2) )*2
j_y = +x/r*(- np.exp(-np.abs(r)) + np.exp(-np.abs(r)/2) )*2

#c = np.arctan2(x, -y)
c = np.sqrt(j_x**2 + j_y**2)
c = (c.ravel() - c.min()) / c.ptp()
c = np.concatenate((c, np.repeat(c, 2)))
c = cm.jet(c)
#c = plt.cm.hsv(c)

fig = plt.figure()
ax = fig.gca(projection='3d')

ax.quiver(x, y, 0, j_x, j_y, 0, colors=c, length=1.2, pivot='middle')


t = np.linspace(-10, 10, 200)

psi = 1 - np.exp(-np.abs(t))
b = np.exp(-t**2)
j_abs = np.abs(t)*np.exp(-t**2)*2
#j_abs = (- np.exp(-np.abs(t)) + np.exp(-np.abs(t)/2) )*2

ax.plot(t, psi, zs=0, zdir='y', label=r"$|\psi|$")
ax.plot(t, b, zs=0, zdir='y', label=r"$|\vec B|$")
ax.plot(t, j_abs, zs=0, zdir='y', label=r"$|\vec j|$")

ax.legend()

ax.set_proj_type('ortho')
ax.set_axis_off()

ax.set_zlim([-0.2, 1.4])
ax.view_init(elev=45, azim=90)
ax.dist=5
fig.savefig("vortex.pdf", bbox_inches="tight")

0 个答案:

没有答案