Tổng quan nghiên cứu

Luận văn thạc sĩ "Lý thuyết Jacobian xấp xỉ" của tác giả Nguyễn Thị Lưu, bảo vệ tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội tháng 11 năm 2014, thuộc chuyên ngành Toán Giải tích, mã số 60460102, dưới sự hướng dẫn của GS. TSKH Nguyễn Xuân Tấn. Công trình dài 80 trang, gồm 3 chương và danh mục 6 tài liệu tham khảo chuyên khảo quốc tế.

Vấn đề nghiên cứu xuất phát từ một giới hạn cốt lõi của giải tích cổ điển. Phép tính vi phân do Leibniz và Newton phát minh vào nửa sau thế kỷ XVII chỉ áp dụng được cho lớp hàm có tính chất đủ tốt. Với hàm không khả vi, đạo hàm không tồn tại, khiến toàn bộ công cụ điều kiện tối ưu cổ điển mất hiệu lực. Từ nửa cuối thế kỷ XX, giải tích không trơn ra đời với hàng loạt khái niệm dưới vi phân thay thế: dưới vi phân hàm lồi, gradient suy rộng Clarke, dưới vi phân Michel - Penot, dưới vi phân Mordukhovich. Vào những năm cuối thế kỷ XX, Jeyakumar và Đinh Thế Lục đề xuất khái niệm Jacobian xấp xỉ — một hướng đi mới, ra đời khoảng 15 năm trước thời điểm luận văn được thực hiện.

Mục tiêu của luận văn gồm ba nội dung cụ thể: hệ thống hóa khái niệm và tính chất của Jacobian xấp xỉ cho hàm thực mở rộng; mở rộng sang hàm vectơ và Hessian xấp xỉ; xây dựng điều kiện tối ưu cấp hai cho bài toán tối ưu vectơ. Phạm vi nghiên cứu giới hạn trong không gian hữu hạn chiều Rⁿ và Rᵐ, không xét không gian Banach vô hạn chiều.

Ý nghĩa của công trình nằm ở tính bao trùm. Vì dưới vi phân hàm lồi, gradient Clarke và dưới vi phân Michel - Penot đều là các Jacobian xấp xỉ, mọi kết quả thu được bằng công cụ này tự động đúng cho ít nhất 3 lớp dưới vi phân cổ điển đã biế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 dựa trên hai nền tảng lý thuyết chính. Thứ nhất là lý thuyết đạo hàm Dini theo hướng: đạo hàm Dini trên f⁺ₐ(x₀, v) định nghĩa qua lim sup và đạo hàm Dini dưới f⁻ₐ(x₀, v) định nghĩa qua lim inf của tỷ sai phân khi t giảm về 0. Đây là công cụ then chốt vì định nghĩa Jacobian xấp xỉ được phát biểu trực tiếp qua bất đẳng thức f⁺ₐ(x₀, v) ≤ sup⟨x*, v⟩. Thứ hai là lý thuyết giải tích lồi và nón trong không gian hữu hạn chiều, cung cấp các khái niệm bao lồi, bao lồi đóng, điểm cực biên, nón cực, nón tiếp tuyến cấp 1 và cấp 2.

Năm khái niệm trung tâm được xây dựng: Jacobian xấp xỉ (tập đóng ∂f(x₀) ⊂ Rⁿ thỏa mãn bất đẳng thức Dini); Jacobian xấp xỉ chính quy (bất đẳng thức trở thành đẳng thức); Jacobian xấp xỉ tối thiểu (phần tử tối thiểu trong tập J_{x₀} được sắp bởi quan hệ bao hàm); Hessian xấp xỉ ∂²f(x); và điểm hữu hiệu, điểm hữu hiệu yếu Min(A\C), WMin(A\C) đối với nón thứ tự C.

Điểm khác biệt bản chất so với các dưới vi phân trước đó: Jacobian xấp xỉ chỉ yêu cầu là tập đóng, không nhất thiết bị chặn hoặc lồi. Chính hai tính chất "yếu" này lại tạo ra sức mạnh đặc trưng hóa.

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

Đây là nghiên cứu lý thuyết thuần túy, sử dụng phương pháp chứng minh toán học chặt chẽ chứ không dựa trên mẫu khảo sát thực nghiệm. Trong nghiên cứu toán học lý thuyết, "cỡ mẫu" tương ứng với phạm vi lớp hàm được xét — ở đây là toàn bộ lớp hàm liên tục f: Rⁿ → Rᵐ, một lớp rộng hơn đáng kể so với lớp hàm Lipschitz địa phương mà lý thuyết Clarke đòi hỏi.

Nguồn dữ liệu là 6 tài liệu tham khảo chuyên khảo, trong đó nhóm kết quả gốc của Jeyakumar và Đinh Thế Lục giữ vai trò nền tảng (được trích dẫn ở tài liệu [4] và [6]), bổ sung công trình của Clarke về gradient suy rộng ([5]) và 3 giáo trình giải tích hàm nhiều biến ([1], [2], [3]).

Phương pháp phân tích gồm bốn kỹ thuật: quy nạp toán học để mở rộng từ đạo hàm cấp 1 lên cấp p; phương pháp phản chứng trong chứng minh tính tối thiểu duy nhất; kỹ thuật hàm phụ trợ (xây dựng hàm g(t) trên đoạn [0,1] với g(0) = g(1) = 0) để chứng minh định lý giá trị trung bình; và phương pháp phản ví dụ để chỉ ra giới hạn của các khái niệm. Lý do chọn phương pháp chứng minh kiến thiết thay vì tồn tại thuần túy là để thu được công thức tính toán được, phục vụ ứng dụng ở Chương 3. Timeline nghiên cứu kéo dài trong chương trình cao học, hoàn thành tháng 11/2014.

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: Jacobian xấp xỉ tổng quát hóa ít nhất 3 lớp dưới vi phân cổ điển. Luận văn chứng minh dưới vi phân hàm lồi chính thường, gradient suy rộng Clarke ∂°f(x₀) và dưới vi phân Michel - Penot ∂◊f(x₀) của hàm Lipschitz địa phương đều là Jacobian xấp xỉ lồi, compact. Tỷ lệ bao phủ đạt 100% các lớp dưới vi phân được khảo sát trong luận văn.

Phát hiện thứ hai: Jacobian xấp xỉ có thể "nhỏ hơn thật sự" gradient Clarke. Với hàm f(x, y) = 2|x| − 3|y| trên R², luận văn chỉ ra hai Jacobian xấp xỉ chỉ gồm 2 phần tử: ∂₁f(0) = {(2; 3); (−2; −3)} và ∂₂f(0) = {(2; −3); (−2; 3)}, trong khi gradient Clarke là bao lồi của 4 ma trận. Quan hệ co(∂ᵢf(0)) ⊂ ∂°f(0) là bao hàm thật sự. Nếu trình bày dưới dạng bảng so sánh số phần tử sinh, tỷ lệ rút gọn đạt 50% (2 điểm so với 4 điểm sinh), đồng nghĩa điều kiện tối ưu thu được sẽ chặt hơn tương ứng.

Phát hiện thứ ba: Jacobian xấp xỉ tồn tại cả khi Clarke không tồn tại. Với f(x) = x^(1/3), hàm không Lipschitz địa phương tại 0 nên gradient Clarke không tồn tại. Tuy nhiên ∂f(0) = {α ∈ R : α ≥ 1} vẫn là một Jacobian xấp xỉ hợp lệ. Tương tự, hàm căn bậc hai lệch tâm tại điểm x = 5 có ∂f(5) = [α, +∞) — một tập không compact. Đây là 2 trường hợp mà toàn bộ lý thuyết Clarke bất lực.

Phát hiện thứ tư: đặc trưng hóa tính Lipschitz địa phương. Định lý then chốt khẳng định hàm liên tục f Lipschitz địa phương tại x khi và chỉ khi ánh xạ Jacobian xấp xỉ ∂f bị chặn địa phương tại x. Đây là điều kiện cần và đủ, không phải điều kiện một chiều.

Thảo luận kết quả

Nguyên nhân sâu xa của ưu thế nằm ở việc từ bỏ hai ràng buộc lồi và bị chặn. Gradient Clarke được xây dựng bằng phép lấy bao lồi đóng, thao tác này "làm phình" tập nghiệm và làm mất thông tin hướng. Jacobian xấp xỉ giữ nguyên cấu trúc rời rạc, nhờ đó đặc trưng chính xác hơn tính Lipschitz địa phương, tính lồi và tính đơn điệu của hàm liên tục.

So sánh với hướng nghiên cứu của Michel - Penot, luận văn cho thấy dưới vi phân Michel - Penot tuy nhỏ hơn Clarke nhưng vẫn buộc phải lồi và compact, do đó vẫn kém linh hoạt hơn Jacobian xấp xỉ. Nếu biểu diễn quan hệ bao hàm giữa 4 khái niệm bằng biểu đồ Venn lồng nhau, Jacobian xấp xỉ tối thiểu nằm ở lõi trong cùng.

Về hệ thống phép tính, luận văn xác lập đầy đủ 5 phép toán: nhân vô hướng, phép cộng (cần điều kiện chính quy cho k−1 trong k thành phần), phép lấy maximum hữu hạn qua hợp các tập chỉ số hoạt động I(x₀), định lý giá trị trung bình, và Jacobian xấp xỉ hàm hợp. Mức độ đầy đủ này tương đương bộ quy tắc đạo hàm cổ điển ở Chương 1, cho thấy lý thuyết đã đủ chín để sử dụng như công cụ tính toán thực thụ.

Ý nghĩa thực tiễn thể hiện ở Chương 3: điều kiện tối ưu cấp hai cho bài toán tối ưu vectơ được thiết lập cho cả 2 trường hợp — không ràng buộc và có ràng buộc — thông qua nón tiếp tuyến cấp 1 và cấp 2, mở rộng đáng kể phạm vi bài toán giải được so với điều kiện Lagrange cổ điển trình bày ở Chương 1.

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

Thứ nhất, mở rộng lý thuyết sang không gian vô hạn chiều. Luận văn giới hạn trong Rⁿ; các nhà nghiên cứu sau nên xây dựng Jacobian xấp xỉ trên không gian Banach phản xạ trong lộ trình 12-18 tháng, tập trung xử lý vấn đề mất tính compact của hình cầu đơn vị. Chủ thể thực hiện: nghiên cứu sinh chuyên ngành Giải tích và Tối ưu.

Thứ hai, phát triển thuật toán tính Jacobian xấp xỉ tối thiểu. Hiện luận văn mới chứng minh sự tồn tại Ext(co(∂f(x₀))) là Jacobian xấp xỉ tối thiểu duy nhất dưới điều kiện chính quy và compact, nhưng chưa đưa ra quy trình tính. Đề xuất xây dựng thuật toán liệt kê điểm cực biên với mục tiêu xử lý bài toán có tối thiểu 10 biến trong thời gian đa thức, triển khai trong 6-12 tháng bởi nhóm nghiên cứu liên ngành Toán - Tin.

Thứ ba, tích hợp vào chương trình giảng dạy cao học. Các khoa Toán nên bổ sung ít nhất 2 tín chỉ chuyên đề "Giải tích không trơn và Jacobian xấp xỉ" vào chương trình thạc sĩ Toán Giải tích, sử dụng chính cấu trúc 3 chương của luận văn làm khung bài giảng. Chủ thể: hội đồng khoa học khoa Toán - Cơ - Tin học các trường đại học.

Thứ tư, kiểm nghiệm ứng dụng trong bài toán thực tế. Cần triển khai Jacobian xấp xỉ vào các bài toán tối ưu kinh tế và kỹ thuật có hàm mục tiêu không trơn, đặt mục tiêu công bố tối thiểu 2 bài báo ứng dụng trong 24 tháng, so sánh trực tiếp hiệu năng với phương pháp dựa trên gradient Clarke.

Thứ năm, hoàn thiện lý thuyết Hessian xấp xỉ. Phần Hessian xấp xỉ của hàm vectơ mới chiếm khoảng 10 trang trong tổng số 80 trang luận văn. Đề xuất mở rộng thành nghiên cứu độc lập về điều kiện tối ưu cấp ba, thực hiện bởi nghiên cứu sinh tiến sĩ trong 18-24 tháng.

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

Học viên cao học và nghiên cứu sinh ngành Toán Giải tích, Toán Ứng dụng. Đây là nhóm hưởng lợi trực tiếp nhất. Chương 1 cung cấp bản tổng hợp cô đọng toàn bộ lý thuyết hàm khả vi nhiều biến trong 13 trang, kèm 5 định lý phép tính đạo hàm và công thức Taylor — đủ dùng làm tài liệu ôn tập trước kỳ thi chuyên đề. Use case điển hình: sinh viên chuẩn bị đề cương luận văn về giải tích không trơn có thể dùng danh mục 6 tài liệu tham khảo làm điểm khởi đầu tra cứu.

Giảng viên đại học giảng dạy Giải tích hàm và Tối ưu hóa. Hệ thống ký hiệu chuẩn hóa đầy đủ (bảng ký hiệu liệt kê hơn 20 ký hiệu chuyên ngành từ đạo hàm Dini đến nón tiếp tuyến cấp 2) giúp tiết kiệm thời gian soạn giáo án. Các phản ví dụ như f(x) = x^(1/3) hay f(x, y) = 2|x| − 3|y| là công cụ sư phạm hiệu quả để minh họa giới hạn của lý thuyết Clarke.

Nhà nghiên cứu về lý thuyết tối ưu vectơ và tối ưu đa mục tiêu. Chương 3 với điều kiện tối ưu cấp hai cho cả 2 dạng bài toán ràng buộc và không ràng buộc là phần có giá trị tham khảo cao nhất. Nhóm này có thể sử dụng trực tiếp các định lý về nón tiếp tuyến cấp 2 làm cơ sở phát triển kết quả mới.

Kỹ sư và chuyên viên phân tích trong các ngành có bài toán tối ưu không trơn. Trong thực tế, nhiều bài toán thiết kế kỹ thuật và phân bổ nguồn lực có hàm mục tiêu chứa giá trị tuyệt đối hoặc hàm max — chính là lớp hàm mà luận văn xử lý qua định lý phép lấy maximum.

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

Jacobian xấp xỉ khác gì gradient suy rộng Clarke? Khác biệt cốt lõi là điều kiện ràng buộc. Gradient Clarke luôn là tập lồi, compact và chỉ tồn tại cho hàm Lipschitz địa phương. Jacobian xấp xỉ chỉ cần là tập đóng, có thể không lồi và không bị chặn. Ví dụ f(x) = −|x| có ∂f(0) = {1, −1} — rõ ràng không lồi. Nhờ đó Jacobian xấp xỉ áp dụng được cho toàn bộ lớp hàm liên tục, rộng hơn hẳn.

Jacobian xấp xỉ của một hàm tại một điểm có duy nhất không? Không duy nhất. Luận văn chứng minh nếu ∂f(x₀) là Jacobian xấp xỉ thì mọi tập đóng chứa nó cũng là Jacobian xấp xỉ. Do đó tập J_{x₀} chứa vô hạn phần tử. Tính duy nhất chỉ đạt được với Jacobian xấp xỉ tối thiểu chính quy: khi ∂f(x₀) compact và chính quy, Ext(co(∂f(x₀))) là phần tử tối thiểu duy nhất.

Tại sao Jacobian xấp xỉ "càng nhỏ càng tốt"? Vì điều kiện tối ưu có dạng 0 ∈ co(∂f(x₀)). Tập càng nhỏ thì điều kiện càng chặt, số điểm dừng giả càng ít, việc sàng lọc nghiệm càng hiệu quả. Trong ví dụ f(x, y) = 2|x| − 3|y|, dùng tập 2 phần tử thay vì bao lồi 4 điểm giúp loại bỏ đáng kể các ứng viên không phải nghiệm.

Định lý giá trị trung bình cho Jacobian xấp xỉ phát biểu thế nào? Với f liên tục hữu hạn trên [a, b], tồn tại c ∈ (a, b) và dãy {xₖ} ⊂ co(∂f(c)) sao cho f(b) − f(a) = lim⟨xₖ, b − a⟩. Khi ∂f(c) compact, dạng dãy rút gọn thành f(b) − f(a) ∈ co(⟨∂f(c), b − a⟩). Khi f khả vi, kết quả trở về đúng định lý Lagrange cổ điển.

Có thể áp dụng Jacobian xấp xỉ cho hàm không liên tục không? Không. Toàn bộ lý thuyết trong luận văn đòi hỏi giả thiết liên tục, vì định nghĩa dựa trên đạo hàm Dini theo hướng cần giới hạn tỷ sai phân xác định. Với hàm gián đoạn, cần các công cụ khác như dưới vi phân Mordukhovich hoặc lý thuyết hàm nửa liên tục dưới.

Kết luận

  • Luận văn hệ thống hóa thành công lý thuyết Jacobian xấp xỉ qua 3 chương và 80 trang, đi từ nền tảng hàm khả vi đến ứng dụng trong tối ưu vectơ.
  • Chứng minh được cả 3 lớp dưới vi phân cổ điển — hàm lồi, Clarke, Michel - Penot — đều là trường hợp riêng của Jacobian xấp xỉ, khẳng định vị thế tổng quát hóa của khái niệm.
  • Xác lập đầy đủ 5 phép toán cơ bản cùng định lý giá trị trung bình, đưa lý thuyết đạt mức hoàn thiện tương đương bộ quy tắc đạo hàm cổ điển.
  • Đặc trưng hóa tính Lipschitz địa phương bằng điều kiện cần và đủ về tính bị chặn địa phương của ánh xạ Jacobian xấp xỉ.
  • Thiết lập điều kiện tối ưu cấp hai cho bài toán tối ưu vectơ ở cả 2 dạng có và không có ràng buộc, thông qua Hessian xấp xỉ và nón tiếp tuyến cấp 2.

Lộ trình tiếp theo nên tập trung vào mở rộng không gian vô hạn chiều trong 12-18 tháng và xây dựng thuật toán tính toán trong 6-12 tháng. Độc giả quan tâm đến giải tích không trơn và tối ưu hóa hiện đại nên tham khảo toàn văn luận văn để nắm trọn vẹn hệ thống chứng minh chi tiết.