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)
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