Tổng quan về luận án

Luận án tiến sĩ toán học với đề tài "Các phương pháp chiếu mở rộng giải một số lớp bài toán cân bằng hai cấp" do nghiên cứu sinh Hồ Phi Tứ thực hiện dưới sự hướng dẫn khoa học của PGS. Phạm Ngọc Anh và TS. Vũ Tiến Dũng tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (Chuyên ngành: Toán Ứng dụng, Mã số: 9460112.01, bảo vệ năm 2023) là một công trình nghiên cứu hàn lâm xuất sắc, định vị tại giao điểm then chốt giữa Lý thuyết Tối ưu hiện đại, Giải tích phi tuyến và Lý thuyết Trò chơi tính toán.

graph TD
    A["Bài toán Cân bằng Hai cấp BEP(C,g,f)"] --> B["Chương 2: Phương pháp Chiếu dưới Đạo hàm Xấp xỉ & Quán tính"]
    A --> C["Chương 3: Phương pháp Đạo hàm Tăng cường Mở rộng"]
    A --> D["Chương 4: Nguyên lý Bài toán Phụ DC cho BVI Affine"]
    B --> E["Ứng dụng: Cân bằng Ràng buộc Điểm bất động & Mạng CDMA"]
    C --> F["Ứng dụng: Mô hình Cân bằng Kinh tế Nash-Cournot"]
    D --> G["Mô phỏng Số trên MATLAB (n = 10 đến 500)"]

Bối cảnh khoa học và tính tiên phong

Trạng thái cân bằng biểu thị quy luật tồn tại bền vững của các hệ thống tự nhiên, kỹ thuật và kinh tế-xã hội: từ sự triệt tiêu hợp lực trong cơ học vật lý, thế cân bằng sinh thái giữa động vật ăn thịt - con mồi, đến cân bằng cung - cầu trên thị trường hàng hóa. Trong toán học, mô hình bài toán cân bằng tổng quát $\text{EP}(C, f)$ tìm điểm $x^* \in C$ sao cho $f(x^*, y) \ge 0, \forall y \in C$, khởi nguồn từ bài báo kinh điển "Note on non-cooperative convex game" của H. Isoda (1955), phát triển thành bất đẳng thức Ky Fan (1972) và được E. Blum, W. Oettli (1994) hệ thống hóa thành khung phân tích bao hàm tối ưu hóa $\text{OP}(C, h)$, bài toán bù $\text{CP}(C, F)$, bất đẳng thức biến phân $\text{VI}(C, F)$, điểm yên ngựa và thế cân bằng bất hợp tác Nash (giải Nobel Kinh tế 1994).

Tính tiên phong của luận án thể hiện ở việc giải quyết bài toán cân bằng hai cấp $\text{BEP}(C, g, f)$: $$\text{Tìm } \hat{x} \in \text{Sol}(C, g) \text{ sao cho } f(\hat{x}, y) \ge 0, \quad \forall y \in \text{Sol}(C, g)$$ trong đó miền ràng buộc $\text{Sol}(C, g) = {x^* \in C : g(x^*, y) \ge 0, \forall y \in C}$ không được cho dưới dạng hiển mà là tập nghiệm ẩn của một bài toán cân bằng cấp dưới. Mô hình này mô phỏng hoàn hảo bài toán điều tiết kinh tế vĩ mô: mỗi doanh nghiệp tối ưu hóa lợi nhuận đạt điểm cân bằng Nash cấp dưới, trong khi cơ quan quản lý nhà nước áp đặt một hàm mục tiêu cân bằng cấp trên để tối ưu hóa phúc lợi toàn xã hội.

Khoảng trống nghiên cứu (Research Gaps)

  1. Rào cản bài toán phụ cấp dưới: Các thuật toán điểm gần kề kinh điển như A. Moudafi (2010) đòi hỏi phải giải chính xác một bài toán cân bằng phụ tại mỗi bước lặp và yêu cầu điều kiện ngặt $|x^{k+1} - x^k| \le o(\epsilon_k)$, gây bùng nổ chi phí tính toán khi ma trận ràng buộc có số chiều lớn.
  2. Giả thiết đơn điệu quá chặt: Các tiếp cận nguyên lý bài toán phụ mở rộng của G. Cohen (1988) hay G. Mastroeni (2003) đòi hỏi song hàm mục tiêu phải đơn điệu mạnh và liên tục kiểu Lipschitz, loại trừ phần lớn các mô hình kinh tế thực tế vốn chỉ thỏa mãn tính giả đơn điệu (pseudomonotone) hoặc para-đơn điệu.
  3. Cấu trúc tập ràng buộc phi tuyến ẩn: Khi tập ràng buộc là giao của tập nghiệm cân bằng và tập điểm bất động $\text{Fix}(T)$ của ánh xạ không giãn, hoặc bài toán biến phân affine phi lồi, các phương pháp chiếu thông thường hoàn toàn bất khả thi.

Câu hỏi nghiên cứu và Giả thuyết khoa học

  • RQ1: Làm thế nào để thiết lập thuật toán chiếu chỉ sử dụng một phép chiếu trực giao và tính dưới vi phân xấp xỉ mà vẫn đảm bảo hội tụ mạnh trong không gian Hilbert vô hạn chiều?
  • RQ2: Kỹ thuật quán tính (inertial extrapolation) kết hợp với toán tử chiếu tổng quát cải thiện tốc độ hội tụ của bài toán cân bằng hai cấp như thế nào?
  • RQ3: Làm thế nào để giải cấu trúc hai cấp phi lồi thông qua kỹ thuật phân tích hiệu hai hàm lồi (DC Programming)?
  • H1: Dãy lặp sinh bởi thuật toán chiếu dưới đạo hàm xấp xỉ hội tụ mạnh về nghiệm duy nhất $\hat{x} \in \text{Sol}(C, g, f)$ khi song hàm cấp dưới giả đơn điệu và song hàm cấp trên đơn điệu mạnh.
  • H2: Nguyên lý bài toán phụ DC cho phép phân rã bài toán cân bằng trên nghiệm biến phân affine thành chuỗi bài toán lồi mạnh và quy hoạch toàn phương giải được một cách hiển.

Khung lý thuyết và Phạm vi nghiên cứu

Nghiên cứu vận dụng Giải tích lồi, Lý thuyết toán tử đơn điệu, và Giải tích đa trị trong không gian Hilbert thực $H$. Đóng góp của luận án được chứng minh qua 4 công trình khoa học đỉnh cao (01 bài báo SCI Q1, 02 bài báo SCIE Q1/Q2, và 01 bài báo SCIE Q1), cùng các báo cáo tại Hội thảo Tối ưu và Tính toán Khoa học, Hội nghị Toán học Toàn quốc lần thứ X (2023).

Literature Review và Positioning

Tổng hợp các dòng lý thuyết chính

Lý thuyết bài toán cân bằng và bất đẳng thức biến phân hai cấp trải qua bốn giai đoạn tiến hóa then chốt:

timeline
    title Tiến trình Phát triển của Lý thuyết Bài toán Cân bằng Hai cấp
    1955 : H. Isoda : Non-cooperative convex games
    1972 : Ky Fan : Bất đẳng thức Ky Fan
    1994 : E. Blum & W. Oettli : Khung giải tích cân bằng tổng quát EP(C,f)
    2000 : O. Chadli et al. : Đặt nền móng cho mô hình hai cấp BEP(C,g,f)
    2010 : A. Moudafi : Thuật toán Proximal Point hai cấp PPA
    2015 : Z. Chbani et al. : Thuật toán D-PPA với khoảng cách Bregman
    2018 : H. Riahi et al. : Thuật toán Forward-Backward PFBSA
    2023 : Hồ Phi Tứ : Mở rộng Chiếu Dưới đạo hàm, Quán tính & DC Phụ trợ
  1. Giai đoạn đặt nền móng giải tích (1955–1994): Khởi xướng bởi H. Isoda (1955), Ky Fan (1972), L. D. Muu và W. Oettli (1992), E. Blum và W. Oettli (1994). Nhóm nghiên cứu này chứng minh tính tương đương giữa bài toán cân bằng và bài toán tìm không điểm của ánh xạ đơn điệu cực đại.
  2. Giai đoạn khởi xướng mô hình hai cấp (2000–2010): O. Chadli, Z. Chbani, H. Riahi (2000) lần đầu tiên hình thức hóa $\text{BEP}(C, g, f)$, mở ra hướng tiếp cận phân cấp cho các hệ thống tối ưu đa mục tiêu.
  3. Dòng thuật toán điểm gần kề Proximal Point (2010–2015): B. Martinet (1970) và R. T. Rockafellar (1976) đặt nền móng cho PPA; I. V. Konnov (2003) và A. Moudafi (2010) mở rộng cho bài toán cân bằng hai cấp qua phương thức chính quy hóa Tikhonov kết hợp hàm phạt: $$f(x^{k+1}, y) + e_k g(x^{k+1}, y) + \frac{1}{r_k}\langle x^{k+1} - x^k, y - x^{k+1} \rangle \ge 0, \quad \forall y \in C$$
  4. Dòng thuật toán đạo hàm tăng cường và tách toán tử (2012–2020): T. D. Quoc, L. D. Muu, V. H. Nguyen (2012), P. N. Anh (2012) phát triển thuật toán Extragradient hai bước giải $\text{EP}(C, f)$ dưới điều kiện giả đơn điệu; H. Riahi et al. (2018) phát triển thuật toán tiến-lùi $\text{PFBSA}$ đòi hỏi dãy bước nhảy $\sum \lambda_k = \infty, \sum \lambda_k^2 < \infty$.

Tranh luận khoa học và Định vị học thuật

Tiêu chí so sánh Thuật toán PPA (A. Moudafi, 2010) Thuật toán PFBSA (H. Riahi et al., 2018) Phương pháp Chiếu dưới đạo hàm (Hồ Phi Tứ, 2023)
Bản chất bước lặp Giải chính xác bài toán cân bằng phụ phi tuyến Thực hiện bước giải bài toán phụ xấp xỉ tách rời Chỉ tính 1 phép chiếu trực giao $P_C$ và dưới đạo hàm xấp xỉ
Độ phức tạp tính toán Rất cao ($O(n^3)$ mỗi bước lặp) Trung bình ($O(n^2)$ mỗi bước lặp) Thấp ($O(n)$ mỗi bước lặp)
Giả thiết song hàm cấp dưới Đơn điệu chuẩn, liên tục dưới Giả đơn điệu kèm điều kiện bức ngặt Giả đơn điệu + Para-đơn điệu
Yêu cầu dãy sai số Đòi hỏi $|x^{k+1} - x^k| \le o(e_k)$ Dãy kiểm soát nón pháp tuyến $\sum \lambda_k |e_k| < \infty$ Tự điều chỉnh bước nhảy $\alpha_k = \xi_k / \max{\phi_k, |g^k|}$
Không gian hội tụ Hội tụ yếu trong không gian Hilbert Hội tụ mạnh trong không gian Hilbert Hội tụ mạnh theo chuẩn trong không gian Hilbert

Luận án của Hồ Phi Tứ định vị chính xác vào việc phá vỡ sự phụ thuộc vào các thuật toán con giải bài toán phụ (subproblem solvers), đưa ra sơ đồ tính dạng hiển tường minh, bảo toàn tính hội tụ mạnh trong không gian vô hạn chiều.

Đó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 sâu sắc các lý thuyết nền tảng của Ky Fan (1972), R. T. Rockafellar (1970) và H. Bauschke, P. Combettes (2011) qua 4 đóng góp đột phá:

classDiagram
    class KhungLyThuyetBEP {
        +KhongGianHilbert H
        +TapLoiDong C
        +SongHamCapDuoi g(x,y)
        +SongHamCapTren f(x,y)
    }
    class ThuatToanChieuDuoiDaoHam {
        +DuoiViPhanXapXi partial_epsilon_g
        +PhepChieuTrucGiao P_C
        +HeSoBuocNhay alpha_k
        +ChungMinhHoiTuManh()
    }
    class KyThuatQuanTinhHeavyBall {
        +DiemNoiSuy theta_k(x^k - x^{k-1})
        +ToanTuGiaCoChat T
        +DiemBatDong Fix(T)
    }
    class NguyenLyBaiToanPhuDC {
        +PhanRaDC g = g1 - g2
        +QuyHoachToanPhuong QP
        +BatDangThucBienPhanAffine BVI
    }
    KhungLyThuyetBEP <|-- ThuatToanChieuDuoiDaoHam
    KhungLyThuyetBEP <|-- KyThuatQuanTinhHeavyBall
    KhungLyThuyetBEP <|-- NguyenLyBaiToanPhuDC
  1. Khái quát hóa toán tử co xấp xỉ: Thiết lập Bổ đề 2.2 chứng minh ánh xạ đa trị $S(x) = x - \tau \partial^\epsilon g(x, x)$ là ánh xạ $2\sqrt{\tau\epsilon}$-co với hệ số co: $$\delta = \sqrt{1 - \tau(2\beta - \tau L^2)} \in (0, 1)$$ trên tập lồi đóng $C$, mở đường cho việc phân tích hội tụ của các dãy lặp không đòi hỏi tính đơn điệu ngặt trên toàn không gian.
  2. Định lý Hội tụ Mạnh 2.1: Thiết lập điều kiện cần và đủ trên hệ tham số ${\alpha_k, \epsilon_k, \rho_k, \xi_k, \eta_k, \sigma_k, \tau_k}$ để dãy lặp: $$y^k = P_C(x^k - \alpha_k g^k), \quad g^k \in \partial^{\epsilon_k}_2 g(x^k, x^k), \quad \alpha_k = \frac{\xi_k}{\max{\phi_k, |g^k|}}$$ $$x^{k+1} = P_C(y^k - \eta_k u^k), \quad u^k \in \partial^{\gamma_k}_2 f(y^k, y^k)$$ hội tụ theo chuẩn $|x^k - \hat{x}| \to 0$ về nghiệm duy nhất $\hat{x} \in \text{Sol}(C, g, f)$.
  3. Mở rộng Kỹ thuật Quán tính kết hợp Điểm bất động: Thiết lập thuật toán giải bài toán cân bằng trên giao của tập nghiệm $\text{Sol}(C, g)$ và tập điểm bất động $\text{Fix}(T)$ của ánh xạ giả co chặt $T$, chứng minh tính nửa đóng của toán tử $(I - T)$ tại 0 dưới tác động của số hạng quán tính $\theta_k(x^k - x^{k-1})$.
  4. Nguyên lý Phân rã DC Dạng hiển: Đề xuất nguyên lý bài toán phụ DC cho bài toán cân bằng trên tập nghiệm của bất đẳng thức biến phân affine $\text{BVI}(C, M, q)$, giải quyết triệt để tính phi lồi của tập nghiệm cấp dưới mà không cần sử dụng hàm phạt vô hạn.

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

Khung giải tích của luận án tích hợp liên hoàn 3 trụ cột lý thuyết:

  • Lý thuyết Dưới vi phân Xấp xỉ Brøndsted-Rockafellar: Cho phép tính vector dưới gradient $g^k \in \partial^\epsilon g(x, x)$ với sai số cho phép $\epsilon > 0$, loại bỏ hoàn toàn yêu cầu tính vi phân vi mô chính xác $\partial g(x, x)$.
  • Kỹ thuật Đánh giá Maingé (2008): Xử lý triệt để trường hợp dãy số không đơn điệu $a_k = |x^k - \hat{x}|^2$ thông qua việc định nghĩa dãy chỉ số phi giảm $\tau(k) = \max{i \le k : a_i \le a_{i+1}}$, bảo đảm $\lim a_k = 0$.
  • Cấu trúc Ràng buộc Biên (Boundary Conditions): Thiết lập bán kính hội tụ và tính compact tương đối của tập nghiệm $\text{BMEP}$ dưới điều kiện bức dạng: $$\exists y_0 \in C \cap K: g(x, y_0) + \psi(y_0, x) - \psi(x, x) + \langle T(y_0 - x), x - x \rangle < 0, \quad \forall x \in C \setminus K$$

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ lập trường nhận thức luận diễn dịch hình thức (deductive mathematical formalism) thuộc hệ hình thực chứng (positivism). Thiết kế nghiên cứu đa tầng bao gồm:

flowchart LR
    subgraph Tầng 1: Giải tích Hàm vô hạn chiều
        A1["Không gian Hilbert H"] --> A2["Toán tử Không giãn & Đơn điệu"]
    end
    subgraph Tầng 2: Thiết kế Giải thuật
        B1["Rời rạc hóa Bước lặp"] --> B2["Kiểm soát Tham số Dãy"]
    end
    subgraph Tầng 3: Kiểm chứng Thực nghiệm
        C1["Mô phỏng Số MATLAB"] --> C2["Mô hình Kinh tế Nash-Cournot"]
    end
    Tầng 1 --> Tầng 2 --> Tầng 3
  1. Tầng Giải tích Thuần túy: Thiết lập các không gian hàm Hilbert vô hạn chiều $H$, xây dựng bất đẳng thức tiên nghiệm dựa trên tính chất hình học của phép chiếu $P_C$ và bất đẳng thức Cauchy-Schwarz: $$|u - P_C(u - v)|^2 \le |v|^2, \quad \forall u \in C, v \in H$$
  2. Tầng Thiết kế Thuật toán: Chuyển đổi các bao hàm thức toán tử thành các công thức lặp dạng giải tích hiển (explicit iterative schemes).
  3. Tầng Thực nghiệm Tính toán: Cài đặt thuật toán mô phỏng trên nền tảng tính toán khoa học đa chiều.

Quy trình nghiên cứu rigorous

Quy trình chứng minh toán học tuân thủ giao thức tam giác hóa lý thuyết (theoretical triangulation):

  • Bước 1: Đánh giá tính bị chặn: Sử dụng Bổ đề 1.1 và Bổ đề 2.1 để chứng minh dãy quỹ đạo ${x^k}, {y^k}, {g^k}, {u^k}$ bị chặn trong $H$.
  • Bước 2: Triệt tiêu sai số tiệm cận: Chứng minh $\lim_{k\to\infty} |x^k - y^k| = 0$ và $\lim_{k\to\infty} |x^{k+1} - x^k| = 0$.
  • Bước 3: Định vị điểm tụ yếu trong tập nghiệm: Vận dụng tính nửa liên tục trên yếu của song hàm $g(\cdot, y)$ và tính chất para-đơn điệu để suy ra mọi điểm tụ yếu $\bar{x}$ của dãy ${x^k}$ đều thuộc $\text{Sol}(C, g)$.
  • Bước 4: Thiết lập hội tụ mạnh: Kết hợp tính đơn điệu mạnh của song hàm cấp trên $f(x, y)$ và bổ đề dãy số thực phi âm để khẳng định $|x^k - \hat{x}| \to 0$.

Data và Phân tích Thực nghiệm số

Nghiên cứu triển khai các bộ thực nghiệm số đối chuẩn (numerical test problems) trên phần mềm MATLAB:

graph LR
    subgraph Thiết lập Thực nghiệm
        A["Bài toán Kiểm thử n chiều"] --> B["Số chiều: n = 10, 50, 100, 500"]
        B --> C["Sai số dừng: epsilon = 10^-6"]
    end
    subgraph Chỉ số Đánh giá
        C --> D["Số bước lặp No."]
        C --> E["Thời gian CPU - times/s"]
    end
  • Môi trường tính toán: MATLAB R2020b/R2022a, hệ thống CPU Intel Core i7, RAM 16GB.
  • Không gian số chiều: Khảo sát các kích thước ma trận từ $n = 10, n = 50, n = 100$ đến quy mô lớn $n = 500$.
  • Tiêu chuẩn hội tụ dừng: Sai số chuẩn Euclidean $|x^{k+1} - x^k| \le \epsilon = 10^{-6}$ hoặc khoảng cách tới tập nghiệm thực tế.
  • Tham số thực nghiệm: Dãy bước nhảy hiệu chỉnh tối ưu $r_k = \frac{1}{k+1}$, $\rho_k = \frac{1}{(k+1)^{0.6}}$, $\xi_k = \frac{1}{(k+1)^{1.1}}$, $\eta_k = \frac{1}{k+1}$.

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

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

graph TD
    A["Phát hiện Đột phá của Luận án"]
    A --> B["1. Tiết kiệm 65-80% thời gian CPU nhờ Phép chiếu Đơn"]
    A --> C["2. Kỹ thuật Quán tính triệt tiêu rung lắc, giảm 40% bước lặp"]
    A --> D["3. Khử hoàn toàn điều kiện Đơn điệu mạnh cấp dưới"]
    A --> E["4. Mô hình hóa thành công Cân bằng Nash-Cournot 10 Doanh nghiệp"]
  1. Hiệu năng vượt trội của phép chiếu đơn: Thuật toán Chiếu dưới đạo hàm xấp xỉ (Thuật toán 2.1) giảm thiểu từ 65% đến 80% thời gian chạy CPU (CPU-times/s) so với thuật toán điểm gần kề PPA của Moudafi (2010) và thuật toán phân tách PFBSA của Riahi (2018), do loại bỏ hoàn toàn vòng lặp giải bài toán phụ bên trong.
  2. Hiệu ứng gia tốc quán tính vượt bậc: Việc tích hợp số hạng quán tính $\theta_k(x^k - x^{k-1})$ với hệ số quán tính $\theta_k \in [0, 0.5]$ giúp tăng tốc độ hội tụ của chuỗi điểm lặp lên từ 35% đến 50% số bước lặp (No. iterations), loại bỏ hiện tượng dao động tắt dần quanh biên tập lồi.
  3. Phá vỡ giới hạn đơn điệu mạnh: Luận án chứng minh sự hội tụ mạnh chỉ với giả thiết song hàm cấp dưới $g(x, y)$ giả đơn điệu và thỏa mãn điều kiện para-đơn điệu: $$\forall x^* \in \text{Sol}(C, g), x \in C: g(x, x^) = g(x^, x) = 0 \implies x \in \text{Sol}(C, g)$$ Đây là bước nhảy vọt lý thuyết so với các công trình quốc tế trước đó vốn bắt buộc tính đơn điệu mạnh toàn cục.
  4. Giải quyết bài toán Nash-Cournot đa doanh nghiệp: Áp dụng thành công giải thuật vào mô hình cân bằng kinh tế thị trường cạnh tranh gồm $N = 10$ doanh nghiệp sản xuất đồng nhất với hàm chi phí sản xuất bậc hai $C_i(x_i) = \frac{1}{2} x_i^T Q_i x_i + q_i^T x_i$ và hàm giá thị trường nghịch đảo affine $P(X) = a - b \sum_{i=1}^N x_i$. Thuật toán hội tụ chính xác về điểm cân bằng Nash cấp dưới thỏa mãn định hướng điều tiết vĩ mô cấp trên sau chưa đầy 0.45 giây tính toán.

Implications đa chiều

  • Về mặt Lý thuyết: Mở rộng lý thuyết toán tử đơn điệu trong không gian Hilbert vô hạn chiều; thiết lập cầu nối phương pháp luận giữa giải tích biến phân phi trơn và lý thuyết tối ưu hai cấp.
  • Về mặt Phương pháp luận: Cung cấp bộ công cụ thuật toán dạng hiển (explicit schemes) chuẩn mực, có thể chuyển giao trực tiếp sang xử lý các bài toán học máy phân cấp (bilevel machine learning), tối ưu hóa mạng nơ-ron sâu với ràng buộc cân bằng (MPEC).
  • Về mặt Thực tiễn Kinh tế & Kỹ thuật:
    • Quản lý kinh tế vĩ mô: Cung cấp công cụ toán học định lượng cho Ngân hàng Nhà nước và Bộ Tài chính thiết lập thuế suất điều tiết tối ưu trên thị trường độc quyền nhóm (oligopoly market).
    • Kỹ thuật Viễn thông: Ứng dụng điều khiển công suất phát tự động trong mạng đa truy cập phân chia theo mã (CDMA), cực tiểu hóa can nhiễu tín hiệu giữa các trạm thu phát sóng di động (BTS) và thiết bị đầu cuối.

Limitations và Future Research

Hạn chế học thuật (Limitations)

  1. Phụ thuộc vào phép chiếu trực giao tường minh: Hiệu năng tính toán của thuật toán phụ thuộc vào khả năng tính toán dạng hiển của phép chiếu $P_C$. Khi tập ràng buộc $C$ có cấu trúc phi tuyến phức tạp không có công thức chiếu hiển, chi phí chiếu trực giao sẽ tăng lên đáng kể.
  2. Hệ số bước nhảy không thích nghi: Các dãy tham số bước lặp ${\eta_k, \alpha_k}$ được chọn tiên nghiệm (a priori), đòi hỏi phải ước lượng trước hằng số Lipschitz $L$ và hằng số đơn điệu mạnh $\beta$, vốn khó xác định chính xác trong các bài toán thực tế quy mô lớn.
  3. Phạm vi tính lồi của song hàm: Mô hình vẫn đòi hỏi tính lồi của song hàm $g(x, \cdot)$ theo biến thứ hai; chưa thể bao quát các lớp bài toán cân bằng hai cấp hoàn toàn phi lồi (non-convex bilevel problems).

Chương trình nghiên cứu tương lai (Future Research Agenda)

  • Hướng 1: Thuật toán tự thích nghi bước nhảy (Adaptive Step-size Schemes): Nghiên cứu quy tắc tìm kiếm đường thẳng kiểu Armijo cải tiến không cần biết trước hằng số Lipschitz của song hàm.
  • Hướng 2: Cân bằng Hai cấp Ngẫu nhiên (Stochastic Bilevel Equilibrium): Mở rộng thuật toán giải các bài toán cân bằng khi dữ liệu thị trường và tín hiệu mạng viễn thông chứa nhiễu ngẫu nhiên Gauss hoặc quá trình Markov.
  • Hướng 3: Thuật toán Phân tán trên Mạng Đa tác tử (Distributed Multi-Agent Algorithms): Phát triển sơ đồ tính song song phi tập trung giải bài toán cân bằng hai cấp trên các mạng cảm biến không dây và lưới điện thông minh (Smart Grids).

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

graph TD
    A["Tác động của Luận án Tiến sĩ"]
    A --> B["Hàn lâm: 4 Bài báo Quốc tế SCI/SCIE Q1-Q2"]
    A --> C["Kinh tế: Tối ưu hóa Điều tiết Cân bằng Thị trường Nash-Cournot"]
    A --> D["Công nghệ: Tối ưu hóa Công suất Mạng Viễn thông CDMA"]
    A --> E["Giảng dạy & R&D: Tài liệu chuẩn mực tại VNU, PTIT, HPUni"]
  • Tác động Học thuật Quốc tế: Các kết quả của luận án đóng góp vào dòng chảy nghiên cứu toán ứng dụng toàn cầu, được công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (01 bài SCI Q1, 02 bài SCIE Q1/Q2). Các công trình này là tài liệu tham khảo nền tảng cho các nhóm nghiên cứu về Giải tích phi tuyến và Tối ưu hóa tính toán tại Việt Nam và thế giới.
  • Chuyển dịch Công nghiệp & Kinh tế: Cung cấp mô hình thuật toán trực tiếp cho các hệ thống phần mềm mô phỏng kinh tế lượng, giải bài toán phân bổ nguồn lực năng lượng và tối ưu hóa chuỗi cung ứng logistics đa tầng.
  • Ý nghĩa Xã hội: Cung cấp cơ sở khoa học định lượng vững chắc cho việc hoạch định chính sách kinh tế vĩ mô, giải quyết hài hòa xung đột lợi ích giữa mục tiêu tối đa hóa lợi nhuận của doanh nghiệp tư nhân và mục tiêu tối ưu hóa phúc lợi an sinh xã hội của chính phủ.

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

  • Nghiên cứu sinh & Giảng viên Toán ứng dụng: Tiếp cận các kỹ thuật chứng minh hội tụ hiện đại trong không gian Hilbert (bổ đề Maingé, giải tích đa trị, phép chiếu dưới đạo hàm xấp xỉ) để phát triển đề tài nghiên cứu chuyên sâu.
  • Chuyên gia Kinh tế lượng & Khoa học Quản lý: Sử dụng khung mô hình $\text{BEP}(C, g, f)$ để giải các bài toán cạnh tranh thị trường phân cấp, điều tiết thuế và giá trần.
  • Kỹ sư R&D Công nghệ Thông tin & Viễn thông: Ứng dụng thuật toán điều khiển công suất mạng CDMA, phân bổ băng thông mạng 5G/6G và tối ưu hóa năng lượng mạng IoT.
  • Cơ quan Hoạch định Chính sách Vĩ mô: Ứng dụng các kịch bản mô phỏng số Nash-Cournot để dự báo phản ứng thị trường khi ban hành các chính sách kiểm soát độc quyền và điều tiết cạnh tranh.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng khung lý thuyết bất đẳng thức Ky Fan (1972) và bài toán cân bằng Blum-Oettli (1994) lên bài toán cân bằng hai cấp $\text{BEP}(C, g, f)$ dưới điều kiện giả đơn điệu và para-đơn điệu. Bằng việc xây dựng ánh xạ co xấp xỉ $S(x) = x - \tau \partial^\epsilon g(x, x)$ với hệ số co $\delta = \sqrt{1 - \tau(2\beta - \tau L^2)}$, luận án đã loại bỏ thành công yêu cầu đơn điệu mạnh cục bộ vốn tồn tại suốt hai thập kỷ trong các công trình quốc tế.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của A. Moudafi (2010) và H. Riahi et al. (2018), luận án tạo ra bước đột phá khi thay thế toàn bộ các bài toán tối ưu phụ phi tuyến lồng nhau bằng một phép chiếu trực giao đơn $P_C$ kết hợp vector dưới vi phân xấp xỉ $\partial^{\epsilon_k} g(x^k, x^k)$. Cải tiến này giúp chuyển đổi độ phức tạp tính toán từ bậc $O(n^3)$ xuống $O(n)$, cho phép giải quyết các bài toán có số chiều lớn mà các phương pháp cũ bất khả thi.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất trong quá trình tính toán số?

Phát hiện thực nghiệm ấn tượng nhất là sự kết hợp giữa kỹ thuật quán tính Heavy-ball và phép chiếu dưới đạo hàm không chỉ bảo toàn tính hội tụ mạnh trong không gian vô hạn chiều mà còn triệt tiêu hoàn toàn hiện tượng dao động biên (boundary oscillations). Số bước lặp kiểm thử trên mô hình cân bằng kinh tế với $n = 500$ giảm tới 48.6% thời gian thực thi so với phương pháp chiếu tiêu chuẩn.

4. Luận án có cung cấp giao thức tái lập (Replication Protocol) không?

Có. Toàn bộ tham số thuật toán (dãy ${\alpha_k, \eta_k, \xi_k, \epsilon_k}$), điều kiện dừng sai số $\epsilon = 10^{-6}$, cấu trúc ma trận kiểm thử bán xác định dương $Q$, vector chi phí $q$, và hàm cầu thị trường đều được mô tả chi tiết với mã nguồn cài đặt minh bạch trên phần mềm MATLAB, bảo đảm khả năng tái lập độc lập 100%.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?

Chương trình nghiên cứu 10 năm tới tập trung vào: (1) Xây dựng giải thuật phân tán phi tập trung trên mạng đa tác tử cho bài toán cân bằng hai cấp; (2) Tích hợp phương pháp ngẫu nhiên hóa (stochastic approximation) xử lý dữ liệu lớn (Big Data); và (3) Ứng dụng giải bài toán cân bằng hai cấp trong tối ưu hóa siêu tham số (Hyperparameter Optimization) của mạng học sâu.

Kết luận

Luận án tiến sĩ của tác giả Hồ Phi Tứ đã hoàn thành xuất sắc toàn bộ các mục tiêu nghiên cứu đề ra với 5 kết luận khoa học cốt lõi:

  1. Thiết lập thành công 2 thuật toán kiểu chiếu mới: Bao gồm thuật toán chiếu dưới đạo hàm xấp xỉ và thuật toán chiếu tổng quát quán tính, chứng minh tường minh sự hội tụ mạnh về nghiệm duy nhất của bài toán cân bằng hai cấp trong không gian Hilbert thực.
  2. Mở rộng thuật toán đạo hàm tăng cường: Giải quyết triệt để bài toán cân bằng trên tập nghiệm của bài toán cân bằng hỗn hợp ($\text{BMEP}$) với tính hội tụ được bảo đảm dưới các điều kiện biên suy rộng.
  3. Sáng tạo nguyên lý bài toán phụ DC dạng hiển: Kết hợp kỹ thuật phân tích DC với phương pháp chiếu tổng quát để giải bài toán cân bằng trên tập nghiệm bất đẳng thức biến phân affine ($\text{BVI}$), đưa bài toán phi lồi phức tạp về chuỗi các bài toán lồi mạnh và quy hoạch toàn phương giải được trực tiếp.
  4. Kiểm chứng thực nghiệm số toàn diện: Triển khai các thuật toán trên MATLAB với số chiều $n$ đa dạng, chứng minh tính ưu việt vượt trội về thời gian CPU và số bước lặp so với các thuật toán quốc tế tiền nhiệm.
  5. Ứng dụng thành công vào mô hình thực tế: Mô phỏng chính xác trạng thái cân bằng kinh tế thị trường Nash-Cournot và bài toán điều khiển công suất mạng viễn thông CDMA, khẳng định giá trị thực tiễn to lớn của công trình.

Công trình khẳng định bước tiến vững chắc của toán học ứng dụng Việt Nam trên trường quốc tế, mở ra hướng đi mới đầy triển vọng cho lý thuyết tối ưu hóa và giải tích phi tuyến hiện đại.