Merge pull request #133 from jaidevd/jd-doc-nonstat

Add documentation / tutorials for nonstationary signals

.gitignore

@@ -52,6 +52,7 @@ coverage.xml
 
 # Sphinx documentation
 doc/_build/
+doc/apiref/*.rst
 
 # PyBuilder
 target/

doc/index.rst

@@ -12,6 +12,7 @@ Contents:
    :maxdepth: 2
 
    introduction
+   nonstationary_signals
    auto_gallery
    api
 

doc/misc_plots/nonstationary_phase_plot.py

@@ -0,0 +1,10 @@
+from tftb.generators import fmlin, amgauss
+import numpy as np
+import matplotlib.pyplot as plt
+
+y_nonstat, _ = fmlin(2048)  # Already analytic, no need of Hilbert transorm
+y_nonstat *= amgauss(2048)
+plt.plot(np.real(y_nonstat), np.imag(y_nonstat))
+plt.xlabel("Real part")
+plt.ylabel("Imaginary part")
+plt.show()

doc/misc_plots/stationary_phase_plot.py

@@ -0,0 +1,17 @@
+import numpy as np
+import matplotlib.pyplot as plt
+from scipy.signal import hilbert
+
+fs = 32768
+ts = np.linspace(0, 1, fs)
+y1 = np.sin(2 * np.pi * 697 * ts)
+y2 = np.sin(2 * np.pi * 1336 * ts)
+y = (y1 + y2) / 2
+
+
+y = y[:(fs / 16)]
+y_analytic = hilbert(y)
+plt.plot(np.real(y_analytic), np.imag(y_analytic))
+plt.xlabel("Real part")
+plt.ylabel("Imaginary part")
+plt.show()

doc/misc_plots/touchtone.py

@@ -0,0 +1,25 @@
+#! /usr/bin/env python
+# -*- coding: utf-8 -*-
+# vim:fenc=utf-8
+#
+# Cube26 product code
+#
+# (C) Copyright 2015 Cube26 Software Pvt Ltd
+# All right reserved.
+#
+# This file is confidential and NOT open source.  Do not distribute.
+#
+
+"""
+
+"""
+import numpy as np
+import matplotlib.pyplot as plt
+fs = 32768
+ts = np.linspace(0, 1, fs)
+y1 = np.sin(2 * np.pi * 697 * ts)
+y2 = np.sin(2 * np.pi * 1336 * ts)
+y = (y1 + y2) / 2
+plt.plot(ts, y)
+plt.xlim(0, 0.1)
+plt.show()