# 7.4 Discrete time circular convolution and the dtfs

 Page 1 / 1
This module describes the circular convolution algorithm and an alternative algorithm

## Introduction

This module relates circular convolution of periodic signals in one domain to multiplication in the other domain.

You should be familiar with Discrete-Time Convolution , which tells us that given two discrete-time signals $x(n)$ , the system's input, and $h(n)$ , the system's response, we define the output of the system as

$y(n)=(x(n), h(n))=\sum_{k=()}$ x k h n k
When we are given two DFTs (finite-length sequences usually oflength $N$ ), we cannot just multiply them together as we do in the above convolutionformula, often referred to as linear convolution . Because the DFTs are periodic, they have nonzero values for $n\ge N$ and thus the multiplication of these two DFTs will be nonzero for $n\ge N$ . We need to define a new type of convolution operation that will result in our convolved signal being zerooutside of the range $n=\{0, 1, \dots , N-1\}$ . This idea led to the development of circular convolution , also called cyclic or periodic convolution.

## Signal circular convolution

Given a signal $f(n)$ with Fourier coefficients ${c}_{k}$ and a signal $g(n)$ with Fourier coefficients ${d}_{k}$ , we can define a new signal, $v(n)$ , where $v(n)=⊛(f(n), g(n))$ We find that the Fourier Series representation of $v(n)$ , ${a}_{k}$ , is such that ${a}_{k}={c}_{k}{d}_{k}$ . $⊛(f(n), g(n))$ is the circular convolution of two periodic signals and is equivalent to the convolution over one interval, i.e. $⊛(f(n), g(n))=\sum_{n=0}^{N} \sum_{\eta =0}^{N} f(\eta )g(n-\eta )$ .

Circular convolution in the time domain is equivalent to multiplication of the Fourier coefficients.
This is proved as follows
${a}_{k}=\frac{1}{N}\sum_{n=0}^{N} v(n)e^{-(j{\omega }_{0}kn)}=\frac{1}{N^{2}}\sum_{n=0}^{N} \sum_{\eta =0}^{N} f(\eta )g(n-\eta )e^{-(\omega {j}_{0}kn)}=\frac{1}{N}\sum_{\eta =0}^{N} f(\eta )\frac{1}{N}\sum_{n=0}^{N} g(n-\eta )e^{-(j{\omega }_{0}kn)}=\forall \nu , \nu =n-\eta \colon \frac{1}{N}\sum_{\eta =0}^{N} f(\eta )\frac{1}{N}\sum_{\nu =-\eta }^{N-\eta } g(\nu )e^{-(j{\omega }_{0}(\nu +\eta ))}=\frac{1}{N}\sum_{\eta =0}^{N} f(\eta )\frac{1}{N}\sum_{\nu =-\eta }^{N-\eta } g(\nu )e^{-(j{\omega }_{0}k\nu )}e^{-(j{\omega }_{0}k\eta )}=\frac{1}{N}\sum_{\eta =0}^{N} f(\eta )d_{k}e^{-(j{\omega }_{0}k\eta )}={d}_{k}\frac{1}{N}\sum_{\eta =0}^{N} f(\eta )e^{-(j{\omega }_{0}k\eta )}={c}_{k}{d}_{k}$

## Circular convolution formula

What happens when we multiply two DFT's together, where $Y(k)$ is the DFT of $y(n)$ ?

$Y(k)=F(k)H(k)$
when $0\le k\le N-1$

Using the DFT synthesis formula for $y(n)$

$y(n)=\frac{1}{N}\sum_{k=0}^{N-1} F(k)H(k)e^{(j\frac{2\pi }{N}kn)}$

And then applying the analysis formula $F(k)=\sum_{m=0}^{N-1} f(m)e^{(-j\frac{2\pi }{N}kn)}$

$y(n)=\frac{1}{N}\sum_{k=0}^{N-1} \sum_{m=0}^{N-1} f(m)e^{(-j\frac{2\pi }{N}kn)}H(k)e^{(j\frac{2\pi }{N}kn)}=\sum_{m=0}^{N-1} f(m)\frac{1}{N}\sum_{k=0}^{N-1} H(k)e^{(j\frac{2\pi }{N}k(n-m))}$
where we can reduce the second summation found in the above equation into $h({\left(\left(n-m\right)\right)()}_{N})=\frac{1}{N}\sum_{k=0}^{N-1} H(k)e^{(j\frac{2\pi }{N}k(n-m))}$ $y(n)=\sum_{m=0}^{N-1} f(m)h({\left(\left(n-m\right)\right)()}_{N})$ which equals circular convolution! When we have $0\le n\le N-1$ in the above, then we get:
$y(n)\equiv ⊛(f(n), h(n))$
The notation $⊛$ represents cyclic convolution "mod N".

## Alternative circular convolution algorithm

• Step 1: Calculate the DFT of $f(n)$ which yields $F(k)$ and calculate the DFT of $h(n)$ which yields $H(k)$ .
• Step 2: Pointwise multiply $Y(k)=F(k)H(k)$
• Step 3: Inverse DFT $Y(k)$ which yields $y(n)$

Seems like a roundabout way of doing things, but it turns out that there are extremely fast ways to calculate the DFT of a sequence.

To circularily convolve $2$ $N$ -point sequences: $y(n)=\sum_{m=0}^{N-1} f(m)h({\left(\left(n-m\right)\right)()}_{N})$ For each $n$ : $N$ multiples, $N-1$ additions

$N$ points implies $N^{2}$ multiplications, $N(N-1)$ additions implies $O(N^{2})$ complexity.

## Steps for circular convolution

We can picture periodic sequences as having discrete points on a circle as the domain

Shifting by $m$ , $f(n+m)$ , corresponds to rotating the cylinder $m$ notches ACW (counter clockwise). For $m=-2$ , we get a shift equal to that in the following illustration:

To cyclic shift we follow these steps:

1) Write $f(n)$ on a cylinder, ACW

2) To cyclic shift by $m$ , spin cylinder m spots ACW $\to (f(n), f({\left(\left(n+m\right)\right)}_{N}))$

## Notes on circular shifting

$f({\left(\left(n+N\right)\right)}_{N})=f(n)$ Spinning $N$ spots is the same as spinning all the way around, or not spinning at all.

$f({\left(\left(n+N\right)\right)}_{N})=f({\left(\left(n-\left(N-m\right)\right)\right)}_{N})$ Shifting ACW $m$ is equivalent to shifting CW $N-m$

$f({\left(\left(-n\right)\right)}_{N})$ The above expression, simply writes the values of $f(n)$ clockwise.

## Convolve (n = 4)

• $h({\left((\left(, -, m, \right), \right))}_{N})$

Multiply $f(m)$ and $\mathrm{sum}$ to yield: $y(0)=3$

• $h({\left((\left(, 1, -, m, \right), \right))}_{N})$

Multiply $f(m)$ and $\mathrm{sum}$ to yield: $y(1)=5$

• $h({\left((\left(, 2, -, m, \right), \right))}_{N})$

Multiply $f(m)$ and $\mathrm{sum}$ to yield: $y(2)=3$

• $h({\left((\left(, 3, -, m, \right), \right))}_{N})$

Multiply $f(m)$ and $\mathrm{sum}$ to yield: $y(3)=1$

## Exercise

Take a look at a square pulse with a period of T.

For this signal ${c}_{k}=\begin{cases}\frac{1}{N} & \text{if k=0}\\ \frac{1}{2}\frac{\sin (\frac{\pi }{2}k)}{\frac{\pi }{2}k} & \text{otherwise}\end{cases}$

Take a look at a triangle pulse train with a period of T.

This signal is created by circularly convolving the square pulse with itself. The Fourier coefficients for this signal are ${a}_{k}={c}_{k}^{2}=\frac{1}{4}\frac{\sin (\frac{\pi }{2}k)^{2}}{(\frac{\pi }{2}k)^{2}}$

Find the Fourier coefficients of the signal that is created when the square pulse and the triangle pulse are convolved.

${a}_{k}=\left\{\begin{array}{cc}\text{undefined}\hfill & k=0\hfill \\ \frac{1}{8}\frac{si{n}^{3}\left[\frac{\pi }{2}k\right]}{{\left[\frac{\pi }{2}k\right]}^{3}}\hfill & \text{otherwise}\hfill \end{array}\right)$

## Circular shifts and dft

If $f(n)\stackrel{\text{DFT}}{↔}F(k)$ then $f({\left(\left(n-m\right)\right)}_{N})\stackrel{\text{DFT}}{↔}e^{-(i\frac{2\pi }{N}km)}F(k)$ ( i.e. circular shift in time domain = phase shift in DFT)

$f(n)=\frac{1}{N}\sum_{k=0}^{N-1} F(k)e^{i\frac{2\pi }{N}kn}$
so phase shifting the DFT
$f(n)=\frac{1}{N}\sum_{k=0}^{N-1} F(k)e^{-(i\frac{2\pi }{N}kn)}e^{i\frac{2\pi }{N}kn}=\frac{1}{N}\sum_{k=0}^{N-1} F(k)e^{i\frac{2\pi }{N}k(n-m)}=f({\left(\left(n-m\right)\right)}_{N})$

## Conclusion

Circular convolution in the time domain is equivalent to multiplication of the Fourier coefficients in the frequency domain.

List and explain four factors of production
capital labour entrepreneur natural resources
Thembi
What is supply
when the supply decreases demand also decreases
Thembi
types of demand and the explanation
what is demand
other things remaining same if demend is increases supply is also decrease and if demend is decrease supply is also increases is called the demand
Mian
if the demand increase supply also increases
Thembi
you are wrong this is the law of demand and not the definition
Tarasum
Demand is the willingness of buy and ability to buy in a specific time period in specific place. Mian you are saying law of demand but not in proper way. you have to keep studying more. because its very basic things in Economics.
Hamza
Demand is the price of Quantity goods and services in which consumer's are willing and able to offer at a price in the market over a period of time
Umar
Demand is the quantity of goods and services that the consumer are willing and able to buy at a alternative prices over a given period of time. But mind you demand is quite different from need and want.
Tarasum
Demand can be defined as the graphical representation between price&demand
alkasim
sorry demand is nt a graphical representation between price and quantity demand but instead that is demand curve.
Ebrima
Demand is the willingness and ability of a consumer to buy a quantity of a good over a given period of time assuming all other things remain constant.
Vedaant
what is commercialization?
How to talk loan for bank?
what is the meaning of gpa?
Answer: GPA stands for Grade Point Average. It is a standard way of measuring academic achievement in the U.S. Basically, it goes as follows: Each course is given a certain number of "units" or "credits", depending on the content of the course.
Yusuf
what is small and Microbuisenes
What is fiscal policy
Dansofo
Who is the funder of Economic
Dansofo
founder , that is Adam Smith
Daniel
what is model
The wealth of Nations
the wealth of nations, is it the first?
Umar
Yes very sure it was released in 1759
Yusuf
thank you Yusuf.
Umar
then when did he died?
Umar
17 July 1790 Born: 16 June 1723, Kirkcaldy, United Kingdom Place of death: Panmure House, Edinburgh, United Kingdom
Yusuf
1790
Yusuf
that's my today questions, thank you Yusuf it's bed time see u after.
Umar
what is fiscal policy
what's mode?
mode is the highest occurring frequency in a distribution
Bola
mode is the most commonly occurring item in a set of data.
Umar
Please, what is the difference between monopoly and monopsony?
is there monopsony word?
Umar
I have no idea though
Umar
Umar
monopsony is when there's only one buyer while monopoly is when there's only one producer.
Bola
who have idea on Banter
Ibrahim
Bola
Monopoly is when there's excessively one seller and there is no entry in the market while monopsony is when there is one buyer
kemigisha
Adam smith was born in 1723
Bola
(uncountable) Good humoured, playful, typically spontaneous conversation. verb (intransitive) To engage in banter or playful conversation. (intransitive) To play or do something amusing. (transitive) To tease mildly.
Umar
which book Adam smith published first? the first book of Adam smith pls.
Umar
wealth on nation, 1776
Daniel
what is market power and how can it affect an economy?
market power:- where a firm is said to be a price setter.market power benefits the powerful at the expense of others.
Umar
Market power refers to the ability of a firm (or group of firms) to raise and maintain price above the level that would prevail under competition is referred to as market or monopoly power. The exercise of market power leads to reduced output and loss of economic welfare
Kartheek
find information about the national budget
Molahlegi
three branches of economics in which tourism is likely to figure
What are those three branches?
IlRegno
Difference between extinct and extici spicies
in a comparison of the stages of meiosis to the stage of mitosis, which stages are unique to meiosis and which stages have the same event in botg meiosis and mitosis
Researchers demonstrated that the hippocampus functions in memory processing by creating lesions in the hippocampi of rats, which resulted in ________.
The formulation of new memories is sometimes called ________, and the process of bringing up old memories is called ________.
Got questions? Join the online conversation and get instant answers!    By By     By Rhodes