Tổng quan về giáo trình

Giáo trình thuộc học phần chuyên ngành trong chương trình đào tạo kỹ sư và cử nhân các ngành Kỹ thuật Điện tử - Viễn thông, Kỹ thuật Điều khiển & Tự động hóa, và Kỹ thuật Máy tính. Trích đoạn tài liệu thuộc Phần II: Ứng dụng MATLAB trong Xử lý tín hiệu và Xử lý ảnh, tập trung trực tiếp vào hai nội dung nền tảng: Chương 9: Tín hiệu và hệ thống và Chương 10: Thiết kế các bộ lọc.

                           KHUNG NỘI DUNG GIÁO TRÌNH
   CHƯƠNG 9: TÍN HIỆU VÀ HỆ THỐNG                              CHƯƠNG 10: THIẾT KẾ CÁC BỘ LỌC

Mục tiêu học tập của giáo trình nhằm trang bị cho người học:

  • Phương pháp biểu diễn, tạo lập và thao tác trên các dạng tín hiệu một chiều, đa kênh, tín hiệu tuần hoàn và không tuần hoàn trong môi trường tính toán số.
  • Kỹ thuật phân tích hệ thống tuyến tính dừng (LTI) thông qua đáp ứng xung, hàm truyền đạt miền $z$, phương trình sai phân, và giản đồ cực - zero.
  • Khả năng chuyển đổi qua lại giữa 7 dạng mô hình hệ thống rời rạc và các mô hình hệ thống liên tục.
  • Nguyên lý và thuật toán thiết kế các bộ lọc số đáp ứng xung vô hạn (IIR) và đáp ứng xung hữu hạn (FIR) đáp ứng chính xác các chỉ tiêu kỹ thuật trong miền tần số.

Giáo trình tiếp cận theo hướng tích hợp chặt chẽ giữa lý thuyết giải tích và công cụ mô phỏng số thông qua thư viện MATLAB Signal Processing Toolbox. Mỗi chủ đề lý thuyết đều gắn liền với thuật toán tính toán cụ thể, cú pháp hàm chuẩn, mã nguồn mẫu và đồ thị kết quả trực quan. Điểm đặc sắc của tài liệu là tính hệ thống cao trong việc phân loại mô hình toán học và bao quát đầy đủ các thuật toán thiết kế từ kinh điển (Butterworth, Chebyshev, Elliptic) đến hiện đại (Parks-McClellan, Constrained Least Squares, Mô hình tham số).

Nội dung kiến thức cốt lõi

Các chương và chủ đề chính

Nội dung trích đoạn được tổ chức thành 2 chương trọng tâm với trình tự logic đi từ biểu diễn, mô hình hóa hệ thống đến tổng hợp và thiết kế bộ lọc:

1. Chương 9: Tín hiệu và hệ thống

  • Biểu diễn và tạo tín hiệu: Biểu diễn tín hiệu rời rạc dưới dạng vector hàng ($1 \times n$), vector cột ($n \times 1$), và tín hiệu $m$ kênh dưới dạng ma trận ($n \times m$). Tạo tín hiệu tuần hoàn bằng hàm square, sawtooth, tín hiệu aperiodic bằng gauspuls, chirp, tín hiệu chuỗi xung pulstran, và hàm Dirichlet diric.
  • Thao tác dữ liệu: Truy xuất dữ liệu qua các hàm cấp thấp (fopen, fread, fscanf), nạp tệp qua load, và tích hợp mã nguồn ngoài qua MEX file.
  • Phân tích hệ thống và cấu trúc lọc: Thực thi tích chập 1D và 2D (conv, conv2). Biểu diễn hàm truyền đạt miền $z$: $$H(z) = \frac{b(1) + b(2)z^{-1} + \dots + b(n+1)z^{-n}}{a(1) + a(2)z^{-1} + \dots + a(m+1)z^{-m}}$$ Thực hiện lọc theo cấu trúc trực tiếp dạng II (Direct-Form II) bằng hàm filter. Các kỹ thuật mở rộng gồm lọc khử méo pha 2 chiều bằng filtfilt, băng lọc đa tốc độ (Multirate polyphase filter bank) bằng upfirdn, lọc miền tần số bằng giải thuật Overlap-and-Add qua fftfilt, và lọc cấu trúc mạng bằng latcfilt.
  • Đáp ứng hệ thống và giản đồ cực - zero: Phân tích đáp ứng tần số miền rời rạc (freqz), miền liên tục (freqs), đáp ứng pha liên tục (unwrap, phasez), độ trễ nhóm $\tau_g(\omega) = -\frac{d\theta(\omega)}{d\omega}$ (grpdelay), độ trễ pha (phasedelay), và giản đồ phân bố cực - zero trên mặt phẳng $z$ (zplane).
  • Hệ thống 7 mô hình toán học rời rạc: Mô hình hàm truyền đạt ($b, a$), mô hình độ lợi - cực - zero ($k, p, z$), mô hình không gian trạng thái ($A, B, C, D$), mô hình khai triển phân thức thặng dư (residuez), mô hình các khâu bậc hai SOS ($L \times 6$ matrix), mô hình mạng Lattice (tf2latc, latc2tf), và mô hình ma trận chập Toeplitz (convmtx).
  • Biến đổi Fourier rời rạc: Cơ chế tính DFT/IDFT 1 chiều và 2 chiều qua giải thuật FFT (fft, ifft, fft2, ifft2), thuật toán Goertzel (goertzel), và dịch chuyển phổ (fftshift).

2. Chương 10: Thiết kế các bộ lọc

  • Chỉ tiêu thiết kế: Xác định dải thông ($\omega_p$), dải chắn ($\omega_s$), gợn sóng dải thông ($\delta_1, R_p$), suy hao dải chắn ($\delta_2, R_s$), và độ rộng dải chuyển tiếp trên thang tần số chuẩn hóa Nyquist ($f_s/2$).
  • Phương pháp thiết kế IIR:
    • Dựa trên nguyên mẫu analog: Thiết kế nguyên mẫu thông thấp chuẩn hóa qua các bộ lọc Butterworth (buttap), Chebyshev loại I (cheb1ap), Chebyshev loại II (cheb2ap), Elliptic (ellipap), Bessel (besselap). Chuyển đổi dải tần analog (lp2lp, lp2hp, lp2bp, lp2bs). Rời rạc hóa sang miền $z$ bằng phương pháp bất biến xung (impinvar) hoặc biến đổi song tuyến tính (bilinear) có hiệu chỉnh tần số prewarping.
    • Hàm thiết kế hoàn chỉnh: Tự động hóa ước lượng bậc tối thiểu và tạo hệ số qua buttord/butter, cheb1ord/cheby1, cheb2ord/cheby2, ellipord/ellip.
    • Thiết kế trực tiếp & mô hình tham số: Thiết kế trực tiếp miền số bằng thuật toán giải hệ phương trình Yule-Walker (yulewalk), bộ lọc Butterworth tổng quát với số cực và zero độc lập (maxflat), mô hình tham số miền thời gian qua mã hóa dự đoán tuyến tính (lpc), giải thuật Prony (prony), giải thuật lặp Steiglitz-McBride (stmcb), và mô hình tối ưu hóa miền tần số (invfreqs, invfreqz).
  • Phương pháp thiết kế FIR:
    • Phân loại 4 kiểu pha tuyến tính: Định nghĩa cấu trúc đối xứng/phản đối xứng chẵn - lẻ của các hệ số (tap) cùng đặc tính triệt tiêu tần số tại $\omega=0$ và $\omega=\pi$, bảo đảm độ trễ nhóm cố định $n/2$.
    • Phương pháp cửa sổ (Windowing): Giới hạn đáp ứng xung lý tưởng sinc bằng các hàm cửa sổ: Rectangular (rectwin), Triangular (triang), Bartlett-Hanning (barthannwin), Blackman (blackman), Blackman-Harris (blackmanharris), Bohman (bohmanwin), Chebyshev (chebwin), Gaussian (gausswin), Parzen (parzenwin), Tukey (tukeywin), Flat top (flattopwin), Nuttall (nuttallwin). Ứng dụng hàm fir1, fir2, và ước lượng tham số cửa sổ Kaiser qua kaiserord.
    • Phương pháp tối ưu hóa: Thiết kế bộ lọc gợn sóng cân bằng (equiripple minimax) qua thuật toán Parks-McClellan (firpm), tối thiểu hóa bình phương sai số tích lũy (firls), thiết kế bộ biến đổi Hilbert và bộ vi phân.
    • Phương pháp bình phương cực tiểu có ràng buộc (CLS): Khống chế sai số gợn sóng tuyệt đối trong toàn bộ dải tần và triệt tiêu hiện tượng Gibbs bằng hàm fircls và fircls1.

Kiến thức nền tảng được xây dựng

Giáo trình củng cố một hệ thống lý thuyết toán học và xử lý tín hiệu hoàn chỉnh:

Khung kiến thức nền tảng Nguyên lý cốt lõi Công cụ thực thi trong MATLAB
Đại số sai phân & Tích chập Phép chập tuyến tính, tích chập 2 chiều, cấu trúc ma trận Toeplitz conv, conv2, convmtx
Không gian trạng thái rời rạc & liên tục Phương trình ma trận hệ thống $\mathbf{x}(n+1) = \mathbf{A}\mathbf{x}(n) + \mathbf{B}\mathbf{u}(n)$, $\mathbf{y}(n) = \mathbf{C}\mathbf{x}(n) + \mathbf{D}\mathbf{u}(n)$ ss2tf, tf2ss, zp2ss, ss2sos
Khai triển thặng dư (Residue) Phân tích phân thức hữu tỷ thành tổng các cực đơn và cực bội residuez
Cấu trúc lọc ghép tầng bậc hai (SOS) Phân tích $H(z) = \prod_{k=1}^L H_k(z)$ giảm nhạy sai số lượng tử sos2tf, tf2sos, zp2sos
Cấu trúc mạng (Lattice) Tính toán hệ số phản xạ $k(i)$ và hệ số bậc thang $v(i)$ tf2latc, latc2tf, latcfilt
Biến đổi Song tuyến tính (Bilinear) Ánh xạ bảo toàn tính ổn định $s = k\frac{z-1}{z+1}$ với kỹ thuật tiền biến dạng tần số bilinear

Kỹ năng phát triển

  1. Kỹ năng kỹ thuật (Technical Skills):
    • Thành thạo lập trình script và hàm chức năng trên MATLAB Signal Processing Toolbox.
    • Kỹ năng xử lý dữ liệu nhị phân, ASCII cấp thấp và liên kết MEX file.
    • Khả năng tính toán phổ DFT/FFT nhiều chiều và phân tích độ trễ nhóm/pha.
  2. Kỹ năng phân tích (Analytical Skills):
    • Đánh giá tính ổn định và tính nhân quả của hệ thống dựa trên vị trí cực - zero và đáp ứng xung.
    • So sánh định lượng giữa các cấu hình bộ lọc IIR (Butterworth phẳng tối đa, Chebyshev tối ưu dải chuyển tiếp, Elliptic tối thiểu hóa bậc, Bessel pha tuyến tính analog).
    • Phân tích hiện tượng méo pha phi tuyến và hiện tượng dao động Gibbs trong thiết kế lọc số.
  3. Kỹ năng thiết kế và thực thi (Practical Competencies):
    • Tổng hợp hệ thống lọc số hoàn chỉnh từ bảng chỉ tiêu kỹ thuật ($R_p, R_s, \omega_p, \omega_s$).
    • Triển khai các cấu trúc lọc phần cứng chuyên dụng như mạng Lattice và chuỗi SOS nhằm hạn chế hiện tượng tràn số trên hệ thống số dấu chấm cố định (Fixed-point).

Phương pháp giảng dạy và học tập

Giáo trình được xây dựng theo phương pháp sư phạm diễn dịch kết hợp thực nghiệm số (Computational Experimentation):

Hệ thống bài tập và nghiên cứu tình huống

Tài liệu cung cấp chuỗi bài tập thực hành cuối chương (từ Bài tập 9-1 đến 9-14) với độ phức tạp tăng dần:

  • Kiểm chứng tính chất đại số: Chứng minh tính giao hoán, kết hợp, phân phối của phép tích chập qua hàm conv với các dãy tín hiệu $x_1(n) = n[u(n+10)-u(n-20)]$, $x_2(n) = \cos(0.1\pi n)[u(n)-u(n-30)]$, $x_3(n) = 1.2^n[u(n+5)-u(n-10)]$.
  • Khảo sát hệ thống sai phân: Đánh giá đáp ứng bộ sai phân số $y(n) = x(n) - x(n-1)$ trên các dạng sóng xung chữ nhật, xung tam giác, xung sin; đối sánh kết quả khi lọc qua filter và filtfilt.
  • Đánh giá hiệu năng truyền thông: Mô phỏng hệ thống với các giá trị tỷ số năng lượng trên bit trên mật độ phổ công suất nhiễu ($E_b/N_0$) từ $0\text{ dB}$ đến $8\text{ dB}$, vẽ đồ thị tỷ lệ lỗi bit (BER).
  • Phát triển thuật toán độc lập: Tự xây dựng hàm tích chập dạng ma trận Toeplitz function [y,H] = conv_tp(h,x).
  • Khảo sát độ phức tạp tính toán: Khảo sát thời gian thực thi của biến đổi FFT với kích thước $N = 4^v$ ($v = 5, 6, \dots, 10$) để chứng minh sự phụ thuộc thời gian tỷ lệ với $N \log_4 N$.

Đánh giá và hướng dẫn tự học

  • Phương pháp đánh giá: Đánh giá dựa trên mã nguồn mô phỏng, độ chính xác của các tham số hàm truyền đạt, tính tối ưu của bậc bộ lọc và các đồ thị phân tích miền thời gian/tần số.
  • Chỉ dẫn tự học: Người học cần cài đặt từng đoạn mã lệnh mẫu trong giáo trình, thay đổi các tham số đầu vào (bậc bộ lọc $n$, độ gợn sóng $R_p$, suy hao $R_s$), quan sát sự dịch chuyển của các điểm cực - zero trên mặt phẳng $z$ bằng lệnh zplane để hiểu rõ cơ chế vật lý của hệ thống.

Điểm nổi bật và cập nhật

  • Tích hợp các thuật toán thiết kế hiện đại: Bên cạnh các phương pháp cổ điển, giáo trình cập nhật giải thuật xấp xỉ tối ưu hóa Parks-McClellan (Remez Exchange Algorithm) qua hàm firpm, thuật toán bình phương cực tiểu có ràng buộc fircls/fircls1, và phương pháp thiết kế bộ lọc nhiều dải tần bất kỳ qua thuật toán Yule-Walker yulewalk.
  • Giải quyết các vấn đề méo phi tuyến trong kỹ thuật lọc:
    • Khử méo pha cho bộ lọc IIR: Trình bày giải thuật lọc 2 chiều (thuận và đảo thời gian) thông qua hàm filtfilt, triệt tiêu hoàn toàn góc pha lệch mà vẫn đảm bảo độ suy hao dải chắn.
    • Khắc phục méo tần số phi tuyến: Phân tích cơ chế biến dạng tần số của phép biến đổi song tuyến tính và thiết lập công thức tiền biến dạng tần số (prewarping): $$k = \frac{2\pi f_p}{\tan(\pi f_p / f_s)}$$
  • Mô hình hóa phần cứng và tối ưu hóa số học: Trình bày chi tiết mô hình khâu bậc hai SOS (Second Order Sections) và cấu trúc mạng Lattice. Đây là các cấu trúc chuẩn để hiện thực hóa các bộ lọc số trên các chip vi xử lý DSP chuyên dụng hoặc FPGA, giúp giảm thiểu sai số lượng tử hóa hệ số và tránh tràn số.
  • Kỹ thuật xử lý tín hiệu đa tốc độ (Multirate): Tích hợp kỹ thuật lọc đa pha (polyphase filtering) thông qua hàm upfirdn, cho phép thay đổi tốc độ lấy mẫu theo tỷ lệ hữu tỷ $P/Q$.

Đối tượng sử dụng giáo trình

  • Sinh viên đại học và học viên cao học: Phù hợp cho sinh viên từ năm thứ 3 trở đi thuộc các ngành Kỹ thuật Điện tử - Viễn thông, Tự động hóa, Kỹ thuật Y sinh, Kỹ thuật Máy tính và Công nghệ Thông tin cần công cụ mô phỏng thực hành môn học Tín hiệu và Hệ thống hoặc Xử lý Tín hiệu Số (DSP).
  • Kiến thức tiên quyết (Prerequisites):
    • Toán học kỹ thuật: Đại số tuyến tính (ma trận, định thức), Giải tích phức, Biến đổi Laplace, Biến đổi Fourier và Biến đổi $z$.
    • Lập trình: Cú pháp cơ bản của MATLAB (ma trận, chỉ mục, vẽ đồ thị hàm số).
  • Giảng viên: Sử dụng làm tài liệu giảng dạy, thiết kế bài tập lớn, xây dựng kịch bản thực hành phòng thí nghiệm (Lab), và biên soạn đề thi môn học.
  • Kỹ sư và nhà nghiên cứu: Tài liệu tham khảo tra cứu nhanh cú pháp, giải thuật chuyển đổi mô hình và quy trình tổng hợp bộ lọc số trong các dự án phát triển hệ thống viễn thông, xử lý âm thanh và xử lý ảnh.

Câu hỏi thường gặp

1. Giáo trình này phù hợp với ai?

Tài liệu phù hợp cho sinh viên đại học từ năm thứ 3, học viên cao học chuyên ngành Kỹ thuật Điện tử, Viễn thông, Tự động hóa và các kỹ sư xử lý tín hiệu cần nắm vững phương pháp mô phỏng số trên MATLAB.

2. Cần kiến thức nền nào để học?

Người học cần có kiến thức nền tảng về Toán kỹ thuật (Giải tích hàm phức, Biến đổi $z$, Chuỗi Fourier, Đại số tuyến tính) và nguyên lý cơ bản của Lý thuyết Mạch hoặc Tín hiệu & Hệ thống.

3. Điểm khác biệt với giáo trình khác?

Giáo trình kết nối trực tiếp lý thuyết giải tích với việc hiện thực hóa mã lệnh thông qua MATLAB Signal Processing Toolbox. Tài liệu không dừng lại ở mức mô tả lý thuyết mà đi sâu vào cấu trúc dữ liệu, các mô hình toán tương đương (SOS, State-space, Lattice) và các giải thuật tối ưu hóa thiết kế bộ lọc.

4. Làm sao để tự học hiệu quả?

Người học nên chạy từng ví dụ mẫu trong giáo trình trên phần mềm MATLAB, khảo sát sự biến thiên của đáp ứng biên độ và pha qua các hàm vẽ trực quan (plot, freqz, zplane), đồng thời tự giải quyết toàn bộ các bài tập lập trình từ 9-1 đến 9-14 ở cuối chương.

5. Có tài liệu bổ trợ nào kèm theo?

Tài liệu cung cấp bảng tra cứu toàn diện hệ thống các hàm xử lý tín hiệu, bảng chuyển đổi đa mô hình hệ thống tuyến tính và các bài toán mô phỏng hệ thống truyền thông thực tế.

Kết luận

Trích đoạn giáo trình trong Phần II: Ứng dụng MATLAB trong Xử lý tín hiệu và Xử lý ảnh cung cấp một hệ thống tri thức hoàn chỉnh, chuẩn xác và mang tính ứng dụng cao về phân tích hệ thống tuyến tính và thiết kế bộ lọc số.

Lộ trình học tập đề xuất:

  1. Nắm vững kỹ thuật tạo lập, biểu diễn tín hiệu và các phép toán cơ bản (tích chập, lọc số Direct-Form II, FFT).
  2. Làm chủ 7 mô hình toán học của hệ thống LTI rời rạc và các hàm chuyển đổi tương đương.
  3. Thực hành thiết kế các bộ lọc số IIR và FIR theo chỉ tiêu kỹ thuật bằng các phương pháp cửa sổ, gợn sóng cân bằng và mô hình tham số.
  4. Triển khai các thuật toán mở rộng (băng lọc đa tốc độ, lọc khử méo pha) vào bài toán kỹ thuật thực tế.

Tài liệu là tài nguyên học thuật nền tảng phục vụ đắc lực cho công tác đào tạo, nghiên cứu và phát triển trong lĩnh vực xử lý tín hiệu số.