blob: 4b3b29de3189a1d742d3575fb458b574549fcf7f (
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
|
//EXAMPLE 3.13
clc;
clear;
a=0.5;
n=0:9;
for i = 0:9
x(i+1) = a^i;
end
//The DTFT of the sequence
K = 4;
k = 0:4/1000:4;
W = k*6*%pi/K;
X1 = (x')*exp(%i*n'*W);
X = %i*diff(X1);
X = [X,0] + X1;
X_mag = abs(X);
[X_phase,db] = phasemag(X);
clf();
a=gca();
figure(0);
plot2d3(mtlb_fliplr(W),X_mag);
xtitle('Magnitude response','W','X_mag');
figure(1);
plot2d3(mtlb_fliplr(W),X_phase);
xtitle('Phase response','W','X_phase');
|
