Tổng quan về luận án

Luận án tiến sĩ Toán học chuyên ngành Giải tích (hướng Phương trình vi phân và tích phân) được thực hiện tại Bộ môn Giải tích, khoa Toán-Tin, Trường Đại học Sư phạm Hà Nội, dưới sự hướng dẫn của PGS.TS Lê Văn Hiện và TS. Trịnh Tuấn Anh, tập trung nghiên cứu tính ổn định của các hệ phương trình vi phân có trễ và ứng dụng trong các mô hình sinh thái. Bối cảnh khoa học của công trình nằm ở giao điểm giữa lý thuyết định tính phương trình vi phân và lý thuyết điều khiển hệ thống — nơi mà, như tác giả nhấn mạnh, "sự xuất hiện của các độ trễ này ảnh hưởng đến dáng điệu tiệm cận nghiệm của hệ nói chung và tính ổn định... nói riêng".

Khoảng trống nghiên cứu được xác định rất cụ thể: các kết quả đã công bố về ổn định và đồng bộ mạng nơron hầu như chỉ áp dụng cho trọng số kết nối hằng và trễ bị chặn, trong khi trễ tỉ lệ (proportional delay) — với dạng qt, 0 < q < 1, thuộc lớp trễ biến thiên không bị chặn — lại phổ biến trong các hệ có cấu trúc mạng nhiều tầng. Tác giả khẳng định một khoảng trống then chốt: "Trước bài báo [1] trong Danh mục công bố của luận án này, chúng tôi chưa tìm thấy một kết quả nghiên cứu nào đề cập đến tính ổn định hữu hạn của mô hình mạng nơron với hệ số biến thiên và trễ tỉ lệ."

Ba câu hỏi nghiên cứu trung tâm được đánh số theo cấu trúc luận án: (1) Điều kiện nào đảm bảo tính ổn định trong thời gian hữu hạn (finite-time stability, FTS) của mạng nơron Hopfield hệ số biến thiên với đa trễ tỉ lệ không đồng nhất? (2) Tính tiêu hao toàn cục của lớp hệ này được thiết lập ra sao trong cả hai trường hợp hệ số phản hồi chính quy và suy biến? (3) Nghiệm tuần hoàn dương của mô hình Nicholson có trễ với hàm suy giảm phi tuyến tồn tại, duy nhất và hút toàn cục dưới điều kiện nào?

Khung lý thuyết dựa trên lý thuyết ổn định Lyapunov, lý thuyết M-ma trận (theo chuyên khảo của Berman–Plemmons về ma trận không âm), bất đẳng thức Halanay cùng các cải biên, và nguyên lí so sánh bằng bất đẳng thức vi-tích phân. Kết quả của luận án được công bố trong 03 bài báo trên các tạp chí quốc tế thuộc danh mục ISI — một chỉ dấu định lượng rõ ràng về chất lượng khoa học của công trình gồm 4 chương với hơn 80 trang nội dung chính.

Literature Review và Positioning

Luận án tổng hợp ba dòng nghiên cứu lớn. Dòng thứ nhất về ổn định thời gian hữu hạn khởi nguồn từ nửa sau thế kỷ XX (tài liệu [20]) và được hệ thống hóa bởi Amato và cộng sự trong lý thuyết điều khiển hữu hạn thời gian [2]; điểm mấu chốt là FTS và ổn định Lyapunov (LS) "là hai khái niệm độc lập theo nghĩa một hệ là FTS nhưng có thể không ổn định theo Lyapunov và ngược lại (xem phản ví dụ trong [16])". Dòng thứ hai về tính tiêu hao dựa trên bất đẳng thức Halanay và các mở rộng gần đây [19, 36], trong đó [19] đưa ra khái niệm hội tụ mũ suy rộng cho phép bỏ điều kiện xác định âm đều của hệ số tiêu hao. Dòng thứ ba về mô hình Nicholson bắt nguồn từ công trình kinh điển của Nicholson (1954) mô tả quần thể ve châu Úc bằng phương trình N′(t) = −αN(t) + βN(t−τ)e^(−γN(t−τ)), với các phát triển gần đây về nghiệm tuần hoàn [17, 22, 25, 26] và yếu tố thu hoạch [10, 23, 27, 35].

Tranh luận học thuật hiện rõ ở hai điểm đối lập. Thứ nhất, về phương pháp xử lý trễ tỉ lệ: trường phái đổi biến x̂(ξ) = x(e^ξ) chuyển về trễ hằng [43, 44, 45] đối lập với cách tiếp cận so sánh trực tiếp — luận án chỉ ra "cách tiếp cận bằng phương pháp đổi biến... không thể mở rộng được cho trường hợp hệ số biến thiên". Thứ hai, về mô hình Nicholson: hầu hết kết quả dùng tốc độ tử vong tuyến tính, trong khi các nhà hải dương học [3] chỉ ra rằng mô hình tuyến tính "thường chỉ đúng đối với các quần thể có mật độ thấp", đòi hỏi các dạng phi tuyến type-I (D(N) = a − be^(−N)) hoặc type-II (D(N) = aN/(b+N)).

So với các nghiên cứu quốc tế, luận án định vị vượt trội trực tiếp: các điều kiện trong [42] (dùng hàm Lyapunov toàn phương) và [44, 45] (bất đẳng thức ma trận với ràng buộc γ > 1 + 1/√q) được chứng minh là ngặt hơn hẳn điều kiện phổ M-ma trận của luận án, với phản ví dụ số cụ thể trong đó ρ(K) > 1 với mọi σ ∈ (1, 3) khiến tiêu chuẩn của [42] mất hiệu lực trong khi tiêu chuẩn của luận án vẫn áp dụng được.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng lý thuyết ổn định Lyapunov–Krasovskii theo ba hướng. Một là, mở rộng khái niệm FTS của Dorato–Amato từ hệ tuyến tính trễ bị chặn [16] sang lớp hệ phi tuyến không dừng với trễ tỉ lệ không đồng nhất q_ij ∈ (0,1) — Định lí 2.1 thiết lập điều kiện (M − γI)ξ ≺ 0 qua vectơ dương ξ và hằng số đồng đều C(ξ) = ξ⁺(ξ₊)⁻¹. Hai là, thay thế khuôn khổ LMI truyền thống bằng điều kiện phổ của M-ma trận, kế thừa Mệnh đề tương đương năm phát biểu (i)–(v) từ lý thuyết Berman–Plemmons [5]. Ba là, với mô hình Nicholson, luận án chuyển từ khung suy giảm tuyến tính sang khung suy giảm phi tuyến D(t, N) = a(t) − b(t)e^(−N) phụ thuộc thời gian — một mở rộng thực chất so với [6], nơi hệ số còn là hằng số.

Khung phân tích độc đáo

Khung phân tích tích hợp ba công cụ lý thuyết có tên tuổi: (i) lý thuyết M-ma trận không suy biến, (ii) bất đẳng thức Halanay cải biên với hội tụ mũ suy rộng [19], và (iii) đạo hàm Dini trên bên phải D⁺v(t) kết hợp hàm dấu suy rộng σ_f. Điểm mới về phương pháp là kỹ thuật so sánh bằng hàm phụ trợ dạng v_i(t) = ξ_l⁻¹‖z⁰‖∞ξ_i·e^(γt) và lập luận phản chứng với tham số λ > 1, tránh hoàn toàn phép đổi biến mũ vốn phá vỡ cấu trúc hệ số biến thiên. Điều kiện biên được phát biểu tường minh: q_ij ∈ (0,1), hàm kích hoạt Lipschitz toàn cục (giả thiết A2.1, A3.1), và hệ số bị chặn (A3.2).

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

Thiết kế nghiên cứu

Là công trình toán học thuần túy, luận án theo hệ hình chứng minh chặt chẽ tiên đề–suy diễn (mathematical rigor), tương ứng với lập trường nhận thức luận duy lý. Thiết kế đa tầng thể hiện qua bốn chương: Chương 1 xây dựng nền tảng (M-ma trận, ổn định Lyapunov, FTS, tính tiêu hao); Chương 2–4 lần lượt giải quyết ba bài toán mở, mỗi chương dựa trên đúng một bài báo ISI trong danh mục công bố.

Quy trình nghiên cứu rigorous

Mỗi kết quả chính tuân theo giao thức: đặt giả thiết định lượng (A2.1–A2.2, A3.1–A3.4), phát biểu định lí, chứng minh bằng phản chứng kết hợp đánh giá vi phân, rồi kiểm chứng bằng ví dụ số và mô phỏng. Tính "triangulation" phương pháp thể hiện ở việc mỗi định lí được đối chiếu chéo với ít nhất hai tiêu chuẩn công bố trước ([8], [42], [44, 45]) để chứng minh tính bao hàm: "Định lí 2.1 suy ra Định lí 2 trong [8]".

Dữ liệu và phân tích

Thay cho dữ liệu thống kê, luận án cung cấp bằng chứng số tường minh. Ví dụ minh họa Chương 2 dùng A = diag{4, 3}, ma trận M có det(M) = 0 (M-ma trận suy biến), với r₁ = 1, r₂ = 1.5, và chứng minh hệ ổn định hữu hạn với mọi T ≤ T_max = 18; mô phỏng thực hiện trên đoạn [0, 15]. Phản ví dụ về phương trình x′(t) = −(1 + e^(−t))x(t) + (1 − e^(−t))x(0.5t) chỉ ra nghiệm thỏa lim inf x(t) ≥ x₀e⁻² > 0 — FTS với mọi tham số nhưng không ổn định tiệm cận Lyapunov. Thuật toán tính tốc độ đồng bộ σ_max được cho tường minh qua nghiệm dương duy nhất của phương trình H_i(σ) = 0.

Phát hiện đột phá và implications

Những phát hiện then chốt

Bốn phát hiện chính, mỗi phát hiện có bằng chứng định lí kèm theo. Thứ nhất, điều kiện phổ M-ma trận (M − γI)ξ ≺ 0 kết hợp ràng buộc r₂/r₁ đảm bảo FTS của mạng Hopfield hệ số biến thiên trễ tỉ lệ (Định lí 2.1) — kết quả đầu tiên thế giới cho lớp hệ này theo khảo cứu của tác giả. Thứ hai, phát hiện phản trực giác: khi −M là M-ma trận không suy biến, nghiệm đồng bộ với tốc độ lũy thừa ‖x(t) − x*(t)‖∞ ≤ β‖x(0) − x*(0)‖∞/(1+t)^σ₀ — tốc độ đa thức chứ không phải mũ, phản ánh bản chất trễ không bị chặn τ_ij(t) = (1−q_ij)t → ∞. Thứ ba, tính tiêu hao toàn cục được chứng minh cả khi hệ số phản hồi suy biến (inf a_i(t) = 0, giả thiết A3.4 với hàm tỉ lệ ϕ(t)) — trường hợp mà toàn bộ hướng tiếp cận Halanay cổ điển thất bại vì đòi hỏi α(t) ≤ −σ < 0. Thứ tư, mô hình Nicholson phi tuyến (4) với p trễ biến thiên có nghiệm tuần hoàn dương duy nhất hút toàn cục; áp dụng cho hệ số hằng thu được điểm cân bằng dương hút toàn cục.

So với nghiên cứu trước, phát hiện thứ hai bao hàm Định lí 2 trong [8] về ổn định đa thức, còn điều kiện tiêu hao (3.15) được chứng minh yếu hơn hẳn điều kiện (3.16) chứa nhân tử 1/(1 − τ_M) của cách tiếp cận trễ bị chặn.

Implications đa chiều

Về lý thuyết, luận án đóng góp đồng thời cho lý thuyết ổn định hữu hạn (mở rộng sang hệ phi tuyến không dừng) và lý thuyết tiêu hao (khái niệm tập hút mũ suy rộng với σ(t) = α̂ln(1+t)). Về phương pháp, kỹ thuật so sánh M-ma trận có thể chuyển giao trực tiếp sang các lớp hệ mạng khác như mạng Cohen–Grossberg hay hệ dương có trễ. Về ứng dụng, kết quả Chương 2–3 phục vụ thiết kế mạng nơron nhân tạo nhiều tầng trong xử lý tín hiệu và nhận dạng mẫu, nơi băng thông hạn chế sinh trễ tỉ lệ; kết quả Chương 4 cung cấp cơ sở toán học cho quản lý khu bảo tồn biển và ngư nghiệp — bối cảnh mà các nhà hải dương học đã chỉ ra tính tất yếu của suy giảm phi tuyến [3]. Điều kiện khả kiểm bằng số (chỉ cần giải bất phương trình tuyến tính Mξ ≺ 0 thay vì LMI) hạ đáng kể chi phí tính toán khi kiểm định ổn định.

Limitations và Future Research

Luận án thừa nhận và thể hiện bốn giới hạn cụ thể. Một, các điều kiện đưa ra là điều kiện đủ, chưa phải điều kiện cần và đủ — khoảng cách với tiêu chuẩn tối ưu chưa được định lượng. Hai, FTS trên [0, T] không kéo theo ổn định tiệm cận Lyapunov, như chính phản ví dụ (2.20) của luận án minh chứng; ứng dụng dài hạn cần kiểm tra bổ sung. Ba, hàm kích hoạt phải Lipschitz toàn cục — loại trừ các lớp kích hoạt không trơn hoặc gián đoạn. Bốn, mô hình Nicholson chỉ xét dạng suy giảm type-I D(t,N) = a(t) − b(t)e^(−N); trường hợp type-II aN/(b+N) chưa được giải quyết trọn vẹn.

Từ đó, năm hướng nghiên cứu tương lai cụ thể: (1) mở rộng sang mô hình Nicholson với suy giảm type-II và yếu tố thu hoạch; (2) phát triển FTS cho mạng nơron trễ tỉ lệ chịu nhiễu ngẫu nhiên; (3) bài toán điều khiển hữu hạn thời gian (finite-time stabilization) cho lớp hệ đã xét; (4) hệ nơron phân thứ (fractional-order) với trễ tỉ lệ; (5) tối ưu hóa tốc độ đồng bộ σ_max thành bài toán quy hoạch trên nón M-ma trận.

Tác động và ảnh hưởng

Về học thuật, ba bài báo ISI trong các lĩnh vực phương trình vi phân và mạng nơron — nơi các công trình về trễ tỉ lệ ([44, 45], Zhou và cộng sự) thu hút trích dẫn cao — có tiềm năng trích dẫn ổn định trong cộng đồng nghiên cứu ổn định hệ trễ vốn được luận án mô tả là "vấn đề nghiên cứu có tính thời sự thu hút sự quan tâm của nhiều tác giả trong và ngoài nước". Về công nghiệp, các tiêu chuẩn kiểm định ổn định khả tính phục vụ thiết kế mạng nơron trong trí tuệ nhân tạo, ước lượng tham số và xử lý tín hiệu số. Về chính sách, mô hình Nicholson phi tuyến cung cấp công cụ dự báo cho quản lý nguồn lợi thủy sản và khu bảo tồn biển ở cấp quốc gia. Tầm quốc tế thể hiện qua việc kết quả trực tiếp cải tiến các công bố của các nhóm nghiên cứu quốc tế [8, 42, 43, 44, 45].

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

Nghiên cứu sinh ngành Giải tích và Toán ứng dụng nhận được ba bài toán mở đã giải cùng năm hướng phát triển tiếp; kỹ thuật so sánh M-ma-trận trong Chương 2–3 là khuôn mẫu chứng minh có thể tái sử dụng. Nhà nghiên cứu lý thuyết điều khiển hưởng lợi từ tiêu chuẩn FTS mới không cần giải LMI. Kỹ sư thiết kế mạng nơron nhận điều kiện ổn định kiểm tra được bằng tính toán ma trận cơ bản với thuật toán tường minh cho σ_max. Nhà quản lý sinh thái và thủy sản có mô hình Nicholson tổng quát hóa với p trễ biến thiên và hệ số tuần hoàn, phù hợp biến động mùa vụ thực tế.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất? Mở rộng khái niệm FTS của Dorato–Amato sang hệ phi tuyến không dừng trễ tỉ lệ bằng lý thuyết M-ma-trận — kết quả đầu tiên cho lớp hệ này trước bài báo [1] của tác giả.

2. Đổi mới phương pháp so với nghiên cứu trước? So với phép đổi biến mũ của [43, 44, 45] (chỉ dùng được cho hệ số hằng) và hàm Lyapunov toàn phương của [42] (cho điều kiện ρ(K) > 1 mất hiệu lực), nguyên lí so sánh trực tiếp của luận án xử lý được hệ số biến thiên và cho điều kiện yếu hơn chứng minh được.

3. Phát hiện bất ngờ nhất? Phản ví dụ (2.20): một hệ ổn định hữu hạn với mọi bộ (r₁, r₂, T) nhưng nghiệm hội tụ về x₀e⁻² > 0, tức không ổn định tiệm cận Lyapunov — khẳng định tính độc lập triệt để của hai khái niệm.

4. Giao thức tái lập có được cung cấp? Có: toàn bộ giả thiết được đánh số, thuật toán tính σ_max ba bước tường minh, dữ liệu ví dụ số (A = diag{4,3}, r₁ = 1, r₂ = 1.5, T_max = 18) đầy đủ để tái lập mô phỏng.

5. Chương trình nghiên cứu dài hạn? Năm hướng nêu trên (Nicholson type-II, nhiễu ngẫu nhiên, điều khiển hữu hạn thời gian, hệ phân thứ, tối ưu tốc độ đồng bộ) tạo thành lộ trình phát triển tự nhiên cho thập kỷ tới của hướng nghiên cứu.

Kết luận

Luận án hoàn thành sáu đóng góp cụ thể: (1) tiêu chuẩn phổ M-ma trận đầu tiên cho FTS của mạng Hopfield hệ số biến thiên trễ tỉ lệ; (2) định lí đồng bộ nghiệm tốc độ lũy thừa kèm thuật toán tính σ_max; (3) chứng minh tiêu hao toàn cục trường hợp hệ số phản hồi chính quy qua tập hút mũ suy rộng; (4) mở rộng sang trường hợp hệ số phản hồi suy biến — vượt ngoài phạm vi mọi biến thể Halanay hiện có; (5) chứng minh tồn tại, bền vững và tiêu hao đều của nghiệm dương mô hình Nicholson suy giảm phi tuyến; (6) định lí tồn tại duy nhất và hút toàn cục của nghiệm tuần hoàn dương kèm hệ quả cho điểm cân bằng. Công trình mở ra ít nhất ba dòng nghiên cứu mới — FTS cho hệ trễ không bị chặn tổng quát, tiêu hao với hệ số suy biến, và động lực học Nicholson phi tuyến không dừng — đồng thời thiết lập chuẩn so sánh mới mà các tiêu chuẩn quốc tế trước đó [8, 42, 44, 45] đều được chứng minh là trường hợp riêng hoặc hệ quả ngặt hơn. Với 03 công bố ISI và các báo cáo tại Viện Toán học thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, luận án khẳng định vị thế của trường phái phương trình vi phân hàm Việt Nam trong dòng nghiên cứu quốc tế về hệ động lực có trễ.