绘制沿y轴旋转的旋转体动画Python 3和Matplotlib

4dc9hkyq 于 2023-03-09 发布在 Python

关注(0)|答案(1)|浏览(158)

我有这个修改过的Python代码来自：
Ploting solid of revolution in Python 3 (matplotlib maybe)
但是，我想把它绕y轴旋转，怎么做？
这是代码，边界在y = 1/（x^2 + 4）和x=1，x=4和y=0之间。

import gif
import numpy as np
import matplotlib.pyplot as plt
import mpl_toolkits.mplot3d.axes3d as axes3d

@gif.frame
def plot_volume(angle):
    fig = plt.figure(figsize = (20, 15))

    ax2 = fig.add_subplot(1, 1, 1, projection = '3d')
    angles = np.linspace(0, 360, 20)
    x = np.linspace(1, 4, 60)
    v = np.linspace(0, 2*angle, 60)

    U, V = np.meshgrid(x, v)
    Y1 = (1 / (U**2 + 4 ))*np.cos(V)
    Z1 = (1 / (U**2 + 4 ))*np.sin(V)
    X = U
    ax2.plot_surface(X, Y1, Z1, alpha = 0.2, color = 'blue', rstride = 6, cstride = 6)
    ax2.set_xlim(-2,5)
    ax2.set_ylim(-0.5,0.5)
    ax2.set_zlim(-1,1)
    ax2.view_init(elev = 50, azim = 30*angle)
    ax2.plot_wireframe(X, Y1, Z1, color = 'black')
    ax2._axis3don = False

frames = []
for i in np.linspace(0, 2*np.pi, 20):
    frame = plot_volume(i)
    frames.append(frame)

gif.save(frames, '/home/browni/LasthrimProjection/Python/solidofrevolution1.gif', duration = 500)

matplotlib

来源：https://stackoverflow.com/questions/75633628/plotting-solid-of-revolution-animation-that-is-revolved-toward-y-axis-python-3-a

1条答案

按热度按时间

jv4diomz1#

这是代码，我终于能够弄清楚：

# https://stackoverflow.com/questions/36464982/ploting-solid-of-revolution-in-python-3-matplotlib-maybe
# R is a region in the xy plane bounded by y = (x^2 + 4)^(-1)
# Create the animation of A solid of revolution that is formed by rotating R around the y axis. 
import gif
import numpy as np
import matplotlib.pyplot as plt
import mpl_toolkits.mplot3d.axes3d as axes3d

@gif.frame
def plot_volume(angle):
    fig = plt.figure(figsize = (20, 15))

    ax2 = fig.add_subplot(1, 1, 1, projection = '3d')
    angles = np.linspace(0, 360, 20)
    x = np.linspace(1/20, 1/5, 60)
    v = np.linspace(0, 2*angle, 60)

    U, V = np.meshgrid(x, v)
    Y1 = np.sqrt((1 / U) - 4)*np.cos(V)
    Z1 = np.sqrt((1 / U) - 4)*np.sin(V)
    X = U
    ax2.plot_surface(X, Y1, Z1, alpha = 0.2, color = 'blue', rstride = 6, cstride = 6)
    ax2.set_xlim(-2,2)
    ax2.set_ylim(-4,4)
    ax2.set_zlim(-2,2)
    ax2.view_init(elev = 50, azim = 30*angle)
    ax2.plot_wireframe(X, Y1, Z1, color = 'black')
    ax2._axis3don = False

frames = []
for i in np.linspace(0, 2*np.pi, 20):
    frame = plot_volume(i)
    frames.append(frame)

gif.save(frames, '/home/browni/LasthrimProjection/Python/solidofrevolutionyaxis.gif', duration = 500)

赞(0）回复(0）举报 2023-03-09

我来回答

绘制沿y轴旋转的旋转体动画Python 3和Matplotlib

1条答案

相关问题

热门标签

最新问答