Implement Theil-Sen's method for fitting a line to noisy data (used in bandwidth estimation).

webrtc/modules/remote_bitrate_estimator/test/plot_dynamics.py

 #!/usr/bin/env python
 # Copyright (c) 2015 The WebRTC project authors. All Rights Reserved.
 #
 # Use of this source code is governed by a BSD-style license
 # that can be found in the LICENSE file in the root of the source
 # tree. An additional intellectual property rights grant can be found
 # in the file PATENTS.  All contributing project authors may
 # be found in the AUTHORS file in the root of the source tree.
 
 # This script is used to plot simulation dynamics. The expected format is
@@ skipping 127 matching lines @@
 
   detector_state = Figure("DetectorState")
   detector_state.addSubplot(['offset_ms'], "Time (s)", "Offset")
   detector_state.addSubplot(['gamma_ms'], "Time (s)", "Gamma")
 
   trendline_state = Figure("TrendlineState")
   trendline_state.addSubplot(["accumulated_delay_ms", "smoothed_delay_ms"],
                              "Time (s)", "Delay (ms)")
   trendline_state.addSubplot(["trendline_slope"], "Time (s)", "Slope")
 
+  target_bitrate = Figure("TargetBitrate")
+  target_bitrate.addSubplot(['target_bitrate_bps'], "Time (s)", "Bitrate (bps)")
+
   # Select which figures to plot here.
-  figures = [receiver, detector_state, trendline_state]
+  figures = [receiver, detector_state, trendline_state, target_bitrate]
 
   # Add samples to the figures.
   for line in sys.stdin:
     if line.startswith("[ RUN      ]"):
       test_name = re.search(r'\.(\w+)', line).group(1)
     if line.startswith("PLOT"):
       try:
         (var_name, ssrc, alg_name, time, value) = parse_plot_line(line)
         for f in figures:
           # The sample will be ignored bv the figures that don't need it.
           f.addSample(var_name, ssrc, alg_name, time, value)
       except ParsePlotLineException as e:
         print e.reason
         print e.line
 
   # Plot figures.
   for f in figures:
     fig = plt.figure(f.name)
     f.plotFigure(fig)
     if save_figure:
       fig.savefig(test_name + f.name + ".png")
   plt.show()
 
 if __name__ == '__main__':
   main()