Tổng quan về luận án

Nghiên cứu về lý thuyết tối ưu và giải tích phi tuyến hiện đại đóng vai trò nền tảng trong việc giải quyết các bài toán cân bằng kinh tế, tối ưu hóa mạng lưới kỹ thuật và xử lý dữ liệu phức tạp. Luận án tiến sĩ toán học của nghiên cứu sinh Nguyễn Thị Thanh Huyền với đề tài "Một số phương pháp giải bài toán chấp nhận tách suy rộng liên quan đến bài toán cân bằng", thực hiện dưới sự hướng dẫn khoa học của GS. Lê Dũng Mưu tại Trường Đại học Sư phạm – Đại học Thái Nguyên, đại diện cho một bước tiến quan trọng trong việc tích hợp cấu trúc của bài toán cân bằng Ky Fan vào bài toán chấp nhận tách (Split Feasibility Problem - SFP).

graph TD
    A["Bài toán chấp nhận tách suy rộng (GSFP)"] --> B["Tập C: Nghiệm bài toán cân bằng EP(K, f)"]
    A --> C["Tập Q: Nghiệm bài toán tối ưu lồi g(u)"]
    B --> D["Toán tử chuyển A: Tuyến tính / Tựa tuyến tính"]
    C --> D
    D --> E["Thuật toán chiếu 1 lần + Phép lặp Mann-Krasnoselskii"]
    D --> F["Thuật toán dưới đạo hàm Clarke cho toán tử phi tuyến"]
    E --> G["Cân bằng Nash-Cournot sản xuất điện & Phí môi trường"]
    F --> H["Mô hình cân bằng Nash có ràng buộc chung"]

Khoảng trống nghiên cứu then chốt (Research Gap) xuất phát từ hạn chế của các phương pháp cổ điển: phần lớn các thuật toán giải bài toán chấp nhận tách truyền thống (như thuật toán $CQ$ của Byrne năm 2002) đòi hỏi tập ràng buộc $C$ và $Q$ phải được cho tường minh dưới dạng hình học đơn giản hoặc là tập điểm bất động của ánh xạ không giãn. Khi $C$ và $Q$ là tập nghiệm ẩn của bài toán cân bằng hoặc bài toán tối ưu lồi, việc áp dụng các phép chiếu mêtric trực tiếp trở nên bất khả thi về mặt tính toán. Hơn nữa, hầu hết các công trình quốc tế trước đây đều giới hạn toán tử chuyển $A: \mathbb{R}^n \to \mathbb{R}^m$ ở dạng tuyến tính; các mở rộng phi tuyến của Li và cộng sự đòi hỏi tính lồi chặt chẽ của hàm khoảng cách hợp thành $p(x) = \frac{1}{2}|P_Q(F(x)) - F(x)|^2$, một giả thiết hiếm khi thỏa mãn trong thực tế.

Luận án đặt ra và giải quyết hai câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Làm thế nào để thiết lập một thuật toán lặp có bước chiếu đơn giản hóa nhằm giải bài toán chấp nhận tách (SEO) khi tập nguồn $C$ là tập nghiệm của bài toán cân bằng para-đơn điệu và tập đích $Q$ là tập nghiệm tối ưu lồi mà không đòi hỏi giải hai bài toán tối ưu phụ ở mỗi bước lặp?
  2. RQ2: Cơ chế toán tử dưới đạo hàm nào cho phép giải bài toán chấp nhận tách phi tuyến (NSEP) khi toán tử chuyển $A$ là phép biến đổi tựa tuyến tính không khả vi?

Khung lý thuyết của công trình được định vị trên sự giao thoa của Lý thuyết bất đẳng thức Ky Fan (1972), Nguyên lý toán tử gần kề Moreau–Rockafellar, Lý thuyết điểm bất động Mann–Krasnoselskii (1953, 1955) và Giải tích biến phân phi trơn của Clarke (1983). Đóng góp mang tính đột phá của luận án là chứng minh sự hội tụ mạnh của dãy lặp dưới điều kiện song hàm para-đơn điệu thông qua phép chiếu duy nhất một lần, loại bỏ hoàn toàn yêu cầu ma trận chuyển vị $A^*$ hoặc ma trận nghịch đảo $A^{-1}$ trong bài toán phi tuyến. Phạm vi nghiên cứu bao quát không gian Euclid hữu hạn chiều $\mathbb{R}^n$ và $\mathbb{R}^m$, đồng thời kiểm chứng bằng mô hình cân bằng bán độc quyền Nash–Cournot trong thị trường điện năng gắn với ràng buộc tối thiểu hóa chi phí phát thải môi trường.

Literature Review và Positioning

Lý thuyết bài toán cân bằng bắt nguồn từ nghiên cứu của Isoda (1955) về trò chơi không hợp tác và được hoàn thiện bởi Ky Fan (1972) dưới dạng bất đẳng thức minimax: "Cho $C$ là một tập lồi, đóng, khác rỗng trong không gian $\mathbb{R}^n$ và $f : C \times C \to \mathbb{R}$ là một song hàm thỏa mãn $f(x, x) = 0, \forall x \in C$ (song hàm có tính chất này thường được gọi là song hàm cân bằng)." Thuật ngữ "Bài toán cân bằng" (Equilibrium Problem - EP) được chính thức định hình bởi Blum và Oettli (1994), tổng quát hóa hàng loạt mô hình như bất đẳng thức biến phân, bài toán tối ưu, bài toán điểm yên ngựa, điểm bất động Kakutani và bài toán bù phi tuyến (Bigi et al., 2019).

graph LR
    A["Ky Fan (1972) / Blum & Oettli (1994): Bài toán cân bằng EP"] --> D["Bài toán chấp nhận tách suy rộng (GSFP)"]
    B["Censor & Elfving (1994) / Byrne (2002): Split Feasibility Problem"] --> D
    C["Moreau (1965) / Combettes (2004): Toán tử gần kề Proximal"] --> D
    D --> E["Luận án Nguyễn Thị Thanh Huyền (2020)"]
    E --> F["Thuật toán chiếu 1 lần Mann-Krasnoselskii (Chương 2)"]
    E --> G["Thuật toán dưới đạo hàm Clarke phi tuyến (Chương 3)"]

Về phương diện bài toán chấp nhận tách, Censor và Elfving (1994) lần đầu tiên giới thiệu SFP cho mô hình bài toán ngược: Tìm $x^* \in C$ sao cho $Ax^* \in Q$. Byrne (2002) phát triển thuật toán $CQ$ dựa trên phép chiếu mêtric $x_{k+1} = P_C(x_k - \gamma A^*(I - P_Q)Ax_k)$, ứng dụng thành công trong phục hồi ảnh y tế và xạ trị ung thư điều biến cường độ (IMRT). Tuy nhiên, một cuộc tranh luận học thuật lớn diễn ra xung quanh chi phí tính toán của các phương pháp giải EP và SFP:

  • Trường phái đạo hàm tăng cường (Extragradient Method): Bắt nguồn từ Antipin (1997) và Langenberg (2012), phương pháp này đảm bảo sự hội tụ cho song hàm đơn điệu nhưng đòi hỏi giải hai bài toán tối ưu phụ (hoặc hai phép chiếu lên tập ràng buộc $C$) tại mỗi bước lặp, gây bùng nổ chi phí tính toán khi số chiều $n$ lớn.
  • Trường phái hàm đánh giá và bài toán phụ: Cohen (1980, 1988), Fukushima (1992) và Quoc et al. (2012) sử dụng hàm khoảng cách Bregman và D-gap functions để chuyển đổi EP thành bài toán tối ưu không ràng buộc, song lại vấp phải bài toán đặt không chỉnh (ill-posed) khi nghiệm không duy nhất.

Luận án của Nguyễn Thị Thanh Huyền định vị chính xác tại điểm giao thoa chưa được giải quyết giữa công trình của Scheimberg, Santos và Santos (2010) về phương pháp chiếu một lần cho bài toán cân bằng và công trình của Thakur (2013) về chuyển đổi SFP thành bài toán tối ưu lồi $\min_{x \in C} \frac{1}{2\lambda}|(I - P_Q)(Ax)|^2$. So với hai nghiên cứu quốc tế tiêu biểu:

  1. Nghiên cứu của Byrne (2002) yêu cầu tập $C, Q$ phải tính được phép chiếu tường minh; luận án mở rộng cho $C = \operatorname{Sol}(EP(K, f))$ và $Q = \arg\min g$ thông qua toán tử gần kề $\operatorname{prox}_{\lambda g}$.
  2. Nghiên cứu của Li et al. (2014) về SFP phi tuyến đòi hỏi hàm dư sai $p(x)$ phải lồi; luận án vượt qua rào cản này bằng cách áp dụng công cụ dưới đạo hàm Clarke cho các toán tử tựa tuyến tính, tạo ra hướng tiếp cận hoàn toàn mới cho bài toán cân bằng Nash có ràng buộc chung.

Đó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 đáng kể lý thuyết đơn điệu hóa trong giải tích lồi bằng việc chứng minh tính khả quy và sự hội tụ của bài toán (SEO): $$\text{Tìm } x^* \in K : f(x^, y) \ge 0, \forall y \in K \quad \text{và} \quad g(Ax^) \le g(u), \forall u \in \mathbb{R}^m$$

Mô hình lý thuyết được thiết lập qua hệ thống mệnh đề và định lý nền tảng:

  • Mệnh đề cấu trúc hàm mục tiêu phụ trợ: Đặt $h(x) := \frac{1}{2}|(I - \operatorname{prox}{\lambda g})Ax|^2$, với $\operatorname{prox}{\lambda g}(u) := \arg\min {g(v) + \frac{1}{2\lambda}|v - u|^2 : v \in \mathbb{R}^m}$, hàm $h(x)$ luôn khả vi và có đạo hàm $\nabla h(x) = A^*(I - \operatorname{prox}_{\lambda g})Ax$, ngay cả khi hàm lồi $g$ không khả vi.
  • Định lý hội tụ Fejér suy rộng: Dãy lặp ${x_k}$ sinh bởi thuật toán thỏa mãn bất đẳng thức co khoảng cách đối với mọi nghiệm $z \in S$: $$|x_{k+1} - z|^2 \le |x_k - z|^2 - (1 - a_k)\rho_k(4 - \rho_k)\frac{h^2(y_k)}{|\nabla h(y_k)|^2} + 2(1 - a_k)\alpha_k f(x_k, z) + A_k$$ trong đó $A_k = 2(1 - a_k)(\alpha_k \varepsilon_k + \beta_k^2)$ với chuỗi $\sum_{k=1}^\infty A_k < +\infty$.
  • Bổ đề tiệm cận triệt tiêu: Chứng minh thành công $\lim_{k\to\infty} h(y_k) = 0$ và $\lim\sup_{k\to\infty} f(x_k, z) = 0$, từ đó suy ra mọi điểm tụ $x^$ đều thỏa mãn đồng thời $x^ \in \operatorname{Sol}(EP)$ và $Ax^* \in \arg\min g$.
flowchart TD
    subgraph IterationStep["Quy trình lặp tại bước k"]
        Xk["Điểm lặp hiện tại: x_k ∈ K"] --> Gk["Lấy dưới đạo hàm: g_k ∈ ∂₂^(ε_k) f(x_k, x_k)"]
        Gk --> AlphaK["Tính bước nhảy: α_k = β_k / max{δ_k, ||g_k||}"]
        AlphaK --> Yk["Phép chiếu mêtric đơn: y_k = P_K(x_k - α_k g_k)"]
        Yk --> Zk["Hiệu chỉnh Proximal: z_k = y_k - µ_k ∇h(y_k)"]
        Zk --> Xk1["Phép lặp Mann-Krasnoselskii: x_(k+1) = a_k x_k + (1 - a_k) z_k"]
    end
    Xk1 --> ConvCheck{"Kiểm tra điều kiện dừng: ||x_(k+1) - x_k|| ≤ ε"}
    ConvCheck -- Chưa đạt --> Xk
    ConvCheck -- Đạt --> Sol["Hội tụ mạnh về nghiệm x* ∈ S = C ∩ A⁻¹ Q"]

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

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết:

  1. Lý thuyết song hàm para-đơn điệu (Paramonotone bifunctions): Một song hàm $f$ được gọi là para-đơn điệu đối với tập nghiệm $S$ nếu với mọi $x \in S, y \in K$ thỏa mãn $f(x, y) = f(y, x) = 0$ thì $y \in S$. Tính chất này đóng vai trò điều kiện biên quan trọng, yếu hơn tính đơn điệu mạnh nhưng đủ để đảm bảo tính đóng và tính đơn trị tiệm cận của quỹ đạo nghiệm.
  2. Ánh xạ gần kề Moreau (Proximal Operator): Thay vì chiếu trực tiếp lên tập nghiệm $Q$, thuật toán sử dụng toán tử $\operatorname{prox}_{\lambda g}$, biến đổi bài toán tối ưu lồi không trơn thành bài toán tìm điểm bất động khả vi của hàm $h(x)$.
  3. Quy tắc bước nhảy tự điều chỉnh (Self-adaptive Stepsize): Việc xác định $\alpha_k = \frac{\beta_k}{\max{\delta_k, |g_k|}}$ với $\sum \beta_k^2 < \infty$ và $\sum \frac{\beta_k}{\delta_k} = \infty$ giải phóng thuật toán khỏi sự phụ thuộc vào hằng số Lipschitz của gradient, vốn rất khó xác định trong thực tế.

Điều kiện biên (Boundary Conditions) được xác định rõ: thuật toán đòi hỏi tập nghiệm $S = C \cap A^{-1}Q \neq \emptyset$. Nếu $S = \emptyset$, dãy lặp sẽ không bị chặn và phân kỳ ra vô hạn, như được chứng minh qua phản ví dụ mẫu với $K = {(u, v) \in \mathbb{R}^2 : u \ge 1, v \ge 1/u}$ và $Q = {(u, v) \in \mathbb{R}^2 : u \ge 1, v = 0}$.

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

Thiết kế nghiên cứu

Luận án quán triệt lập trường triết học nhận thức luận Diễn dịch - Thực chứng (Deductive Positivism) kết hợp với Kiến tạo Toán học (Mathematical Constructivism). Thiết kế nghiên cứu đa tầng (Multi-level Design) chuyển hóa từ không gian giải tích trừu tượng $\mathbb{R}^n \times \mathbb{R}^m$ xuống cấu trúc thuật toán số và mô hình kinh tế kỹ thuật vi mô.

Cấp độ phân tích Đối tượng nghiên cứu Công cụ toán học & Phương pháp Tiêu chuẩn kiểm định
Cấp độ 1: Lý thuyết giải tích Cấu trúc không gian, tập lồi, song hàm para-đơn điệu Giải tích lồi, phép chiếu mêtric $P_K$, nón pháp tuyến $N_K$ Bổ đề Fejér, Bổ đề Opial, dãy Cauchy trong $\mathbb{R}^n$
Cấp độ 2: Thiết kế thuật toán Sơ đồ lặp chiếu kết hợp Mann–Krasnoselskii & Clarke Dưới vi phân xấp xỉ $\partial_2^\varepsilon f$, toán tử gần kề $\operatorname{prox}_{\lambda g}$ Chuỗi hội tụ $\sum \beta_k^2 < \infty$, $\sum \frac{\beta_k}{\delta_k} = \infty$
Cấp độ 3: Mô phỏng số Thị trường điện Nash–Cournot & Cân bằng ràng buộc chung Lập trình số trị, ma trận phát thải, hàm chi phí bậc hai Sai số nghiệm $|x_{k+1} - x_k| \le 10^{-6}$, thời gian CPU

Quy trình nghiên cứu rigorous

Quy trình chứng minh toán học tuân thủ tiêu chuẩn chặt chẽ thông qua tam giác đạc phương pháp luận (Methodological Triangulation):

  • Tính hợp thức cấu trúc (Construct Validity): Sử dụng $\varepsilon$-dưới vi phân $\partial_2^{\varepsilon_k} f(x_k, x_k)$ nhằm khắc phục hiện tượng kỳ dị khi gradient không tồn tại, cho phép vector dưới đạo hàm $g_k$ thỏa mãn bất đẳng thức: $$f(x_k, z) + \varepsilon_k \ge \langle g_k, z - x_k \rangle, \quad \forall z \in K$$
  • Tính hợp thức nội tại (Internal Validity): Thiết lập dãy Cauchy cho dãy hình chiếu ${P_S(x_{k_j})}$: $$|P_S(x_{k_m}) - P_S(x_{k_j})|^2 \le 2|x_{k_j} - P_S(x_{k_j})|^2 + 2\sum_{i=k_j}^{k_m - 1} A_i - 2|x_{k_m} - P_S(x_{k_m})|^2 \xrightarrow{j, m \to \infty} 0$$ chứng minh tính hội tụ duy nhất về điểm $x^* \in S$.
  • Độ tin cậy thuật toán (Algorithmic Reliability): Dãy tham số nội suy ${a_k}$ và ${\rho_k}$ bị chặn chặt chẽ: $0 < a \le a_k \le b < 1$ và $0 < \xi \le \rho_k \le 4 - \xi$, bảo đảm bán kính co của phép lặp Mann luôn nằm trong khoảng ổn định tuyệt đối.

Data và phân tích

Trong phần thử nghiệm số trị, thuật toán được hiện thực hóa bằng ngôn ngữ khoa học trên nền tảng MATLAB. Dữ liệu mô phỏng bài toán cân bằng bán độc quyền Nash–Cournot bao gồm $n$ nhà máy điện, mỗi nhà máy $i$ sở hữu $I_i$ tổ máy phát ($N = \sum |I_i|$ tổng số máy phát). Hàm giá thị trường được mô hình hóa dưới dạng hàm affin nghịch đảo theo tổng sản lượng $s = \sum_{j=1}^N x_j$: $$p_i(s) = \alpha - \beta_i s, \quad (\alpha > 0, \beta_i > 0)$$ Hàm chi phí sản xuất $c_j(x_j)$ của tổ máy $j$ là hàm lồi thực sự: $c_j(x_j) = \frac{1}{2} d_j x_j^2 + e_j x_j + f_j$.

Ma trận công nghệ phát thải $A = (a_{l,j})_{m \times N}$ quy định mức tiêu hao $m$ loại tài nguyên và phát thải môi trường. Hàm chi phí môi trường tổng thể $g(Ax) = \frac{1}{2}(Ax)^T D (Ax) + q^T (Ax)$ với $D$ là ma trận nửa xác định dương. Các kiểm tra độ vững (Robustness Checks) được thực hiện với các kích thước ma trận biến thiên từ số chiều nhỏ ($N=10, m=3$) đến quy mô lớn ($N=1000, m=50$), khẳng định tốc độ hội tụ siêu tuyến tính và tính ổn định vượt trội so với các thuật toán gradient hai bước truyền thống.

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

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

  1. Hội tụ mạnh với một phép chiếu duy nhất: Luận án chứng minh rằng không cần đến hai phép chiếu như phương pháp Antipin, thuật toán kết hợp chiếu Mann–Krasnoselskii vẫn đạt được sự hội tụ mạnh tuyệt đối: $$\lim_{k\to+\infty} x_k = \lim_{k\to+\infty} y_k = \lim_{k\to+\infty} P_S(x_k) = x^*$$
  2. Cơ chế tiệm cận triệt tiêu của hàm dư sai Proximal: Xác lập tốc độ suy giảm của chuỗi khoảng cách $\sum_{k=1}^\infty \frac{h^2(y_k)}{|\nabla h(y_k)|^2} < +\infty$, bảo đảm rằng điều kiện chấp nhận tách $Ax^* \in \arg\min g$ được thỏa mãn chính xác tuyệt đối ở giới hạn vô hạn.
  3. Phát hiện về động thái phân kỳ khi tập nghiệm rỗng: Khi $K \cap A^{-1}Q = \emptyset$, nghiệm xấp xỉ thỏa mãn bất đẳng thức tăng trưởng biên: $$u_{k+1} \ge u_k + \frac{1}{32 [u_k]^3} \implies \lim_{k\to+\infty} u_k = +\infty$$ Đây là bằng chứng toán học quan trọng khẳng định tính nhạy cảm của thuật toán trong việc nhận biết tính khả thi của bài toán.
  4. Khả năng giải quyết bài toán phi tuyến tựa tuyến tính: Ứng dụng thành công dưới đạo hàm Clarke $\partial^C \phi(x)$ cho các toán tử chuyển tựa tuyến tính dạng phân thức $A(x) = \frac{a^T x + b}{c^T x + d}$, mở rộng biên giới giải tích phi tuyến sang các bài toán không lồi cục bộ.

Implications đa chiều

graph TD
    A["Phát hiện cốt lõi của luận án"] --> B["Ý nghĩa Lý thuyết"]
    A --> C["Ý nghĩa Phương pháp"]
    A --> D["Ý nghĩa Thực tiễn & Chính sách"]
    
    B --> B1["Mở rộng lý thuyết Ky Fan & Moreau-Rockafellar"]
    B --> B2["Chuẩn hóa điều kiện para-đơn điệu trong GSFP"]
    
    C --> C1["Giảm 50% số phép chiếu mêtric tại mỗi bước lặp"]
    C --> C2["Thuật toán tự thích nghi không cần hằng số Lipschitz"]
    
    D --> D1["Tối ưu hóa thị trường điện lực cạnh tranh"]
    D --> D2["Cân bằng giữa tối đa hóa lợi nhuận & Giảm phát thải CO₂"]
  • Ý nghĩa lý thuyết: Hợp nhất lý thuyết bài toán cân bằng đơn điệu và bài toán chấp nhận tách thành một khung giải tích thống nhất, tạo cơ sở toán học vững chắc cho lớp bài toán tối ưu hai cấp (Bilevel Optimization).
  • Ý nghĩa phương pháp luận: Cung cấp một sơ đồ thuật toán tiết kiệm tài nguyên tính toán (giảm 50% số phép chiếu mêtric trên mỗi bước lặp so với phương pháp Extragradient), dễ dàng cài đặt trong các hệ thống tính toán hiệu năng cao (HPC).
  • Ý nghĩa thực tiễn và chính sách năng lượng: Mô hình cung cấp công cụ định lượng chính xác cho các cơ quan điều tiết thị trường điện lực nhằm xác định mức thuế môi trường tối ưu, buộc các tập đoàn phát điện độc quyền nhóm (Oligopoly) phải tự động điều chỉnh sản lượng về mức cân bằng sinh thái mà không làm suy giảm tính ổn định của lưới điện quốc gia.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nội tại:

  1. Không gian hữu hạn chiều: Toàn bộ các định lý hội tụ mạnh được thiết lập trong không gian Euclid $\mathbb{R}^n$. Việc mở rộng sang không gian Hilbert vô hạn chiều $L^2(\Omega)$ hoặc không gian Banach lồi đều đòi hỏi phải bổ sung các điều kiện về tính compact yếu.
  2. Rào cản para-đơn điệu: Thuật toán chưa bao quát được các song hàm không đơn điệu hoặc hoàn toàn không lồi (Non-convex/Non-monotone Equilibrium Problems).
  3. Độ nhạy tham số: Dù thuật toán có tính tự điều chỉnh, việc lựa chọn dãy hạt nhân xấp xỉ ${\varepsilon_k}$ và ${\beta_k}$ ban đầu vẫn ảnh hưởng đáng kể đến tốc độ hội tụ trong các bài toán quy mô siêu lớn.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra các hướng phát triển:

  • Hướng 1: Mở rộng thuật toán chiếu Mann–Krasnoselskii sang không gian Hilbert vô hạn chiều ứng dụng trong phương trình đạo hàm riêng (PDE-constrained optimization).
  • Hướng 2: Phát triển các biến thể ngẫu nhiên (Stochastic Split Equilibrium Algorithms) nhằm xử lý dữ liệu lớn chứa nhiễu đo lường trong thời gian thực.
  • Hướng 3: Thiết kế thuật toán phân tán (Decentralized & Parallel Algorithms) cho mạng lưới các tác tử độc lập (Multi-agent Systems) trong điều khiển lưới điện thông minh (Smart Grid).

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

Các kết quả nghiên cứu trong luận án đã được công bố trên các tạp chí quốc tế uy tín hàng đầu thuộc danh mục ISI/Scopus (Q1/Q2):

  • Mathematical Methods of Operations Research (Springer)
  • Journal of Global Optimization (Springer)

Công trình đã được báo cáo và nhận được sự đánh giá cao tại các diễn đàn khoa học chuyên sâu:

  • Hội nghị Toán ứng dụng và Tin học toàn quốc (Trường Đại học Bách Khoa Hà Nội, 2016).
  • The 7th International Conference on High Performance Scientific Computing (Hà Nội, 2018).
  • Hệ thống Xêmina liên cơ quan: Viện Toán học – Viện Hàn lâm KH&CN Việt Nam, Đại học Khoa học Tự nhiên và Đại học Thăng Long.

Về mặt ứng dụng công nghiệp, phương pháp chấp nhận tách của luận án cung cấp giải pháp toán học cốt lõi cho các thuật toán tối ưu hóa xạ trị ung thư điều biến cường độ chùm tia (IMRT) và tái tạo ảnh chụp cắt lớp y tế (CT/MRI), giúp giảm thiểu thời gian xử lý ảnh y khoa và tối ưu hóa liều lượng bức xạ bảo vệ mô lành.

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

mindmap
  root((Hệ sinh thái Hưởng lợi))
    Nghiên cứu sinh & Nhà toán học
      Giải tích lồi & Tối ưu hóa
      Lý thuyết điểm bất động
      Bất đẳng thức Ky Fan
    Kỹ sư R&D & Hệ thống Năng lượng
      Điều độ kinh tế hệ thống điện
      Tối ưu hóa mạng lưới Smart Grid
      Mô hình hóa chuỗi cung ứng
    Chuyên gia Y sinh & Xử lý Tín hiệu
      Tái tạo hình ảnh CT / MRI
      Tối ưu hóa liều lượng xạ trị IMRT
      Xử lý tín hiệu nén
    Nhà hoạch định Chính sách Môi trường
      Cơ chế định giá hạn ngạch Carbon
      Thiết kế thị trường điện cạnh tranh
      Kiểm soát phát thải công nghiệp
  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung chứng minh chặt chẽ về tính hội tụ mạnh của phương pháp chiếu đơn kết hợp toán tử gần kề, giải quyết triệt để các bài toán cân bằng suy rộng.
  • Kỹ sư R&D và Điều độ Hệ thống Điện: Ứng dụng trực tiếp thuật toán để giải bài toán phân bổ công suất phát kinh tế (Economic Dispatch) có xét đến chi phí xử lý khí thải và hạn ngạch carbon.
  • Nhà hoạch định Chính sách Năng lượng & Môi trường: Sử dụng mô hình mô phỏng cân bằng Nash–Cournot để dự báo hành vi giá và sản lượng của các doanh nghiệp phát điện khi áp dụng các sắc thuế môi trường mới.

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 tích hợp thành công toán tử gần kề Moreau $\operatorname{prox}_{\lambda g}$ vào thuật toán chiếu một lần của Scheimberg (2010), giải quyết trọn vẹn bài toán chấp nhận tách (SEO) cho song hàm para-đơn điệu mà không cần giả thiết đơn điệu mạnh. Công trình đã mở rộng trực tiếp Lý thuyết bất đẳng thức Ky Fan (1972) và Lý thuyết điểm bất động Mann–Krasnoselskii sang cấu trúc không gian tách lồi kép.

2. Sự đổi mới về phương pháp luận so với các công trình quốc tế tiền nhiệm?

So với thuật toán $CQ$ của Byrne (2002) và thuật toán đạo hàm tăng cường của Antipin (1997), phương pháp của luận án cắt giảm chính xác 50% số lượng bài toán tối ưu phụ cần giải tại mỗi vòng lặp. Đồng thời, cơ chế bước nhảy tự thích nghi $\alpha_k = \beta_k / \max{\delta_k, |g_k|}$ loại bỏ hoàn toàn yêu cầu tính toán trước hằng số Lipschitz hay chuẩn ma trận $|A|$, tạo ưu thế vượt trội về hiệu năng tính toán.

3. Phát hiện thực nghiệm/toán học nào gây bất ngờ nhất trong nghiên cứu?

Phát hiện toán học ấn tượng nhất là đặc tính phân kỳ biên $u_{k+1} \ge u_k + \frac{1}{32[u_k]^3}$ khi $S = \emptyset$. Trong khi đa số các thuật toán khác rơi vào trạng thái lặp vô hạn hoặc dao động cục bộ không thể kiểm soát khi không có nghiệm, thuật toán đề xuất tạo ra quỹ đạo chuẩn tiến thẳng ra vô hạn, đóng vai trò như một tiêu chuẩn tự động phát hiện tính vô nghiệm của bài toán.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp đầy đủ và minh bạch 100% các điều kiện tham số: dãy số thực $0 < a < a_k < b < 1$, chuỗi độ lệch $\sum \varepsilon_k < \infty$, $\sum \beta_k^2 < \infty$, $\sum \beta_k / \delta_k = \infty$, cùng cấu trúc ma trận chi phí $c_j(x_j)$ và ma trận phát thải $A$, cho phép tái lập chính xác toàn bộ kết quả mô phỏng trên bất kỳ nền tảng tính toán số nào.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình phát triển tập trung vào ba trục chính: (1) Số hóa và song song hóa thuật toán trên kiến trúc GPU giải quyết các bài toán tối ưu hàng triệu chiều; (2) Tích hợp học máy tăng cường (Reinforcement Learning) để tự động tối ưu hóa việc chọn tham số bước nhảy; (3) Mở rộng giải quyết các bài toán cân bằng Nash tổng quát ngẫu nhiên trong mạng lưới giao thông thông minh và kinh tế số.

Kết luận

  1. Xây dựng thành công thuật toán chiếu đơn kết hợp phép lặp Mann–Krasnoselskii giải bài toán chấp nhận tách liên quan đến bài toán cân bằng para-đơn điệu và tối ưu lồi (SEO).
  2. Chứng minh toán học nghiêm ngặt sự hội tụ mạnh của dãy lặp về nghiệm chính xác của bài toán mà không đòi hỏi tính đơn điệu mạnh hay khả vi của hàm mục tiêu.
  3. Phát triển thuật toán dưới đạo hàm Clarke xử lý thành công bài toán chấp nhận tách phi tuyến (NSEP) với toán tử chuyển tựa tuyến tính, vượt qua rào cản tính lồi của các phương pháp truyền thống.
  4. Xây dựng mô hình toán kinh tế hoàn chỉnh cho bài toán cân bằng bán độc quyền Nash–Cournot sản xuất điện gắn với tối thiểu hóa phát thải môi trường.
  5. Cung cấp phản ví dụ mẫu và chứng minh giải tích chuẩn xác về động thái phân kỳ của thuật toán khi tập nghiệm khả thi rỗng.
  6. Mở ra 3 nhánh nghiên cứu liên ngành mới: tối ưu hóa phân tán trên lưới điện thông minh, giải thuật toán chấp nhận tách ngẫu nhiên và giải tích biến phân phi lồi trong xử lý dữ liệu lớn. Công trình xác lập vị thế học thuật vững chắc với 2 bài báo quốc tế thuộc danh mục ISI, tạo nền tảng lý thuyết giá trị cho chuyên ngành Toán giải tích và Tối ưu hóa tại Việt Nam.