blob: 35c37dc004c036fd72eacbc6056cde55524d6564 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
//EXAMPLE 3.12
//x[n]=((-1)^n)*(a^n)*u[n].....given a=0.5;
clc;
clear;
a=0.5;
n=0:9;
for i = 0:9
x(i+1) = (a*exp(-%i*%pi))^i;
end
//The DTFT of the sequence
K = 4;
k = 0:4/1000:4;
W = k*6*%pi/K;
X = (x')*exp(%i*n'*W);
X_mag = abs(X);
X_phase = phasemag(X);
//PLOTTING GRAPHS FOR THE INTERVAL OF 0 TO 6*%pi
clf();
a=gca();
figure(0);
plot2d3(mtlb_fliplr(W),X_mag);
xtitle('Magnitude response','W','Amplitude');
figure(1);
plot2d3(mtlb_fliplr(W),X_phase);
xtitle('Phase response','W','X_phase,degrees');
|
