Bộ Môn Viễn Thông-ĐH Bách
Khoa TpHCM
TÀI LIỆU THÍ NGHIỆM XỬ LÝ SỐ TÍN HIỆU
68
%multibandfir63.m: Multiband FIR filter with 63 coefficients
f = [0 0.1 0.12 0.18 0.2 0.3 0.32 0.38 0.4 1];
m = [0 0 1 1 0 0 1 1 0 0];
n = 63;
cof = remez(n-1,f,m);
% frequency response with 256
points
[h w] = freqz(cof,1,256);
% plot magnitude of
the filter
plot(5000*f,m);
figure;
plot(w/pi,abs(h));
Trong đoạn chương trình trên, có một số lưu ý:
- Dòng lệnh
cof = remez(n-1,f,m) trả về vector hệ số của bộ lọc FIR bậc n – 1, với f và
m xác định các dải tần số theo bảng ở trên.
- Lệnh freqz để tính đáp ứng tần số của bộ lọc
- Lệnh plot thứ nhất vẽ đáp ứng tần số
mong muốn dựa trên f và m.
- Lệnh figure tạo ra một cửa sổ mới và lệnh plot thứ hai vẽ đáp ứng tần số của bộ lọc
đã
thiết kế được lên cửa sổ mới này.
Ở đây, kết quả của quá trình thiết kế mà ta cần nhận được chính là các hệ số chứa trong biến
cof. Chúng được sử dụng khi thực hiện bộ lọc lên trên kit DSP (lưu ý đưa về định dạng 16 bit
có dấu khi thực hiện trên kit DSP). Tiếp theo thực hiện tương tự như phần II.1.1.
1. Đáp ứng tần số của bộ lọc thiết kế dùng MATLAB:
Bộ Môn Viễn Thông-ĐH Bách Khoa TpHCM
TÀI LIỆU THÍ NGHIỆM XỬ LÝ SỐ TÍN HIỆU
72
Chia sẻ với bạn bè của bạn: