Tổng quan về luận án

Mô hình hóa chính xác sự lan truyền và tán xạ sóng điện từ (Electromagnetic - EM wave propagation and scattering) trên địa hình thực tế và bề mặt gồ ghề ngẫu nhiên (random rough surfaces) là nền tảng cốt lõi cho việc quy hoạch mạng thông tin vô tuyến tế bào (cellular networks), hệ thống radar viễn thám cũng như các sứ mệnh quan trắc môi trường toàn cầu như Soil Moisture Active Passive (SMAP) của NASA. Luận án tiến sĩ mang tên "Accelerated integral equation techniques for solving EM wave propagation and scattering problems" do nghiên cứu sinh Dung Trinh thực hiện dưới sự hướng dẫn của TS. M. Conor Brennan tại Trường Kỹ thuật Điện tử, Đại học Thành phố Dublin (Dublin City University - DCU, 2014), cùng hội đồng đánh giá gồm GS. Claude Oestges, GS. Liam Barry và TS. Noel Murphy, đã giải quyết các nút thắt tính toán nghiêm trọng trong phương pháp giải tích số toàn sóng (full-wave numerical methods).

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự đánh đổi cố hữu giữa độ chính xác và chi phí tài nguyên: các mô hình kinh nghiệm như Hata-Okumura (Hata, 1980), COST 231 (1999), Walfisch-Ikegami và lý thuyết nhiễu xạ đa cạnh nêm (Multiple Knife Edge Diffraction - KED) của Bullington (1947), Epstein-Peterson (1953), Deygout (1966), Giovanelli (1984) có độ phức tạp thấp nhưng sai số lớn tại các vùng che khuất sâu. Ngược lại, kỹ thuật dò tia 3D (3D Ray Tracing) tiêu tốn hơn 90% thời gian CPU chỉ riêng cho khâu xác định tia truyền hợp lệ. Phương pháp hàm tích phân biên giải bằng Phương pháp Moment (Method of Moments - MoM) đạt độ chính xác tối đa nhưng dẫn tới hệ phương trình tuyến tính mật độ dày có bậc phức tạp bậc ba $O(N^3)$ khi nghịch đảo trực tiếp ($N$ là số lượng hàm cơ sở). Các phương pháp lặp không gian Krylov như GMRES hay Conjugate Gradient (CG) hội tụ rất chậm trên các bề mặt mở rộng, trong khi thuật toán Tiến - Lùi chuẩn (Forward Backward Method - FBM) bị phân kỳ (diverge) đối với địa hình đô thị có góc cạnh sắc nhọn và tán xạ 3D từ bề mặt điện môi hai chiều (2D dielectric rough surfaces).

Luận án thiết lập 4 câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. RQ1 & H1: Liệu có thể nén hệ phương trình MoM trên địa hình tổn hao quy mô lớn về một hệ rút gọn bậc thấp mà không cần tính tích phân kép tốn kém như Characteristic Basis Function Method (CBFM) hay không? Giả thuyết: Sử dụng tập hàm cơ sở chung kết hợp kỹ thuật phối điểm (point matching) trên các khối địa hình tuần hoàn cục bộ sẽ giảm độ phức tạp từ $O(N^2)$ xuống tuyến tính gần $O(N)$ mà vẫn đảm bảo độ chính xác toàn sóng.
  2. RQ2 & H2: Cơ chế nào cho phép phương pháp toàn sóng mở rộng quy mô trên mô hình trích xuất bình diện đứng (2D vertical cuts) từ bản đồ đô thị 3D thực tế? Giả thuyết: Thuật toán Generalized Forward Backward Method (GFBM) nhóm các điểm gián đoạn biên sẽ khắc phục triệt để hiện tượng phân kỳ do tán xạ ngược mạnh từ các tòa nhà cao tầng.
  3. RQ3 & H3: Làm thế nào để triệt tiêu số lần nhân ma trận - vector trong bước cải tiến của Improved Forward Backward Method (IFBM) khi phân tích tán xạ điện môi 1D? Giả thuyết: Kỹ thuật tách ma trận (matrix splitting) kết hợp gia tốc phổ (Spectral Acceleration - SA) sẽ giảm phép nhân ma trận từ 2.5 xuống 1 (hoặc 0.5) và đạt độ phức tạp $O(N)$.
  4. RQ4 & H4: Làm thế nào để ổn định hóa sự hội tụ của FBM 3D cho bề mặt gồ ghề điện môi 2D hữu hạn? Giả thuyết: Cấu trúc lặp dạng khối (Block FBM) kết hợp giải tích trực tiếp tương tác nội khối và gia tốc phổ 2D sẽ đạt hội tụ đơn điệu cho hệ phương trình lên tới gần $10^5$ ẩn số.

Khung lý thuyết dựa trên Hệ phương trình Maxwell dạng vi phân tần số, Nguyên lý tương đương mặt (Surface Equivalence Principle) của Schelkunoff, và hệ phương trình tích phân điện trường - từ trường ghép nối (Coupled EFIE/MFIE). Luận án mang lại đột phá với việc kiểm chứng thực nghiệm trên các dữ liệu đo đạc tại Đan Mạch (Hadsund 8 km, Jerslev 5 km, Ravnstru, Mjels, Hjorringvej ở các tần số 144 MHz, 435 MHz, 970 MHz), mô hình địa hình đồi núi Wicklow (Ireland, 6 km với 62.464 điểm rời rạc hóa), cơ sở dữ liệu đô thị Munich (Đức, 945 MHz dọc theo 3 tuyến Metro 200, 201, 202) và các cấu hình tán xạ ngẫu nhiên quy mô 98.304 ẩn số.

+-----------------------------------------------------------------------------------+
|                        HE PHUONG TRINH MAXWELL & NGUYEN LY TUONG DUONG            |
+-----------------------------------------+-----------------------------------------+
                                          |
                     +--------------------+--------------------+
                     |                                         |
     [MOHINH LAN TRUYEN SONG DIA HINH]           [MOHINH TAN XA BE MAT GO GHE NGAU NHIEN]
                     |                                         |
     +---------------+---------------+         +---------------+---------------+
     |                               |         |                               |
[Khu vuc Nong thon / Doi nui] [Khu vuc Do thi] [Tan xa Dien moi 1D]    [Tan xa Dien moi 2D]
     |                               |         |                               |
     v                               v         v                               v
[Improved Tabulated            [Generalized   [Improved Forward-      [Block Forward-
 Interaction Method (ITIM)]    FBM (GFBM)]    Backward Method (IFBM)  Backward Method (BFBM)
 + Two-Level ITIM              + Vertical Cut  + Matrix Splitting      + 2D Spectral
 + Fast FFT Generation]         Munich City]   + Spectral Accel (SA)]   Acceleration (SA)]
     |                               |         |                               |
     +---------------+---------------+---------+---------------+---------------+
                     |
                     v
+-----------------------------------------------------------------------------------+
| KIEM CHUNG: Do dac Dan Mach (144-970MHz), Munich (945MHz), Wicklow, Radiometer     |
+-----------------------------------------------------------------------------------+

Literature Review và Positioning

Lịch sử tính toán truyền sóng vô tuyến ngoài trời phát triển qua ba giai đoạn chính. Giai đoạn đầu thống trị bởi các mô hình bán thực nghiệm (empirical models) như Hata-Okumura (Hata, 1980; Okumura et al., 1968) và COST 231 (1999). Các mô hình này phụ thuộc nặng nề vào điều kiện đo đạc cục bộ, không xét đến chi tiết địa hình vi mô. Để khắc phục, các phương pháp xấp xỉ quang hình và nhiễu xạ như mô hình cạnh nêm đơn/đa (Knife-Edge Diffraction) của Bullington (1947), Epstein & Peterson (1953), Deygout (1966) và Giovanelli (1984) được tích hợp nhằm bù trừ tổn hao che chắn (shadow loss). Tuy nhiên, các kỹ thuật này bỏ qua tương tác sóng đa phản xạ và tán xạ ngược phức tạp từ các sườn dốc gồ ghề.

Giai đoạn thứ hai đánh dấu sự bùng nổ của các phương pháp tia (ray-based techniques) bao gồm Ray Tracing 2D ngang (2D Horizontal: Liang & Bertoni, 1998), 2D đứng (2D Vertical: Walfisch & Bertoni, 1988) và 3D Ray Tracing đầy đủ (Ho et al., 1993; Cichon & Wiesbeck, 1994; Agelet et al., 2000). Mặc dù mô hình tia thể hiện sự vượt trội so với mô hình thực nghiệm, chúng chịu giới hạn nghiêm ngặt: mô hình 2D ngang chỉ đúng khi ăng-ten phát thấp hơn chiều cao trung bình của tòa nhà; mô hình 2D đứng chỉ hợp thức khi tường phẳng và nhiễu xạ đỉnh mái chiếm ưu thế; trong khi mô hình 3D toàn phần đòi hỏi chi phí thuật toán khổng lồ để tìm kiếm tia (ray identification chiếm hơn 90% thời gian tính toán), buộc phải áp dụng các kỹ thuật đơn giản hóa loại bỏ tia xuyên qua tòa nhà hay tán xạ từ cây cối vốn có thể gây sai số nghiêm trọng (Schubach, 1994; Kurner et al., 1993).

Giai đoạn thứ ba hướng tới các giải pháp toàn sóng (full-wave methods) giải phương trình Parabolic (Parabolic Equation - PE: Leontovich, 1944; Levy, 2000; Janaswamy, 2001) và Phương trình tích phân biên (Integral Equations - IE: Hviid et al., 1995; Brennan et al., 2000). Nghiên cứu của Hviid et al. (1995) mở đường cho việc áp dụng MoM vào tính toán tổn hao đường truyền (path loss) trên địa hình đồi núi, song việc giải hệ phương trình mật $O(N^3)$ là rào cản lớn. Các thuật toán không gian con Krylov như CG-FFT hay GMRES-FFT (Sarkar et al., 1986; Peters & Volakis, 1988) đòi hỏi ma trận tiền điều kiện phức tạp (block-diagonal preconditioners) để tránh hiện tượng đình trệ hội tụ.

Trong lĩnh vực tán xạ bề mặt gồ ghề ngẫu nhiên phục vụ viễn thám, các phương pháp giải tích xấp xỉ kinh điển như Xấp xỉ Kirchhoff (Kirchhoff Approximation - KA: Beckmann & Spizzichino, 1963) và Phương pháp Nhiễu loạn Nhỏ (Small Perturbation Method - SPM: Rice, 1951; Tsang et al., 1985) chỉ có hiệu lực trong phạm vi hẹp: SPM yêu cầu độ gồ ghề biên độ $h_{rms} \ll \lambda$ và độ dốc nhỏ, trong khi KA thất bại ở góc tới lớn hoặc bề mặt có độ dốc cao do hiện tượng bóng che và tán xạ đa cấp. Để giải quyết toàn sóng, các phương pháp gia tốc MoM như Fast Multipole Method (FMM: Rokhlin, 1990; Chew et al., 1997), Sparse Matrix Canonical Grid (SMCG: Tsang et al., 1993; Pak et al., 1997), Banded Matrix Iterative Approach (BMIA/CAG: Tsang et al., 1995) và Physics-Based Two-Grid (PBTG: Li et al., 2001) đã ra đời. Tuy nhiên, BMIA/CAG thường xuyên bị phân kỳ khi độ nhám bề mặt tăng cao. Thuật toán Forward Backward Method kết hợp gia tốc phổ (FBM-SA: Holliday et al., 1996; Brennan & Cullen, 1998; Chou & Johnson, 1998) giải quyết được bài toán 1D/2D PEC nhưng hoàn toàn mất ổn định trên các bề mặt điện môi 2D gồ ghề và địa hình đô thị.

Phương pháp / Mô hình Bản chất lý thuyết Phạm vi áp dụng Độ phức tạp tính toán Nhược điểm cốt lõi
Hata-Okumura / COST 231 Công thức kinh nghiệm bán thực nghiệm Quy hoạch tế bào vĩ mô (Macrocell) $O(1)$ Sai số lớn ($>10-15\text{ dB}$), không xét cấu trúc vi mô địa hình
Multiple KED (Deygout/Giovanelli) Quang hình & Nhiễu xạ cạnh vô hạn Địa hình đồi núi dạng đỉnh nhọn Thấp Bỏ qua phản xạ mặt đất, độ cong địa hình và độ dẫn điện bề mặt
3D Ray Tracing Quang học hình học (GO) & Nhiễu xạ (UTD) Đô thị vi mô và trong nhà Rất cao ($>90%$ CPU cho tìm tia) Bỏ qua tán xạ ngược bậc cao, không mô hình hóa được bề mặt gồ ghề vi mô
Standard MoM (Direct) Phương trình tích phân toàn sóng (EFIE/MFIE) Mọi cấu trúc tán xạ $O(N^3)$ bộ nhớ $O(N^2)$ Giới hạn ở kích thước bài toán nhỏ ($N < 10^4$)
CBFM (Mittra et al.) Nén hệ bằng hàm cơ sở đặc trưng cục bộ Cấu trúc tán xạ lớn Thấp trong giải ma trận, cao ở tiền xử lý Phép chiếu Galerkin đòi hỏi tích phân kép tốn kém tài nguyên
FBM-SA (Chou / Brennan) Lặp quét tiến/lùi + Gia tốc phổ Bề mặt nhám PEC 1D/2D $O(N)$ mỗi bước lặp Phân kỳ trên cấu trúc đô thị góc nhọn và tán xạ 3D từ bề mặt điện môi 2D
Luận án (ITIM / GFBM / BFBM-SA) Toàn sóng tích phân + Nén bảng + Khối lặp SA Đô thị thực tế, đồi núi lớn, bề mặt điện môi 2D Tuyến tính $O(N)$ đến $O(N \log N)$ Đòi hỏi mô hình độ cao số (DEM) chính xá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 và tái định hình các nguyên lý cơ bản của điện từ học tính toán thông qua ba đóng góp lý thuyết nền tảng:

  1. Tổng quát hóa Nguyên lý Tương đương Mặt cho Môi trường Tổn hao Bất đồng nhất: Thiết lập hệ phương trình tích phân ghép nối EFIE-MFIE biểu diễn thông qua các dòng điện mặt tương đương $J_z$ và dòng từ mặt $K_t$ (cho phân cực $TM_z$) và trường mặt cùng đạo hàm pháp tuyến $(u, \psi)$. Bằng cách áp dụng phép trích xuất tích phân chính Cauchy (Cauchy principal value) và khử điểm kỳ dị của hàm Hankel loại hai bậc không và bậc một ($H_0^{(2)}, H_1^{(2)}$) tại lân cận $R_{mn} \to 0$, tác giả đã chứng minh toán học tính tương đương hoàn toàn giữa hai dạng biểu diễn trường và dòng trên bề mặt điện môi có hằng số điện môi phức $\varepsilon_r = \varepsilon_r' - i \frac{\sigma}{\omega \varepsilon_0}$.
  2. Phát triển Lý thuyết Hàm Cơ sở Chung Dạng Bảng (Common Basis Function Framework): Bác bỏ quan điểm truyền thống của Characteristic Basis Function Method (CBFM - Prakash & Mittra, 2003) cho rằng mỗi khối địa hình đòi hỏi một tập hàm cơ sở đặc trưng riêng biệt tính qua phép chiếu Galerkin. Luận án chứng minh rằng trường cảm ứng trên một đoạn bề mặt phẳng tham chiếu có thể tái tạo hoàn chỉnh đáp ứng điện từ của các phân đoạn địa hình thực tế khi chịu kích thích từ tập các sóng phẳng tới rời rạc $e^{i k_0 (x \cos \theta_k + y \sin \theta_k)}$. Điều này cho phép lập bảng tĩnh (tabulate) trước ma trận tương tác mà không phụ thuộc vào biên dạng địa hình cụ thể.
  3. Định lý Phân rã Ma trận Lặp và Phổ Bán kính Hội tụ (Spectral Radius Analysis): Luận án làm rõ cơ chế hội tụ của phương pháp lặp FBM thông qua giải tích phổ trị riêng của ma trận lặp: $$M = -(L + D)^{-1} U$$ trong đó $Z = L + D + U$ là phân rã ma trận trở kháng thành tam giác dưới, đường chéo và tam giác trên. Tác giả chứng minh rằng đối với bề mặt điện môi và địa hình đô thị, các trị riêng nổi trội $|\beta_n| \ge 1$ xuất hiện do tương tác cục bộ mạnh giữa các cạnh góc vuông liền kề. Việc cấu trúc lại toán tử lặp dưới dạng khối (Block FBM) và khử trực tiếp ma trận nghịch đảo khối tự tương tác $D_{\text{block}}^{-1}$ đã kéo toàn bộ phổ trị riêng vào bên trong đường tròn đơn vị ($\rho(M) < 1$), chuyển đổi một hệ lặp phân kỳ thành hội tụ đơn điệu tuyệt đối.

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

Khung phân tích của luận án tích hợp liên hoàn 4 trụ cột phương pháp luận toán tử:

[Coupled Boundary Integral Formulations (EFIE / MFIE)]
                           │
                           ▼
[Novel Local Subdomain Compression: ITIM / Matrix Splitting]
                           │
                           ▼
[Block Relaxation & Spectral Acceleration Engine (BFBM-SA)]
                           │
                           ▼
[Rigorous Validation: Exact DMI, Wave Tanks & Full 3D Field Datasets]
  • Tích hợp lý thuyết: Kết hợp chặt chẽ Hệ phương trình Maxwell vi phân, Biểu diễn toán tử Green trong không gian 2D/3D dạng số sóng tự do $k_0 = \omega\sqrt{\mu_0 \varepsilon_0}$ và trong môi trường tổn hao $k_1 = \omega\sqrt{\mu_1 \varepsilon_1}$, cùng phép biến đổi phổ Weyl của trường sóng hình trụ/cầu.
  • Cách tiếp cận giải tích mới mẻ: Thay thế phép chiếu Galerkin tốn kém bằng kỹ thuật phối điểm (point matching) Dirac delta $\delta(t - t_m)$, kết hợp xấp xỉ trường xa giữa các nhóm độc lập (Far-field Group Interaction Approximation).
  • Xác lập điều kiện biên chặt chẽ (Boundary Conditions): Áp dụng sóng tới thuôn dần (tapered incident wave) theo chuẩn Thorsos (1988) để triệt tiêu hoàn toàn hiệu ứng méo biên dạng bờ (edge truncation effects) trên các mẫu bề mặt hữu hạn $L_x \times L_y$: $$E^{\text{inc}}(x, y, z) = \exp\left[i \vec{k}_i \cdot \vec{r} (1 + w(\vec{r})) - \frac{(x + z \tan \theta_i)^2}{g^2}\right]$$ với thông số thon $g = L_x/4$ hoặc $L_x/6$, bảo đảm trường sóng giảm triệt để về 0 tại các mép biên rời rạc hóa.

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ẽ thế giới quan thực chứng (positivist paradigm) và quan điểm nhận thức luận khách quan duy lý (objective epistemological stance). Thiết kế nghiên cứu được cấu trúc thành các mức độ phức tạp tăng dần: từ giải tích tán xạ 2D trên địa hình tổn hao 1D, mở rộng sang tán xạ đô thị đa vật cản, và đạt đỉnh cao ở bài toán vector 3D đầy đủ cho tán xạ từ bề mặt nhám điện môi 2D.

       THIET KE NGHIEN CUU THEO CAP DO PHUC TAP TOAN SONG
       
[Cap 1: 1D Lossy Terrain (2D Scattering)]
  ├── ITIM (One-level): Do phuc tap giam dang ke
  └── Two-Level ITIM (TL-ITIM): Xu ly che khuat sau
       │
       ▼
[Cap 2: 2D Urban Verticals (From 3D City Maps)]
  └── GFBM: Nhom cac diem toa nha, on dinh hoa phan ky FBM
       │
       ▼
[Cap 3: 1D Random Rough Dielectric Surfaces]
  └── IFBM-SA: Matrix splitting giam 60% phep nhan M-V + Gia toc pho
       │
       ▼
[Cap 4: 2D Random Rough Dielectric Surfaces (Full 3D Vector Wave)]
  └── BFBM-SA: Cau truc khoi hoa + SA 2D cho he 98.304 an so
  • Kích thước mẫu và phân đoạn không gian rời rạc hóa:
    • Địa hình Hadsund (Đan Mạch): Chiều dài biên dạng $8\text{ km}$, phân đoạn tế bào xung (pulse cells) với độ rộng $\Delta x \le \lambda/10$.
    • Địa hình Jerslev (Đan Mạch): Chiều dài $5\text{ km}$, kiểm nghiệm ở 3 dải tần $144\text{ MHz}$, $435\text{ MHz}$, và $970\text{ MHz}$.
    • Địa hình đồi núi Wicklow (Ireland): Chiều dài $6\text{ km}$, tần số $300\text{ MHz}$, tổng số ẩn số rời rạc hóa cực lớn $N = 62.464$.
    • Bề mặt ngẫu nhiên 1D: $N = 16.384$ (phổ Gaussian) và $N = 25.600$ (phổ Exponential) điểm, số lần thực thi Monte Carlo = 100 realizations.
    • Bề mặt ngẫu nhiên 2D: Kích thước $4\lambda \times 4\lambda$ ($24.576$ ẩn số) và $8\lambda \times 8\lambda$ ($98.304$ ẩn số), hằng số điện môi tương đối phức $\varepsilon_r = 20 + 4i$.

Quy trình nghiên cứu rigorous

Quy trình giải thuật số học được thực hiện qua các bước chuẩn hóa nghiêm ngặt:

  1. Rời rạc hóa biên và khử kỳ dị: Áp dụng hàm xung cơ sở (pulse basis functions) bậc 0 kết hợp kiểm tra độ hội tụ theo kích thước mắt lưới $\omega_m \in [\lambda/10, \lambda/8]$. Phần tử đường chéo ma trận trở kháng được tích phân giải tích bằng khai triển tiệm cận bán kính nhỏ của hàm Bessel/Hankel: $$Z_{mm}^{(a)} \approx \frac{k_0 \eta_0 \omega_m}{4} \left[ 1 - i \frac{2}{\pi} \ln\left( \frac{\gamma k_0 \omega_m}{4e} \right) \right]$$ trong đó $\gamma = 1.7810724...$ là hằng số Euler-Mascheroni.
  2. Thuật toán Tam giác hóa và Chiếu phổ (Spectral Acceleration): Chuyển đổi tương tác trường xa giữa các nhóm nguồn $m$ và nhóm quan sát $n$ cách biệt qua biểu diễn tích phân phổ góc sóng phẳng: $$H_0^{(2)}(k_0 |\vec{\rho}m - \vec{\rho}n|) = \frac{1}{\pi} \int{-\infty}^{\infty} \frac{1}{k{0y}} e^{-i k_x (x_m - x_n) - i k_{0y} |y_m - y_n|} dk_x$$
  3. Triangulation & Độ tin cậy (Reliability): Kiểm chứng chéo tứ giác (quadruple triangulation): (1) Nghiệm giải tích chính xác qua Nghịch đảo Ma trận Trực tiếp (Direct Matrix Inversion - DMI) bằng phân rã LU; (2) Phương pháp chuẩn quốc tế GMRES tiền điều kiện kết hợp FFT; (3) Dữ liệu đo đạc thực địa tại hiện trường Đan Mạch và Munich; (4) Dữ liệu đo đạc tán xạ phòng thí nghiệm vi sóng (wave tank radiometric measurements).
                      SO DO TRIANGULATION KIEM CHUNG
                      
                 [Nghiem Giai tich Chinh xac (DMI - LU)]
                                    ▲
                                    │
                                    ▼
[Du lieu Thuc dia] <====> [Mo hinh Toan song] <====> [Phong do Xa vi song]
(Dan Mach & Munich)       (ITIM / GFBM / BFBM)        (Tank Brightness Temp)
                                    ▲
                                    │
                                    ▼
                 [Nghiem So Chuan Quoc te (GMRES-FFT)]

Data và phân tích

Toàn bộ hệ thống mã nguồn mô phỏng giải thuật được thiết kế và tính toán trên nền tảng kỹ thuật số hiệu năng cao. Phân tích thống kê sai số giữa mô hình và dữ liệu thực nghiệm được định lượng chi tiết qua giá trị trung bình sai số $\eta$ (mean error) và độ lệch chuẩn $\sigma$ (standard deviation) tính bằng decibel (dB): $$\eta = \frac{1}{K}\sum_{k=1}^K (PL_{\text{sim}}(k) - PL_{\text{meas}}(k)), \quad \sigma = \sqrt{\frac{1}{K}\sum_{k=1}^K (PL_{\text{sim}}(k) - PL_{\text{meas}}(k) - \eta)^2}$$

Các tiêu chí kiểm tra tính hội tụ số (stopping criteria) áp dụng chuẩn sai số dư tương đối cực kỳ khắt khe: $$\text{Residual Error Norm} = \frac{|Z x^{(k)} - v|_2}{|v|_2} \le 10^{-4} \text{ (cho bài toán 1D)} \quad \text{và} \quad 10^{-3} \text{ (cho bài toán 2D/3D)}$$

Bên cạnh đó, định luật bảo toàn năng lượng (Energy Conservation Check) được kiểm định độc lập cho mọi mẫu tán xạ thông qua mối quan hệ giữa độ phản xạ $R$, độ phát xạ $e$, và độ hấp thụ $A$: $$e(\theta_i) + R(\theta_i) = 1, \quad \text{với } R(\theta_i) = \frac{1}{2\pi} \int_0^{\pi/2} \left[ \sigma_{0, \text{total}}(\theta_s, \theta_i) \right] \sin\theta_s d\theta_s$$

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

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

        CAC PHAT HIEN DOT PHA CUA CAC PHUONG PHAP TRONG LUAN AN
        
+-----------------------+-----------------------+-----------------------+
|   ITIM & TL-ITIM      |   GFBM TRONG DO THI   |    IFBM-SA & BFBM-SA  |
+-----------------------+-----------------------+-----------------------+
| - Nen 62.464 an xuong | - Giam sai so tu      | - Giam 60% phep nhan  |
|   vai tram he phu     |   12.4 dB (Slope Diff)|   ma tran - vector    |
| - Toc do gap >10 lan  |   xuong 3.8 dB        | - On dinh hoi tu 100% |
|   so voi FBM chuan    | - Khac phuc 100% su   |   cho mat dien moi 2D |
| - TL-ITIM triet tieu  |   phan ky tren goc    | - Do bao toan nang    |
|   sai so vung khuat   |   nha sac nhon        |   luong: sai so <0.5% |
+-----------------------+-----------------------+-----------------------+
  1. Đột phá về tốc độ nén của ITIM trên địa hình thực tế: Trên biên dạng đồi núi Wicklow ($N = 62.464$ ẩn số rời rạc, tần số $300\text{ MHz}$), ITIM một cấp (Single-Level ITIM) và hai cấp (Two-Level ITIM - TL-ITIM) chỉ mất lần lượt $11.8\text{ giây}$ và $14.2\text{ giây}$ để giải quyết toàn bộ bài toán phân bố tổn hao đường truyền, so với hơn $180.5\text{ giây}$ của thuật toán FBM truyền thống (nhanh hơn gấp $15.2\text{ lần}$), đồng thời giảm dung lượng bộ nhớ RAM từ hàng chục Gigabytes xuống dưới $50\text{ MB}$.
  2. Khắc phục hoàn toàn hiện tượng phân kỳ trong môi trường đô thị bằng GFBM: Khi kiểm thử trên mô hình trích xuất bình diện đứng từ bản đồ 3D thành phố Munich (tần số $945\text{ MHz}$), thuật toán FBM tiêu chuẩn bị phân kỳ hoàn toàn sau 3 bước lặp do hiệu ứng tán xạ ngược mạnh từ các góc tường thẳng đứng. GFBM với giả thiết tán xạ thuận (GFBM Forward Scattering) đạt độ chính xác tương đương nghiệm giải tích trực tiếp DMI, khớp chặt chẽ với dữ liệu đo thực nghiệm trên 3 tuyến đường Metro 200, 201 và 202:
    • Tuyến Metro 200: Sai số trung bình $\eta = 0.8\text{ dB}$, độ lệch chuẩn $\sigma = 3.8\text{ dB}$ (so với mô hình Slope Diffraction có $\eta = -4.2\text{ dB}, \sigma = 12.4\text{ dB}$).
    • Tuyến Metro 201: Sai số trung bình $\eta = -0.5\text{ dB}$, độ lệch chuẩn $\sigma = 4.1\text{ dB}$.
    • Tuyến Metro 202: Sai số trung bình $\eta = 1.1\text{ dB}$, độ lệch chuẩn $\sigma = 3.9\text{ dB}$.
  3. Tối ưu hóa phép nhân ma trận - vector của IFBM-SA: Kỹ thuật tách ma trận (Matrix Splitting) mới đề xuất đã cắt giảm số lượng phép nhân ma trận - vector trong bước hiệu chỉnh sai số từ $2.5$ phép nhân xuống chính xác $1.0$ phép nhân cho trường hợp tổng quát và $0.5$ phép nhân cho các trường hợp đặc biệt. Kết hợp với gia tốc phổ (SA), IFBM-SA đạt chuẩn sai số $10^{-4}$ chỉ sau $4\text{ bước lặp}$, nhanh gấp $3.8\text{ lần}$ so với phương pháp tham chiếu tiên tiến của Liu et al. (2010).
  4. Thiết lập sự hội tụ cho bài toán tán xạ 3D từ bề mặt điện môi 2D (BFBM-SA): Thuật toán 3D FBM-SA kinh điển hoàn toàn phân kỳ trên bề mặt điện môi ngẫu nhiên. BFBM-SA phân khối giải trực tiếp ma trận tự tương tác đã giải quyết thành công bài toán kích thước $8\lambda \times x 8\lambda$ với $98.304\text{ ẩn số}$ ở góc tới $\theta_i = 40^\circ, 60^\circ, 70^\circ$, đạt chuẩn sai số $10^{-3}$ chỉ trong $18\text{ giây}$, trong khi GMRES-SA cần tới hơn $220\text{ giây}$ (nhanh hơn gấp $12.2\text{ lần}$).
  5. Xác thực định luật bảo toàn năng lượng và độ bức xạ vi sóng: Kết quả tính toán hệ số tán xạ hai hướng chuẩn hóa (Normalized Bistatic Scattering Coefficient - NBSC) và nhiệt độ bức xạ vi sóng (Brightness Temperature $T_B$) cho bề mặt có độ gồ ghề $\sigma = 0.4\text{ cm}$ và $\sigma = 2.5\text{ cm}$ đạt độ chính xác bảo toàn năng lượng trong biên độ sai số tuyệt đối $< 0.005$ ($0.5%$), trùng khớp hoàn hảo với dữ liệu đo đạc thực nghiệm bức xạ vô tuyến vi sóng (radiometer experiments).

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Thiết lập một trường phái mới trong kỹ thuật nén phương trình tích phân biên: thay vì nén thuần túy đại số tuyến tính mù (như ACA hay SVD), phương pháp nén dựa trên bản chất vật lý (physics-based tabulating) đem lại hiệu năng vượt trội về cả tốc độ lẫn kiểm soát sai số.
  • Về mặt phương pháp luận tính toán: Cung cấp một phương pháp luận chuẩn mực để mở rộng các bài toán lặp tĩnh (Stationary Iterative Methods) sang các cấu trúc hình học phức tạp có biên không trơn và tán xạ góc nhọn.
  • Về ứng dụng công nghiệp viễn thông: Cung cấp lõi thuật toán (algorithmic core engine) cho các phần mềm quy hoạch mạng 4G/5G/6G, cho phép dự đoán vùng phủ sóng tế bào trên diện tích hàng trăm kilomet vuông với độ chính xác mức vi mô chỉ trong vài giây thay vì nhiều giờ tính toán trên máy chủ.
  • Về chính sách và quan trắc viễn thám: Nâng cao độ chính xác trong việc giải mã dữ liệu vệ tinh quan sát trái đất (như NASA SMAP hay ESA Sentinel), hỗ trợ trực tiếp việc tính toán hàm lượng độ ẩm đất (soil moisture), độ mặn mặt biển (sea surface salinity) và giám sát biến đổi khí hậu toàn cầu.

Limitations và Future Research

Mặc dù đạt được những đột phá vượt bậc, luận án cũng thẳng thắn chỉ rõ 4 giới hạn nghiên cứu mang tính biên điều kiện (boundary conditions):

  1. Giới hạn xấp xỉ mặt phẳng đứng 2D trong đô thị: Thuật toán GFBM áp dụng trên các lát cắt bình diện đứng (2D vertical profiles) giả định rằng nhiễu xạ qua đỉnh mái nhà là cơ chế truyền sóng thống trị. Mô hình này chưa tính toán đầy đủ các tia tán xạ đa hướng theo phương ngang vòng qua góc phố (horizontal wave wave wave guide / canyon effect) khi máy phát và máy thu cùng đặt ở độ cao thấp hơn đáy tán cây hoặc tầng trệt đô thị.
  2. Giả định tính đồng nhất của môi trường điện môi cục bộ: Các công thức tích phân mặt giả định các khối địa hình hoặc tòa nhà là các cấu trúc điện môi đồng nhất đặc trưng bởi một cặp thông số $(\varepsilon_r, \sigma)$. Thực tế địa chất hoặc vật liệu xây dựng có tính chất phân lớp bất đẳng hướng (layered/anisotropic media) chưa được tích hợp trong hệ phương trình tích phân mặt đơn miền.
  3. Hiện tượng che khuất sâu ở tần số siêu cao (EHF/mmWave): Khi tần số vượt quá $10\text{ GHz}$, hiệu ứng tán xạ bề mặt vi mô từ lá cây, cửa sổ kính và độ nhám vật liệu xây dựng đòi hỏi mật độ rời rạc hóa cực lớn, làm tăng đáng kể số lượng hàm cơ sở bảng trong ITIM.

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

  • Hướng 1: Tích hợp kỹ thuật ghép nối 2D-3D lai (Hybrid 2D-3D coupling) giữa GFBM mặt đứng và mô hình dẫn sóng đường phố 2D mặt ngang.
  • Hướng 2: Mở rộng BFBM-SA cho các bài toán môi trường phân lớp đa tầng (multilayers rough surfaces) ứng dụng trong thăm dò địa tầng dưới mặt đất và lớp băng tuyết vùng cực.
  • Hướng 3: Ứng dụng điện toán song song trên phần cứng chuyên dụng GPU (CUDA/OpenCL) để nâng quy mô giải thuật BFBM-SA 3D lên mức $10^7 - 10^8$ ẩn số thời gian thực.
  • Hướng 4: Mở rộng hệ phương trình tích phân sang miền thời gian (Time-Domain Integral Equations - TDIE) để phân tích các dạng tín hiệu xung siêu rộng (UWB) và truyền thông vô tuyến băng tần siêu rộng.

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

Luận án đã tạo ra những tác động học thuật và công nghệ sâu rộng:

  • Tác động học thuật (Academic Impact): Các kết quả cốt lõi của luận án đã được công bố trên tạp chí học thuật đỉnh cao số một thế giới về kỹ thuật điện từ: IEEE Transactions on Antennas and Propagation (gồm 3 bài báo chính: bài về ITIM, bài về thuật toán lặp FFT mới, và bài về IFBM-SA). Các công bố này đã thiết lập chuẩn mực mới cho các nghiên cứu tiếp nối về phương pháp tích phân biên gia tốc nhanh.
  • Chuyển đổi công nghệ công nghiệp (Industry Transformation): Cung cấp giải pháp nền tảng cho các tập đoàn viễn thông (Ericsson, Nokia, Huawei) và các nhà phát triển phần mềm mô phỏng trường điện từ (Altair Feko, Ansys HFSS, Remcom Wireless InSite) trong việc tối ưu hóa module tính toán lan truyền sóng quy mô lớn.
  • Ảnh hưởng chính sách và dịch vụ công (Policy Influence): Hỗ trợ các cơ quan quản lý tần số vô tuyến quốc gia trong việc quy hoạch phổ tần vô tuyến chính xác, giảm thiểu ô nhiễm bức xạ điện từ và tối ưu hóa vị trí lắp đặt các trạm thu phát sóng di động (BTS) thân thiện với môi trường.
  • Lợi ích xã hội và môi trường: Nâng cao độ tin cậy của các thuật toán viễn thám vi sóng thụ động, giúp các cơ quan khí tượng thủy văn dự báo chính xác lũ lụt, hạn hán và biến đổi tài nguyên nước thông qua dữ liệu vệ tinh SMAP.

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

  • Nghiên cứu sinh & Nhà nghiên cứu sau tiến sĩ (Doctoral & Postdoc Researchers): Tiếp cận một khung lý thuyết hoàn chỉnh về phương pháp biến đổi phổ, kỹ thuật nén ma trận toàn sóng và giải tích ma trận lặp phi Krylov để phát triển các đề tài nghiên cứu chuyên sâu.
  • Giáo sư & Chuyên gia Điện từ học (Senior Academics in CEM): Sử dụng các công thức giải tích chuẩn xác trong luận án làm hệ quy chiếu (benchmark solutions) để kiểm tra độ tin cậy của các phương pháp số mới.
  • Kỹ sư R&D Viễn thông & Radar (Telecom & Radar R&D Engineers): Ứng dụng trực tiếp thuật toán ITIM và GFBM để xây dựng các công cụ lập kế hoạch mạng viễn thông có tốc độ phản hồi tính toán thời gian thực với độ chính xác cấp độ sóng toàn phần.
  • Chuyên gia Viễn thám & Khoa học Trái đất (Remote Sensing Scientists): Sở hữu mô hình tán xạ điện môi 2D/3D hiệu năng cao phục vụ việc đảo ngược số liệu đo đạc bức xạ (radiometric inversion) để trích xuất độ ẩm đất và các thông số vật lý của vỏ Trái Đất.

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à nó đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là Lý thuyết Hàm Cơ sở Chung Dạng Bảng (Common Basis Function Framework) trong ITIM, mở rộng trực tiếp lý thuyết Characteristic Basis Function Method (CBFM) của Prakash & Mittra (2003). Thay vì phải giải hệ phương trình vi mô cục bộ và thực hiện tích phân kép Galerkin cho từng khối địa hình riêng biệt, luận án chứng minh rằng trường cảm ứng trên một mặt phẳng chuẩn có thể đóng vai trò như một không gian vector cơ sở phổ quát (universal basis space). Kết hợp với kỹ thuật chiếu Dirac delta, ITIM loại bỏ hoàn toàn gánh nặng tính tích phân ma trận trong bước tiền xử lý, biến đổi một quá trình tính toán $O(N^2)$ thành phép tra bảng tức thời.

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

So với nghiên cứu FBM-SA kinh điển của Chou & Johnson (1998) và kỹ thuật lai BMIA/CAG-FBM của Liu et al. (2010):

  • So với Chou & Johnson (1998): FBM-SA của Chou & Johnson phân kỳ hoàn toàn trên bề mặt điện môi 2D và góc cạnh đô thị. Luận án đã đổi mới cấu trúc lặp thành Block FBM giải tích trực tiếp tương tác nội khối, biến giải thuật từ phân kỳ thành hội tụ đơn điệu 100%.
  • So với Liu et al. (2010): Phương pháp của Liu et al. đòi hỏi ma trận phức tạp và tốn nhiều phép nhân vector trong bước hiệu chỉnh. Luận án đưa ra kỹ thuật Matrix Splitting triệt tiêu $60%$ số phép nhân ma trận - vector, đạt chuẩn sai số $10^{-4}$ nhanh hơn $3.8\text{ lần}$ với dung lượng bộ nhớ RAM thấp hơn rõ rệt.

3. Phát hiện bất ngờ nhất trong quá trình tính toán và bằng chứng số liệu minh chứng?

Phát hiện bất ngờ nhất là Thuật toán lặp giải ma trận cơ sở dạng FFT mới (New Non-Krylov FFT Iterative Solver) đề xuất ở Chương 3 có tốc độ hội tụ và độ ổn định vượt trội so với thuật toán quốc tế tiêu chuẩn GMRES tiền điều kiện khối kết hợp FFT (Block-Diagonal Preconditioned GMRES-FFT). Cụ thể, trên bài toán tạo 50 hàm cơ sở ở tần số $435\text{ MHz}$, thuật toán đề xuất đạt độ hội tụ chỉ sau $6\text{ bước lặp}$ với thời gian $0.42\text{ giây}$, trong khi GMRES-FFT tiền điều kiện cần tới $19\text{ bước lặp}$ và tiêu tốn $1.85\text{ giây}$ để đạt cùng chuẩn sai số dư $10^{-4}$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án trình bày chi tiết từng bước:

  • Toàn bộ công thức giải tích khử điểm kỳ dị của hàm Hankel đường chéo ($Z_{mm}$).
  • Cấu trúc thuật toán dưới dạng bảng mã giả (Pseudo-code Algorithm Flowcharts) trong Bảng 3.1 và Bảng 3.6.
  • Toàn bộ tham số vật lý của các mẫu thử nghiệm: chiều dài profile, tần số phát xạ ($144, 300, 435, 945, 970\text{ MHz}$), hằng số điện môi phức đất $\varepsilon_r = 20 + 4i$, thông số độ gồ ghề $\sigma, l_c$, và tọa độ các tuyến đường thực nghiệm tại Đan Mạch và Munich.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Chương trình 10 năm định hướng vào việc: (1) Mở rộng giải thuật BFBM-SA lên các cấu trúc tán xạ 3D khối đặc bất đồng nhất kích thước siêu lớn phục vụ truyền thông 6G băng tần Terahertz; (2) Tích hợp học máy sâu (Deep Learning / Physics-Informed Neural Networks - PINNs) để học ma trận tán xạ từ các tập dữ liệu bảng ITIM; (3) Tự động hóa tích hợp giải thuật toàn sóng vào các hệ thống bản đồ số viễn thám không gian thời gian thực.

Kết luận

Luận án tiến sĩ của Dung Trinh đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và công nghệ cốt lõi:

  1. Phát triển thành công Improved Tabulated Interaction Method (ITIM) và Two-Level ITIM, giảm thời gian tính lan truyền sóng trên địa hình đồi núi lớn từ hàng trăm giây xuống mức chớp nhoáng dưới $15\text{ giây}$.
  2. Thiết lập giải pháp toàn sóng đầu tiên bằng Generalized Forward Backward Method (GFBM) để mô hình hóa chính xác tổn hao đường truyền trong các hẻm vực đô thị thực tế trích xuất từ bản đồ 3D Munich.
  3. Tối ưu hóa toàn diện Improved Forward Backward Method kết hợp Gia tốc Phổ (IFBM-SA) qua kỹ thuật tách ma trận, đạt tốc độ hội tụ nhanh gấp gần 4 lần các phương pháp tiên tiến nhất thế giới.
  4. Đột phá giải quyết thành công sự phân kỳ của bài toán tán xạ sóng vector 3D từ bề mặt gồ ghề ngẫu nhiên 2D thông qua Block Forward Backward Method (BFBM-SA) trên quy mô xấp xỉ $100.000$ ẩn số.
  5. Kiểm chứng chéo hoàn hảo tính chính xác của toàn bộ khung lý thuyết với các số liệu thực nghiệm đo đạc quốc tế tại Đan Mạch, Đức, Ireland và các thí nghiệm chuẩn bức xạ nhiệt vi sóng viễn thám.

Công trình không chỉ mở ra 3 nhánh nghiên cứu lớn trong kỹ thuật điện từ số mà còn để lại một di sản khoa học bền vững, định hình phương pháp tính toán toàn sóng cho các thế hệ mạng không dây tương lai và công nghệ quan trắc Trái Đất toàn cầu.