blob: b03bac3ddaecfa5031f115ae1f51d207e77a855a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
//compute circular convolution of x1(n)={-1 2 -2 1} & x2(n)={1 -2 -1 2}
clc;
clear;
n=0:1:3;//for four point DFT
x1=[-1 2 -2 1];
disp(x1,'Given signal=> x1(n)=');
X1=dft(x1,-1);
disp(X1,'four point DFT of x(n)=>');
disp(abs(X1),'magnitude of four point DFT x(n)=>');
disp(atan(imag(X1),real(X1)),'phase of four point DFT x(n)=>');
x2=[1 -2 -1 2];
disp(x2,'Given signal=> x1(n)=');
X2=dft(x2,-1);
disp(X2,'four point DFT of x(n)=>');
disp(abs(X2),'magnitude of four point DFT x(n)=>');
disp(atan(imag(X2),real(X2)),'phase of four point DFT x(n)=>');
X3=X1.*X2;
x3=dft(X3,1);
disp(x3,'by circular convolution x3(n)=>');
|
