Bao hàm vi phân khác phương trình vi phân thông thường ở điểm nào?

Phương trình vi phân thông thường có đạo hàm bằng một vector đơn trị duy nhất, trong khi bao hàm vi phân có đạo hàm thuộc về một tập hợp vector xác định bởi ánh xạ đa trị, cho phép mô tả các quá trình đa nghiệm và hệ thống không trơn.

Nghiệm của bao hàm vi phân được định nghĩa như thế nào?

Nghiệm của bao hàm vi phân trên một khoảng thời gian được định nghĩa là một hàm vector liên tục tuyệt đối, thỏa mãn điều kiện đạo hàm tức thời thuộc vào tập giá trị của ánh xạ đa trị tại hầu khắp nơi theo độ đo Lebesgue.

Định lý khả sinh Aubin có ý nghĩa gì trong lý thuyết bao hàm vi phân?

Định lý khả sinh do Aubin thiết lập năm 1984 xác định điều kiện cần và đủ về mặt hình học để quỹ đạo trạng thái luôn nằm trong một miền an toàn cho trước, đòi hỏi giao của tập vận tốc và nón tiếp tuyến Bouligand phải khác rỗng tại mọi điểm trên biên.

Bao hàm vi phân được ứng dụng trong những lĩnh vực nào?

Bao hàm vi phân là công cụ then chốt trong lý thuyết điều khiển tối ưu phi trơn, mô hình hóa hệ động lực lai kết hợp chuyển động liên tục và chuyển mạch rời rạc, cũng như phân tích các hệ cơ học có ma sát khô và va chạm.

Bao hàm vi phân (differential inclusion) là gì?

Bao hàm vi phân (differential inclusion) là một dạng tổng quát hóa của phương trình vi phân thông thường, trong đó đạo hàm theo thời gian của biến trạng thái được yêu cầu thuộc về một tập hợp các vector vận tốc khả dĩ xác định bởi một ánh xạ đa trị thay vì một giá trị vector đơn lẻ duy nhất. Khái niệm này đóng vai trò nền tảng trong giải tích biến phân hiện đại, cho phép mô tả chính xác các quá trình động lực có vận tốc không xác định duy nhất hoặc chịu tác động của các ràng buộc gián đoạn. Bài viết này trình bày cơ sở toán học của ánh xạ đa trị, các định lý tồn tại nghiệm kinh điển, lý thuyết khả sinh, cấu trúc tập nghiệm và các ứng dụng then chốt trong lý thuyết điều khiển tối ưu cũng như hệ động lực lai.

Bản chất toán học và cấu trúc của bao hàm vi phân

Trong lý thuyết phương trình vi phân cổ điển, sự tiến triển của một hệ thống vật lý hoặc kỹ thuật được mô tả qua phương trình vi phân thường có dạng đạo hàm bậc nhất của biến trạng thái theo biến thời gian bằng một hàm vector đơn trị. Cụ thể, với biến trạng thái $x$ thuộc không gian Euclid và biến thời gian $t$ , phương trình vi phân cổ điển được viết dưới dạng:

\dot{x}(t) = f(t, x(t))

Trong biểu thức trên, với mỗi thời điểm $t$ và trạng thái $x(t)$ , hàm số $f$ chỉ cung cấp một vector vận tốc duy nhất. Tuy nhiên, trong thực tế kỹ thuật và tự nhiên, rất nhiều hệ thống vận hành với các quy luật không trơn hoặc có sự tham gia của các tham số điều khiển chưa xác định trước. Để vượt qua giới hạn của phương trình vi phân đơn trị, lý thuyết bao hàm vi phân thay thế hàm số $f$ bằng một ánh xạ đa trị $F$ , đưa đến hệ thức tổng quát có dạng:

\dot{x}(t) \in F(t, x(t))

Trong hệ thức này, với biến trạng thái $x(t)$ tại thời điểm $t$ , biểu thức $F(t, x(t))$ trả về một tập con chứa các vector khả dĩ trong không gian trạng thái. Như Aubin và Cellina (1984) đã chỉ ra, bao hàm vi phân không chỉ đơn thuần là việc thay đổi ký hiệu toán học mà là một bước chuyển biến về chất trong việc nghiên cứu các hệ động lực đa nghiệm, không tất định và có ràng buộc pha.

Khái niệm nghiệm trong bao hàm vi phân

Do vế phải là một tập hợp, khái niệm nghiệm cổ điển khả vi liên tục khắp nơi không còn phù hợp. Theo Filippov (1988), nghiệm của một bao hàm vi phân trên một khoảng thời gian được định nghĩa là một hàm vector liên tục tuyệt đối, ký hiệu là $x(t)$ , thỏa mãn điều kiện vi phân hầu khắp nơi theo độ đo Lebesgue. Cụ thể, đạo hàm thông thường $\dot{x}(t)$ tồn tại tại hầu hết các điểm thời gian và vector vận tốc tức thời này phải thuộc về tập hợp $F(t, x(t))$ .

Việc sử dụng lớp hàm liên tục tuyệt đối cho phép nghiệm có thể thay đổi vận tốc đột ngột tại các tập hợp điểm có độ đo Lebesgue bằng không, phản ánh đúng tính chất vật lý của các va chạm cơ học, hiện tượng ma sát khô hoặc hành vi chuyển mạch tức thời của các linh kiện bán dẫn trong mạch điện tử.

Các tính chất tô pô của ánh xạ đa trị

Để nghiên cứu tính đúng đắn về mặt toán học của bao hàm vi phân, ánh xạ đa trị $F$ cần thỏa mãn các điều kiện giải tích và tô pô nhất định. Các khái niệm tính liên tục đối với ánh xạ đa trị phức tạp hơn ánh xạ đơn trị và được phân chia thành hai khái niệm cơ bản:

Nửa liên tục trên (Upper Semicontinuity): Ánh xạ đa trị $F$ được gọi là nửa liên tục trên tại điểm trạng thái nếu với mọi tập mở chứa ảnh của điểm đó, luôn tồn tại một lân cận của trạng thái sao cho ảnh của mọi điểm trong lân cận đều nằm trọn vẹn trong tập mở đã cho. Tính chất này đảm bảo tập giá trị không bị nổ đột ngột khi trạng thái biến thiên nhỏ.
Nửa liên tục dưới (Lower Semicontinuity): Ánh xạ đa trị $F$ được gọi là nửa liên tục dưới nếu với mọi tập mở có giao khác rỗng với ảnh của trạng thái, luôn tồn tại lân cận sao cho ảnh của mọi điểm trong lân cận cũng giao khác rỗng với tập mở đó. Tính chất này đảm bảo các hướng vận tốc không bị biến mất đột ngột.
Tính chất compact và lồi: Trong hầu hết các bài toán tồn tại nghiệm, tập giá trị $F(t, x)$ được giả thiết là các tập lồi và compact khác rỗng. Tính lồi cho phép áp dụng các định lý điểm bất động trong không gian vô hạn chiều và đảm bảo tính đóng của tập các quỹ đạo nghiệm.

Các định lý tồn tại nghiệm cơ bản

Bài toán Cauchy đối với bao hàm vi phân yêu cầu tìm một quỹ đạo nghiệm $x(t)$ xuất phát từ một điểm trạng thái ban đầu cho trước tại thời điểm bắt đầu:

\dot{x}(t) \in F(t, x(t)), \quad x(0) = x_0

Khác với phương trình vi phân Lipschitz thông thường có nghiệm duy nhất theo định lý Picard - Lindelöf, bao hàm vi phân thường sinh ra một họ các quỹ đạo nghiệm xuất phát từ cùng một trạng thái ban đầu.

Định lý tồn tại nghiệm của Smirnov và Deimling

Đối với ánh xạ đa trị có giá trị lồi và compact, định lý tồn tại nghiệm được thiết lập dựa trên nguyên lý nửa liên tục trên. Smirnov (2001) đã tổng kết rằng nếu ánh xạ đa trị $F$ thỏa mãn tính chất nửa liên tục trên theo biến trạng thái, đo được theo biến thời gian, và có giá trị là các tập con lồi compact khác rỗng, đồng thời bị chặn bởi một hàm khả tích, thì bài toán Cauchy luôn có ít nhất một nghiệm liên tục tuyệt đối cục bộ.

Hơn thế nữa, Deimling (1992) đã chứng minh rằng tập hợp tất cả các nghiệm xuất phát từ cùng một điểm trạng thái ban đầu là một tập con compact và liên thông trong không gian Banach của các hàm liên tục trang bị chuẩn hội tụ đều. Tính chất liên thông của tập nghiệm ngụ ý rằng tập các điểm có thể đạt tới tại bất kỳ thời điểm tương lai nào cũng là một tập liên thông, không bị chia cắt thành các thành phần rời rạc.

Định lý làm giãn Filippov - Ważewski cho ánh xạ không lồi

Một thách thức lớn trong giải tích biến phân xuất hiện khi ánh xạ đa trị $F$ không có giá trị lồi. Trường hợp này xuất hiện phổ biến trong bài toán điều khiển tối ưu với tập điều khiển rời rạc hoặc dạng bang-bang. Khi giả thiết tính lồi bị vi phạm, bài toán Cauchy có thể không có nghiệm cổ điển hoặc tập nghiệm có thể không đóng.

Để giải quyết bài toán này, định lý làm giãn Filippov - Ważewski được Aubin và Cellina (1984) trình bày chi tiết. Định lý khẳng định rằng nếu ánh xạ đa trị $F$ liên tục Lipschitz theo biến trạng thái $x$ (theo khoảng cách Hausdorff) và có giá trị compact khác rỗng, thì tập nghiệm của bao hàm vi phân ban đầu là một tập con trù mật trong tập nghiệm của bao hàm vi phân đã được bao lồi hóa:

\dot{x}(t) \in \overline{\operatorname{conv}} \, F(t, x(t))

Kết quả này mang ý nghĩa thực tiễn sâu sắc: mọi quỹ đạo chuyển động lý thuyết tạo ra bởi việc làm lồi hóa vận tốc đều có thể được xấp xỉ với độ chính xác tùy ý bằng một quỹ đạo thực của hệ thống không lồi ban đầu thông qua kỹ thuật chuyển mạch vận tốc với tần số cao.

Lý thuyết khả sinh và tập bất biến

Trong nhiều bài toán thực tế, trạng thái của hệ thống không được phép di chuyển tự do trong toàn bộ không gian mà phải luôn nằm trong một miền an toàn cho trước, ví dụ giới hạn nhiệt độ của lò phản ứng hạt nhân, điện áp cho phép của lưới điện hoặc hành lang bay của máy bay không người lái. Lý thuyết khả sinh (viability theory) do Aubin phát triển năm 1984 cung cấp công cụ toán học giải quyết trọn vẹn bài toán này.

Miền khả sinh và nón tiếp tuyến Bouligand

Cho một tập đóng $K$ trong không gian trạng thái. Tập $K$ được gọi là một miền khả sinh đối với bao hàm vi phân nếu với mọi trạng thái ban đầu xuất phát bên trong tập này, luôn tồn tại ít nhất một quỹ đạo nghiệm của hệ nằm trọn vẹn bên trong tập $K$ trong suốt khoảng thời gian chuyển động.

Để đặc trưng hóa điều kiện khả sinh tại biên của miền trạng thái, Aubin (1984) sử dụng khái niệm nón tiếp tuyến Bouligand (contingent cone), ký hiệu là $T_K(x)$ . Với điểm trạng thái $x$ thuộc tập $K$ , nón tiếp tuyến Bouligand tập hợp tất cả các vector vận tốc giới hạn mà từ đó một chuyển động nhỏ sẽ không đưa trạng thái ra ngoài tập $K$ :

T_K(x) = \left\{ v \in \mathbb{R}^n \;\middle|\; \liminf_{h \to 0^+} \frac{\operatorname{dist}(x + hv, K)}{h} = 0 \right\}

Trong công thức toán học trên, ký hiệu $v$ biểu diễn vector vận tốc thử nghiệm, ký hiệu $h$ là bước thời gian tiến dần về không từ phía dương, và hàm số khoảng cách biểu thị khoảng cách Euclid từ điểm trạng thái dịch chuyển tới tập đóng đang xét.

Định lý khả sinh Aubin

Định lý khả sinh nổi tiếng của Aubin (1984) thiết lập điều kiện hình học cần và đủ để một tập đóng trở thành miền khả sinh. Cụ thể, định lý phát biểu rằng một tập đóng $K$ là miền khả sinh của bao hàm vi phân với ánh xạ đa trị lồi compact nếu và chỉ nếu tại mọi điểm trạng thái $x$ thuộc tập $K$ , giao của tập vận tốc khả dĩ và nón tiếp tuyến Bouligand luôn khác rỗng:

F(x) \cap T_K(x) \neq \emptyset, \quad \forall x \in K

Điều kiện này mang ý nghĩa hình học trực quan: tại bất kỳ vị trí nào trên biên của miền an toàn, hệ thống phải luôn sở hữu ít nhất một hướng vận tốc khả thi chỉ vào phía trong hoặc trượt dọc theo bề mặt tiếp tuyến của miền an toàn, ngăn không cho quỹ đạo trạng thái bị ép buộc văng ra ngoài phạm vi cho phép.

Các lớp bao hàm vi phân tiêu biểu

Tùy thuộc vào bản chất vật lý và cấu trúc giải tích của ánh xạ đa trị, bao hàm vi phân được phân chia thành nhiều lớp chuyên biệt với các phương pháp nghiên cứu tương ứng.

Bao hàm vi phân Filippov

Khi một hệ thống động lực được mô tả bởi phương trình vi phân có vế phải gián đoạn trên một mặt cong trơn, các nghiệm cổ điển không thể tiếp tục kéo dài khi chạm vào mặt gián đoạn. Filippov (1988) đã đề xuất phương pháp chính quy hóa bằng cách biến phương trình vi phân gián đoạn thành một bao hàm vi phân. Ánh xạ đa trị Filippov tại điểm gián đoạn được xây dựng bằng cách lấy bao lồi đóng của giới hạn các giá trị hàm số tại các điểm lân cận ngoài mặt gián đoạn:

F_F(t, x) = \bigcap_{\delta > 0} \bigcap_{\mu(N) = 0} \overline{\operatorname{conv}} \, f(t, (x + \delta B) \setminus N)

Trong đó, ký hiệu $B$ là hình cầu đơn vị đóng, $\delta$ là bán kính lân cận, và $N$ là tập con có độ đo Lebesgue $\mu(N)$ bằng không. Quy tắc này cho phép hệ thống trượt dọc theo mặt gián đoạn, giải thích hoàn hảo hiện tượng điều khiển trượt (sliding mode control) trong ngành kỹ thuật tự động hóa.

Quá trình quét Moreau

Được nhà cơ học Jean Jacques Moreau đề xuất, quá trình quét (Moreau sweeping process) là một lớp bao hàm vi phân đặc biệt mô tả chuyển động của trạng thái bị giới hạn trong một tập lồi biến thiên theo thời gian. Phương trình có dạng:

-\dot{x}(t) \in N_{C(t)}(x(t))

Trong đó, biểu thức $C(t)$ là tập lồi đóng biến thiên, và $N_{C(t)}(x(t))$ là nón pháp tuyến ngoài của tập lồi tại điểm trạng thái. Quá trình quét là công cụ toán học tiêu chuẩn để mô hình hóa hiện tượng va chạm không đàn hồi hoàn toàn, cơ cấu cam cơ học và các mạng lưới điện phi tuyến tính có diode lý tưởng.

Bao hàm vi phân dưới vi phân

Trong lý thuyết tối ưu hóa và giải tích lồi, khi một hàm thế năng không khả vi nhưng lồi, vector vận tốc của hệ tiêu tán năng lượng được điều khiển bởi toán tử dưới vi phân (subdifferential). Hệ thức vi phân có dạng:

\dot{x}(t) \in -\partial \varphi(x(t))

Với $\varphi$ là một hàm lồi thực sự nửa liên tục dưới và $\partial \varphi$ là dưới vi phân theo nghĩa phân tích lồi. Lớp bao hàm vi phân này đảm bảo tính đơn điệu cực đại của toán tử vế phải, từ đó suy ra tính duy nhất nghiệm của quỹ đạo chuyển động mặc dù vế phải là ánh xạ đa trị.

Bảng so sánh các cấu trúc toán học của hệ động lực

Để làm rõ vị trí của bao hàm vi phân trong bức tranh toàn cảnh của lý thuyết vi phân, bảng dưới đây đối chiếu các đặc tính cơ bản giữa phương trình vi phân thông thường, phương trình vi phân gián đoạn và bao hàm vi phân:

Đặc tính toán học	Phương trình vi phân cổ điển	Phương trình vi phân gián đoạn	Bao hàm vi phân
Dạng biểu thức vế phải	Vector đơn trị liên tục	Vector đơn trị có bước nhảy	Ánh xạ đa trị gán tập hợp
Không gian nghiệm	Khả vi liên tục khắp nơi	Liên tục tuyệt đối (đạo hàm liên tục từng khúc)	Liên tục tuyệt đối hầu khắp nơi
Tính duy nhất nghiệm	Được đảm bảo khi thỏa mãn Lipschitz	Thường không bảo đảm	Tập nghiệm đa trị, liên thông
Xử lý hiện tượng trượt	Không áp dụng được	Dễ dẫn đến bế tắc tại mặt gián đoạn	Mô tả tự nhiên qua bao lồi Filippov
Khung lý thuyết phân tích	Giải tích cổ điển và không gian hàm	Giải tích không trơn cục bộ	Giải tích đa trị và lý thuyết biến phân

Ứng dụng trong điều khiển tối ưu và hệ động lực lai

Lý thuyết bao hàm vi phân phát triển mạnh mẽ nhờ những đóng góp quyết định trong việc giải quyết các bài toán công nghệ kỹ thuật và lý thuyết điều khiển tự động hiện đại.

Mô hình hóa bài toán điều khiển tối ưu

Trong lý thuyết điều khiển tối ưu truyền thống, một hệ thống điều khiển được mô tả bởi phương trình vi phân phụ thuộc vào biến trạng thái $x$ và biến điều khiển $u$ , trong đó biến điều khiển bị ràng buộc phải thuộc một tập hợp cho trước:

\dot{x}(t) = g(t, x(t), u(t)), \quad u(t) \in U(t, x(t))

Bằng cách định nghĩa ánh xạ đa trị vận tốc thông qua ảnh của tập điều khiển qua hàm động lực học:

F(t, x) = \left\{ g(t, x, u) \;\middle|\; u \in U(t, x) \right\}

Hệ điều khiển có ràng buộc được chuyển hóa hoàn toàn thành một bao hàm vi phân thuần túy không còn chứa biến điều khiển hiển nhiên. Như Ioffe (2020) đã phân tích trên tạp chí Vietnam Journal of Mathematics, cách tiếp cận qua bao hàm vi phân cho phép thiết lập nguyên lý cực đại Pontryagin và các điều kiện tối ưu cần bậc nhất mà không đòi hỏi tính khả vi trơn của hàm số, đồng thời xử lý trọn vẹn các bài toán điều khiển có ràng buộc pha phức tạp.

Hệ động lực lai

Hệ động lực lai (hybrid dynamical systems) là các hệ thống trong đó chuyển động liên tục của trạng thái xen kẽ với các bước nhảy rời rạc tại những thời điểm biến cố. Điển hình là các hệ thống điều khiển số nhúng, bộ chuyển đổi nguồn điện tử công suất và hệ thống điều hòa nhiệt độ thông minh.

Goebel, Sanfelice và Teel (2009) đã xây dựng một khung lý thuyết tổng quát cho hệ động lực lai dựa trên nền tảng bao hàm vi phân và bao hàm sai phân. Quá trình biến thiên trạng thái được mô tả qua hai quan hệ:

Dòng chảy liên tục (Flow): Diễn ra khi trạng thái thuộc tập dòng chảy, thỏa mãn bao hàm vi phân $\dot{x} \in F(x)$ .
Bước nhảy rời rạc (Jump): Diễn ra khi trạng thái chạm vào tập bước nhảy, thỏa mãn bao hàm sai phân $x^+ \in G(x)$ .

Nhờ sử dụng bao hàm vi phân cho pha dòng chảy, khung giải tích của Goebel, Sanfelice và Teel (2009) đảm bảo tính ổn định tiệm cận của hệ thống có tính bền vững tự nhiên trước các nhiễu đo lường nhỏ và sai số mô hình, khắc phục triệt để nhược điểm mất ổn định của các mô hình hệ lai giải tích cổ điển.

Hạn chế lý thuyết và các hướng nghiên cứu hiện đại

Mặc dù bao hàm vi phân cung cấp năng lực mô hình hóa mạnh mẽ, việc ứng dụng và phân tích lý thuyết này vẫn đối mặt với những thách thức nội tại đáng kể:

Sự mất mát tính duy nhất nghiệm: Do nghiệm của bao hàm vi phân là một họ các quỹ đạo, việc thiết kế các bộ quan sát trạng thái (state observers) và bộ lọc ước lượng đòi hỏi các công cụ xác định tập hợp (set-membership estimation) thay vì các bộ lọc điểm như bộ lọc Kalman truyền thống.
Độ phức tạp tính toán xấp xỉ số: Việc giải số bao hàm vi phân đòi hỏi xấp xỉ các tập hợp đa diện hoặc ellipsoid trong không gian nhiều chiều. Khối lượng tính toán tăng theo hàm mũ đối với số chiều trạng thái, gây trở ngại khi triển khai các thuật toán điều khiển dự báo mô hình theo thời gian thực.
Bao hàm vi phân ngẫu nhiên và vô hạn chiều: Nghiên cứu mở rộng bao hàm vi phân sang không gian Banach vô hạn chiều để mô tả các phương trình đạo hàm riêng phi trơn, cũng như việc kết hợp toán tử Ito ngẫu nhiên vào ánh xạ đa trị đang là các chủ đề nghiên cứu mũi nhọn của giải tích biến phân đương đại.