使用python 的matplotlib 畫(huà)軌道實(shí)例
更新時(shí)間：2020年01月19日 09:29:42 作者：cq_pf
今天小編就為大家分享一篇使用python 的matplotlib 畫(huà)軌道實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
如下所示：
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.patches as mpatches
from scipy import stats
fig = plt.figure()
ax = fig.add_subplot(111, xlim=(0, 10), ylim=(-4, 4))
sx=0;sy=0;r=1.5 ;
circle = mpatches.Circle((sx,sy),r,ec='b',fc = 'b',alpha=0.6)
ax.add_patch(circle)
n = 5
alpha = np.linspace(-np.pi/2,np.pi/2,100)
for i in range(1,n+1):
a = (i+1)*2
b = (i+1)
plt.plot(a*np.cos(alpha),b*np.sin(alpha))
if i <= 3:
ax.text(a+1,0,'layer.%d'%(i-1),ha='center',va='center')
rate = 3
Jnum = 100
for i in range(Jnum):
data = stats.poisson.rvs(mu = rate, loc = 0 , size = 1)
while data == 0:
data = stats.poisson.rvs(mu = rate, loc = 0 , size = 1)
data = data[0]
r = np.random.random()
beta = (np.random.random()-0.5)*np.pi
print beta
data += r
a = (data+1)*2*np.cos(beta)
b = (data+1)*np.sin(beta)
ax.scatter(a,b, c='y',edgecolors = 'y')
xx = [3.5 , 6.8,6.8]
yy = [-1.8,0.9,3.1]
pp = [1.00,0.85,0.75]
rota = [-45,0,0]
for i in range(3):
plt.plot((1.5,xx[i]) , (0,yy[i]) , 'r')
ax.scatter(xx[i],yy[i], c='y',edgecolors = 'b')
ax.text(xx[i] , yy[i] , 'probability %lf'%pp[i] , ha = 'center' )
tx = np.linspace(0,10,10)
ty = np.linspace(-4,-4,10)
ax.set_xticks(tx)
ax.set_xticklabels((''))
ax.set_yticks(ty)
ax.set_yticklabels((''))
#plt.axis('off')
plt.show()
以上這篇使用python 的matplotlib 畫(huà)軌道實(shí)例就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
相關(guān)文章
深入理解Python的import機(jī)制有助于更好地組織代碼、提高代碼復(fù)用性,本文將深入研究Python的Import機(jī)制,包括模塊的導(dǎo)入過(guò)程、命名空間與作用域、相對(duì)導(dǎo)入以及包的結(jié)構(gòu)和導(dǎo)入等方面,通過(guò)豐富的示例代碼,助你更全面地理解和應(yīng)用這2024-01-01
python構(gòu)建深度神經(jīng)網(wǎng)絡(luò)（DNN）
這篇文章主要為大家詳細(xì)介紹了python構(gòu)建深度神經(jīng)網(wǎng)絡(luò)DNN，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-03-03
python爬蟲(chóng)之selenium庫(kù)的安裝及使用教程
今天帶大家來(lái)學(xué)習(xí)怎么安裝使用selenium庫(kù),文中有非常詳細(xì)的圖文介紹,對(duì)正在學(xué)習(xí)python爬蟲(chóng)的小伙伴們很有幫助,需要的朋友可以參考下2021-05-05
python統(tǒng)計(jì)mysql數(shù)據(jù)量變化并調(diào)用接口告警的示例代碼
這篇文章主要介紹了python統(tǒng)計(jì)mysql數(shù)據(jù)量變化并調(diào)用接口告警的示例代碼，幫助大家更好的利用python操作數(shù)據(jù)庫(kù)，感興趣的朋友可以了解下2020-09-09
分享7個(gè) Python 實(shí)戰(zhàn)項(xiàng)目練習(xí)
這篇文章主要介紹了分享7個(gè) Python 實(shí)戰(zhàn)項(xiàng)目代碼,經(jīng)過(guò)Python3.6.4調(diào)試通過(guò)的代碼，就具一點(diǎn)的參考價(jià)值，需要的小伙伴可以參考一下2022-03-03
Python實(shí)現(xiàn)常見(jiàn)坐標(biāo)系的相互轉(zhuǎn)換
WGS84坐標(biāo)系、GCJ02坐標(biāo)系、BD09坐標(biāo)系和Web?墨卡托投影坐標(biāo)系是我們常見(jiàn)的四個(gè)坐標(biāo)系。這篇文章為大家整理了這四個(gè)坐標(biāo)系之間相互轉(zhuǎn)換的方法，需要的可以參考一下2023-02-02
python代碼打印100-999之間的回文數(shù)示例
今天小編就為大家分享一篇python代碼打印100-999之間的回文數(shù)示例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-11-11
Django 大文件下載實(shí)現(xiàn)過(guò)程解析
這篇文章主要介紹了Django 大文件下載實(shí)現(xiàn)過(guò)程解析,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-08-08
Python實(shí)現(xiàn)爬取騰訊招聘網(wǎng)崗位信息
這篇文章主要介紹了如何用python爬取騰訊招聘網(wǎng)崗位信息保存到表格，并做成簡(jiǎn)單可視化。文中的示例代碼對(duì)學(xué)習(xí)Python有一定的幫助，感興趣的可以了解一下2022-01-01
Python3.5 Json與pickle實(shí)現(xiàn)數(shù)據(jù)序列化與反序列化操作示例
這篇文章主要介紹了Python3.5 Json與pickle實(shí)現(xiàn)數(shù)據(jù)序列化與反序列化操作,結(jié)合實(shí)例形式分析了Python3.5使用Json與pickle模塊實(shí)現(xiàn)json格式數(shù)據(jù)的序列化及反序列化操作相關(guān)步驟與注意事項(xiàng),需要的朋友可以參考下2019-04-04