最新国产好看的视频,伊人天堂AV在线,国产Aaaaaa视频,蜜臀视频在线观看一区,人妻av色图,密臀久久久精品影片,青青视频免费观看毛片,久草在线观看视,国产三级精品色情在线

使用python 的matplotlib 畫軌道實(shí)例

 更新時間:2020年01月19日 09:29:42   作者:cq_pf  
今天小編就為大家分享一篇使用python 的matplotlib 畫軌道實(shí)例,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧

如下所示:

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 畫軌道實(shí)例就是小編分享給大家的全部內(nèi)容了,希望能給大家一個參考,也希望大家多多支持腳本之家。

相關(guān)文章

最新評論

兴义市| 临泽县| 大荔县| 潞城市| 高碑店市| 轮台县| 温宿县| 甘谷县| 中西区| 敦煌市| 宜章县| 临猗县| 始兴县| 滦平县| 南部县| 内黄县| 岳池县| 三原县| 四平市| 蒙山县| 金昌市| 商洛市| 沈阳市| 翁牛特旗| 广州市| 安陆市| 丽江市| 旬邑县| 恩施市| 女性| 丽江市| 呼玛县| 岳西县| 河东区| 云霄县| 平泉县| 鸡西市| 龙陵县| 任丘市| 图们市| 安多县|