How not to fall victim to the star notation of the convolution?
Please
accept marketing cookies
to access the video player.
The Convolution Series
- Definition of convolution and intuition behind it
- Mathematical properties of convolution
- Convolution property of Fourier, Laplace, and z-transforms
- Identity element of the convolution
- Star notation of the convolution
- Circular vs. linear convolution
- Fast convolution
- Convolution vs. correlation
- Convolution in MATLAB, NumPy, and SciPy
- Deconvolution: Inverse convolution
- Convolution in probability: Sum of independent random variables
Introduction
Taking advantage of the introduction of delays in the previous article, I wanted to warn you against a common pitfall when talking about the convolution [1].
The star notation x[n]∗h[n] is very convenient. It must, however, be used with caution.
The following notation should be clear to you by now
y[n]=x[n]∗h[n],n∈Z.(1)What if we wanted to obtain a delayed version of y[n], i. e., y[n−n0]? A natural move would be to substitute n←n−n0
y[n−n0]=?x[n−n0]∗h[n−n0],(2)but…
Let’s evaluate the right hand side of “Equation” 2
x[n−n0]∗h[n−n0]=k=−∞∑∞x[k−n0]h[n−n0−k]=k=−∞∑∞x[k]h[n−2n0−k]=y[n−2n0].(3)By blindly substituting n←n−n0, we overshot the desired delay by a factor of two.
The correct way to write this is [1, Eq. 2.52]
y[n−n0]=x[n]∗h[n−n0]=k=−∞∑∞x[k]h[n−n0−k].(4)This is just one of many problems that arise when using the star notation.
Useful notational tip
When the convolution looks any way different from the typical x[n]∗h[n], I try to bring it back to that basic form by defining “helper functions”. Then I use the definition of the convolution and substitute the original functions again, inserting the correct argument.
How this works is best explained through an example.
Example 1: Both convolved signals delayed
Let’s say the operands of the convolution we need to perform are both delayed by different amounts, i. e., we want to calculate
x[n−nx]∗h[n−nh]=…n,nx,nh∈Z.(5)Let’s define two “helper functions” x1[n],h1[n]
x1[n]=x[n−nx],(6)h1[n]=h[n−nh].(7)Now we can insert these functions into Equation 5
x[n−nx]∗h[n−nh]=x1[n]∗h1[n],(8)use the definition of the convolution
x1[n]∗h1[n]=k=−∞∑∞x1[k]h1[n−k],(9)and finally substitute the original functions x[n] and h[n] according to Equations 6 and 7 respectively
k=−∞∑∞x1[k]h1[n−k]=k=−∞∑∞x[k−nx]h[n−k−nh].(10)This approach always works for me. At the same time, any shortcuts in an attempt not to use it inevitably led me to an error in calculations.
Example 2: One of the convolved signal is time-reversed
This final example should make clear why “helper functions” ensure us that we correctly evaluate the star notation. In this example, one of the operands is time-reversed.
x[n]∗h[−n]=h2[n]=h[−n]x[n]∗h2[n]=k=−∞∑∞x[k]h2[n−k]=h2[n−k]=h[−(n−k)]k=−∞∑∞x[k]h[k−n].(11)Is it possible to guess the correct answer right away? Yes, definitely. But it is not easy, especially if the arguments get even more complicated and the convolution is a part of a much larger body of derivations.
Summary
In this article, we discussed notational issues concerning the discrete convolution and how to avoid common pitfalls when using the star notation. It all boils down to the idea of using “helper functions”.
Bibliography
[1] A. V. Oppenheim, R. W. Schafer Discrete-Time Signal Processing, 3rd Edition, Pearson 2010.
Comments powered by Talkyard.