|
程序如下:
while position_r[1] > 0:
while position_r[0] <4:
rss_test=[]
position_t=[]
for i in range(MAX):
dl=abs(link[i][1]*position_r[0]+link[i][2]*position_r[1]+link[i][3])/math.sqrt(link[i][1]**2+link[i][2]**2)
if dl<radius:
print(link[i][0])
for j in link[i][4]:
rss_test.append(j)
if len(rss_test)!=0:
fig, (ax0, ax1) = plt.subplots(ncols=2, figsize=(8, 4))
bins=np.arange(-20.,20.,1.)
(n1,bins,patches)=ax0.hist(rss1,bins, normed=1, facecolor='green', alpha=0.5)
plt.xlabel('RSS')
plt.ylabel('Probability')
plt.title(r'Histogram of RSS')
plt.subplots_adjust(left=0.15)
bins=np.arange(-20.,20.,1.)
(n2,bins,patches)=ax1.hist(rss_test,bins, normed=1, facecolor='green', alpha=0.5)
plt.xlabel('RSS')
plt.ylabel('Probability')
plt.title(r'Histogram of RSS')
# Tweak spacing to prevent clipping of ylabel
plt.subplots_adjust(left=0.15)
#plt.show()
print(n2)
sum=0.0
for p in range(len(n1)):
sum=math.sqrt(n1[p]*n2[p])+sum
dis=math.sqrt(1-sum)
print(dis)
print(d)
d.append(dis)
position_t.append(position_r[0])
position_t.append(position_r[1])
position_t.append(dis)
print(position_t)
target.append(position_t)
print(target)
print(position_r)
position_r[0]=position_r[0]+2*radius
position_r[0]=0.5#置回为初始x值
position_r[1]=position_r[1]-2*radius
print(d)
print(target)
程序运行可以输出想要的结果,但最后一行会出现alloc: invalid block: 034F24A8: 88 3,请问问题出在哪呢?谢谢高手解答~
|
|