示例#1
0
from pylab import *
from sim import simulate
from curve import curvature

l = 10.0

#Destination
dests = array(curvature(10))

#initial point
initial = [0.0, 0.0]

#stud mag without grid
[x, y, gx, gy, corr, err] = simulate(0, l, initial, pi/2, dests, 1.0,1,100, c_mag = 1e-2,c_align=0)
print err
plot(x,y, 'black', label='High-res sensor', linewidth=2)

#stud mag with grid
[x, y, gx, gy, corr, err] = simulate(9, l, initial, pi/2, dests, 1.0,1,100, c_mag = 1e-2,c_align=0)
print err
plot(x,y, 'black', label='High-res sensor with Grid', linestyle='-.', linewidth=3)

#hagga mag without grid
[x, y, gx, gy, corr, err] = simulate(0, l, initial, pi/2, dests, 1.0, 1,100,c_mag = 1.5e-1,c_align=0)
print err
plot(x,y, 'gray', label='Low-res sensor', linewidth=2)

#hagga mag with grid
[x, y, gx, gy, corr, err] = simulate(9, l, initial, pi/2, dests, 1.0, 1,100,c_mag = 1.5e-1,c_align=0)
print err
plot(x,y, 'gray',label='Low-res sensor with Grid', linestyle='-.', linewidth=3)
示例#2
0
from pylab import *
from sim import simulate
from curve import curvature
p = []
avg_errors = []
max_errors = []
corrs = []
initial = [0.0,0.0]
l = 10.0

for n in arange(l, l*10, 3):
  errs = []
  dests=array(curvature(n))
  
  corrs.append(0)
  for j in range(0,100):
    print n, j
    tmp = simulate(9, l, initial, pi/2, dests, 1.0,1,100)
    errs.append(tmp[-1])
    corrs[-1] = tmp[-2]
    
    
  p.append(n)
  avg_errors.append(average(errs))  
  max_errors.append(max(errs))

  print n, avg_errors[-1], max_errors[-1]
print p
print avg_errors
print max_errors
print corrs