Digital Signal Processing Reference
In-Depth Information
Fig. 1.32. M
ATLAB
plot for
Example 1.23. (a)
x
1
(
t
);
(b)
x
2
(
t
); (c) Re
x
3
(
t
);
(d) Im
x
3
(
t
).
5
40
4
30
3
20
2
10
1
0
0
−1
−10
−2
−20
−3
−30
−4
−5
−40
−5
−4
−3
−2
−1
−10
−8
−6
−4
−2
1 0
(a)
(b)
time (
t
)
time (
t
)
real component
imaginary component
1
1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
0
−0.2
−0.2
−0.4
−0.4
−0.6
−0.6
−0.8
−0.8
−1
−1
−5
1 0
1 5
−5
1 0
1 5
(c)
(d)
time (
t
)
time (
t
)
>> % plot function 1
>> subplot(2,2,1), stem(k, f1, 'filled'), grid
>> xlabel('k')
>> ylabel('-9.2sin(0.1
\
pi k-0.75
\
pi')
>> title('Part (a)')
>> %%%%%%%%%%%%
>> % Part(b) %
>> %%%%%%%%%%%%
>> k = [-5:25];
>>f2=2*1.1.ˆ(-1.8*k) - 2.1 * 0.9.ˆ(0.7*k);
>> subplot(2,2,2), stem(k, f2, 'filled'), grid
>> xlabel('k')
>> ylabel('2(1.1)ˆ
{
-1.8k
}
- 2.1(0.9)ˆ0.7k')
>> title('Part (b)')
>> %%%%%%%%%%%%
>> % Part(c) %
>> %%%%%%%%%%%%
>> k = [0:50];
>> f3 = (-0.93).ˆk .* exp(j*pi*k/sqrt(350));
>> subplot(2,2,3), stem(k, real(f3), 'filled'), grid
Search WWH ::
Custom Search