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

Tài liệu bài giảng Thiết kế thí nghiệm và Xử lý số liệu (tên tiếng Anh: Data Analysis and Design of Experiment, mã học phần: CH3309) do PGS. Trịnh Văn Dũng biên soạn tại Bộ môn Quá trình và Thiết bị, Khoa Kỹ thuật Hóa học, Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh (phiên bản cập nhật năm 2017). Môn học giữ vị trí cốt lõi trong khối kiến thức kỹ thuật chuyên ngành, cung cấp cơ sở phương pháp luận toán học – thống kê nhằm phục vụ công tác nghiên cứu khoa học, tối ưu hóa công nghệ và phân tích thực nghiệm trong các lĩnh vực Công nghệ Hóa học, Công nghệ Sinh học và Công nghệ Thực phẩm.

Mục tiêu học tập của tài liệu tập trung vào việc trang bị cho người học phương pháp tiếp cận đối tượng nghiên cứu thông qua mô hình "hộp đen" (black-box model), kỹ thuật lập kế hoạch thực nghiệm để cực tiểu hóa số lần thí nghiệm mà vẫn bảo đảm độ tin cậy thống kê, cùng các thuật toán ước lượng hệ số hồi quy đa biến. Cấu trúc bài giảng đi từ các khái niệm cơ bản về thực nghiệm, phân loại mục tiêu nghiên cứu, đến các mô hình quy hoạch trực giao cấp 1 (mô hình tuyến tính và tương tác) và quy hoạch trực giao cấp 2 (mô hình phi tuyến bậc hai mô tả bề mặt đáp ứng).

Điểm đặc thù của tài liệu là tính chặt chẽ về mặt toán học thông qua việc giải thích nguồn gốc của phương trình hồi quy từ chuỗi Taylor, kết hợp phương pháp bình phương cực tiểu (Ordinary Least Squares - OLS) với các điều kiện đại số ma trận nhằm phân tách độc lập các hệ số hồi quy.

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

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

Nội dung tài liệu được tổ chức theo 3 phần trọng tâm với tiến trình logic chặt chẽ:

  1. Khái niệm chung về quy hoạch thực nghiệm và mô hình hộp đen:

    • Định nghĩa bản chất của thực nghiệm: Tác động có chủ đích lên đối tượng nghiên cứu để thu thập dữ liệu tin cậy.
    • Phân loại thực nghiệm theo 3 giai đoạn: Thí nghiệm phòng thí nghiệm (nghiên cứu quy luật chung, kiểm định giả thuyết), thí nghiệm ứng dụng (nghiên cứu quy trình cụ thể, xác định tính chất lý – hóa), và thí nghiệm sản xuất (kiểm tra, đánh giá kết quả tính toán thiết kế).
    • Quy trình 3 bước thực hiện nghiên cứu tích cực: (1) Đặt bài toán và ước tính chi phí, (2) Thu thập thông tin ưu tiên về đối tượng, (3) Lựa chọn phương pháp giải, biến đầu vào, thông số và hàm mục tiêu.
    • Giới thiệu tổng quan hệ thống các phương pháp: Quy hoạch trực giao cấp I, cấp II, quy hoạch riêng phần, tâm xoay, đơn hình, tối ưu và đa mục tiêu.

  2. Cơ sở lý thuyết hồi quy và Quy hoạch trực giao cấp 1:

    • Cơ sở toán học: Khai triển chuỗi Taylor hàm đáp ứng $y = f(x_1, x_2, \dots, x_k)$ tại lân cận điểm cơ sở $x_{i0}$ để xác định dạng phương trình lý thuyết và phương trình thực nghiệm: $$y = b_0 + \sum_{i=1}^k b_i x_i + \sum_{i < j}^k b_{ij} x_i x_j + \dots$$
    • Kỹ thuật chuyển đổi thang đo: Mã hóa biến thực $u_i \in [a_i, b_i]$ sang biến quy ước (biến không thứ nguyên) $x_i \in [-1, +1]$ thông qua các công thức: $$u_{i0} = \frac{a_i + b_i}{2}, \quad \Delta u_i = \frac{b_i - a_i}{2}, \quad x_i = \frac{u_i - u_{i0}}{\Delta u_i}$$
    • Tính chất của ma trận trực giao cấp 1: Điều kiện trực giao giữa các cột biến và cột tương tác ($\sum_{j=1}^N x_{ij} = 0$; $\sum_{j=1}^N x_{ij} x_{mj} = 0$ khi $i \neq m$).
    • Thuật toán tính nhanh hệ số hồi quy không phụ thuộc lẫn nhau bằng tích vô hướng: $$b_j = \frac{1}{N}\sum_{i=1}^N x_{ji} y_i, \quad b_{ij} = \frac{1}{N}\sum_{m=1}^N (x_i x_j)_m y_m$$
    • Ví dụ áp dụng đầy đủ cho hệ 2 biến ($2^2 = 4$ thí nghiệm) và 3 biến ($2^3 = 8$ thí nghiệm) có lặp thí nghiệm để xử lý số liệu trung bình và phương sai.

  3. Quy hoạch trực giao cấp 2 (Second-Order Orthogonal Designs):

    • Sự cần thiết của đa thức bậc hai khi mô tả bề mặt đáp ứng phi tuyến có cực trị: $$y = b_0 + \sum_{i=1}^k b_i x_i + \sum_{i=1}^k b_{ii} x_i^2 + \sum_{i < j}^k b_{ij} x_i x_j$$
    • Cấu tạo không gian điểm thí nghiệm: Bổ sung các điểm sao (axial points) với cánh tay đòn $\alpha$ ($(\pm \alpha, 0, \dots, 0)$) và các điểm lặp ở tâm $n_0$ $(0, 0, \dots, 0)$. Tổng số thí nghiệm xác định theo công thức: $$N = 2^k + 2k + n_0 \quad (\text{với } k < 5) \quad \text{hoặc} \quad N = 2^{k-l} + 2k + n_0 \quad (\text{với } k \ge 5)$$
    • Phương pháp trực giao hóa ma trận: Do $\sum x_{ij}^2 \neq 0$, tài liệu xây dựng phép biến đổi biến bậc hai $x'{ij} = x{ij}^2 - \overline{x_i^2}$ với $\overline{x_i^2} = \frac{1}{N}\sum_{j=1}^N x_{ij}^2 = \frac{2^k + 2\alpha^2}{N}$, giúp ma trận chuyển về dạng trực giao hoàn toàn.
    • Bảng tra cứu giá trị cánh tay đòn $\alpha$ tối ưu phụ thuộc vào số yếu tố $k$ ($k = 2, 3, 4$) và số thí nghiệm lặp ở tâm $n_0$ ($n_0 = 1, 2, 3, 4$).
    • Công thức tính toán các hệ số hồi quy hiệu chỉnh $b_0', b_i, b_{ij}, b_{ii}'$ và chuyển đổi về hệ số thực nghiệm ban đầu.

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

Tài liệu củng cố và tích hợp ba mảng kiến thức nền tảng:

  • Giải tích đa biến: Ứng dụng chuỗi Taylor để xấp xỉ hóa các hàm phi tuyến phức tạp trong kỹ thuật thành các đa thức bậc 1 và bậc 2.
  • Đại số tuyến tính: Ứng dụng tích vô hướng vectơ, tính trực giao của không gian cột trong ma trận thực nghiệm $X$ nhằm giải hệ phương trình chuẩn tắc Gauss trong phương pháp bình phương cực tiểu một cách độc lập.
  • Xác suất thống kê thực nghiệm: Nguyên lý phân tích phương sai, tính toán giá trị trung bình từ các thí nghiệm lặp tại tâm $n_0$, đánh giá ảnh hưởng riêng phần và ảnh hưởng tương tác giữa các yếu tố công nghệ.

Kỹ năng phát triển

  • Kỹ năng thiết kế ma trận thực nghiệm: Khả năng xác định miền khảo sát, chọn mức cơ sở $u_{i0}$, bước nhảy $\Delta u_i$, và lập bảng ma trận thực nghiệm $2^k$ hoặc ma trận cấp 2 mở rộng.
  • Kỹ năng xử lý dữ liệu định lượng: Kỹ thuật mã hóa/giải mã số liệu, tính toán ma trận thực nghiệm và ước lượng trực tiếp các hệ số hồi quy không cần qua phép nghịch đảo ma trận phức tạp.
  • Kỹ năng mô hình hóa công nghệ: Chuyển đổi dữ liệu thực nghiệm đo đạc từ phòng thí nghiệm thành mô hình toán học giải tích phục vụ công tác tối ưu hóa chế độ công nghệ trong CN Hóa học, Sinh học và Thực phẩm.

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

Phương pháp tiếp cận sư phạm

Tài liệu bài giảng áp dụng phương pháp tiếp cận suy diễn logic từ trừu tượng toán học đến cụ thể hóa kỹ thuật:

  • Bắt đầu từ bài toán tổng quát mô hình hóa "hộp đen" $\xi \to [X] \to y$ để người học nắm rõ bản chất tương tác giữa đầu vào, nhiễu và đầu ra.
  • Dẫn dắt cơ sở toán học từ chuỗi Taylor và phương pháp bình phương cực tiểu trước khi rút ra các công thức rút gọn dựa trên tính trực giao.
  • Minh họa không gian thực nghiệm bằng sơ đồ hình học 2D/3D (các đỉnh khối lập phương, các điểm sao trên trục tọa độ và điểm tâm).
Không gian bố trí thí nghiệm trực giao cấp 2 (k = 2):
          x2
[Điểm sao 6]      [Tâm 9]   [Điểm sao 5]

Bài tập và bảng tính mẫu

Tài liệu tích hợp các bảng dữ liệu số thực nghiệm đầy đủ (như ví dụ ma trận thực nghiệm $2^3$ với 8 thí nghiệm, đo lặp 2 lần $y_1, y_2$, tính trung bình $\overline{y}$ và hàm mục tiêu $f(x)$). Qua đó, sinh viên có thể tự thực hành từng bước tính tích vô hướng để tìm hệ số hồi quy: $$b_0 = 0.119, \quad b_2 = -0.131, \quad b_5 = -0.256 \dots$$

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

  • Sinh viên cần nắm vững quy trình biến đổi đại số khi trực giao hóa các cột bậc hai $x'{ij} = x{ij}^2 - \overline{x_i^2}$.
  • Tự kiểm tra tính trực giao của ma trận bằng cách nhân vô hướng các cặp cột bất kỳ trước khi áp dụng công thức tính hệ số $b_j$.
  • Học phần thường được đánh giá qua bài tập lớn xử lý số liệu thực tế tại phòng thí nghiệm chuyên ngành và bài thi tự luận/tính toán cuối kỳ.

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

  1. Hệ thống hóa thuật toán đại số trực giao: Tài liệu chỉ rõ bản chất vì sao ma trận cấp 1 trực giao nhưng ma trận cấp 2 đầy đủ lại mất tính trực giao ở cột bậc hai ($\sum x_{0j} x_{uj}^2 \neq 0$), từ đó đưa ra phương pháp trừ số hạng trung bình $\overline{x_i^2}$ để bảo toàn tính trực giao toàn phần.
  2. Cung cấp bảng tra cánh tay đòn $\alpha$ chuẩn hóa: Xây dựng sẵn bảng tra số học cho giá trị $\alpha$ tương ứng với các trường hợp số biến $k = 2, 3, 4$ và số thí nghiệm tâm $n_0 = 1, 2, 3, 4$, giảm thiểu thao tác giải phương trình phi tuyến bậc 4 phức tạp cho người học: $$2^k(1 - \overline{x^2})^2 - 4(\alpha^2 - \overline{x^2}) + (2^k + n_0 - 4)\overline{x^2}^2 = 0$$
  3. Liên hệ mật thiết với ngành CN Hóa học – Sinh học – Thực phẩm: Các ví dụ và biến số ($Z_1, Z_2, u_1, u_2$) được xây dựng bám sát các thông số công nghệ như nhiệt độ, nồng độ, thời gian phản ứng, áp suất hay tỷ lệ phối liệu.

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

  • Sinh viên đại học: Sinh viên năm thứ 3 và năm thứ 4 các ngành Kỹ thuật Hóa học, Kỹ thuật Sinh học, Công nghệ Thực phẩm, Kỹ thuật Môi trường chuẩn bị làm thí nghiệm chuyên ngành, đồ án môn học và khóa luận tốt nghiệp.
  • Học viên cao học và nghiên cứu sinh: Sử dụng làm tài liệu tra cứu phương pháp luận quy hoạch thực nghiệm khi thực hiện luận văn thạc sĩ và luận án tiến sĩ kỹ thuật.
  • Giảng viên và cán bộ nghiên cứu: Làm đề cương bài giảng học phần CH3309 hoặc tài liệu tham khảo cho các nghiên cứu R&D tại các phòng thí nghiệm phân tích, viện nghiên cứu công nghệ.
  • Điều kiện tiên quyết: Người học cần hoàn thành các học phần Toán cao cấp (Giải tích, Đại số tuyến tính) và Xác suất Thống kê ứng dụng.

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

1. Giáo trình này phù hợp với đối tượng nào?

Tài liệu được thiết kế chuyên biệt cho sinh viên đại học, học viên sau đại học và kỹ sư khối ngành Kỹ thuật Hóa học, Sinh học, Thực phẩm có nhu cầu tối ưu hóa quá trình công nghệ thông qua thực nghiệm khoa học.

2. Cần chuẩn bị kiến thức nền tảng nào trước khi học?

Người học cần có kiến thức vững về đại số tuyến tính (phép tính ma trận, tích vô hướng), giải tích (chuỗi Taylor, đạo hàm riêng) và xác suất thống kê cơ bản (ước lượng kỳ vọng, phương sai, phân phối xác suất).

3. Quy hoạch trực giao cấp 1 khác quy hoạch trực giao cấp 2 ở điểm nào?

  • Cấp 1: Mô tả mối quan hệ tuyến tính giữa biến vào và hàm đáp ứng (kèm tương tác bậc 1), chỉ sử dụng $2^k$ thí nghiệm ở 2 mức $(-1, +1)$.
  • Cấp 2: Mô tả bề mặt đáp ứng phi tuyến có cực trị (chứa các số hạng $x_i^2$), đòi hỏi khảo sát ở ít nhất 3 đến 5 mức, bổ sung thêm các điểm sao $\pm\alpha$ và điểm lặp ở tâm $n_0$.

4. Vì sao cần biến đổi từ biến thực $u_i$ sang biến mã hóa $x_i$?

Việc mã hóa đưa các biến có thứ nguyên và khoảng biến thiên khác nhau về cùng một miền không thứ nguyên đối xứng $[-1, +1]$, giúp chuẩn hóa ma trận, triệt tiêu sự phụ thuộc đơn vị đo và tạo điều kiện trực giao để tính toán hệ số dễ dàng.

5. Làm thế nào để tự học giáo trình này hiệu quả?

Người học nên bám sát trình tự bài giảng: Đọc hiểu lý thuyết mô hình hộp đen $\to$ Thực hành lập ma trận mã hóa $\to$ Tự giải tay các bảng số liệu mẫu trong slide $\to$ Đối chiếu kết quả phương trình hồi quy trước khi ứng dụng vào phần mềm chuyên dụng.

Kết luận

Tài liệu bài giảng Thiết kế thí nghiệm và Xử lý số liệu (CH3309) của PGS. Trịnh Văn Dũng cung cấp hệ thống lý thuyết chuẩn xác và phương pháp luận thực nghiệm chặt chẽ cho khối ngành Kỹ thuật Hóa học, Sinh học và Thực phẩm. Với tiến trình từ mô hình tuyến tính cấp 1 đến mô hình bề mặt đáp ứng cấp 2, tài liệu giúp người học nắm vững bản chất toán học của quy hoạch trực giao và nâng cao năng lực phân tích số liệu thực nghiệm khoa học.