Demonstrates the convolution theorem with the use of the Complex FFT, Complex-by-Complex Multiplication, and Support Functions.
The convolution theorem states that convolution in the time domain corresponds to multiplication in the frequency domain. Therefore, the Fourier transform of the convoution of two signals is equal to the product of their individual Fourier transforms. The Fourier transform of a signal can be evaluated efficiently using the Fast Fourier Transform (FFT).
Two input signals,
b[n], with lengths
n2respectively, are zero padded so that their lengths become
N, which is greater than or equal to
(n1+n2-1)and is a power of 4 as FFT implementation is radix-4. The convolution of
b[n]is obtained by taking the FFT of the input signals, multiplying the Fourier transforms of the two signals, and taking the inverse FFT of the multiplied result.
This is denoted by the following equations: where
B[k]are the N-point FFTs of the signals
b[n]respectively. The length of the convolved signal is
- Block Diagram:
- Variables Description:
testInputA_f32points to the first input sequence
srcALenlength of the first input sequence
testInputB_f32points to the second input sequence
srcBLenlength of the second input sequence
outLenlength of convolution output sequence,
(srcALen + srcBLen - 1)
AxBpoints to the output array where the product of individual FFTs of inputs is stored.
- NMSIS DSP Software Library Functions Used: