Tổng quan luận án

Luận án tiến sĩ toán học "Approximation Problems for Dynamic Equations on Time Scales" (tên tiếng Việt: "Bài toán xấp xỉ cho phương trình động lực trên thang thời gian") do nghiên cứu sinh Nguyễn Thu Hà thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Hữu Dư tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội (bảo vệ năm 2017), thuộc chuyên ngành Phương trình Vi phân và Tích phân (mã số: 62 46 01 03).

  • Tính cấp thiết và khoảng trống nghiên cứu: Lý thuyết giải tích trên thang thời gian (Analysis on Time Scales) do Stefan Hilger đề xuất nhằm thống nhất và mở rộng hai nhánh giải tích: giải tích liên tục ($\mathbb{R}$) và giải tích rời rạc ($\mathbb{Z}$). Trong các thập niên qua, các nghiên cứu tập trung chủ yếu vào việc tổng quát hóa các tính chất định tính (tính ổn định, dao động, bài toán biên) từ phương trình vi phân và phương trình sai phân sang phương trình động lực trên thang thời gian nói chung. Tuy nhiên, một hướng tiếp cận chưa được khảo sát đầy đủ là xem bản thân lý thuyết thang thời gian như một công cụ và bài toán xấp xỉ (approximation problem). Cụ thể, khi xét một dãy các thang thời gian ${\mathbb{T}n}{n=1}^\infty$ hội tụ về thang thời gian $\mathbb{T}$ theo khoảng cách Hausdorff, các câu hỏi đặt ra là: nghiệm của phương trình động lực trên $\mathbb{T}_n$ có hội tụ về nghiệm trên $\mathbb{T}$ hay không; tốc độ hội tụ được đánh giá như thế nào; và các đặc trưng định tính như miền ổn định mũ, phổ ma trận và bán kính ổn định của phương trình động lực ẩn tuyến tính phụ thuộc liên tục ra sao vào dữ liệu phương trình khi cả hệ số ma trận lẫn miền thang thời gian cùng biến thiên.

  • Mục tiêu nghiên cứu:

    1. Thiết lập điều kiện cho sự tồn tại, duy nhất và sự hội tụ của nghiệm đối với phương trình động lực Delta ($\Delta$) và phương trình động lực Nabla ($\nabla$) trên dãy thang thời gian biến thiên; đồng thời đánh giá tốc độ hội tụ của nghiệm thông qua khoảng cách Hausdorff giữa các thang thời gian.
    2. Nghiên cứu sự phụ thuộc dữ liệu (data-dependence) của các đặc trưng định tính trong phương trình động lực ẩn tuyến tính (Linear Implicit Dynamic Equations - LIDEs), bao gồm: sự phụ thuộc liên tục của miền ổn định mũ đều theo thang thời gian, tính liên tục của phổ chùm ma trận và tính nửa liên tục dưới của bán kính ổn định dưới tác động của nhiễu cấu trúc ở cả hai vế phương trình.

  • Đối tượng và phạm vi nghiên cứu:

    • Đối tượng: Phương trình động lực tường minh cấp một $x^\Delta(t) = f(t, x(t))$ và $x^\nabla(t) = f(t, x(t))$; phương trình động lực ẩn tuyến tính hệ số hằng $A x^\Delta(t) = B x(t)$ và $A_n x^{\Delta_n}(t) = B_n x(t)$ với $A, B, A_n, B_n \in \mathbb{C}^{m \times m}$ ($A, A_n$ có thể suy biến); chùm ma trận ${A, B}$ với chỉ số bất kỳ; bán kính ổn định phức dưới nhiễu cấu trúc hai vế $[A + D\Sigma E_1, B + D\Sigma E_2]$.
    • Phạm vi không gian và thời gian: Thang thời gian $\mathbb{T}$ là tập con đóng, khác rỗng bất kỳ của $\mathbb{R}$ có độ hạt bị chặn (bounded graininess) chứa điểm ban đầu $t_0$; không gian trạng thái hữu hạn chiều $\mathbb{R}^m$ hoặc $\mathbb{C}^m$; khoảng thời gian xét nghiệm là khoảng compact $[t_0, T]\mathbb{T}$ hoặc nửa đường thẳng thực $[t_0, \infty)\mathbb{T}$.

Tổng quan tài liệu và vị trí của luận án

Lý thuyết giải tích trên thang thời gian được khởi xướng bởi Stefan Hilger trong luận án tiến sĩ năm 1988 (dưới sự hướng dẫn của Bernd Aulbach) và công bố năm 1990, đặt nền móng cho việc mô hình hóa đồng thời các hiện tượng liên tục và rời rạc. Về sau, hệ thống lý thuyết giải tích và phương trình động lực trên thang thời gian được hệ thống hóa và phát triển mạnh mẽ qua các chuyên khảo của Bohner và Peterson (2001, 2003), Agarwal, O'Regan và Lakshmikantham (2002), Hilscher và Zeidan (2004), Kaymakcalan (1996).

Trong phân tích số trị cho phương trình vi phân thường, các phương pháp sai phân kinh điển như Euler tường minh (Euler explicit/forward) và Euler ẩn (Euler implicit/backward) đã được nghiên cứu sâu sắc bởi Hairer, Wanner (1996), Butcher (2008), Deuflhard (2012), Ascher và Petzold (1998). Về bản chất, việc rời rạc hóa phương trình vi phân trên lưới điểm chính là việc chuyển bài toán sang một thang thời gian rời rạc $\mathbb{T}n = {t_0^{(n)}, t_1^{(n)}, \dots, t{k_n}^{(n)}}$.

Đối với hệ vi phân - đại số (DAEs) và phương trình động lực ẩn tuyến tính (LIDEs), cấu trúc nghiệm và tính giải được phụ thuộc chặt chẽ vào chỉ số (index) của chùm ma trận ${A, B}$, được khảo sát bởi Campbell (1980, 1982), Griepentrog và März (1986), Kunkel và Mehrmann (2006), Dai (1989), Gantmacher (1959).

Về lý thuyết ổn định và bán kính ổn định:

  • Khái niệm bán kính ổn định cho hệ phương trình vi phân bất biến theo thời gian được Hinrichsen và Pritchard đưa ra năm 1986.
  • Son và Hinrichsen (1989) mở rộng khái niệm bán kính ổn định cho phương trình sai phân.
  • Pötzsche, Siegmund và Wirth (2003) nghiên cứu lý thuyết phổ và tính ổn định mũ đều trên thang thời gian.
  • Doan, Kalauch và Siegmund (2011) thiết lập công thức bán kính ổn định cho phương trình động lực tuyến tính tường minh $x^\Delta = Bx$.
  • Du và Loi (2012), Du và Thuan (2014) khảo sát bán kính ổn định cho phương trình vi phân - đại số và phương trình động lực ẩn tuyến tính.
  • Du, Lien, Linh (2003, 2004) nghiên cứu sự phụ thuộc liên tục của bán kính ổn định đối với hệ vi phân đại số khi tham số nhiễu tiến về 0.

Vị trí của luận án: Luận án giải quyết khoảng trống học thuật bằng cách kết nối lý thuyết giải tích thang thời gian với bài toán xấp xỉ giải tích thông qua metric Hausdorff. Luận án tổng quát hóa sự hội tụ của các lược đồ giải số (như Euler) thành định lý hội tụ nghiệm của phương trình vi phân trên dãy thang thời gian tùy ý; đồng thời mở rộng bài toán phụ thuộc dữ liệu của bán kính ổn định cho hệ ẩn từ trường hợp suy biến tham số nhỏ sang trường hợp biến thiên đồng thời cả ma trận hệ số và cấu trúc hình học của thang thời gian.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và các khái niệm sử dụng

  • Thang thời gian và các toán tử bước nhảy: Thang thời gian $\mathbb{T}$ là tập con đóng, khác rỗng của $\mathbb{R}$. Toán tử nhảy tiến $\sigma(t) := \inf{s \in \mathbb{T}: s > t}$, toán tử nhảy lùi $\rho(t) := \sup{s \in \mathbb{T}: s < t}$. Độ hạt tiến $\mu(t) := \sigma(t) - t$ và độ hạt lùi $\nu(t) := t - \rho(t)$. Phân loại điểm: điểm phân tán phải ($\sigma(t) > t$), điểm dày phải ($\sigma(t) = t$), điểm phân tán trái ($\rho(t) < t$), điểm dày trái ($\rho(t) = t$).
  • Không gian hàm: Hàm liên tục dày phải (rd-continuous, ký hiệu $C_{rd}$), hàm liên tục dày trái (ld-continuous, ký hiệu $C_{ld}$) và hàm chính quy (regulated function).
  • Đạo hàm Delta và Nabla: Đạo hàm $\Delta$ định nghĩa bởi $f^\Delta(t) = \lim_{s \to t, \sigma(t)} \frac{f(\sigma(t)) - f(s)}{\sigma(t) - s}$; đạo hàm $\nabla$ định nghĩa qua bước lùi $\rho(t)$. Tại điểm phân tán phải, $f^\Delta(t) = \frac{f(\sigma(t)) - f(t)}{\mu(t)}$.
  • Độ đo và tích phân Lebesgue trên thang thời gian: Độ đo Carathéodory $\Delta$-measure $m_\Delta$ và $\nabla$-measure $m_\nabla$. Với $a, b \in \mathbb{T}$ ($a \le b$), $m_\Delta([a, b)) = b - a$, $m_\Delta({t_0}) = \sigma(t_0) - t_0$; $m_\nabla((a, b]) = b - a$, $m_\nabla({t_0}) = t_0 - \rho(t_0)$.
  • Biến đổi trụ (Cylinder transformation) và hàm mũ: Với $h > 0$, phép biến đổi trụ $\xi_h: \mathbb{C}h \to \mathbb{Z}h$ định nghĩa bởi $\xi_h(z) = \frac{\text{Ln}(1 + hz)}{h}$ (với $\text{Ln}$ là nhánh logarit chính có miền giá trị $[-\pi i, \pi i)$), và $\xi_0(z) = z$. Hàm mũ suy rộng trên $\mathbb{T}$ với $p \in \mathcal{R}$ (nhóm thụt lùi): $$e_p(t, t_0) := \exp\left(\int{t_0}^t \xi{\mu(s)}(p(s))\Delta s\right)$$
  • Đa thức suy rộng: Dãy đa thức ${h_k(t, \tau)}{k=0}^\infty$ xác định bởi $h_0(t, \tau) = 1$ và $h_k(t, \tau) = \int\tau^t h_{k-1}(s, \tau)\Delta s$. Ta có bất đẳng thức $h_k(t, t_0) \le \frac{(t - t_0)^k}{k!}$ với mọi $t \ge t_0$.
  • Khoảng cách Hausdorff giữa các thang thời gian: Với $\mathbb{T}_1, \mathbb{T}_2 \in \mathcal{T}(t_0)$ (tập các thang thời gian có độ hạt bị chặn chứa $t_0$): $$d_H(\mathbb{T}_1, \mathbb{T}2) := \max\left{\sup{t_1 \in \mathbb{T}_1} d(t_1, \mathbb{T}2), \sup{t_2 \in \mathbb{T}_2} d(t_2, \mathbb{T}_1)\right}$$

Phương pháp nghiên cứu thực tế

Luận án sử dụng thuần túy phương pháp nghiên cứu toán học lý thuyết và giải tích định lượng:

  1. Phương pháp xấp xỉ liên tiếp (Picard iteration) và giải tích hàm: Thiết lập toán tử tích phân tương đương trên thang thời gian, áp dụng nguyên lý ánh xạ co Banach và bổ đề Gronwall-Bellman trên thang thời gian để chứng minh sự tồn tại, duy nhất và tốc độ hội tụ nghiệm.
  2. Phương pháp tô-pô không gian metric: Sử dụng metric Hausdorff để định lượng mức độ sai khác giữa các tập rời rạc/liên tục và dãy thang thời gian.
  3. Đại số tuyến tính và lý thuyết phổ chùm ma trận: Sử dụng dạng chuẩn tắc Weierstrass-Kronecker, phép chiếu phổ, chỉ số chùm ma trận $\text{Ind}(A, B)$ để giải tách hệ động lực ẩn thành phần vi phân và phần đại số.
  4. Phương pháp phân tích hàm phức và tối ưu hóa giải tích: Xác định các miền ổn định đều $\mathcal{U}_\mathbb{T}$ và cực trị chuẩn toán tử truyền đạt trên biên phổ để đánh giá bán kính ổn định $r(A, B; D, E; \mathbb{T})$.

Nội dung chính theo từng chương

Chapter 1: Preliminary

Chương 1 hệ thống hóa các khái niệm, bổ đề và công cụ cơ bản của giải tích trên thang thời gian và lý thuyết ổn định mũ.

  • Toán tử và độ đo: Thiết lập chi tiết cấu trúc tập hợp các điểm của thang thời gian $\mathbb{T}$ (xem Bảng 1), tính chất của các hàm $C_{rd}$, $C_{ld}$, nhóm hàm thụt lùi $\mathcal{R}(\mathbb{T}, \mathbb{K})$, nhóm thụt lùi dương $\mathcal{R}^+(\mathbb{T}, \mathbb{R})$. Định nghĩa độ đo Lebesgue $\Delta$ và $\nabla$, chỉ ra mối liên hệ giữa tích phân $\Delta$ trên thang thời gian với tích phân Lebesgue thông thường: $$\int_a^b f(t)\Delta t = \int_{[a, b]\mathbb{T}} f(t)\text{mes}(dt) + \sum{a \le t < b} f(t)(\sigma(t) - t)$$

Bảng 1: So sánh các đặc trưng giải tích trên các thang thời gian điển hình

Thang thời gian $\mathbb{T}$ Toán tử bước nhảy $\sigma(t)$ Độ hạt $\mu(t)$ Đạo hàm Delta $f^\Delta(t)$ Tích phân $\int_a^b f(t)\Delta t$
$\mathbb{R}$ (Liên tục) $t$ $0$ $f'(t)$ (Đạo hàm cổ điển) $\int_a^b f(t)dt$ (Tích phân Riemann/Lebesgue)
$\mathbb{Z}$ (Rời rạc nguyên) $t + 1$ $1$ $\Delta f(t) = f(t+1) - f(t)$ $\sum_{t=a}^{b-1} f(t)$ (Tổng sai phân)
$h\mathbb{Z}$ ($h > 0$) $t + h$ $h$ $\frac{f(t+h)-f(t)}{h}$ $\sum_{k=a/h}^{b/h - 1} f(kh)h$
$P_{a, b} = \bigcup_{k=0}^\infty [k(a+b), k(a+b)+a]$ $t$ nếu $t \in [k(a+b), k(a+b)+a)$
$t+b$ nếu $t = k(a+b)+a$		 $0$ hoặc $b$ Hỗn hợp vi phân và sai phân Tích phân từng đoạn cộng tổng bước nhảy
  • Toán tử Cauchy và ổn định mũ: Xét hệ phương trình vi phân thuần nhất $x^\Delta = A(t)x$. Nghiệm duy nhất của $X^\Delta = A(t)X, X(t_0) = I$ là toán tử Cauchy $\Phi_A(t, t_0)$. Hệ được gọi là ổn định mũ đều khi và chỉ khi tồn tại $N \ge 1, \alpha > 0$ với $-\alpha \in \mathcal{R}^+$ sao cho $|\Phi_A(t, t_0)| \le N e_{-\alpha}(t, t_0)$ với mọi $t \ge t_0$.
  • Miền ổn định mũ: Đối với hệ hệ số hằng $x^\Delta = Ax$, miền ổn định mũ được định nghĩa: $$\mathcal{S}(\mathbb{T}) := \left{\lambda \in \mathbb{C} : \limsup_{T \to \infty} \frac{1}{T - t_0} \int_{t_0}^T \lim_{s \downarrow \mu(t)} \frac{\ln|1 + s\lambda|}{s}\Delta t < 0\right}$$ Miền ổn định mũ đều là $\mathcal{U}\mathbb{T} := {\lambda \in \mathbb{C} : x^\Delta = \lambda x \text{ ổn định mũ đều}}$. Định lý Pötzsche-Siegmund-Wirth chứng minh rằng hệ $x^\Delta = Ax$ ổn định mũ đều khi và chỉ khi phổ $\sigma(A) \subset \mathcal{U}\mathbb{T}$.

Chapter 2: On the convergence of solutions for dynamic equations on time scales

Chương 2 trình bày hướng tiếp cận coi thang thời gian là bài toán xấp xỉ, chứng minh tính hội tụ của nghiệm phương trình động lực khi dãy thang thời gian hội tụ.

  • Sự tồn tại và duy nhất nghiệm: Xét bài toán giá trị ban đầu (IVP): $$\begin{cases} x^\Delta(t) = f(t, x(t)), \quad t \in \mathbb{T} \ x(t_0) = x_0 \end{cases}$$ Với $f: \mathbb{T} \times \mathbb{R}^m \to \mathbb{R}^m$ liên tục $rd$ và thỏa mãn điều kiện Lipschitz toàn cục theo $x$ với hằng số $k$ ($|f(t, x_1) - f(t, x_2)| \le k|x_1 - x_2|$). Định lý 2.3 chứng minh bài toán luôn tồn tại duy nhất nghiệm $x(t)$ trên $[t_0, T]\mathbb{T}$ thông qua dãy xấp xỉ Picard $x_0(t) \equiv x_0$, $x{n+1}(t) = x_0 + \int_{t_0}^t f(s, x_n(s))\Delta s$, với đánh giá sai số: $$\sup_{t_0 \le s \le T} |x_{n+1}(s) - x_n(s)| \le C k^n h_n(T, t_0) \le C \frac{k^n (T - t_0)^n}{n!}$$
  • Sự hội tụ của nghiệm cho phương trình $\Delta$: Giả sử dãy thang thời gian ${\mathbb{T}n}{n=1}^\infty \subset \mathcal{T}(t_0)$ hội tụ về $\mathbb{T}$ theo khoảng cách Hausdorff ($d_H(\mathbb{T}_n, \mathbb{T}) \to 0$ khi $n \to \infty$). Định lý 2.7 khẳng định dãy nghiệm $x_n(t)$ của phương trình $x_n^\Delta(t) = f(t, x_n(t)), x_n(t_0) = x_0$ trên $\mathbb{T}_n$ hội tụ từng điểm về nghiệm $x(t)$ trên $\mathbb{T}$.
  • Đánh giá tốc độ hội tụ: Khi hàm $f(t, x)$ thỏa mãn điều kiện Lipschitz theo cả hai biến $t$ và $x$, tồn tại hằng số $C > 0$ sao cho sai số giữa hai nghiệm bị chặn bởi khoảng cách Hausdorff: $$|x_n(t) - x(t)| \le C d_H(\mathbb{T}, \mathbb{T}_n), \quad \forall t \in \mathbb{T} \cap \mathbb{T}_n : t_0 \le t \le T$$
  • Mở rộng cho phương trình Nabla và phương pháp Euler: Các kết quả tương tự được thiết lập đầy đủ cho phương trình vi phân lùi $x^\nabla(t) = f(t, x(t))$. Luận án chứng minh rằng phương pháp Euler tường minh và Euler ẩn cổ điển chỉ là các trường hợp riêng của định lý hội tụ này khi chọn dãy thang thời gian rời rạc $\mathbb{T}n = {t_0^{(n)} < t_1^{(n)} < \dots < t{k_n}^{(n)}}$ có bước lưới tiến về 0.

Chapter 3: On data-dependence of implicit dynamic equations on time scales

Chương 3 tập trung vào tính chất định tính của phương trình động lực ẩn tuyến tính (LIDEs): $$A_n x^{\Delta_n}(t) = B_n x(t), \quad t \in \mathbb{T}_n, \quad n \in \mathbb{N}$$ với các ma trận hệ số $A_n, B_n \in \mathbb{C}^{m \times m}$ và ma trận dẫn đầu $A_n$ có thể suy biến ($\det A_n = 0$).

  • Sự phụ thuộc liên tục của miền ổn định mũ đều: Định lý 3.7 thiết lập mối liên hệ topo giữa các miền ổn định khi $\mathbb{T}n \to \mathbb{T}$ trong tô-pô Hausdorff: $$\mathcal{U}\mathbb{T} \subset \bigcup_{n=1}^\infty \bigcap_{m \ge n} \mathcal{U}{\mathbb{T}m} \quad \text{và} \quad \bigcup{n=1}^\infty \bigcap{m \ge n} \mathcal{U}_{\mathbb{T}m} \setminus \mathbb{R} \subset \mathcal{U}\mathbb{T} \setminus \mathbb{R}$$
  • Sự phụ thuộc dữ liệu của phổ chùm ma trận: Phổ $\sigma(A, B) := {\lambda \in \mathbb{C} : \det(\lambda A - B) = 0}$. Khi các hệ số biến thiên $(A_n, B_n) \to (A, B)$, phổ có thể biến đổi gián đoạn nếu chỉ số chùm ma trận thay đổi. Luận án đưa ra các điều kiện đại số:
    • Trường hợp $\text{Ind}(A, B) = 1$: Mệnh đề 3.11 chỉ ra nếu $\lim_{n \to \infty} (A_n, B_n) = (A, B)$ và $(A_n - A)Q = 0$ với mọi $n \in \mathbb{N}$ (trong đó $Q$ là toán tử chiếu tương ứng), thì $\lim_{n \to \infty} \sigma(A_n, B_n) = \sigma(A, B)$ theo khoảng cách Hausdorff.
    • Trường hợp $\text{Ind}(A, B) > 1$: Mệnh đề 3.13 chứng minh nếu $(A_n - A)\widehat{Q} = (B_n - B)\widehat{Q} = 0$ với mọi $n \in \mathbb{N}$, thì tính liên tục của phổ $\lim_{n \to \infty} \sigma(A_n, B_n) = \sigma(A, B)$ tiếp tục được bảo toàn.
  • Bán kính ổn định dưới nhiễu cấu trúc hai vế: Xét hệ chịu nhiễu: $$[\widetilde{A}, \widetilde{B}] = [A, B] + D\Sigma E = [A + D\Sigma E_1, B + D\Sigma E_2]$$ trong đó $D \in \mathbb{C}^{m \times l}, E = [E_1, E_2] \in \mathbb{C}^{q \times 2m}$, và khối nhiễu không xác định $\Sigma \in \mathbb{C}^{l \times q}$. Bán kính ổn định phức được xác định bởi: $$r(A, B; D, E; \mathbb{T}) = \left(\sup_{\lambda \in \mathcal{U}_\mathbb{T}^c} |G(\lambda)|\right)^{-1}$$ với hàm truyền $G(\lambda) = (E_1 \lambda - E_2)(\lambda A - B)^{-1} D$.
  • Tính nửa liên tục dưới của bán kính ổn định: Định lý 3.18 chứng minh rằng nếu dãy hệ số và thang thời gian cùng hội tụ $\lim_{n \to \infty} (A_n, B_n, D_n, E_n; \mathbb{T}n) = (A, B, D, E; \mathbb{T})$ kèm theo các điều kiện triệt tiêu trên toán tử chiếu, thì bán kính ổn định đạt tính nửa liên tục dưới: $$r(A, B; D, E; \mathbb{T}) = \liminf{n \to \infty} r(A_n, B_n, D_n, E_n; \mathbb{T}_n)$$

Kết quả và những đóng góp mới

Các đóng góp mới về mặt khoa học của luận án bao gồm:

  1. Xây dựng khung lý thuyết xấp xỉ nghiệm qua khoảng cách Hausdorff: Đề xuất cách tiếp cận đưa cấu trúc tô-pô metric Hausdorff lên họ các thang thời gian $\mathcal{T}(t_0)$, chứng minh định lý hội tụ nghiệm tổng quát cho phương trình động lực Delta và Nabla khi miền thời gian biến thiên.
  2. Đánh giá định lượng tốc độ hội tụ nghiệm: Thiết lập đánh giá sai số nghiệm cấp 1 theo khoảng cách Hausdorff $|x_n(t) - x(t)| \le C d_H(\mathbb{T}, \mathbb{T}_n)$ dưới điều kiện Lipschitz hai biến, giải thích bản chất giải tích của các thuật toán số trị như Euler tường minh và Euler ẩn.
  3. Đặc trưng hóa quan hệ giữa các miền ổn định đều: Thiết lập các bao hàm thức topo giữa miền ổn định mũ đều $\mathcal{U}\mathbb{T}$ và giới hạn của dãy miền ổn định ${\mathcal{U}{\mathbb{T}_n}}$ khi thang thời gian biến thiên.
  4. Xác định điều kiện liên tục của phổ chùm ma trận suy biến: Tìm ra các điều kiện triệt tiêu đại số trên ma trận chiếu $((A_n - A)Q = 0$ hoặc $(A_n - A)\widehat{Q} = (B_n - B)\widehat{Q} = 0)$ bảo đảm phổ $\sigma(A_n, B_n)$ hội tụ về $\sigma(A, B)$ đối với chùm ma trận có chỉ số $\text{Ind}(A, B) = 1$ và $\text{Ind}(A, B) > 1$.
  5. Chứng minh tính nửa liên tục dưới của bán kính ổn định hai vế: Mở rộng bài toán phân tích độ bền vững ổn định cho hệ động lực ẩn tuyến tính chịu nhiễu cấu trúc đồng thời ở cả hai vế khi cả ma trận hệ số và miền thang thời gian cùng biến thiên.

Các kết quả này đã được công bố trên 3 bài báo khoa học quốc tế uy tín:

  • Qualitative Theory of Dynamical Systems (2016)
  • Dynamics of Continuous, Discrete and Impulsive Systems Series A: Mathematical Analysis (2015)
  • Mathematical Control and Signals Systems (2016)

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế được đề cập:

    • Luận án chỉ khảo sát không gian trạng thái hữu hạn chiều $\mathbb{R}^m, \mathbb{C}^m$ và các thang thời gian có độ hạt bị chặn (bounded graininess).
    • Phương trình động lực ẩn tuyến tính mới chỉ xét trường hợp hệ số hằng $A, B$, chưa bao quát trường hợp ma trận hệ số biến thiên theo thời gian $A(t)x^\Delta(t) = B(t)x(t)$.
    • Điều kiện liên tục của phổ và bán kính ổn định đòi hỏi các giả thiết đại số tương đối chặt chẽ trên các ma trận chiếu ($(A_n - A)Q = 0$).

  • Hướng nghiên cứu tiếp theo:

    • Mở rộng bài toán xấp xỉ và khảo sát sự phụ thuộc dữ liệu cho phương trình động lực ẩn phi tuyến hoặc hệ số biến thiên trên thang thời gian.
    • Áp dụng quan điểm xấp xỉ Hausdorff để đánh giá các phương pháp giải số bậc cao (như họ phương pháp Runge-Kutta trên thang thời gian).
    • Nghiên cứu bài toán ổn định và bán kính ổn định cho hệ chuyển mạch (switched systems) và hệ điều khiển vô hạn chiều trên miền thời gian không đồng nhất.

Giá trị tham khảo

  • Đối tượng thụ hưởng:

    • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Toán giải tích, Phương trình vi phân và tích phân, Giải tích số, và Lý thuyết điều khiển tự động.
    • Kỹ sư và nhà mô hình hóa hệ thống: Làm việc trong lĩnh vực mạng điện tử, cơ hệ robot, kinh tế lượng hoặc sinh thái học quần thể - nơi xuất hiện các hệ động lực vừa liên tục vừa rời rạc hoặc hệ vi phân - đại số suy biến.

  • Nội dung có giá trị tham khảo chuyên sâu:

    • Chương 1: Tài liệu tổng hợp chuẩn xác về giải tích cơ bản trên thang thời gian, biến đổi trụ $\xi_h(z)$, toán tử giải và các bất đẳng thức giải tích.
    • Chương 2: Kỹ thuật đánh giá xấp xỉ định lượng sử dụng đa thức suy rộng $h_k(t, \tau)$ và khoảng cách Hausdorff.
    • Chương 3: Phương pháp đại số xử lý chùm ma trận suy biến với chỉ số $\text{Ind} \ge 1$ và công thức tính bán kính ổn định chịu nhiễu hai vế.

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

1. Độ đo Lebesgue $\Delta$ và $\nabla$ của một điểm đơn lẻ ${t_0}$ trên thang thời gian được xác định như thế nào?
Theo Mục 1.3 của luận án, độ đo của một điểm đơn lẻ phụ thuộc vào vị trí và tính chất của điểm đó:

  • Đối với độ đo $m_\Delta$: Nếu $t_0 \in \mathbb{T} \setminus {T_{\max}}$, thì $m_\Delta({t_0}) = \sigma(t_0) - t_0 = \mu(t_0)$. Điểm dày phải có độ đo bằng 0, điểm phân tán phải có độ đo dương. Nếu $\mathbb{T}$ có điểm cực đại hữu hạn $T_{\max}$, thì $m_\Delta({T_{\max}}) = \infty$.
  • Đối với độ đo $m_\nabla$: Nếu $t_0 \in \mathbb{T} \setminus {T_{\min}}$, thì $m_\nabla({t_0}) = t_0 - \rho(t_0) = \nu(t_0)$. Điểm dày trái có độ đo bằng 0, điểm phân tán trái có độ đo dương.

2. Biến đổi trụ (Cylinder transformation) $\xi_h(z)$ được định nghĩa như thế nào và giữ vai trò gì?
Biến đổi trụ $\xi_h: \mathbb{C}h \to \mathbb{Z}h$ được định nghĩa bởi: $$\xi_h(z) = \begin{cases} \frac{\text{Ln}(1 + hz)}{h} & \text{khi } h > 0 \ z & \text{khi } h = 0 \end{cases}$$ trong đó $\text{Ln}$ là hàm logarit chính với miền giá trị $[-\pi i, \pi i)$. Biến đổi trụ đóng vai trò then chốt giúp xây dựng tường minh hàm mũ $e_p(t, t_0)$ trên thang thời gian bất kỳ thông qua tích phân $\int{t_0}^t \xi{\mu(s)}(p(s))\Delta s$, bảo đảm hàm mũ là nghiệm duy nhất của phương trình động lực vô hướng tuyến tính $x^\Delta = p(t)x, x(t_0) = 1$.

3. Khi hàm $f(t, x)$ thỏa mãn điều kiện Lipschitz theo cả hai biến, tốc độ hội tụ nghiệm được đánh giá ra sao?
Theo kết quả tại Chương 2 (công thức 0.24), nếu hàm vế phải $f(t, x)$ thỏa mãn điều kiện Lipschitz đối với cả biến thời gian $t$ và biến không gian $x$, thì sai số giữa nghiệm $x_n(t)$ trên thang thời gian $\mathbb{T}_n$ và nghiệm $x(t)$ trên thang thời gian $\mathbb{T}$ bị chặn tuyến tính theo khoảng cách Hausdorff giữa hai thang thời gian: $$|x_n(t) - x(t)| \le C d_H(\mathbb{T}, \mathbb{T}_n), \quad \forall t \in \mathbb{T} \cap \mathbb{T}_n : t_0 \le t \le T$$ với $C$ là một hằng số dương không phụ thuộc vào $n$.

4. Khái niệm chỉ số (index) của chùm ma trận ${A, B}$ ảnh hưởng như thế nào đến tính liên tục của phổ $\sigma(A, B)$?
Đối với phương trình ẩn $A x^\Delta = B x$ với $A$ suy biến, nghiệm chứa các thành phần bị ràng buộc bởi quan hệ đại số xác định bởi chỉ số $\text{Ind}(A, B)$. Dưới tác động của nhiễu hệ số $(A_n, B_n) \to (A, B)$, nếu chỉ số bị thay đổi đột ngột thì phổ $\sigma(A_n, B_n)$ sẽ không hội tụ về $\sigma(A, B)$. Do đó, luận án phải đưa ra điều kiện triệt tiêu ma trận chiếu $(A_n - A)Q = 0$ (với $\text{Ind} = 1$) hoặc $(A_n - A)\widehat{Q} = (B_n - B)\widehat{Q} = 0$ (với $\text{Ind} > 1$) để bảo đảm phổ hội tụ liên tục theo metric Hausdorff.

5. Công thức xác định bán kính ổn định phức cho phương trình động lực ẩn tuyến tính chịu nhiễu cấu trúc ở cả hai vế là gì?
Theo công thức (0.22) tại Chương 1 và Chương 3, đối với hệ $A x^\Delta(t) = B x(t)$ chịu nhiễu dạng $[A + D\Sigma E_1, B + D\Sigma E_2]$, bán kính ổn định phức $r(A, B; D, E; \mathbb{T})$ được tính bằng nghịch đảo của giá trị supremum chuẩn hàm truyền trên phần bù của miền ổn định đều $\mathcal{U}\mathbb{T}^c$: $$r(A, B; D, E; \mathbb{T}) = \left(\sup{\lambda \in \mathcal{U}_\mathbb{T}^c} |(E_1 \lambda - E_2)(\lambda A - B)^{-1} D|\right)^{-1}$$

Kết luận

Luận án của nghiên cứu sinh Nguyễn Thu Hà đã nghiên cứu thành công lý thuyết giải tích trên thang thời gian dưới giác độ một bài toán xấp xỉ định lượng và khảo sát sự phụ thuộc dữ liệu của phương trình động lực. Tác giả đã thiết lập các định lý về sự hội tụ nghiệm cùng tốc độ hội tụ cấp một theo khoảng cách Hausdorff cho phương trình động lực Delta và Nabla, thống nhất bản chất hội tụ của các phương pháp số trị sai phân kinh điển. Đồng thời, công trình giải quyết trọn vẹn bài toán về sự phụ thuộc liên tục của miền ổn định, phổ chùm ma trận và tính nửa liên tục dưới của bán kính ổn định cho phương trình động lực ẩn tuyến tính chịu nhiễu cấu trúc hai vế khi cả ma trận hệ số và miền thời gian cùng biến thiên.