我知道如何得到一個(gè)圓柱體，但我想描繪一些東西，比如衛(wèi)生紙卷，用圓柱體在表面繪制阿基米德螺線。但我需要的是像情節(jié)一樣的衛(wèi)生紙卷。我想出了這背后的數(shù)學(xué)問題有人可以幫助我使用 python 我需要在 3D 中繪制它以獲得以下等式 在實(shí)踐中公式 我想用 L 作為參數(shù)，我的方程變成這里 h 是金屬厚度，r 是軋輥的內(nèi)半徑。該公式使用螺旋輥的圓形近似值。我也知道長度 L =50 有人可以幫我寫 matplotlib 代碼嗎這正是我需要的 http://pgfplots.net/tikz/examples/ cylinder-spiral/ 有人可以幫我把它放在 Matplotlib 中嗎
1 回答
慕虎7371278
TA貢獻(xiàn)1802條經(jīng)驗(yàn) 獲得超4個(gè)贊
以原點(diǎn)和半徑r為圓心的圓的參數(shù)方程為， x = r \times sin(\theta) y = r \times cos(\theta) 其中\(zhòng)theta \in [0,2\pi]。
對于螺旋，半徑隨著 \$\theta\$ 增加。假設(shè) \$r\$ 依賴，\theta因?yàn)閞 = (a+b\theta)它可以， x = (a+b\theta) sin(\theta) y = (a+b\theta) cos(\theta)
要使其成為具有垂直軸的 3D 圖形，您可以添加圓柱體的長度z在linspace(0, L)哪里L(fēng)。
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import math
import numpy as np
L = 50
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
xpoints=[]
ypoints=[]
a = 0.1
b = 0.1
for theta in np.linspace(0, 2*math.pi, 20):
xpoints.append((a+b*theta)*math.cos(theta))
ypoints.append((a+b*theta)*theta*math.sin(theta))
z = np.linspace(0,L)
theta, z = np.meshgrid(theta, z)
ax.plot_surface(xpoints,ypoints,z)
plt.show()
既然你有一個(gè)工作代碼，你可以把它發(fā)布在代碼審查堆棧交換中，我可以在那里用排版數(shù)學(xué)來解釋。
以下是我以某種方式想出的解決方案，如果有人幫助我改進(jìn)我的解決方案，我會(huì)很高興
L = 50
h= 0.5
r= 5.0[![plot][1]][1]
R = np.sqrt((L*h)/(np.pi)+r**r)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
xpoints=[]
ypoints=[]
for theta in np.linspace(0,20,R):
xpoints.append(1.1*theta*cos(theta))
ypoints.append(1.1*theta*sin(theta))
z = np.linspace(0,R)
theta, z = np.meshgrid(t, z)
ax.plot_surface(xpoints,ypoints,z)
plt.show()
添加回答
舉報(bào)
0/150
提交
取消