Questions tagged [scipy]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1
vote
1answer
1k views

How to design bandpass filter in python when centre frequency is greater than sampling frequency?

I implemented basic bandpass filter using Scipy resource. But how to implement the same when signal is centred at non-zero centre frequency with centre frequency > sampling frequency. I trying to ...
1
vote
1answer
4k views

How to normalize output of $\tt scipy.signal.correlate$

I have 2 different signals and I'm trying to cross-correlate then using Python 2.7 and scipy.signal.correlate. How do I normalize my results (such that the max ...
3
votes
2answers
657 views

Do DCT results contain phase spectrum?

I've been working much with numpy.fft lately where just like the documentation says: When the input a is a time-domain signal and ...
4
votes
3answers
20k views

Calculating values of frequency bins in Python

I have to use FFT to determine the period of waves inside a signal. After applying FFT on a window of 10000 point from a signal, I get something like this: What I don't understand is that FFT is ...
3
votes
2answers
2k views

Python butter filter: choosing between analog and digital filter types

I have a system that performs wireless sampling of data about every 7.5ms (133Hz). Due to it being wireless, I get occasional data drop out. I want to construct a LP butter filter with cut-off ...
1
vote
2answers
198 views

problem with decimate in scipy 0.18

I use scipy.signal.decimate to downsample a large spectroscopic data-set. This has worked superbly so far. ...
1
vote
2answers
1k views

Why does python's scipy.signal.dimpulse introduce delay in impulse response?

Consider a simple linear, time invariant system of the form: $y_k = cy_{k-1} + (1-c)x_k$ The impulse response of this system can be computed by either dimpulse or by applying lfilter to a vector ...
1
vote
2answers
229 views

Need of $\tt abs()$ method when plotting a power spectral density for a given dataset

I am a newbie in signal processing and would like to know the significance of using the abs() function and squaring the values received as an output of ...
1
vote
1answer
70 views

$n\textrm{ Hz}$ waveform sampled at $m\textrm{ Hz}$ per second

Here is an example of plotting a square wave given in SciPy Documentation A $5\textrm{ Hz}$ waveform sampled at $500\textrm{ Hz}$ for 1 second: ...
3
votes
2answers
120 views

How to discard coefficients in the frequency domain without losing too much quality?

I was reading a text that shows that one can throw away 66% of the coefficients and gets an image that is still acceptable. I tried to replicate it using python (last example) and selectively ...
1
vote
2answers
4k views

Is there a function in python which returns the Amplitude/Sound Pressure Level of a sound file(like .wav) file?

I found a function called scipy.io.wavefile.read() which returns the sampling rate and an array. But in the documentation, it is not specified that whether the ...
0
votes
3answers
334 views

Problems designing high pass FIR filter with $f_\textrm{cutoff} \ll f_\textrm{s}$

Is there any inherant problems with high-passing a signal with a very low cutoff frequency compared to the sampling rate? I'm not sure what I am missing as I cant find any info on it. Examples ...
-1
votes
1answer
1k views

Converting from radians/samples to normalized frequency range $[0.0, 1.0]$?

I'm trying to understand the following code, specifically how to use the results of scipy.signal.freqz where w in ...
2
votes
0answers
225 views

Scipy.signal.dimpulse(system): how to translate a StateSpace to a “system”? [closed]

The command scipy.signal.dimpulse(system, x0=None, t=None, n=None) does not seem to accept a scipy StateSpace as an input to its ...
0
votes
1answer
637 views

How does scipy.signal.deconvolve work?

The deconvolve(signal, divisor) function of scipy ...
3
votes
1answer
3k views

Butterworth band pass filter

I am trying to build a Butterworth bandpass filter. My intent is to have a pass band between $250\textrm{ Hz}$ and $1000\textrm{ Hz}$. However the plot doesn't look right, as can be seen from the ...
0
votes
0answers
147 views

Translating Bessel filter from python to C - odd behavior past the cutoff frequency

I recently translated the code from scipy.signal.bessel to C in order to implement a bessel filter in a piece of analysis ...
2
votes
1answer
456 views

discrete digital filtering in python

I am not from DSP background and also not sure if I am posting this in right community. I have basic understanding of resistors, capacitors and filters. My requirement is as follows- The supply ...
1
vote
0answers
2k views

How to design an interpolated notch filter in Python to filter powerline harmonics

I have a a time domain data set that records the magnetic field vs time, which must be processed to reveal an embedded signal. This data also contains power line harmonics (i.e. multiples of 60 Hz) ...
2
votes
2answers
91 views

How to know when a sound in a given frequency starts and ends?

I need to detect two sounds (913 Hz and 1370 Hz), knowing the time in which each starts and ends. At this time, I am using a band pass filter (butter), and this is the result: Here is the first time ...
0
votes
2answers
207 views

Beginner Question: Extracting Features for Image Classification

I'm a Python programmer, but a beginner to image processing :) Apologies if this isn't a suitable question for the forum, happy to rewrite or move. I want to create a supervised classifier to ...
2
votes
0answers
178 views

Alternative to scipy maximum_filter to find peaks?

I'm using scipy maximum_filter together with binary_erosion to find peaks in my spectrum. I just wondering if there are any ...