One Paragraph of project description goes here
These instructions will get you a copy of the project up and running on your local machine for development and testing purposes. See deployment for notes on how to deploy the project on a live system.
| #!/usr/bin/env python | |
| # Simple example of Wiener deconvolution in Python. | |
| # We use a fixed SNR across all frequencies in this example. | |
| # | |
| # Written 2015 by Dan Stowell. Public domain. | |
| import numpy as np | |
| from numpy.fft import fft, ifft, ifftshift |
| import numpy as np | |
| ## Tri Diagonal Matrix Algorithm(a.k.a Thomas algorithm) solver | |
| def TDMAsolver(a, b, c, d): | |
| ''' | |
| TDMA solver, a b c d can be NumPy array type or Python list type. | |
| refer to http://en.wikipedia.org/wiki/Tridiagonal_matrix_algorithm | |
| and to http://www.cfd-online.com/Wiki/Tridiagonal_matrix_algorithm_-_TDMA_(Thomas_algorithm) | |
| ''' | |
| nf = len(d) # number of equations |