12
FFT -> zeroing coefficients -> IFFT is not the correct way of doing filtering - the actual filter realized by doing so has poor characteristics.
The correct way of filtering signals is to compute the coefficients of a digital filter, a process known as filter design and for which a large body of software tools/documentation is available, and apply it to ...
8
The best way to apply frequency domain filtering for signal streams is overlap add (or related flavors overlap save, or block convolvers, etc.).
You basically take in one frame at a time (say 1024 samples). Zero pad to twice the length (2048), do an FFT, multiply with (also zero padded) transfer function of the filter, do an inverse FFT. Save the last 1024 ...
8
In theory you can do this, but in practice it is difficult to do because the time and phase alignment must be pretty good for it to work. If the alignment is good you will get the destructive interference that you are seeking. If they aren't, you will get some constructive interference. Even worse, whether they are destructively or constructively ...
8
The result will indeed be a high pass filter. From your difference equation, the transfer function of the low pass filter is
$$H_l(z)=\frac{\beta}{1-(1-\beta)z^{-1}}\tag{1}$$
with $\beta=1/\alpha$.
Note that this is actually a leaky integrator, not a classic low pass filter, because its frequency response does not have a zero at Nyquist.
The high pass ...
8
In general you can't simply subtract a low-pass filtered version of a signal from the original one to obtain a high-pass filtered signal. The reason is as follows. What you're actually doing is implement a system with frequency response
$$H(\omega)=1-H_{LP}(\omega)\tag{1}$$
where $H_{LP}(\omega)$ is the frequency response of the low-pass filter. Note that $...
7
A high pass filter will remove a large portion to all of the 0 Hz component or DC offset from a signal. That's similar to subtracting the mean average, which for a normally distributed signal will make around half of the signal negative (since it was below the original average mean).
If you want the mean to be the same after a high pass filter, you could ...
7
In the $s$-domain, the LPF-to-BPF transformation doubles the order of the filter. that is because the LPF has one transition from passband to stopband, but the BPF has two such transitions.
Remember that when we do frequency response, we substitute $s=j\omega$. That is, we evaluate the $s$-plane transfer function $H(s)$ as $H(j\omega)$. The magnitude of $...
answered Mar 12 '15 at 15:42
robert bristow-johnson
13.9k3 gold badges24 silver badges58 bronze badges
6
A high pass filter has some properties that are similar to taking a derivative. The derivative of a function that is always positive will be negative wherever the slope is negative. There's nothing you can do about it. For example, take the function $f(t) = 400 + \sin \omega t$. This wiggles around between 399 and 401 (very positive), yet the derivative ...
6
Attenuation in this context refers to a reduction in amplitude of some signal component as it passes through a system (in your case, the system being the highpass filter. Linear filters are designed to pass certain frequency components while stopping (or attenuating) others. As the name implies, a highpass filter is designed to pass high-frequency content ...
6
The input to your filter is essentially a step function (since it has value 0 prior to t = 0, and a positive value for t > 0), so you see the step response of the filter, hence the initial ringing. This is expected behaviour, and after a suitable amount of time the step response will have settled to zero.
5
This is a first order, low pass infinite impulse response filter (IIR). It is some times called the exponentially weighted moving average. You can find the derivation on Wikipedia. Orfanidis' Introduction to Signal Processing gives a nice description as well.
Don't think of it as high passing y and low passing x, this will lead to confusion. It is ...
5
You can design the different filter types directly with the scipy.signal functions.
There are three main functions for creating
finite impulse response filters with the scipy.signal
package.
signal.remez
signal.firwin
signal.firwin2
The remez function, as arguments, takes the number of taps
(order+1), the "bands", and the "desired" gain. The
"bands" is ...
5
That depends on the definition of high-pass filter. If you define a high-pass filter as a filter that has high response in the high frequencies in frequency domain, then the easiest way is to take a look at the magnitude of Fourier transform, (by definition).
Applying Fourier transform (in Matlab)
A = fftshift(abs(fft2(padarray([-1 -1 -1; 0 0 0; 1 1 1],[...
5
They are both highpass type filters, but used with very different intentions.
One should immediately observe the fundamental difference that the output of unsharp masking filter is an enhanced image to be viewed by humans, whereas the output of the Sobel (edge detector) filter is not an image to be viewed by humans, but rather a description of the edges to ...
5
In general there is no straightforward analytical solution. As you know, you need to solve
$$\left|H(e^{j\omega_c})\right|=\frac{1}{\sqrt{2}}\tag{1}$$
for $\omega_c$, where it is assumed that the maximum filter gain equals $1$.
For Butterworth filters, the specified cut-off frequency always equals the $3\textrm{ dB}$ frequency. This is not the case for other ...
4
Consider this filter in 1D (for the sake of simplicity):
[-1 1 0]
This filter responds to pixel differences, i.e. edges (left pixel * (-1) + center pixel * 1 + right pixel * 0) which is (center pixel - left pixel).
Response of this filter is positive on rising edges and negative on falling edges.
The filter is not symmetric and hence it will "shift" the ...
4
You indicated that you're having trouble figuring out how to design a suitable highpass filter. One method is to first design a lowpass filter prototype, then apply a transform that warps the filter's response into a filter of another type (such as a highpass or bandpass filter). This is done by substituting an expression for $z^{-1}$ into the transfer ...
4
You're probably best off asking your professor for clarification. As you reasoned, it sounds like you're being asked to multiply the image in the frequency domain by either a rectangular or circular mask that eliminates all frequencies outside the mask. This is not typically how you would apply a filter, however, as "ideal" frequency-domain filter masks have ...
4
For an LTI system another method to generate a highpass filter impulse response $h_{hp}[n]$, from an existing lowpass filter impulse reponse $h_{lp}[n]$ is the following:
$$h_{high}[n]=(-1)^n h_{low}[n]$$
From the DTFT properties we can see that:
$$h_{high}[n]=e^{j\pi n} h_{low}[n]$$
$$H_{high}(e^{j\omega})=H_{low}(e^{j(\omega - \pi)})$$
which would ...
4
This is a first-order discrete-time Butterworth filter, which can be obtained via the bilinear transform from the following analog prototype filter (with normalized cut-off frequency):
$$H(s)=\frac{s}{1+s}\tag{1}$$
The bilinear transform replaces $s$ by $k\frac{z-1}{z+1}$, where $k$ is chosen such that the discrete-time filter has the specified cut-off ...
4
Note that the pole locations would be the same for all $4$ classic types of frequency-selective filters (low pass, high pass, band pass, band stop). It's the location of the zeros that determines which filter type it is. Since we don't see any zeros, we can assume that they are at infinity, hence, the filter is a low pass filter. For a high pass filter you ...
4
A discrete-time first-order high pass filter with unity gain at Nyquist and a zero at DC is described by the following difference equation:
$$y[n]=\frac{1+\alpha}{2}\big(x[n]-x[n-1]\big)+\alpha y[n-1],\qquad -1<\alpha<1\tag{1}$$
Its transfer function is given by
$$H(z)=\frac{1+\alpha}{2}\frac{1-z^{-1}}{1-\alpha z^{-1}}\tag{2}$$
Evaluating the squared ...
3
If the DFT is the Uniform sampling from $ 0 $ to $ 2 \pi $ then the first bin is given by:
$$ x[k] = \sum_{n = 0}^{N - 1} x[n] $$
Namely it is the sum of all the samples.
Hence in order to remove the DC (Mean) all you need is a filter which has zero in its DC bin.
Since, the filtered signal, which is a convolution (Circular) of the DFT of the input signal ...
3
If you want to detect the peak value of the time-domain signal, it is important not to change phase of each frequency signals. Otherwise, the peak value will be affected by it and will slightly change.
I think you already have recorded signal and you can process it after the recording.
If so, it is better to use zero-phased filters, like this function in ...
3
Regarding scaling:
If you sum up the coefficients, you get the magnitude for DC. Hence, it makes sense that you get those numbers ($\approx 1$ for LP, $\approx 0$ for HP).
In addition to Matt L.'s excellent answer one can just point out that what he is using is referred to as magnitude complementary filters, which is the common case for linear-phase FIR ...
3
I agree that a large-order (12 taps+) high-pass Butterworth filter will be sufficient.
Why? Well, Butterworth filters are commonly used in professional audio-editing applications to do filtering, because their frequency response is designed to be as flat as possible in the passband. This means that the filtered audio is of higher quality/fidelity, as ...
3
The simplest option is to remove the low-pass signal from the data. Just subtract the outcome of your code from the original image. This amounts to computing:
$$I_\textrm{high} = \delta \ast I - g\ast I = (\delta -g) \ast I\,.$$
The term $\delta \ast I$ corresponds to the convolution of a Dirac delta to the image, which is an invariant operation.
As the ...
3
out[i] = out[i-1] + (alpha*(in[i] - out[i-1]));
From Wikipedia's article on Low pass filters:
That is, this discrete-time implementation of a simple RC low-pass filter is the exponentially weighted moving average
[…]
The loop that calculates each of the n outputs can be refactored into the equivalent:
for i from 1 to n
y[i] := y[i-1] + α *...
3
In addition to @MarcusMüller answer, you can find other names and historical details on this type of filter in Is there a technical term for this simple method of smoothing out a signal?.
I will not rewrite the other answer, just give here a few aspects:
This filter is so useful and common that its has several names, for instance first order exponential ...
3
Today I received an answer from the authors of the book. It turns out that figure 1.7a was wrong, and they sent me the updated one.
The corrected diagram will be added in this Summer's update!
Only top voted, non community-wiki answers of a minimum length are eligible
