Wei Cui shineboycw
danstowell
/ wiener_deconvolution_example.py
Last active
October 1, 2025 12:43
Simple example of Wiener deconvolution in Python
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #!/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 |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| function [maxtab, mintab]=peakdet(v, delta, x) | |
| %PEAKDET Detect peaks in a vector | |
| % [MAXTAB, MINTAB] = PEAKDET(V, DELTA) finds the local | |
| % maxima and minima ("peaks") in the vector V. | |
| % MAXTAB and MINTAB consists of two columns. Column 1 | |
| % contains indices in V, and column 2 the found values. | |
| % | |
| % With [MAXTAB, MINTAB] = PEAKDET(V, DELTA, X) the indices | |
| % in MAXTAB and MINTAB are replaced with the corresponding | |
| % X-values. |