通過(guò)Python繪制九種二次曲面的示例詳解
二次曲面
python中繪制三維圖需要將坐標(biāo)系聲明為3d。
球面方程為
x2+y2+z2=R2
寫(xiě)為極坐標(biāo)形式為
x=Rsinθcosφ
y=Rsinθsinφ
z=Rcosθ
令R = 1 R=1R=1，則畫(huà)圖為
代碼如下
>>> import matplotlib.pyplot as plt >>> import numpy as np >>> theta = np.arange(0,6.4,0.1).reshape(64,1) >>> phi = np.arange(0,3.2,0.1).reshape(1,32) >>> x = np.sin(theta)*np.cos(phi) >>> y = np.sin(theta)*np.sin(phi) >>> z = np.cos(theta) >>> ax = plt.gca(projection='3d') >>> ax.plot_surface(x,y,z) <mpl_toolkits.mplot3d.art3d.Poly3DCollection object at 0x000001CECF13A730> >>> plt.show()
二次曲面共有九種，代碼均與橢球曲面類(lèi)似，為了加強(qiáng)立體感，可在畫(huà)圖的時(shí)候設(shè)置顏色映射，下列各圖部分用到
from matplotlib import cm #... ax.plot_surface(x,y,z,cmap=cm.coolwarm)
a,b,c均為1時(shí)的曲面
1.橢圓錐面 x2/a2+y2/b2-z2/c2=0
2.橢球面 x2/a2+y2/b2+z2/c2=1
3.單葉雙曲面 x2/a2+y2/b2-z2/c2=1
4.雙葉雙曲面 x2/a2+y2/b2-z2/c2=-1
5.橢圓拋物面 z=x2/a2+y2/b2
6.雙曲拋物面 z=x2/a2-y2/b2
7.橢圓柱面 x2/a2+y2/b2=1
8.雙曲柱面 x2/a2-y2/b2=1
9.拋物柱面 y2=2px
在上面各式中，橢圓錐面、單葉雙曲面、雙葉雙曲面具有極為相似的表達(dá)式
故可繪制動(dòng)態(tài)圖來(lái)表示這一過(guò)程，由于animation中無(wú)法繪制plot_surface，所以采用將單張圖片生成gif的方式。
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm
import imageio
theta = np.arange(0,6.4,0.1)
z = np.arange(-2,2,0.02).reshape(200,1)
gifImgs = []
fig = plt.figure()
for i in np.arange(-1,1,0.02):
theta = np.arange(0,6.4,0.1).reshape(1,64)
Z = np.repeat(z,64).reshape(200,64)
x = np.sqrt(z**2+i)*np.cos(theta)
y = np.sqrt(z**2+i)*np.sin(theta)
ax = plt.gca(projection='3d')
ax.plot_surface(x,y,Z,cmap=cm.coolwarm)
plt.savefig("%.2f.jpg" % i)
gifImgs.append(imageio.imread("%.2f.jpg" % i))
imageio.mimsave("test.gif",gifImgs,fps=5)
到此這篇關(guān)于通過(guò)Python繪制九種二次曲面的示例詳解的文章就介紹到這了,更多相關(guān)Python曲面內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
Elasticsearch py客戶(hù)端庫(kù)安裝及使用方法解析
這篇文章主要介紹了Elasticsearch py客戶(hù)端庫(kù)安裝及使用方法解析,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-09-09
這篇文章主要介紹了Python?list?append方法如何給列表追加元素，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。如有錯(cuò)誤或未考慮完全的地方，望不吝賜教2022-11-11
下面小編就為大家?guī)?lái)一篇使用Turtle畫(huà)正螺旋線(xiàn)的方法。小編覺(jué)得挺不錯(cuò)的，現(xiàn)在就分享給大家，也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧2017-09-09
這篇文章主要介紹了Matlab求解數(shù)組中的最大值及它所在的具體位置，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2021-04-04
Python中關(guān)于元組 集合 字符串 函數(shù) 異常處理的全面詳解
本篇文章介紹了我在學(xué)習(xí)python過(guò)程中對(duì)元組、集合、字符串、函數(shù)、異常處理的總結(jié)，通讀本篇對(duì)大家的學(xué)習(xí)或工作具有一定的價(jià)值，需要的朋友可以參考下2021-10-10
這篇文章主要介紹了Django視圖、傳參和forms驗(yàn)證操作，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-07-07
使用Python開(kāi)發(fā)一個(gè)功能完善的媒體文件截幀工具
在視頻處理和分析工作中,我們經(jīng)常需要從視頻中提取特定幀作為圖片,今天,我將帶大家使用?Python?的?wxPython?和?OpenCV?庫(kù),開(kāi)發(fā)一個(gè)功能完善的媒體文件截幀工具,感興趣的小伙伴跟著小編一起來(lái)看看吧2025-10-10
matlab、python中矩陣的互相導(dǎo)入導(dǎo)出方式
這篇文章主要介紹了matlab、python中矩陣的互相導(dǎo)入導(dǎo)出方式，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2020-06-06
python實(shí)現(xiàn)Flappy Bird源碼
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)Flappy Bird源碼，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-12-12
Python下載網(wǎng)絡(luò)文本數(shù)據(jù)到本地內(nèi)存的四種實(shí)現(xiàn)方法示例
這篇文章主要介紹了Python下載網(wǎng)絡(luò)文本數(shù)據(jù)到本地內(nèi)存的四種實(shí)現(xiàn)方法,涉及Python網(wǎng)絡(luò)傳輸、文本讀寫(xiě)、內(nèi)存I/O、矩陣運(yùn)算等相關(guān)操作技巧,代碼中包含了較為詳盡的注釋說(shuō)明便于理解,需要的朋友可以參考下2018-02-02