Tín hiệu và hệ thống
93
dài chuỗi xung là 10ms, tốc độ lặp lại của chuỗi xung Gauss là 1kHz, suy hao theo hàm mũ cơ
số 0.8. Xung Gauss có tần số trung tâm là 10kHz, băng thông 50%.
T = 0:1/50E3:10E-3; % vector thời gian của chuỗi xung (0 – 10ms)
D = [0:1/1E3:10E-3;0.8.^(0:10)]'; % cột 1 của D xác định các thời điểm lặp
% cột 2 của D xác định biên độ tương ứng của xung Gauss (bị suy hao)
Y = pulstran(T,D,'gauspuls',10E3,0.5); % gọi hàm pulstran, hai thông số
cuối
% là tần số trung tâm và tỷ lệ băng thông của xung Gauss
plot(T,Y) % vẽ tín hiệu
Kết quả như sau:
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
0.009
0.01
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
Hình 9.2
Hàm sinc(x): là biến đổi Fourier ngược của xung chữ nhật có chiều rộng băng 2π và
chiều cao bằng 1:
x
x
x
c
π
π )
sin(
)
(
sin
=
(9.1)
x có thể là một vector hay một ma trận.
Hàm Dirichlet: còn gọi là hàm sinc tuần hoàn hay hàm sinc chồng lấn (aliased sinc), được
định nghĩa như sau:
)
(x
D
n
=
⎪
⎩
⎪
⎨
⎧
∈
=
−
−
otherwise
x
n
nx
Z
k
k
x
n
k
)
2
/
sin(
)
2
/
sin(
,
2
)
1
(
)
1
(
π
(9.2)
Trong MATLAB, để thực hiện hàm Dirichlet ta gọi hàm diric(x,n).