0
$\begingroup$

I'am pretty new to this topic so i dont have much experience with DSP. I want to filter(highpass) a WAVE file, i have already programmed the FFT, and now want to filter the FFT vector via convolution (which is a multiplication in frequency domain ?!).

1) Can anyone give me an easy equation for such a filter?

2) If i have an equation i compute it, also Fourier-Transform it, an multiply it with my FFT vector and inverse Transform it an i have applied the filter?!

Improve this question
$\endgroup$
3
  • $\begingroup$ Google "overlap add", that's the algorithm you want to use. Highpass it typically much easier done with an IIR filter than with an FFT. $\endgroup$
    – Hilmar
    Commented Jan 25, 2014 at 14:25
  • $\begingroup$ dspguide.com/ch18.htm $\endgroup$
    – sellibitze
    Commented Jan 25, 2014 at 21:25
  • $\begingroup$ IF you don't already know the impulse response you want your signal convolved with then "filter design" is the topic you should research. try "windowed sinc" or use a tool like Gnu/Octave which comes with various nice functions for designing these things. Apart from the "windowed since" approach most of the others are kind of a black art and many people just rely on tools to compute the coefficients for them. $\endgroup$
    – sellibitze
    Commented Jan 25, 2014 at 21:30

1 Answer 1

Reset to default
0
$\begingroup$

Actually, yes. As Hilmar pointed out, it might be cheaper to use an ordinary filter, but that depends on the circumstances.

If you have two signals $A$ and $B$ with Fourier transform being $\mathcal{F}$, their convolution $A\ast B$ will be

\begin{equation} A\ast B = \tilde{\mathcal{F}}(\mathcal{F}(A)\cdot\mathcal{F}(B)) \end{equation}

With $\tilde{\mathcal{F}}$ being the inverse transform. This is the well known "convolution lemma". There is not much of a theory to see. In the Fourier domain, you just replace convolution by (point-wise) multiplication. Your question already contained the answer, to me it just appears that you were not sure about the details.

Improve this answer
$\endgroup$
5
  • $\begingroup$ Well thank you =) But can you give me an equation to campute the signal to convolute with? $\endgroup$
    – ShuftY
    Commented Jan 25, 2014 at 15:10
  • $\begingroup$ I didn't notice you are asking for that. The signal you mean depends on the situation. Eg. if you want to have a highpass, you should set the "low" Fourier coefficients (The ones with low indices and their symmetry partners with high indices) to zero. You would then multiply with a "box" function. The signal you would be convolving with then is the inverse transform of that box function. $\endgroup$
    – user7358
    Commented Jan 25, 2014 at 15:17
  • $\begingroup$ Okay just to get it right. Lets say i want to boost all above 100Hz. I got 44100 samplepoints. I make a new vector with vector[0] to[99] and [44099] to[43999] = 0 an the rest is set to 1. This is my "Box" function. Then i inverse Transform this vector and then multiply it with my Fouriertransformed vector of samplingpoints and denn inverse transform it? $\endgroup$
    – ShuftY
    Commented Jan 25, 2014 at 15:53
  • $\begingroup$ No. You would multiply the transformed of the signal you want to filter with that box function. You could also convolve the signal with the inverse transform of the box function instead. $\endgroup$
    – user7358
    Commented Jan 25, 2014 at 15:54
  • $\begingroup$ Note that you have to use an FFT that is longer than the sum of the length of your samples and the length of the impulse response of your filter. Just zeroing bins can cause ringing artifacts. See: dsp.stackexchange.com/questions/6220/… If you have more than 1 second of data, you will have to use overlap-add or overlap save fast convolution to combine your 1 second frames. $\endgroup$
    – hotpaw2
    Commented Jan 25, 2014 at 17:19

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.