Question

我想想象VM相对于时间t的大小。在以下程序中，t设置为零。如果有人能告诉我如何制作双线性插值等高线图的动画，我将不胜感激。谢谢！请在这里输入代码：

import numpy as np
from numpy import cos, sin, pi, exp, sqrt
from matplotlib import pyplot as plt
import matplotlib.cm as cm

nx=30
ny=30

nu=0.0025
p=2*pi

x = np.linspace(0,p,nx)
y = np.linspace(0,p,ny)

fig,ax = plt.subplots(1,1)
ax.set_xlim(0,p)
ax.set_ylim(0,p)

t=0
X,Y = np.meshgrid(x,y)
U=-cos(X)*sin(Y)*exp(-2*nu*t)
V=sin(X)*cos(Y)*exp(-2*nu*t)
VM=sqrt(U**2+V**2)

im = plt.imshow(VM, interpolation='bilinear', origin='lower',
            cmap=cm.gray, extent=(0, p, 0, p))
levels = np.arange(-1.2, 1.6, 0.2)

# add a vertical colorbar.
CBI = plt.colorbar(im, orientation='vertical', shrink=0.8)
plt.show()

Answer 1

主要是从here复制的。在那里查看深入的解释。

import numpy as np
from numpy import cos, sin, pi, exp, sqrt
from matplotlib import pyplot as plt
import matplotlib.cm as cm
from matplotlib import animation

nx=30
ny=30

nu=0.0025
p=2*pi

x = np.linspace(0,p,nx)
y = np.linspace(0,p,ny)

fig,ax = plt.subplots(1,1)
ax.set_xlim(0,p)
ax.set_ylim(0,p)

t=0
X,Y = np.meshgrid(x,y)
U=-cos(X)*sin(Y)*exp(-2*nu*t)
V=sin(X)*cos(Y)*exp(-2*nu*t)
VM=sqrt(U**2+V**2)

im = plt.imshow(VM, interpolation='bilinear', origin='lower',
            cmap=cm.gray, extent=(0, p, 0, p))
levels = np.arange(-1.2, 1.6, 0.2)

# add a vertical colorbar.
CBI = plt.colorbar(im, orientation='vertical', shrink=0.8)


def init():
    im.set_data((np.empty_like(VM)))
    return im,

def animate(t):
    U=-cos(X)*sin(Y)*exp(-2*nu*t)
    V=sin(X)*cos(Y)*exp(-2*nu*t)
    VM=sqrt(U**2+V**2)
    im.set_data(VM)
    return im,

# # call the animator.  blit=True means only re-draw the parts that have changed.
anim = animation.FuncAnimation(fig, animate, init_func=init,
                               frames=200, interval=20, blit=True)

plt.show()

如何为双线性插值产生的等高线图设置动画？

1 个答案: