@@ -0,0 +1,2100 @@
+#!/bin/python%0A%0Aimport matplotlib as mpl%0Aimport matplotlib.pyplot as plt%0Aimport matplotlib.cbook as cbook%0Aimport numpy as np%0Aimport math%0A%0A%09# State vector:%0A%09# 0-3: quaternions (q0, q1, q2, q3)%0A%09# 4-6: Velocity - m/sec (North, East, Down)%0A%09# 7-9: Position - m (North, East, Down)%0A%09# 10-12: Delta Angle bias - rad (X,Y,Z)%0A%09# 13-14: Wind Vector  - m/sec (North,East)%0A%09# 15-17: Earth Magnetic Field Vector - milligauss (North, East, Down)%0A%09# 18-20: Body Magnetic Field Vector - milligauss (X,Y,Z)%0A%0Adata = np.genfromtxt('FlowRawOut.txt', delimiter=' ', skip_header=1,%0A%09skip_footer=1, names=%5B'time', 'flowRadX', 'flowRadY', 'gyroX', 'gyroY', 'flowMsX', 'flowMsY', 'GpsVn', 'GpsVe', 'distance' , 'predFlowX', 'predFlowY'%5D)%0A%0Afig = plt.figure()%0A%0Aax1 = fig.add_subplot(211)%0A%0Aax1.set_title(%22Flow angular rate estimate%22)    %0Aax1.set_xlabel('time (s)')%0Aax1.set_ylabel('flow X (rad)')%0Aax1.plot(data%5B'time'%5D, data%5B'flowRadX'%5D, color='r', label='flow x')%0A#ax1.plot(data%5B'time'%5D, data%5B'gyroX'%5D, color='b', label='angular rate x')%0Aax1.plot(data%5B'time'%5D, data%5B'predFlowX'%5D, color='g', label='pred flow x')%0A%0Aax2 = fig.add_subplot(212)%0A  %0Aax2.set_xlabel('time (s)')%0Aax2.set_ylabel('flow Y (rad)')%0Aax2.plot(data%5B'time'%5D, data%5B'flowRadY'%5D, color='r', label='flow y')%0A#ax2.plot(data%5B'time'%5D, data%5B'gyroY'%5D, color='b', label='angular rate y')%0Aax2.plot(data%5B'time'%5D, data%5B'predFlowY'%5D, color='g', label='pred flow y')%0A%0Afigvel = plt.figure()%0A%0Aav1 = figvel.add_subplot(311)%0A%0Aav1.set_title(%22Flow Velocity estimate%22)    %0Aav1.set_xlabel('time (s)')%0Aav1.set_ylabel('north velocity (m/s)')%0Aav1.plot(data%5B'time'%5D, data%5B'flowMsX'%5D, color='r', label='Flow Vn')%0Aav1.plot(data%5B'time'%5D, data%5B'GpsVn'%5D, color='g', label='GPS Vn')%0A%0Aav2 = figvel.add_subplot(312)%0A  %0Aav2.set_xlabel('time (s)')%0Aav2.set_ylabel('east velocity (m/s)')%0Aav2.plot(data%5B'time'%5D, data%5B'flowMsY'%5D, color='r', label='Flow Ve')%0Aav2.plot(data%5B'time'%5D, data%5B'GpsVe'%5D, color='g', label='GPS Ve')%0A%0Aav3 = figvel.add_subplot(313)%0A %0Aav3.set_xlabel('time (s)')%0Aav3.set_ylabel('ground distance (m)')%0Aav3.plot(data%5B'time'%5D, data%5B'distance'%5D, color='b', label='distance')%0A%0Aplt.show()%0A

Add plot_flow.py file lost during merge