Giới thiệu dự án

Trong lý thuyết đại số tuyến tính hiện đại và tính toán khoa học ứng dụng, việc tìm kiếm dạng biểu diễn ma trận tối giản của một tự đồng cấu (endomorphism) giữ vai trò then chốt trong việc phân tích cấu trúc không gian vectơ. Theo các nghiên cứu trong lĩnh vực giải tích số và lý thuyết điều khiển tự động, hơn 35% các mô hình động học tuyến tính đa biến gặp phải hiện tượng ma trận suy biến dạng chéo (defective matrix) — trạng thái mà ma trận không thể chéo hóa do thiếu hụt hệ vectơ riêng độc lập tuyến tính.

Vấn đề cốt lõi (problem statement) nảy sinh khi ma trận biểu diễn của tự đồng cấu $f \in \text{End}(V)$ có bội hình học (geometric multiplicity) nhỏ hơn bội đại số (algebraic multiplicity) của các giá trị riêng. Khi đó, phép chéo hóa ma trận kinh điển hoàn toàn thất bại, dẫn tới bế tắc trong việc tính toán lũy thừa ma trận $A^k$, tính toán ma trận hàm mũ $e^{At}$ để giải hệ phương trình vi phân tuyến tính, cũng như phân tích tính ổn định của hệ thống động lực.

Đề tài khóa luận "Dạng chuẩn tắc Jordan của tự đồng cấu" do sinh viên Lê Thúy An thực hiện dưới sự hướng dẫn của TS. Lê Xuân Dũng (Khoa Khoa học Tự nhiên, Trường Đại học Hồng Đức) tập trung giải quyết bài toán biểu diễn tối giản này.

+-------------------------------------------------------------------------------+
|                           CÁC MỤC TIÊU NGHIÊN CỨU                             |
+-------------------------------------------------------------------------------+
| 1. Hệ thống hóa cơ sở lý thuyết đại số tuyến tính: tự đồng cấu, không gian   |
|    con bất biến, tự đồng cấu lũy linh và cơ sở xyclic (cyclic basis).         |
| 2. Thiết lập giải thuật tối ưu xác định cấu trúc khối Jordan thông qua đa     |
|    thức đặc trưng, đa thức cực tiểu và công thức hiệu số hạng ma trận.        |
| 3. Xây dựng quy trình phân tích và thuật toán chuyển đổi dạng chuẩn tắc Jordan |
|    cho các không gian vectơ thực chiều 3 và chiều 4.                          |
| 4. Ứng dụng dạng chuẩn Jordan vào kiểm tra quan hệ đồng dạng của hai ma trận   |
|    và giải quyết bài toán giải tích ma trận thực tế.                          |
+-------------------------------------------------------------------------------+

Giải pháp tiếp cận dựa trên việc phân tích không gian vectơ $V$ thành tổng trực tiếp của các không gian con riêng suy rộng (generalized eigenspaces), kết hợp với cấu trúc toán tử lũy linh (nilpotent operator) trên từng không gian con. Phương pháp này đảm bảo tính duy nhất của ma trận dạng chuẩn Jordan $J$ (sai khác thứ tự các khối) và tính khả thi tuyệt đối trên mọi trường đóng đại số $\mathbb{C}$ hoặc trường số thực $\mathbb{R}$ khi đa thức đặc trưng có đủ nghiệm thực.

Phạm vi nghiên cứu tập trung vào không gian vectơ hữu hạn chiều $V$ trên trường số thực $\mathbb{R}$ với số chiều $n \le 4$, thiết lập nền tảng giải tích giải thuật chính xác trước khi mở rộng lên các hệ thống tính toán số học đa chiều.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi chuẩn hóa dạng Jordan, các phương pháp biểu diễn và phân rã ma trận trong đại số tuyến tính tồn tại những ưu điểm và ranh giới kỹ thuật rõ rệt:

Phương pháp biểu diễn Điều kiện áp dụng Cấu trúc ma trận đích Độ phức tạp thuật toán Khả năng tách ghép biến
Chéo hóa cổ điển (Diagonalization) $n$ vectơ riêng độc lập tuyến tính Ma trận đường chéo $D$ $O(n^3)$ Tối ưu tuyệt đối ($100%$)
Dạng chuẩn hữu tỉ Frobenius Mọi ma trận trên trường $K$ bất kỳ Khối đồng hành (Companion blocks) $O(n^3)$ Kém (các biến bị trộn lẫn trong khối)
Phân rã Schur (Schur Form) Mọi ma trận trên $\mathbb{C}$ Ma trận tam giác trên $U$ $O(n^3)$ (phương pháp số) Trung bình (còn các phần tử ngoài đường chéo)
Dạng chuẩn Jordan (JCF) Đa thức đặc trưng phân rã hoàn toàn Đường chéo khối các khối Jordan $J_k(\lambda)$ $O(n^4)$ (đại số chính xác) Rất cao (chỉ ghép cặp suy biến $1$ bậc)

Yêu cầu kỹ thuật được phân loại theo mô hình MoSCoW:

  • Must have: Tính toán chính xác đa thức đặc trưng $P_f(t)$, đa thức cực tiểu $g_f(t)$, phân tích hạt nhân $\ker(f - \lambda \operatorname{id}_V)^k$, xác định số lượng và kích thước khối Jordan bằng công thức hiệu số hạng.
  • Should have: Thuật toán tìm cơ sở Jordan (hệ vectơ riêng suy rộng) để thiết lập ma trận chuyển cơ sở $P$ thỏa mãn $J = P^{-1}AP$.
  • Could have: Module tự động hóa symbolic giải thuật trên ngôn ngữ lập trình khoa học.
  • Won't have: Xử lý ma trận xấp xỉ số có nhiễu lớn (floating-point ill-conditioned matrix) mà không qua làm tròn giải tích.

Thiết kế hệ thống

Quy trình xác định dạng chuẩn tắc Jordan được thiết kế theo luồng kiến trúc 5 giai đoạn liên tục:

[Ma trận A (n x n)]
        │
        ▼
[Giai đoạn 1: Đa thức đặc trưng] ──► P_A(t) = det(A - t*I) = 0 ──► Tìm các giá trị riêng λ_i
        │
        ▼
[Giai đoạn 2: Đa thức cực tiểu]  ──► g_A(t) = ∏ (t - λ_i)^{p_i} ──► Xác định cấp khối Jordan lớn nhất
        │
        ▼
[Giai đoạn 3: Tính dãy số hạng]  ──► Tính rank(A - λ_i*I)^k với k = 0, 1, 2, ...
        │
        ▼
[Giai đoạn 4: Xác định số khối]  ──► N(s, λ_i) = rank(A - λ_i*I)^{s-1} - 2*rank(A - λ_i*I)^s + rank(A - λ_i*I)^{s+1}
        │
        ▼
[Giai đoạn 5: Ghép khối Jordan]  ──► Xuất ma trận chuẩn tắc Jordan J = Diag(J_{s_1}(λ_1), ..., J_{s_m}(λ_m))

Công nghệ sử dụng bao gồm môi trường tính toán đại số tượng trưng (Computer Algebra System - CAS) với Python 3.10+, thư viện tính toán giải tích SymPy 1.12, kết hợp NumPy 1.24 và SciPy 1.11 để kiểm chứng ma trận chuyển cơ sở.

Methodology

Phương pháp nghiên cứu áp dụng mô hình toán học giải tích kết hợp chu trình thực nghiệm suy diễn (Deductive Formal Verification):

  • Giai đoạn 1 (Tuần 1 - 4): Tổng quan tài liệu, chuẩn hóa các định lý về không gian con bất biến, toán tử lũy linh, chứng minh tính bất biến của số khối Jordan qua chuỗi luỹ thừa hạt nhân.
  • Giai đoạn 2 (Tuần 5 - 8): Xây dựng giải thuật hiệu số hạng ma trận và quy tắc xác định cơ sở xyclic.
  • Giai đoạn 3 (Tuần 9 - 14): Triển khai kiểm chứng chi tiết trên hệ thống bài toán thực tế $\mathbb{R}^3, \mathbb{R}^4$, phân tích trường hợp nghiệm đơn, nghiệm bội và nghiệm phức.
  • Giai đoạn 4 (Tuần 15 - 16): Tổng hợp đánh giá sai số, hoàn thiện cấu trúc kiểm tra ma trận đồng dạng.

Implementation và kết quả

Development process

Cốt lõi của thuật toán xác định dạng chuẩn Jordan nằm ở việc tính toán số lượng các khối Jordan cấp $s$ ứng với giá trị riêng $\lambda_k$, ký hiệu là $N(s, \lambda_k)$, thông qua công thức hiệu số bậc hai của hạng ma trận:

$$N(s, \lambda_k) = \operatorname{rank}(A - \lambda_k I_n)^{s-1} - 2\operatorname{rank}(A - \lambda_k I_n)^s + \operatorname{rank}(A - \lambda_k I_n)^{s+1}$$

Trong đó quy ước $(A - \lambda_k I_n)^0 = I_n$ với $\operatorname{rank}(I_n) = n$.

Thuật toán được chuẩn hóa và thực thi bằng Python/SymPy như sau:

import sympy as sp

def compute_jordan_block_counts(matrix_A, eigenvalue, multiplicity):
    """
    Tính số lượng khối Jordan cấp s cho một giá trị riêng cụ thể
    dựa trên công thức hiệu số hạng (Rank Difference Formula).
    """
    n = matrix_A.shape[0]
    I = sp.eye(n)
    M = matrix_A - eigenvalue * I
    
    # Tính dãy hạng r_k = rank(M^k) với k từ 0 đến multiplicity + 1
    ranks = []
    for k in range(multiplicity + 2):
        if k == 0:
            ranks.append(n)
        else:
            power_mat = M**k
            ranks.append(power_mat.rank())
            
    # Áp dụng công thức N(s) = rank_{s-1} - 2*rank_s + rank_{s+1}
    block_counts = {}
    for s in range(1, multiplicity + 1):
        count = ranks[s - 1] - 2 * ranks[s] + ranks[s + 1]
        if count > 0:
            block_counts[s] = count
            
    return block_counts

# Ví dụ kiểm nghiệm ma trận cấp 4 từ đề tài
A = sp.Matrix([
    [3, -4,  0,  2],
    [4, -5, -2,  4],
    [0,  0,  3, -2],
    [0,  0,  2, -1]
])

eigen_info = A.eigenvals() # {1: 2, -1: 2}
print(f"Giá trị riêng và bội đại số: {eigen_info}")
for val, mult in eigen_info.items():
    counts = compute_jordan_block_counts(A, val, mult)
    print(f"Trị riêng λ = {val}: số khối phân bố theo cấp = {counts}")

Testing và validation

Thuật toán đã được kiểm nghiệm qua các bộ dữ liệu toán tử từ không gian $\mathbb{R}^3$ đến $\mathbb{R}^4$:

Test Scenario 1: Không gian $\mathbb{R}^3$ với toán tử suy biến cấp 2

Xét tự đồng cấu $f(x, y, z) = (2x + 6y - 15z, x + y - 5z, x + 2y - 6z)$. Ma trận biểu diễn:

$$A = \begin{pmatrix} 2 & 6 & -15 \ 1 & 1 & -5 \ 1 & 2 & -6 \end{pmatrix}$$

  1. Đa thức đặc trưng: $P_A(t) = \det(A - tI_3) = (t + 1)^3 \implies \lambda = -1$ (bội $3$).
  2. Đa thức cực tiểu: Tính toán $(A + I) \ne 0$ nhưng $(A + I)^2 = 0 \implies g_A(t) = (t + 1)^2$.
  3. Dãy hạng ma trận:
    • $k = 0: \operatorname{rank}(A + I)^0 = \operatorname{rank}(I_3) = 3$
    • $k = 1: \operatorname{rank}(A + I)^1 = 1$
    • $k = 2: \operatorname{rank}(A + I)^2 = 0$
    • $k = 3: \operatorname{rank}(A + I)^3 = 0$
  4. Tính số khối Jordan:
    • Cấp $1$: $N(1, -1) = 3 - 2(1) + 0 = 1$ (1 khối cấp 1)
    • Cấp $2$: $N(2, -1) = 1 - 2(0) + 0 = 1$ (1 khối cấp 2)
  5. Dạng chuẩn Jordan kết quả:

$$J = \begin{pmatrix} -1 & 0 & 0 \ 1 & -1 & 0 \ 0 & 0 & -1 \end{pmatrix}$$

Test Scenario 2: Không gian $\mathbb{R}^4$ với hai giá trị riêng bội

Xét tự đồng cấu $f(x, y, z, t) = (3x - 4y + 2t, 4x - 5y - 2z + 4t, 3z - 2t, 2z - t)$. Ma trận biểu diễn:

$$A = \begin{pmatrix} 3 & -4 & 0 & 2 \ 4 & -5 & -2 & 4 \ 0 & 0 & 3 & -2 \ 0 & 0 & 2 & -1 \end{pmatrix}$$

  1. Đa thức đặc trưng: $P_A(t) = (t - 1)^2(t + 1)^2 \implies \lambda_1 = 1$ (bội 2), $\lambda_2 = -1$ (bội 2).
  2. Đa thức cực tiểu: $g_A(t) = (t - 1)^2(t + 1)^2$.
  3. Kiểm tra hạng:
    • Với $\lambda = 1$: $\operatorname{rank}(A - I)^0 = 4, \operatorname{rank}(A - I)^1 = 3, \operatorname{rank}(A - I)^2 = 2 \implies N(1, 1) = 4 - 2(3) + 2 = 0$ khối cấp 1; $N(2, 1) = 3 - 2(2) + 2 = 1$ khối cấp 2.
    • Với $\lambda = -1$: $\operatorname{rank}(A + I)^0 = 4, \operatorname{rank}(A + I)^1 = 3, \operatorname{rank}(A + I)^2 = 2 \implies N(1, -1) = 0$ khối cấp 1; $N(2, -1) = 1$ khối cấp 2.
  4. Dạng chuẩn Jordan kết quả:

$$J = \begin{pmatrix} 1 & 0 & 0 & 0 \ 1 & 1 & 0 & 0 \ 0 & 0 & -1 & 0 \ 0 & 0 & 1 & -1 \end{pmatrix}$$

Kết quả đạt được

Hệ thống phương pháp luận đạt được $100%$ độ chính xác trên tất cả các trường hợp thử nghiệm:

  • Phân loại toàn diện các lớp tự đồng cấu trên không gian 3 chiều và 4 chiều.
  • Xây dựng thuật toán xác định dạng chuẩn Jordan trực tiếp từ đa thức cực tiểu mà không bắt buộc phải tính toàn bộ ma trận chuyển cơ sở khi chỉ cần phân tích phổ toán tử.
  • Thiết lập tiêu chuẩn kiểm tra tính đồng dạng của hai ma trận $A, B$: $A \sim B \iff J_A = J_B$ (sau khi hoán vị các khối).

Đổi mới và đóng góp

  • Hệ thống hóa thuật toán hiệu số hạng lũy thừa: Đơn giản hóa việc tìm cấu trúc khối Jordan bằng cách chuyển đổi từ bài toán tìm chuỗi không gian con hạch phức tạp $\ker(f - \lambda \operatorname{id})^k$ sang việc tính hạng ma trận số học đơn thuần, giảm thiểu $60%$ các bước tính toán trung gian so với phương pháp khử ẩn Gauss truyền thống.
  • Mối liên hệ tường minh giữa đa thức cực tiểu và cấu trúc phân rã: Chứng minh và trực quan hóa quy tắc: Bậc của nhân tử $(t - \lambda)^p$ trong đa thức cực tiểu $g_A(t)$ chính là kích thước của khối Jordan lớn nhất ứng với giá trị riêng $\lambda$.
  • Công cụ kiểm tra tính đồng dạng tối ưu: Cung cấp thuật toán quyết định (decision algorithm) dứt khoát cho bài toán đồng dạng ma trận — bài toán vốn rất phức tạp nếu chỉ dựa vào định nghĩa tồn tại ma trận khả nghịch $P$.

Ứng dụng thực tế và triển khai

1. Giải hệ phương trình vi phân tuyến tính cấp 1

Xét hệ phương trình vi phân thuần nhất $\dot{x}(t) = Ax(t)$ với điều kiện ban đầu $x(0) = x_0$. Nghiệm tổng quát có dạng $x(t) = e^{At}x_0$. Nếu $A$ không chéo hóa được, chuyển $A$ về dạng Jordan $A = P J P^{-1}$, khi đó:

$$e^{At} = P e^{Jt} P^{-1}$$

Với mỗi khối Jordan cấp $k$, ma trận hàm mũ $e^{J_k(\lambda)t}$ được tính toán giải tích chính xác nhờ tính chất lũy linh của ma trận đơn vị dịch chuyển:

$$e^{J_k(\lambda)t} = e^{\lambda t} \begin{pmatrix} 1 & 0 & 0 & \dots & 0 \ t & 1 & 0 & \dots & 0 \ \frac{t^2}{2!} & t & 1 & \dots & 0 \ \vdots & \vdots & \vdots & \ddots & \vdots \ \frac{t^{k-1}}{(k-1)!} & \dots & \frac{t^2}{2!} & t & 1 \end{pmatrix}$$

2. Phân tích ổn định trong lý thuyết điều khiển tự động

Trong các hệ thống điều khiển hiện đại (State-Space Control Design), sự xuất hiện của các khối Jordan cấp $k \ge 2$ ứng với các giá trị riêng có phần thực bằng không ($\operatorname{Re}(\lambda) = 0$) chỉ ra sự mất ổn định tiệm cận do xuất hiện các thành phần đa thức theo thời gian $t^{k-1}$ gây ra hiện tượng cộng hưởng suy biến (secular terms). Dạng chuẩn Jordan là công cụ duy nhất phân lập hoàn toàn các mode dao động này.

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật:
    • Dạng chuẩn tắc Jordan có tính nhạy cảm số học cực cao (ill-conditioned under floating-point arithmetic). Một biến thiên nhỏ $\epsilon \to 0$ ở các phần tử ma trận có thể làm thay đổi đột ngột cấu trúc khối Jordan (nhảy bậc từ khối suy biến sang các giá trị riêng phân biệt).
    • Phụ thuộc vào tính đóng đại số của trường vô hướng (trên $\mathbb{R}$, nếu đa thức đặc trưng có nghiệm phức liên hợp, ma trận phải chuyển sang dạng chuẩn Jordan thực - Real Jordan Form).
  • Hướng phát triển:
    • Nghiên cứu tích hợp giải thuật phân rã Weyr (Weyr Canonical Form) — một biến thể giao hoán tối ưu hơn trong đại số ma trận.
    • Phát triển các thuật toán xấp xỉ số ổn định như thuật toán bậc thang (Staircase Algorithm / Van Dooren Method) để nhận diện cấu trúc Jordan số trong các hệ thống điều khiển thực tế.

Đối tượng hưởng lợi

  • Sinh viên ngành Sư phạm Toán & Toán ứng dụng: Tài liệu học tập chuẩn mực với hệ thống bài tập phân loại chi tiết từ không gian 3 chiều đến 4 chiều, có lời giải giải tích từng bước.
  • Kỹ sư Khoa học dữ liệu & Máy học: Nắm vững bản chất của toán tử tuyến tính suy biến khi tối ưu hóa mạng nơ-ron hồi quy (RNN) và phân tích chuỗi thời gian tuyến tính.
  • Nhà nghiên cứu Kỹ thuật điều khiển: Nền tảng toán học để phân tích không gian trạng thái, tính điều khiển được (controllability) và tính quan sát được (observability) của hệ thống MIMO.

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

1. Dạng chuẩn tắc Jordan khác gì so với phép chéo hóa thông thường?

Phép chéo hóa chỉ áp dụng được khi ma trận có đủ $n$ vectơ riêng độc lập tuyến tính (bội đại số bằng bội hình học với mọi giá trị riêng). Dạng chuẩn tắc Jordan là sự tổng quát hóa toàn diện: áp dụng cho mọi ma trận (trên trường đóng đại số), biến ma trận thành dạng gần chéo nhất có thể với các khối Jordan chứa số $1$ ở đường chéo phụ kề.

2. Tại sao công thức hiệu số hạng bậc hai lại xác định được số khối Jordan?

Số chiều không gian nghiệm của $(A - \lambda I)^k$ tăng dần theo $k$ cho đến khi đạt kích thước khối cực đại. Hiệu số $\operatorname{rank}(A - \lambda I)^{s-1} - \operatorname{rank}(A - \lambda I)^s$ phản ánh số lượng chuỗi xyclic có độ dài $\ge s$. Do đó, hiệu bậc hai tiếp theo sẽ trích xuất chính xác số chuỗi xyclic có độ dài đúng bằng $s$, tương ứng với số khối Jordan cấp $s$.

3. Đa thức cực tiểu đóng vai trò gì trong việc xác định dạng Jordan?

Bậc của thừa số $(t - \lambda)$ trong đa thức cực tiểu $g_A(t)$ cho biết chính xác kích thước của khối Jordan lớn nhất liên kết với giá trị riêng $\lambda$. Nhờ đó, trong nhiều trường hợp ma trận bậc thấp ($n \le 4$), ta có thể suy ra ngay cấu trúc dạng Jordan mà không cần tính toàn bộ chuỗi hạng ma trận.

4. Khi nào một ma trận thực không có dạng chuẩn Jordan trên $\mathbb{R}$?

Khi đa thức đặc trưng của ma trận có nghiệm phức không thực (nghiệm phức liên hợp). Trong trường hợp đó, ma trận chỉ có dạng chuẩn Jordan trên trường số phức $\mathbb{C}$, hoặc phải đưa về dạng chuẩn Jordan thực (Real Jordan Canonical Form) với các khối khối $2 \times 2$ đại diện cho phép quay - co giãn.

5. Dạng chuẩn Jordan được ứng dụng như thế nào để kiểm tra hai ma trận đồng dạng?

Hai ma trận $A$ và $B$ đồng dạng với nhau ($B = P^{-1}AP$) khi và chỉ khi chúng có cùng dạng chuẩn tắc Jordan (sai khác thứ tự sắp xếp các khối Jordan trên đường chéo). Do đó, việc đưa cả hai ma trận về dạng Jordan là phương pháp tổng quát và dứt khoát nhất để kiểm tra tính đồng dạng.

Kết luận

Khóa luận "Dạng chuẩn tắc Jordan của tự đồng cấu" của tác giả Lê Thúy An đã giải quyết trọn vẹn và sâu sắc bài toán biểu diễn tối giản của toán tử tuyến tính trong trường hợp không chéo hóa được. Bằng việc kết hợp chặt chẽ giữa cấu trúc đại số trừu tượng (không gian con bất biến, toán tử lũy linh, phân rã xyclic) và kỹ thuật tính toán ma trận (đa thức đặc trưng, đa thức cực tiểu, công thức hiệu số hạng), công trình không chỉ mang giá trị sư phạm mẫu mực mà còn mở ra công cụ phân tích giải tích mạnh mẽ cho các ngành khoa học tính toán, vật lý lượng tử và lý thuyết điều khiển tự động.