docs: def plot_prh_filtered(p, r, h, p_lf, r_lf, h_lf):
      '''Plot original and low-pass filtered PRH data
  
      Args
      ----
      p: ndarray
          Derived pitch data
      r: ndarray
          Derived roll data
      h: ndarray
          Derived heading data
      p_lf: ndarray
          Low-pass filtered pitch data
      r_lf: ndarray
          Low-pass filtered roll data
      h_lf: ndarray
          Low-pass filtered heading data
      '''
      import numpy
  
      fig, (ax1, ax2, ax3) = plt.subplots(3, 1, sharex='col')
  
      #rad2deg = lambda x: x*180/numpy.pi
  
      ax1.title.set_text('Pitch')
      ax1.plot(range(len(p)), p, color=_colors[0], linewidth=_linewidth,
              label='original')
      ax1.plot(range(len(p_lf)), p_lf, color=_colors[1], linewidth=_linewidth,
              label='filtered')
  
      ax2.title.set_text('Roll')
      ax2.plot(range(len(r)), r, color=_colors[2], linewidth=_linewidth,
              label='original')
      ax2.plot(range(len(r_lf)), r_lf, color=_colors[3], linewidth=_linewidth,
              label='filtered')
  
      ax3.title.set_text('Heading')
      ax3.plot(range(len(h)), h, color=_colors[4], linewidth=_linewidth,
              label='original')
      ax3.plot(range(len(h_lf)), h_lf, color=_colors[5], linewidth=_linewidth,
              label='filtered')
  
      plt.ylabel('Radians')
      plt.xlabel('Samples')
      for ax in [ax1, ax2, ax3]:
          ax.legend(loc="upper right")
  
      plt.show()
  
      return None