Bối cảnh và vấn đề nghiên cứu

Khóa luận tốt nghiệp chuyên ngành Toán Ứng Dụng với đề tài "Một số ứng dụng của đại số tuyến tính vào các lĩnh vực khác" được thực hiện bởi sinh viên Trần Đình Chiến dưới sự hướng dẫn của cô Phạm Bích Như tại Bộ môn Toán học, Trường Đại học Cần Thơ (tháng 12 năm 2022).

Trong bối cảnh toàn cầu hóa kinh tế và sự phát triển mạnh mẽ của cuộc Cách mạng công nghiệp lần thứ tư (CMCN 4.0), toán học giữ vai trò nền tảng cho nhiều lĩnh vực khoa học công nghệ mũi nhọn như công nghệ thông tin, viễn thông, vật liệu mới và công nghệ sinh học. Khi giải quyết các bài toán thực tiễn phức tạp, các phương pháp tính toán thủ công truyền thống thường bộc lộ hạn chế về mặt thời gian và độ chính xác. Đại số tuyến tính—đặc biệt là lý thuyết hệ phương trình tuyến tính và đại số ma trận—đã trở thành công cụ toán học thiết yếu để mô hình hóa và tối ưu hóa quy trình tính toán trong nhiều ngành nghề.

Khóa luận tập trung thực hiện các nhiệm vụ nghiên cứu sau:

  1. Hệ thống hóa các cơ sở lý thuyết nền tảng về hệ phương trình đại số tuyến tính, định lý về nghiệm, thuật toán khử Gauss và các phép toán trên ma trận.
  2. Khảo sát phương pháp ứng dụng hệ phương trình tuyến tính trong hình học giải tích (xác định đường cong, mặt cong) và giải hệ phương trình vi phân tuyến tính (mô hình động học quần thể).
  3. Xây dựng và phân tích các mô hình toán học ứng dụng ma trận vào 7 lĩnh vực thực tế: Di truyền học, Mật mã học, Hóa học, Kinh tế học Leontief, Chuỗi Markov, Mạch điện và Phân luồng giao thông.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Hệ phương trình đại số tuyến tính, lý thuyết đại số ma trận, định thức, hạng ma trận, ma trận nghịch đảo, bài toán giá trị riêng - vectơ riêng và kỹ thuật chéo hóa ma trận.
  • Phạm vi nội dung: Các bài toán ứng dụng liên ngành trong sinh học quần thể, hệ mật mã cổ điển/Hill, cân bằng phản ứng hóa học, mô hình kinh tế vĩ mô cân đối liên ngành (Input-Output).
  • Phạm vi thời gian và địa bàn: Khóa luận được hoàn thành và nghiệm thu vào tháng 12 năm 2022 tại Trường Đại học Cần Thơ.

Cơ sở lý thuyết và phương pháp

Khóa luận sử dụng khung lý thuyết giải tích và đại số tuyến tính chuẩn tắc kết hợp với phương pháp mô hình hóa toán học liên ngành:

  • Khung lý thuyết và mô hình toán học:

    • Dạng phương trình ma trận tổng quát: $Ax = b$, trong đó $A \in M_{m \times n}(\mathbb{R})$ là ma trận hệ số, $x \in \mathbb{R}^n$ là vectơ ẩn, $b \in \mathbb{R}^m$ là vectơ hệ số tự do.
    • Quy tắc Cramer: Nghiệm duy nhất $x_i = \frac{\Delta_i}{\Delta}$ áp dụng cho hệ có ma trận vuông $A$ không suy biến ($\det(A) \neq 0$).
    • Định lý Kronecker – Capelli: Điều kiện cần và đủ để hệ phương trình có nghiệm là $\text{rank}(A) = \text{rank}(B)$, với $B = (A|b)$ là ma trận mở rộng. Hệ có nghiệm duy nhất khi $\text{rank}(A) = \text{rank}(B) = n$ và vô số nghiệm khi $\text{rank}(A) = \text{rank}(B) < n$.
    • Kỹ thuật chéo hóa ma trận: Biểu diễn ma trận $A$ dưới dạng $A = PDP^{-1}$ (với $D$ là ma trận đường chéo chứa các giá trị riêng $\lambda_i$, $P$ là ma trận chuyển cơ sở tạo bởi các vectơ riêng độc lập tuyến tính) nhằm tối ưu phép tính lũy thừa $A^n = PD^n P^{-1}$.
    • Mô hình kinh tế đầu vào - đầu ra mở của Wassily Leontief (Mô hình Input - Output, Nobel Kinh tế 1973) thông qua hệ phương trình $(I - C)p = d$.
  • Phương pháp nghiên cứu:

    • Phương pháp phân tích - tổng hợp lý thuyết: Thu thập, hệ thống hóa các định lý, hệ quả và thuật toán đại số ma trận từ tài liệu chuẩn.
    • Phương pháp mô hình hóa toán học: Chuyển đổi các bài toán từ ngôn ngữ thực tế (sinh học, hóa học, mật mã, kinh tế) thành các bài toán đại số ma trận tương đương.
    • Phương pháp biến đổi sơ cấp dòng (Gauss - Jordan): Khử ẩn tuần tự đưa ma trận về dạng bậc thang hoặc bậc thang rút gọn để giải hệ phương trình, tính hạng, tìm ma trận nghịch đảo và giải mã thông tin.
    • Phương pháp phân tích phổ: Sử dụng giá trị riêng và vectơ riêng để đánh giá trạng thái dừng dài hạn trong các mô hình tiến hóa di truyền và cân bằng kinh tế.
  • Nguồn dữ liệu:

    • Dữ liệu giả định và dữ liệu thực nghiệm tiêu chuẩn trong các bài toán hóa học (phản ứng cháy, phản ứng phức tạp), các bảng xác suất lai di truyền Mendel, bảng hoán vị mật mã ký tự chữ số và bảng thông số hệ số kỹ thuật của nền kinh tế giả định.

Nội dung chính theo từng chương

Chương 1: Kiến thức chuẩn bị

Chương 1 thiết lập hệ thống khái niệm và công cụ tính toán cơ bản của đại số tuyến tính:

  • Hệ phương trình tuyến tính: Phân loại hệ thuần nhất ($b_i = 0, \forall i$) và không thuần nhất; hệ tương thích (có nghiệm) và không tương thích (vô nghiệm); hệ xác định (nghiệm duy nhất) và bất định (vô số nghiệm). Trình bày chứng minh quy tắc Cramer và Định lý Kronecker – Capelli.
  • Phương pháp khử Gauss: Mô tả thuật toán biến đổi sơ cấp theo hàng đưa hệ phương trình về dạng tam giác để truy hồi ngược tìm nghiệm. Văn bản minh họa bằng ví dụ giải hệ 3 phương trình 3 ẩn: $$\begin{cases} 4x_1 + 6x_2 + 10x_3 = 4 \ 6x_1 - 5x_2 + 8x_3 = 20 \ -8x_1 + 6x_2 + 4x_3 = 4 \end{cases} \implies (x_1, x_2, x_3) = (-1, -2, 2)$$
  • Đại số ma trận: Định nghĩa các dạng ma trận đặc biệt (ma trận dòng, ma trận không, ma trận vuông, ma trận tam giác trên/dưới, ma trận chéo, ma trận đơn vị $I_n$, ma trận chuyển vị $A^t$, ma trận đối xứng). Định nghĩa các phép toán: cộng ma trận cùng cấp, nhân vô hướng, nhân ma trận.
  • Biến đổi sơ cấp và ma trận nghịch đảo: Định nghĩa ma trận bậc thang, khái niệm hạng ma trận $\text{rank}(A)$ và thuật toán tìm ma trận nghịch đảo $A^{-1}$ bằng cách biến đổi sơ cấp trên ma trận ghép $(A | I_n) \to (I_n | A^{-1})$.
Khái niệm / Công cụ Điều kiện áp dụng Kết quả / Ý nghĩa toán học
Quy tắc Cramer Hệ có số phương trình bằng số ẩn ($m=n$) và $\det(A) \neq 0$ Nghiệm duy nhất $x_i = \Delta_i / \Delta$
Định lý Kronecker – Capelli Áp dụng cho mọi hệ $m \times n$ Hệ có nghiệm $\iff \text{rank}(A) = \text{rank}(A|b)$
Khử Gauss - Jordan Mọi ma trận đại số Đưa về ma trận bậc thang rút gọn để giải hệ và tìm $A^{-1}$
Chéo hóa ma trận Ma trận vuông cấp $n$ có $n$ vectơ riêng độc lập tuyến tính $A = PDP^{-1}$, tính nhanh $A^k = PD^k P^{-1}$

Chương 2: Ứng dụng hệ phương trình tuyến tính trong các bài toán thực tế

Chương này trình bày hai hướng ứng dụng hình học và giải tích của hệ phương trình tuyến tính:

  • Xây dựng đường cong và mặt cong: Xác định phương trình đường thẳng $ax + by + c = 0$, đường tròn $ax^2 + by^2 + cx + dy + e = 0$, hoặc các đường conic qua các điểm chỉ định bằng cách thiết lập hệ phương trình tuyến tính thuần nhất $AX = 0$. Hệ có nghiệm không tầm thường khi và chỉ khi $\det(A) = 0$.
    • Ví dụ: Tìm đường tròn đi qua ba điểm $A(1,1)$, $B(-2,1)$, $C(3,2)$ thông qua giải hệ định thức cấp 4 với các ẩn hệ số $a, b, c, d$.
  • Giải hệ phương trình vi phân: Mô hình hóa các bài toán sinh thái học về mối quan hệ giữa quần thể con mồi và động vật ăn thịt (predator - prey). Kỹ thuật chuyển đổi phương trình vi phân tuyến tính bậc cao (cấp 2, cấp 4) về hệ phương trình vi phân tuyến tính cấp 1 dạng ma trận: $\vec{x}'(t) = A\vec{x}(t) + \vec{g}(t)$.

Chương 3: Ứng dụng hệ số ma trận trong các bài toán thực tế

Chương 3 triển khai chi tiết các ứng dụng của ma trận trong các ngành khoa học chuyên biệt:

  • Ứng dụng trong Di truyền học:

    • Mô hình hóa sự phân bố kiểu gen ($AA, Aa, aa$) qua $n$ thế hệ bằng phương trình trạng thái $X_n = A^n X_0$, với $A$ là ma trận chuyển tiếp xác suất.
    • Áp dụng kỹ thuật chéo hóa ma trận $A = PDP^{-1}$ để tính phân bố giới hạn khi $n \to \infty$. Trong bài toán lai cá thể với kiểu gen $AA$, ma trận chuyển tiếp có các giá trị riêng $\lambda_1 = 1, \lambda_2 = 1/2, \lambda_3 = 0$; phân phối kiểu gen dài hạn hội tụ về vector $[1, 0, 0]^T$ (100% cá thể mang kiểu gen $AA$).
    • Phân tích di truyền liên kết giới tính trên nhiễm sắc thể X.
    • Bài toán chuyển đổi nguyện vọng giữa 200 sinh viên ngành Toán và 350 sinh viên ngành Sinh với ma trận chuyển tiếp $A = \begin{pmatrix} 0.92 & 0.04 \ 0.08 & 0.96 \end{pmatrix}$, tìm ra số sinh viên toán học sau 15 năm là 186 người và sau 20 năm có 365 người theo ngành Sinh học.
  • Ứng dụng trong Mật mã học:

    • Nghiên cứu mật mã ma trận (mật mã Hill): Biểu diễn bản rõ thành chuỗi số nguyên, gom cụm thành các vectơ kích thước $n \times 1$ và mã hóa bằng ma trận khả nghịch $A$ cấp $n$ ($C = AP$). Quá trình giải mã thực hiện bằng tích $P = A^{-1}C$.
    • Thám mã (phá mã) khi biết một phần bản rõ và bản mã tương ứng: Thiết lập phương trình $Q = AP \implies A^{-1} = PQ^{-1}$, tìm ma trận khóa giải mã $(A^{-1})^T$ bằng cách khử hàng Gauss trên ma trận ghép $[Q^T | P^T] \to [I | (A^{-1})^T]$. Minh họa qua việc giải mã thành công thông điệp ẩn "IAMAGOODSTUDENTT" từ chuỗi mật mã cho trước.
  • Ứng dụng trong Hóa học:

    • Thiết lập phương trình cân bằng phản ứng hóa học dựa trên định luật bảo toàn nguyên tố. Chuyển đổi số lượng nguyên tử từng nguyên tố thành các vectơ thành phần, ghép thành ma trận hệ số và đưa về dạng bậc thang rút gọn (RREF).
    • Nghiệm của hệ là không gian nghiệm chứa biến tự do; chọn giá trị nguyên dương nhỏ nhất cho biến tự do để thu được các hệ số nguyên tối giản của phương trình hóa học.
  • Ứng dụng trong Mô hình kinh tế Leontief:

    • Mô hình Input - Output mở: Xét nền kinh tế gồm $n$ ngành sản xuất, ma trận tiêu thụ $C = (c_{ij}){n \times n}$ (trong đó $c{ij}$ là giá trị sản lượng ngành $i$ cần để ngành $j$ tạo ra một đơn vị giá trị), vectơ sản xuất $p$, vectơ nhu cầu mở bên ngoài $d$.
    • Mối quan hệ cân bằng: $(I - C)p = d \implies p = (I - C)^{-1}d$.
    • Điều kiện sinh lời và tính khả dụng của mô hình: Ma trận tiêu thụ $C$ là hiệu quả khi $(I - C)^{-1}$ tồn tại và không âm. Điều kiện đủ là tổng các phần tử trên mỗi cột (hoặc mỗi hàng) của ma trận $C$ đều nhỏ hơn 1.

Kết quả và đóng góp

  • Kết quả tính toán cụ thể trong văn bản:

    • Giải quyết chính xác các bài toán nghiệm đại số bằng phương pháp ma trận nghịch đảo, quy tắc Cramer và phương pháp khử Gauss.
    • Xác định trạng thái dừng phân bố di truyền sau $n$ thế hệ: chứng minh bằng giải tích phổ rằng khi cho lai liên tục với dòng thuần chủng $AA$, vector trạng thái luôn tiệm cận về $[1, 0, 0]^T$.
    • Xác định nghiệm của bài toán dịch chuyển chuyên ngành với giá trị riêng $\lambda_1 = 22/25, \lambda_2 = 1$, đưa ra dự báo chính xác số lượng 186 sinh viên Toán sau 15 năm và 365 sinh viên Sinh sau 20 năm.
    • Thiết lập quy trình thám mã tìm ma trận nghịch đảo $A^{-1} = \begin{pmatrix} -2/13 & 5/13 \ 3/13 & -1/13 \end{pmatrix}$ để khôi phục bản rõ từ bản mã mà không cần biết trước khóa mã hóa ban đầu.
    • Đưa ra nghiệm tổng quát cho bài toán cân bằng kinh tế Leontief với ma trận tiêu thụ $C$ có các phần tử $k$ ($0 < k < 1/3$), xác định vectơ sản xuất $p = \frac{1000}{1-3k}[1, 1, 1]^T$ tương ứng với nhu cầu tiêu dùng $d = [1000, 1000, 1000]^T$.
  • Đóng góp của khóa luận:

    • Hệ thống hóa quy trình ứng dụng công cụ đại số tuyến tính vào các bài toán liên ngành từ khoa học tự nhiên (Hóa học, Sinh học), kỹ thuật tính toán (Mật mã học, Phương trình vi phân) đến khoa học xã hội (Kinh tế học vĩ mô).
    • Minh họa tính ưu việt của phương pháp ma trận (như thuật toán Gauss-Jordan, kỹ thuật chéo hóa ma trận) so với các phép biến đổi đại số sơ cấp rời rạc khi xử lý bài toán nhiều biến số.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế được ghi nhận:

    • Tác giả tự nhận định kinh nghiệm làm đề tài nghiên cứu còn hạn chế nên luận văn khó tránh khỏi một số thiếu sót về mặt học thuật.
    • Các ví dụ tính toán minh họa trong khóa luận chủ yếu dừng lại ở quy mô nhỏ (ma trận cấp $2 \times 2, 3 \times 3, 4 \times 4$), mang tính minh họa phương pháp lý thuyết chứ chưa xử lý các tập dữ liệu lớn trong thực tế.
    • Khóa luận chưa kết hợp việc lập trình tự động hóa các thuật toán đại số tuyến tính trên các phần mềm tính toán chuyên dụng (như MATLAB, Python NumPy, R).
  • Hướng nghiên cứu tiếp theo:

    • Triển khai chi tiết các mục ứng dụng đã nêu trong đề cương mục lục bao gồm: Mô hình chuỗi Markov trong dự báo xác suất, phân tích mạng lưới mạch điện theo định luật Kirchhoff, và mô hình hóa bài toán phân luồng giao thông đô thị.
    • Mở rộng quy mô nghiên cứu sang các hệ thống dữ liệu lớn bằng cách tích hợp thuật toán số trị và lập trình tính toán ma trận.

Giá trị tham khảo

Khóa luận là tài liệu học liệu phù hợp cho sinh viên, giảng viên và người học thuộc các nhóm chuyên ngành sau:

  • Sinh viên ngành Toán học, Toán Ứng dụng, Sư phạm Toán: Tham khảo các ví dụ thực tiễn sinh động để làm tài liệu học tập, giảng dạy môn Đại số tuyến tính, giúp kết nối lý thuyết trừu tượng với ứng dụng thực tế.
  • Sinh viên ngành Kinh tế, Quản trị kinh doanh: Tham khảo phương pháp thiết lập và giải ma trận đầu vào - đầu ra Leontief (Chương 3, Mục 3.4), hiểu rõ điều kiện khả nghịch và tính sinh lời của các ngành kinh tế.
  • Sinh viên ngành Khoa học máy tính, An toàn thông tin: Tham khảo nguyên lý hoạt động và kỹ thuật thám mã của hệ mật mã Hill dựa trên phép nhân ma trận và khử Gauss-Jordan (Chương 3, Mục 3.2).
  • Sinh viên ngành Công nghệ sinh học, Hóa học: Tham khảo quy trình vector hóa để cân bằng phản ứng hóa học phức tạp và mô hình ma trận chuyển tiếp xác suất kiểu gen qua các thế hệ.

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

1. Theo Định lý Kronecker – Capelli trong khóa luận, điều kiện nào quyết định số lượng nghiệm của một hệ phương trình tuyến tính?
Điều kiện cần và đủ để hệ phương trình $Ax=b$ có nghiệm là hạng của ma trận hệ số bằng hạng của ma trận mở rộng, tức $\text{rank}(A) = \text{rank}(A|b)$. Nếu $\text{rank}(A) = \text{rank}(A|b) = n$ (số ẩn), hệ có nghiệm duy nhất. Nếu $\text{rank}(A) = \text{rank}(A|b) < n$, hệ có vô số nghiệm phụ thuộc vào $n - \text{rank}(A)$ ẩn tự do.

2. Tại sao trong bài toán di truyền học, tác giả lại sử dụng kỹ thuật chéo hóa ma trận $A = PDP^{-1}$ thay vì nhân ma trận $A^n$ trực tiếp?
Khi tính toán phân bố kiểu gen sau $n$ thế hệ lớn (ví dụ $n = 20$), việc nhân ma trận liên tiếp $A^n$ rất phức tạp và dễ dẫn đến sai số tính toán tích lũy. Kỹ thuật chéo hóa $A = PDP^{-1}$ cho phép tính $A^n = PD^n P^{-1}$ dễ dàng thông qua việc lũy thừa các phần tử trên đường chéo chính của ma trận $D$, giúp giảm khối lượng tính toán và hạn chế sai số.

3. Phương pháp thám mã hệ mật mã ma trận Hill được thực hiện như thế nào khi biết một phần cặp bản rõ - bản mã?
Khi biết trước $n$ vectơ bản rõ $P = [P_1 | P_2 | \dots | P_n]$ và các vectơ bản mã tương ứng $Q = [AP_1 | AP_2 | \dots | AP_n]$, ta có hệ thức $Q = AP \implies A^{-1} = PQ^{-1}$. Để tìm ma trận giải mã $(A^{-1})^T$, tác giả lập ma trận ghép $[Q^T | P^T]$ rồi thực hiện các phép biến đổi sơ cấp dòng của phương pháp Gauss để đưa về dạng $[I | (A^{-1})^T]$, từ đó suy ra ma trận giải mã $A^{-1}$.

4. Đại số tuyến tính hỗ trợ cân bằng phương trình phản ứng hóa học bằng nguyên lý nào?
Phương pháp dựa trên định luật bảo toàn nguyên tố để lập hệ phương trình tuyến tính thuần nhất cho các hệ số phản ứng. Mỗi chất được biểu diễn bằng một vectơ số lượng nguyên tử của từng nguyên tố cấu thành. Bằng cách đưa ma trận hệ số về dạng bậc thang rút gọn (RREF), ta xác định được biến tự do và tìm nghiệm nguyên dương nhỏ nhất thỏa mãn phương trình.

5. Trong mô hình kinh tế Leontief, điều kiện nào đảm bảo ma trận tiêu thụ $C$ là hiệu quả và có thể đáp ứng mọi vectơ nhu cầu bên ngoài?
Ma trận tiêu thụ $C$ được gọi là hiệu quả khi ma trận nghịch đảo $(I - C)^{-1}$ tồn tại và không âm. Điều kiện đủ để đạt được điều này là tổng các phần tử trên mỗi cột (hoặc mỗi hàng) của ma trận $C$ phải nhỏ hơn 1, biểu thị rằng chi phí đầu vào của mỗi ngành để tạo ra một đơn vị sản lượng luôn nhỏ hơn giá trị của đơn vị sản lượng đó (ngành có lãi).


Kết luận

Khóa luận tốt nghiệp của sinh viên Trần Đình Chiến đã hệ thống hóa toàn diện các kiến thức nền tảng của đại số tuyến tính và chứng minh rõ nét tính ứng dụng liên ngành của công cụ ma trận. Bằng việc giải quyết cụ thể các bài toán từ di truyền học, mật mã học, hóa học đến mô hình kinh tế Leontief, công trình làm sáng tỏ vai trò cầu nối của toán học ứng dụng trong thực tiễn. Đề tài là nguồn tài liệu tham khảo học thuật có giá trị cho sinh viên chuyên ngành Toán cũng như các ngành khoa học tự nhiên, kỹ thuật và kinh tế.