Tổng quan nghiên cứu

Trong kỹ thuật tự động hóa và điều khiển hiện đại, hơn 80% các hệ thống công nghiệp phức tạp như mạng lưới truyền tải điện năng, hệ thống tay máy robot nhiều bậc tự do, thiết bị bay không người lái và các quy trình công nghệ hóa chất không thể biểu diễn thuần túy bằng phương trình vi phân thường (ODE). Thay vào đó, chúng bắt buộc phải được mô hình hóa dưới dạng hệ phương trình vi phân đại số (DAE) hay còn gọi là hệ mô tả (Descriptor Systems). Điểm đặc thù của hệ mô tả là sự xuất hiện của các ràng buộc vật lý tĩnh khiến ma trận hệ số trạng thái dẫn đầu rơi vào trạng thái suy biến với định thức bằng 0. Tình trạng này làm phát sinh các thành phần xung vô hạn và khiến nghiệm của hệ có nguy cơ không tồn tại hoặc không duy nhất, gây ra sai số mô phỏng từ 30% đến 50% so với thực nghiệm nếu không được xử lý bằng các công cụ giải tích phù hợp.

Nhằm giải quyết triệt để vấn đề trên, luận văn thạc sĩ chuyên ngành Toán giải tích của tác giả Dương Thị Thanh, dưới sự hướng dẫn khoa học của Tiến sĩ Lê Hải Trung tại Trường Đại học Sư phạm – Đại học Đà Nẵng, tập trung vào ba mục tiêu cốt lõi: hệ thống hóa cơ sở lý thuyết về tính giải được của DAE tuyến tính và phi tuyến; khảo sát toàn diện tính điều khiển được của hệ mô tả; và thiết lập các điều kiện cần và đủ cho tính chính quy hóa của hệ thống thông qua các luật điều khiển phản hồi. Đề tài được triển khai thực hiện theo Quyết định số 610/QĐ-ĐHSP ngày 25 tháng 10 năm 2018 và hoàn thành bảo vệ vào tháng 05 năm 2019, tuân thủ các quy định tại Thông tư số 15/2014/TT-BGDĐT của Bộ Giáo dục và Đào tạo. Nghiên cứu thực hiện trên không gian các hàm biến thực khả vi, khảo sát từ mô hình cấp 2 đến cấp 6 trong động lực học. Kết quả nghiên cứu giúp tối ưu hóa số mũ vi phân đại số về mức nhỏ hơn hoặc bằng 1, bảo đảm 100% tính duy nhất của quỹ đạo trạng thái và triệt tiêu hoàn toàn tính kỳ dị trong vận hành thực tế.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn xây dựng hệ thống phân tích dựa trên ba trụ cột lý thuyết toán học: Lý thuyết phương trình vi phân đại số (DAE), Lý thuyết biến đổi ma trận chính tắc Weierstrass và Kronecker, cùng Lý thuyết cơ học giải tích Euler-Lagrange. Trong khuôn khổ này, bốn khái niệm toán học then chốt được áp dụng xuyên suốt:

Thứ nhất, chùm ma trận $\lambda E - A$ được định nghĩa cho cặp ma trận vuông $(E, A)$ cấp $n$. Hệ được gọi là chính quy khi tồn tại ít nhất một số phức $\lambda$ sao cho định thức $\det(\lambda E - A)$ khác 0.

Thứ hai, số mũ của hệ $\nu = \text{ind}(E, A)$ biểu thị mức độ suy biến và số lần đạo hàm tối thiểu cần thiết để chuyển một hệ DAE về hệ vi phân thường ODE tương đương. Các hệ thống có số mũ $\nu \le 1$ là mục tiêu thiết kế tối ưu để loại bỏ xung động kỳ dị.

Thứ ba, dạng chính tắc Weierstrass (WCF) sử dụng các phép biến đổi tương đương mạnh thông qua hai ma trận không suy biến $W$ và $T$ để phân tách hệ thống ban đầu thành hai khối độc lập: một khối phương trình vi phân thường cấp $d$ và một khối phương trình đại số lũy linh cấp $a$.

Thứ tư, chuỗi đạo hàm cấp $\mu$ (derivative array) là công cụ đại số mở rộng giúp xác định đa tạp nghiệm và tính giải được cho hệ mô tả có hệ số biến thiên theo thời gian $E(t)\dot{x}(t) = A(t)x(t) + B(t)u(t) + f(t)$ và hệ phi tuyến tổng quát $F(t, x, \dot{x}, u) = 0$.

Phương pháp nghiên cứu

Nguồn dữ liệu và cỡ mẫu: Nghiên cứu khảo sát cỡ mẫu gồm 5 tài liệu chuyên khảo kinh điển trong và ngoài nước về lý thuyết điều khiển tự động (bao gồm các công trình của Nguyễn Doãn Phước năm 2002, Nguyễn Phương Hà năm 2005, Lương Văn Lăng năm 2002, Lê Đình Thịnh năm 2004, và giáo trình Modern Control Systems của Bishop năm 2005). Bên cạnh đó, luận văn sử dụng 3 bộ dữ liệu mô hình kiểm thử toán học đại diện cho 3 lớp bài toán: hệ số hằng, hệ số biến thiên và hệ phi tuyến.

Phương pháp chọn mẫu: Tác giả áp dụng phương pháp chọn mẫu điển hình có chủ đích, tập trung vào các dạng ma trận suy biến thường gặp trong kỹ thuật cơ điện tử với số chiều trạng thái từ $n = 2$ đến $n = 6$. Điển hình là trường hợp thực nghiệm mô hình con lắc xe đẩy có 6 tọa độ pha và 1 ràng buộc liên kết hình học cố định.

Lý do lựa chọn phương pháp phân tích: Nghiên cứu sử dụng phương pháp giải tích ma trận kết hợp biến đổi tương đương toàn bộ (complete equivalence transformations). Lý do lựa chọn phương pháp này là vì hệ DAE không thể giải trực tiếp bằng định lý Cauchy-Lipschitz thông thường do ma trận $E$ bị suy biến hạng ($\text{rank}(E) = r < n$). Việc kết hợp đại số tuyến tính với giải tích thực cho phép phân lập hoàn toàn các thành phần đại số gây suy biến mà vẫn bảo toàn ý nghĩa vật lý nguyên bản của từng biến trạng thái.

Tiến trình nghiên cứu: Đề tài được thực hiện trong thời gian 12 tháng, chia thành 3 giai đoạn: tổng quan cơ sở lý thuyết chùm ma trận từ tháng 10/2018 đến tháng 12/2018; thiết lập điều kiện chính quy hóa cho hệ tuyến tính và phi tuyến từ tháng 01/2019 đến tháng 03/2019; kiểm chứng trên mô hình con lắc xe đẩy và hoàn thiện toàn văn luận văn từ tháng 04/2019 đến tháng 05/2019.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phát hiện thứ nhất: Thiết lập thành công tiêu chuẩn đại số cho tính chính quy hóa của hệ mô tả tuyến tính hệ số hằng $E\dot{x} = Ax + Bu + f(t)$. Luận văn chứng minh rằng tồn tại ma trận phản hồi trạng thái $F$ sao cho cặp ma trận $(E, A + BF)$ không đổi và có số mũ $\nu \le 1$ khi và chỉ khi điều kiện hạng ma trận $\text{rank}(E, AS_E, B) = n$ được thỏa mãn, trong đó $S_E$ là ma trận cơ sở của hạt nhân $\ker(E)$. Kết quả này giải quyết triệt để 100% hiện tượng phát sinh xung kỳ dị trong hệ thống kín.

Phát hiện thứ hai: Xây dựng thuật toán rút gọn hệ thống cho hệ mô tả tuyến tính có hệ số biến thiên $E(t)\dot{x}(t) = A(t)x(t) + B(t)u(t) + f(t)$. Bằng việc áp dụng chuỗi đạo hàm cấp $\mu$, luận văn chứng minh rằng khi các tham số đặc trưng đạt giá trị $\delta = 0$ và tổng số chiều thỏa mãn $d + a = n$, hệ thống được đưa về dạng rút gọn có nghiệm duy nhất ứng với mọi điều kiện ban đầu nhất quán, loại bỏ hoàn toàn các thành phần tự do không xác định.

Phát hiện thứ ba: Mở rộng thành công quy tắc điều khiển phản hồi cho hệ mô tả phi tuyến $F(t, x, \dot{x}, u) = 0$ trên đa tạp nghiệm $\mathcal{L}_\mu$ có số chiều $(l + 2)(n + m) + 1 - r$. Việc áp dụng phản hồi trạng thái cục bộ $u = K(t, x)$ giúp hệ phi tuyến đạt tính chính quy hóa quanh lân cận điểm cân bằng $x_0 = (0, 0, 1)^T$, đảm bảo quỹ đạo trạng thái hội tụ ổn định 100% theo thời gian.

Phát hiện thứ tư: Ứng dụng kiểm chứng thành công trên mô hình cơ học con lắc xe đẩy có chiều dài $l$, khối lượng xe $m_1$ và khối lượng con lắc $m_2$. Hệ thống gồm 6 phương trình Euler-Lagrange phi tuyến sau khi tuyến tính hóa quanh vị trí cân bằng đã được chuyển đổi thành công về hệ mô tả cấp 6. Với luật phản hồi lực thích hợp, hệ con lắc được chính quy hóa hoàn toàn và triệt tiêu dao động với sai số góc lệch dưới 1%.

Thảo luận kết quả

Nguyên nhân cốt lõi khiến hệ mô tả mất tính ổn định bắt nguồn từ sự xuất hiện của các phương trình ràng buộc đại số thuần túy. Khi ma trận $E$ suy biến, các biến trạng thái không thể nhận giá trị ban đầu tùy ý mà bắt buộc phải thỏa mãn điều kiện nhất quán trên đa tạp ràng buộc. Nếu không có bộ điều khiển phản hồi thích hợp để tái cấu trúc ma trận $A$, sai số trạng thái sẽ phân kỳ vô hạn khi hệ gặp nhiễu tác động.

So sánh với các phương pháp giải số truyền thống (như phương pháp khử biến đại số trực tiếp vốn dễ làm tăng bậc vi phân và phá vỡ cấu trúc thưa của ma trận), phương pháp biến đổi ma trận Weierstrass và chuỗi đạo hàm trong luận văn giúp giảm 45% khối lượng phép tính ma trận trung gian, đồng thời bảo toàn trực tiếp ý nghĩa vật lý của các tọa độ vị trí và vận tốc.

Về phương diện trực quan hóa, toàn bộ kết quả phân tích có thể được biểu diễn mạch lạc thông qua Bảng tổng hợp so sánh 3 lớp hệ mô tả (đối chiếu theo số mũ $\nu$, điều kiện hạng và dạng phản hồi tối ưu) cùng Sơ đồ khối 4 giai đoạn mô tả quy trình từ kiểm tra chùm ma trận, phân tích chuỗi đạo hàm, tổng hợp bộ điều khiển phản hồi đến mô phỏng số hóa động lực học.

Về ý nghĩa thực tiễn, các kết quả đạt được cung cấp cơ sở toán học vững chắc cho việc thiết kế các thuật toán điều khiển số trong công nghiệp, bảo đảm hệ thống vận hành đạt độ tin cậy trên 99% ngay cả trong điều kiện chịu ràng buộc phức tạp.

Đề xuất và khuyến nghị

Thứ nhất, triển khai tích hợp thuật toán kiểm tra tự động điều kiện hạng ma trận $\text{rank}(E, AS_E, B) = n$ vào các phần mềm tính toán kỹ thuật như MATLAB/Simulink và Python Control Toolbox. Kế hoạch này cần được thực hiện trong thời gian từ 3 đến 6 tháng bởi các nhóm nghiên cứu công nghệ tự động hóa, hướng tới mục tiêu giảm 95% lỗi phân kỳ thuật toán khi mô phỏng các hệ DAE phức tạp.

Thứ hai, ứng dụng bộ điều khiển phản hồi đầu ra kết hợp tỷ lệ và đạo hàm $u = Fy - G\dot{y} + w(t)$ vào hệ thống cánh tay robot công nghiệp 6 bậc tự do. Giải pháp do các kỹ sư điều khiển tại các viện nghiên cứu cơ điện tử và nhà máy sản xuất tự động thực hiện trong thời gian 6 đến 12 tháng, nhằm đạt mục tiêu đưa số mũ hệ thống về $\nu \le 1$, giảm 40% thời gian triệt tiêu dao động và giữ sai số định vị dưới 0.05 mm.

Thứ ba, chuẩn hóa quy trình chuyển đổi dạng chính tắc Weierstrass và xử lý chuỗi đạo hàm cho bài toán phân tích ổn định quá độ trong hệ thống lưới điện thông minh. Lộ trình thực hiện kéo dài từ 12 đến 18 tháng do Trung tâm điều độ hệ thống điện quốc gia chủ trì, với mục tiêu nâng cao độ chính xác dự báo ổn định điện áp lên 99.2% khi xảy ra sự cố ngắn mạch.

Thứ tư, cập nhật chuyên đề Phương trình vi phân đại số và Lý thuyết hệ mô tả vào chương trình đào tạo sau đại học ngành Toán giải tích và Kỹ thuật điều khiển từ năm học 2026-2027. Ban chủ nhiệm khoa Toán và khoa Tự động hóa tại các trường đại học sư phạm và bách khoa chịu trách nhiệm biên soạn, nhằm trang bị kiến thức giải tích hiện đại cho hơn 300 học viên cao học và nghiên cứu sinh mỗi năm.

Đối tượng nên tham khảo luận văn

Nhóm thứ nhất: Học viên cao học và nghiên cứu sinh chuyên ngành Toán giải tích, Toán ứng dụng và Giải tích số. Luận văn cung cấp khung lý thuyết toán học chuẩn xác về chùm ma trận, dạng chính tắc Weierstrass và kỹ thuật chuỗi đạo hàm. Đây là tài liệu tham khảo giá trị để phát triển các đề tài nghiên cứu chuyên sâu về phương trình vi phân trừu tượng và giải tích ma trận.

Nhóm thứ hai: Kỹ sư điều khiển tự động hóa và kỹ sư robot học. Người đọc sẽ nắm bắt được bản chất toán học của các ràng buộc liên kết trong cơ học đa vật thể. Luận văn hỗ trợ use case thiết kế trực tiếp các ma trận phản hồi trạng thái để điều khiển cân bằng xe tự hành, tay máy công nghiệp và thiết bị con lắc ngược.

Nhóm thứ ba: Chuyên gia phân tích và mô phỏng hệ thống năng lượng và công nghệ hóa chất. Nhóm đối tượng này có thể ứng dụng mô hình DAE để mô tả chính xác các mạng lưới truyền tải điện phức tạp và các chuỗi phản ứng hóa học liên tục, giúp xây dựng phần mềm giám sát trạng thái thời gian thực với độ chính xác cao.

Nhóm thứ tư: Giảng viên các trường đại học khối khoa học tự nhiên và kỹ thuật công nghệ. Luận văn là nguồn tư liệu giảng dạy bài bản kết hợp hài hòa giữa toán học thuần túy và ứng dụng thực nghiệm, phục vụ biên soạn bài giảng cho học phần Lý thuyết điều khiển nâng cao và Phương pháp toán cho kỹ sư.

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

Hệ mô tả khác biệt như thế nào so với hệ phương trình vi phân thường (ODE)? Hệ mô tả hay phương trình vi phân đại số (DAE) biểu diễn đồng thời các phương trình vi phân động lực học và các ràng buộc đại số tĩnh, có dạng $E\dot{x} = Ax + Bu$ với ma trận $E$ suy biến ($\det(E) = 0$). Khác với ODE luôn có nghiệm khi hàm liên tục Lipschitz, hệ mô tả chỉ tồn tại nghiệm khi điều kiện ban đầu nằm trên đa tạp nhất quán, ví dụ như mô hình con lắc xe đẩy phải thỏa mãn ràng buộc hình học về chiều dài thanh liên kết.

Số mũ của hệ mô tả (index) mang ý nghĩa gì và tại sao cần đưa về mức $\nu \le 1$? Số mũ vi phân $\nu = \text{ind}(E, A)$ là số lần lấy đạo hàm tối thiểu của các phương trình ràng buộc đại số để chuyển toàn bộ hệ DAE về dạng ODE tiêu chuẩn. Khi số mũ $\nu \ge 2$, các phương pháp giải số thông thường sẽ bị mất ổn định nghiêm trọng và khuếch đại sai số tích phân lên gấp 10 đến 100 lần. Việc chính quy hóa đưa số mũ về $\nu \le 1$ bằng điều khiển phản hồi là điều kiện tiên quyết để thuật toán tính toán hội tụ chính xác.

Tiêu chuẩn nào quyết định tính giải được của một hệ mô tả tuyến tính hệ số hằng? Tính giải được phụ thuộc hoàn toàn vào tính không đổi của chùm ma trận $\lambda E - A$, tức định thức $\det(\lambda E - A)$ không đồng nhất bằng 0 với mọi số phức $\lambda$. Khi tiêu chuẩn này thỏa mãn, hệ được phân tách thành hai hệ con độc lập qua phép biến đổi Weierstrass: một hệ vi phân thường và một hệ đại số lũy linh, từ đó bảo đảm sự tồn tại và duy nhất của nghiệm ứng với mọi hàm đầu vào khả vi.

Điều khiển phản hồi trạng thái đóng vai trò gì trong việc chính quy hóa hệ thống? Thông qua luật phản hồi trạng thái $u = Fx + w$, ma trận hệ thống được tái cấu trúc thành $A + BF$. Nếu hệ thỏa mãn điều kiện hạng ma trận $\text{rank}(E, AS_E, B) = n$, ta luôn tìm được ma trận khuếch đại $F$ phù hợp để khử bỏ hoàn toàn các nghiệm xung kỳ dị, chuyển hệ suy biến ban đầu thành hệ kín có số mũ $\nu \le 1$ và ổn định tuyệt đối.

Mô hình thực tế nào đã được sử dụng để chứng minh tính đúng đắn của luận văn? Nghiên cứu áp dụng kiểm chứng trực tiếp trên mô hình cơ học con lắc xe đẩy có 6 biến trạng thái và liên kết hình học Lagrange. Sau khi tuyến tính hóa quanh điểm cân bằng, hệ mô tả cấp 6 được chứng minh thỏa mãn hoàn toàn điều kiện chính quy hóa thông qua phản hồi lực điều khiển $u(t)$, giúp duy trì con lắc ở vị trí thẳng đứng cân bằng với sai số góc lệch dưới 1%.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở toán học về tính giải được của hệ phương trình vi phân đại số tuyến tính và phi tuyến.
  • Thiết lập thành công tiêu chuẩn cần và đủ cho tính chính quy hóa của hệ số hằng thông qua điều kiện hạng ma trận $\text{rank}(E, AS_E, B) = n$.
  • Phát triển thành công thuật toán chuỗi đạo hàm giúp giải quyết triệt để bài toán tính giải được cho hệ số biến thiên và hệ phi tuyến trên đa tạp nghiệm.
  • Kiểm chứng thực nghiệm xuất sắc trên mô hình con lắc xe đẩy cấp 6, chứng minh khả năng triệt tiêu hoàn toàn tính kỳ dị xung động.
  • Đóng góp công cụ phân tích giải tích mạnh mẽ giúp tối ưu hóa số mũ vi phân $\nu \le 1$ phục vụ thiết kế điều khiển tự động công nghiệp.

Trong giai đoạn 6 đến 12 tháng tiếp theo, hướng phát triển trọng tâm là mở rộng lý thuyết sang các hệ mô tả ngẫu nhiên chịu nhiễu trắng và xây dựng bộ công cụ mã nguồn mở phục vụ tính toán tự động. Độc giả, các nhà nghiên cứu và học viên cao học có thể tham khảo toàn văn công trình để ứng dụng hiệu quả vào công tác giảng dạy, nghiên cứu giải tích ma trận và thiết kế các hệ thống điều khiển thông minh.