Tổng quan về luận án
Hệ thống vệ tinh định vị toàn cầu (Global Navigation Satellite Systems - GNSS), điển hình là GPS, GLONASS, Galileo và BeiDou, đóng vai trò nền tảng trong việc cung cấp thông tin Vị trí, Vận tốc và Thời gian (Position, Velocity, and Time - PVT) cho các hạ tầng trọng yếu toàn cầu như hàng không dân dụng, phương tiện tự hành (UAV/UGV), lưới điện thông minh, nông nghiệp chính xác và đồng bộ hóa giao dịch tài chính liên ngân hàng. Tuy nhiên, do quỹ đạo truyền dẫn không gian sâu với khoảng cách hơn 20.000 km, cường độ tín hiệu Signal-in-Space (SIS) khi đến bề mặt Trái Đất cực kỳ yếu, thường dao động ở mức danh định $-158.5\text{ dBW}$ đối với GPS L1 C/A và $-157.0\text{ dBW}$ đối với Galileo E1. Cấu trúc tín hiệu dân dụng mở cùng mức công suất cực thấp này khiến máy thu GNSS đối mặt với nguy cơ tê liệt trước các dạng can nhiễu vô tuyến có chủ đích (Intentional Radio Frequency Interference - RFI), đặc biệt là kỹ thuật tấn công giả mạo tín hiệu (Spoofing Attack).
Khác với can nhiễu phá sóng (Jamming) chỉ gây mất tín hiệu cục bộ, tấn công giả mạo phát ra các cấu trúc mã giả mạo tinh vi nhằm chiếm quyền điều khiển vòng bám (tracking loops) của máy thu, từ đó thao túng dữ liệu PVT mà hệ thống nạn nhân không hề phát hiện. Trước thực trạng các nền tảng vô tuyến điều khiển bằng phần mềm (Software-Defined Radio - SDR) chi phí thấp ngày càng phổ biến, các hình thức tấn công đã chuyển dịch từ mức độ sơ cấp (Simplistic) sang trung cấp (Intermediate) và tinh vi (Sophisticated).
Khoảng trống nghiên cứu (Research Gap) cốt lõi mà luận án chỉ ra là: Các thuật toán phòng thủ góc tới (Angle of Arrival - AoA) chi phí thấp hiện nay, tiêu biểu là bộ dò Tổng bình phương vi phân pha (Sum of Squares - SoS detector) phát triển bởi Falletti và Nicola, dựa trên giả định tuyệt đối rằng máy thu bị thao túng toàn bộ chùm vệ tinh từ một nguồn phát duy nhất. Khi rơi vào trạng thái bám hỗn hợp (Mixed Tracking) — tức máy thu chỉ bị lừa bởi một tập con vệ tinh giả mạo ($S$) trong khi vẫn bám bắt các vệ tinh xác thực ($A$) — chỉ số SoS bị vô hiệu hóa hoàn toàn, dẫn đến việc chấp nhận giả thuyết sai lầm $H_1$ (bỏ sót tấn công). Đồng thời, hiện tượng trượt chu kỳ sóng mang (Cycle Slips) và sự giao cắt quỹ đạo vi phân pha chưa được mô hình hóa chặt chẽ về mặt toán học, dẫn đến tỷ lệ cảnh báo sai (False Alarm Rate) và bỏ sót cảnh báo (Missed Detection Rate) cao trong môi trường vận hành thực tế.
Để giải quyết triệt để vấn đề này, luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:
- Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để phát hiện và cô lập chính xác một tập con tín hiệu GNSS giả mạo trong điều kiện bám hỗn hợp (Mixed Tracking) mà không yêu cầu đồng bộ hóa xung nhịp phức tạp giữa hai máy thu?
- Câu hỏi nghiên cứu 2 ($RQ_2$): Làm cách nào để triệt tiêu ảnh hưởng của hiện tượng trượt chu kỳ (Cycle Slips) và phân tách hiện tượng giao cắt quỹ đạo hình học giữa vệ tinh thật và vệ tinh giả trên mặt phẳng pha vi phân?
- Câu hỏi nghiên cứu 3 ($RQ_3$): Cơ chế toán học nào cho phép tự động phân cụm và nhận dạng tín hiệu giả mạo đa hướng (Sophisticated Attacks) mà không cần cấu hình ngưỡng kinh nghiệm tĩnh?
Hệ thống giả thuyết tương ứng được kiểm định:
- Giả thuyết 1 ($H_1$): Sự phân tán của các giá trị vi phân bậc hai sóng mang (Double Differences - DD) giữa các ăng-ten phân cách cố định sẽ co cụm cục bộ quanh một hằng số hình học $M_{cnt}$ đối với chùm tín hiệu phát từ cùng một bộ phát giả mạo, tạo thành một vùng tương đồng $\Sigma_j^o$ có thể phân tách thống kê.
- Giả thuyết 2 ($H_2$): Mô hình hóa xu hướng vi phân pha thập phân theo dạng tuyến tính từng đoạn (Piecewise Linear Model) sẽ cung cấp thêm tham số độ dốc (Slope) và hệ số chặn (Intercept), loại bỏ hoàn toàn các điểm giao cắt tức thời gây cảnh báo sai.
- Giả thuyết 3 ($H_3$): Mô hình hỗn hợp Gauss (Gaussian Mixture Model - GMM) kết hợp thuật toán Kỳ vọng - Cực đại hóa (Expectation-Maximization - EM) có khả năng tự động học không gian phân phối xác suất của góc tới để phân loại tín hiệu thực/giả trong kịch bản tấn công đa hướng.
Khung lý thuyết của luận án tích hợp chặt chẽ giữa xử lý tín hiệu thống kê (Statistical Signal Processing), lý thuyết kiểm định tỷ số hợp lý tổng quát (Generalized Likelihood Ratio Test - GLRT), lý thuyết vi phân sóng mang kép (Differential Carrier Phase) và học máy thống kê (Unsupervised Machine Learning). Luận án được thực hiện trên cấu hình 2 ăng-ten với đường cơ sở ngắn (Baseline $D \approx 2\text{ m}$), sử dụng máy thu thương mại Septentrio AsteRx4 OEM, kiểm thử trên dải tỷ số sóng mang trên mật độ tạp âm ($C/N_0$) từ $32\text{ dB-Hz}$ đến $52\text{ dB-Hz}$, trên các tập dữ liệu mô phỏng tĩnh và động (quỹ đạo bay TRJ1, TRJ2) kéo dài hàng giờ với $100%$ độ chính xác trong các thử nghiệm đề xuất.
Literature Review và Positioning
Các nghiên cứu chống giả mạo GNSS trong y văn quốc tế được phân chia thành ba nhóm phương pháp tiếp cận chính:
CÁC HƯỚNG PHÒNG THỦ GIẢ MẠO GNSS
|
+----------------------------------------+---------------------------------------+
| | |
MẬT MÃ HỌC (Cryptographic) XÁC THỰC NGOẠI VI (External) ĐẶC TRƯNG TÍN HIỆU (Signal Features)
- Humphrey et al., Wesson et al. - Lo et al., Blanch et al. - Psiaki et al. (VSD)
- Mã hóa tin điện & lan truyền - Tích hợp INS/IMU, eLORAN, - Akos et al. (AGC Monitoring)
- Hạn chế: Cần sửa đổi chuẩn SIS, Cellular, Wi-Fi - Broumandan et al. (Moving Antenna)
không bảo vệ được tín hiệu mở - Hạn chế: Tăng chi phí phần cứng, - Falletti & Nicola (Dual-Antenna SoS)
sai số tích lũy của cảm biến phụ - Hien et al. (D3, LR-D3 & GMM-AoA)
- Nhóm phương pháp mật mã học (Cryptographic Defenses): Các công trình của Humphreys et al. và Wesson et al. đề xuất xác thực bản tin hàng hải (Navigation Message Authentication - NMA) và mã hóa mã trải phổ (Spreading Code Encryption). Mặc dù cung cấp mức độ bảo mật toán học vững chắc, phương pháp này đòi hỏi sự thay đổi toàn diện từ phân khúc điều khiển vệ tinh, không thể triển khai trên các chùm vệ tinh hiện hữu và hoàn toàn vô hiệu trước các cuộc tấn công tái phát lặp trễ (Meaconing).
- Nhóm phương pháp xác thực ngoại vi (External Verification): Nghiên cứu của Lo et al. và Blanch et al. sử dụng hệ thống dẫn đường quán tính (Inertial Navigation Systems - INS/IMU), trạm mặt đất eLORAN hoặc định vị tế bào sóng viễn thông để đối soát PVT. Nhược điểm chí mạng là chi phí cảm biến cấp chiến thuật rất cao, trọng lượng lớn và hiện tượng sai số trôi dạt (Drift error) tích lũy theo thời gian.
- Nhóm phương pháp dựa trên đặc trưng tín hiệu (Signal Features): Bao gồm kỹ thuật giám sát méo hàm tương quan (Vestigial Signal Defense - VSD) của Psiaki et al. (2013), kiểm định tính phù hợp Chi-square Goodness-of-Fit (GoF) của Broumandan et al. (2016), và kỹ thuật giám sát độ lợi điều khiển khuếch đại tự động (AGC Monitoring) của Akos et al. (2012). Mặc dù VSD và GoF có khả năng phát hiện tương đối tốt, chúng đòi hỏi quyền can thiệp trực tiếp vào miền tương quan phức tạp bên trong máy thu phần mềm (Software Receiver), không thể ứng dụng trên các module phần cứng OEM thương mại đại trà.
Trong bối cảnh đó, phương pháp phòng thủ không gian dựa trên góc tới (Angle of Arrival - AoA) được công nhận là cơ chế phòng thủ độc lập và vững chắc nhất (Montgomery et al., 2011). Tuy nhiên, các hệ thống ước lượng hướng sóng truyền thống đòi hỏi mảng ăng-ten đồng bộ pha dao động phức tạp (Coherent Multi-Antenna Arrays), cồng kềnh và đắt đỏ. Bước đột phá gần đây của Falletti và Nicola (2018, 2021) đã giới thiệu bộ dò Tổng bình phương vi phân pha (Sum of Squares - SoS detector) sử dụng cặp ăng-ten độc lập. Dẫu vậy, SoS bộc lộ hai điểm nghẽn nghiêm trọng:
- Thứ nhất, SoS chỉ thiết lập cho bài toán phân loại nhị phân toàn phần: toàn bộ vệ tinh là thật ($H_1$) hoặc toàn bộ vệ tinh bị giả mạo ($H_0$). Khi kẻ tấn công thực hiện chiến thuật giả mạo trung cấp (Intermediate Attack), chỉ phát sóng 3 hoặc 4 vệ tinh giả mạo nhằm lái dần vị trí của máy thu trong khi các vệ tinh còn lại vẫn là tín hiệu thật (hiện tượng Mixed Tracking), đại lượng kiểm định $\Lambda_{\text{SoS}}$ bị chi phối bởi các vệ tinh thật, dẫn đến việc bỏ sót hoàn toàn cuộc tấn công.
- Thứ hai, SoS cực kỳ nhạy cảm với hiện tượng trượt chu kỳ (Cycle Slips) của vòng khóa pha (PLL). Một bước nhảy pha nguyên lần bước sóng làm $\Lambda_{\text{SoS}}$ tăng vọt, gây ra báo động giả liên tục.
Vị trí học thuật của luận án được xác lập bằng việc mở rộng và vượt qua các giới hạn của Falletti & Nicola (2021) và Broumandan et al. (2016): Luận án phát triển thuật toán Độ phân tán của vi phân kép (Dispersion of Double Differences - D3), cải tiến thành mô hình hồi quy tuyến tính LR-D3 và mở rộng sang mô hình học máy GMM. Luận án thiết lập một hệ phương pháp hoàn chỉnh, có khả năng phát hiện tập con giả mạo, triệt tiêu trượt chu kỳ, xử lý giao cắt quỹ đạo pha và thích ứng với tấn công đa hướng, đồng thời giữ nguyên ưu điểm không cần đồng bộ xung nhịp và tương thích trực tiếp với các máy thu thương mại chi phí thấp.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đóng góp quan trọng vào lý thuyết xử lý tín hiệu định vị vệ tinh vi phân bằng cách tái cấu trúc và hoàn thiện mô hình vi phân pha sóng mang kép trong bài toán an ninh không gian:
- Mở rộng mô hình vi phân sóng mang kép (Differential Carrier-Phase Model):
Phương trình đo pha sóng mang của máy thu đơn lẻ cho vệ tinh thứ $i$ được biểu diễn:
$$\phi_i = r_i + N_i \lambda + c(\delta t_i - \delta t_u) - \varepsilon_{I,i} + \varepsilon_{T,i} + \varepsilon_i$$
Trong đó $r_i$ là khoảng cách hình học, $N_i$ là số nguyên bất định vòng pha, $\lambda$ là bước sóng tín hiệu, $\delta t_i$ và $\delta t_u$ là sai số đồng hồ vệ tinh và máy thu, $\varepsilon_{I,i}$ và $\varepsilon_{T,i}$ là trễ tầng điện ly và đối lưu, $\varepsilon_i$ là nhiễu trắng và đa đường.
Khi sử dụng hai máy thu tạo thành đường cơ sở ngắn ($D \approx 2\text{ m}$), sai số tầng khí quyển bị triệt tiêu hoàn toàn. Hiệu vi phân đơn ($\Delta\phi_i$) và vi phân kép ($\Delta\nabla\varphi_i$) đối với vệ tinh tham chiếu $r$ được rút gọn thành:
$$\Delta\nabla\varphi_i = \frac{D}{\lambda}[\cos(\alpha_i) - \cos(\alpha_r)] + \Delta\nabla N_i + \Delta\nabla\varepsilon_i$$
Bằng toán tử lấy phần lẻ (Fractional DD), thành phần bất định nguyên $\Delta\nabla N_i$ bị loại bỏ:
$$\mu_i = \Delta\nabla\varphi_i - \text{round}(\Delta\nabla\varphi_i) = \beta[\cos(\alpha_i) - \cos(\alpha_r)] + \Delta\nabla\varepsilon_i$$
với $\beta = D/\lambda$.
KHUNG LÝ THUYẾT VI PHÂN PHA D3 & LR-D3
|
+---------------------------------------+---------------------------------------+
| |
[TÍN HIỆU XÁC THỰC: Tập A] [TÍN HIỆU GIẢ MẠO: Tập S]
- Các góc tới độc lập: \alpha_l \neq \alpha_k - Góc tới chung duy nhất: \alpha_{cnt}
- Giá trị hình học phân tán ngẫu nhiên: - Giá trị hình học hội tụ đồng nhất:
m_l = \beta[\cos(\alpha_l) - \cos(\alpha_r)] M_{cnt} = \beta[\cos(\alpha_{cnt}) - \cos(\alpha_r)]
- Khoảng cách cặp: |\mu_j - \mu_k|^2 > \xi_k^2 - Khoảng cách cặp: |\mu_i - \mu_k|^2 \le \xi_k^2
- Phân phối Chi-square phi tâm \chi^2(\lambda) - Phân phối Chi-square trung tâm \chi^2(0)
-
Chứng minh toán học về sự đồng nhất góc tới của tập con giả mạo:
Luận án chứng minh rằng đối với mọi vệ tinh giả mạo $i \in S$ phát ra từ cùng một nguồn, góc tới vật lý là duy nhất ($\alpha_i = \alpha_{cnt}$). Do đó:
$$M_{cnt} = \beta[\cos(\alpha_{cnt}) - \cos(\alpha_r)], \quad \forall i \in S$$
Đặc biệt, nếu vệ tinh tham chiếu $r$ cũng thuộc tập giả mạo $S$, giá trị hình học $M_{cnt} \equiv 0$. Ngược lại, đối với các vệ tinh xác thực $l \in A$, các giá trị $m_l = \beta[\cos(\alpha_l) - \cos(\alpha_r)]$ là hoàn toàn độc lập và biến thiên liên tục theo quỹ đạo chuyển động của vệ tinh trên bầu trời.
-
Mô hình phân phối thống kê của đại lượng kiểm định $\Lambda_{D3}$:
Luận án chỉ ra rằng đại lượng kiểm định cặp $\Lambda_{D3}(j,k) = |\mu_j - \mu_k|^2$, khi chuẩn hóa theo phương sai nhiễu $(\sigma_j^2 + \sigma_k^2)$, tuân theo phân phối Chi-square ($\chi^2$) với 1 bậc tự do:
- Dưới giả thuyết $H_0$ (cặp tín hiệu giả mạo cùng nguồn), thống kê tuân theo phân phối $\chi^2$ trung tâm (Central $\chi^2$).
- Dưới giả thuyết $H_1$ (cặp tín hiệu thật hoặc một thật - một giả), thống kê tuân theo phân phối $\chi^2$ phi tâm (Non-central $\chi^2$) với tham số phi tâm $\lambda$:
$$\lambda = \frac{(m_j - m_k)^2}{\sigma_j^2 + \sigma_k^2}$$
Điều này đặt cơ sở lý thuyết chuẩn xác để xác định giải tích xác suất bỏ sót cảnh báo ($P_{md}$) và xác suất cảnh báo sai ($P_{fa}$), thay thế hoàn toàn việc gán ngưỡng mò mẫm theo kinh nghiệm.
Khung phân tích độc đáo
Khung phân tích của luận án kết hợp 3 trụ cột toán học vững chắc:
+---------------------------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH BA TRỤ CỘT ĐỘC ĐÁO |
+------------------------------------+----------------------------------+---------------------------+
| TRỤ CỘT 1: THUẬT TOÁN D3 | TRỤ CỘT 2: HỒI QUY TUYẾN TÍNH | TRỤ CỘT 3: HỌC MÁY GMM |
| (Dispersion of Double Differences) | (LR-D3 Piecewise Linear Model) | (Gaussian Mixture Models) |
| - Không gian tương đồng \Sigma_j^o | - Bổ sung tham số độ dốc (Slope) | - Tự động học mật độ xác |
| - Kiểm định cụm \ge 3 phần tử | và hệ số chặn (Intercept) | suất phân phối pha |
| - Giám sát trượt chu kỳ Doppler | - Triệt tiêu lỗi giao cắt quỹ đạo| - Phân loại đa hướng |
| - Lọc trung bình trượt thời gian | trên mặt phẳng pha vi phân | - Tối ưu bằng thuật toán |
| | - Đạt P_fa = 0% và P_md = 0% | Kỳ vọng - Cực đại (EM) |
+------------------------------------+----------------------------------+---------------------------+
- Vùng tương đồng (Similarity Region): Xác định không gian quan trắc $\Sigma_j^o = [\mu_j - \xi_j, \mu_j + \xi_j]$. Tiêu chuẩn kích hoạt báo động giả mạo được lượng hóa: Phát hiện cuộc tấn công khi tồn tại ít nhất 3 tín hiệu vi phân pha rơi vào cùng một vùng tương đồng:
$$\exists (i,j,k) \in (S \cup A): |\mu_i - \mu_k|^2 \le \xi_k^2 \quad \text{và} \quad |\mu_j - \mu_k|^2 \le \xi_k^2$$
- Điều kiện biên (Boundary Conditions): Chiều dài đường cơ sở $D \approx 2\text{ m}$ (đảm bảo triệt tiêu trễ khí quyển mà không gây tương quan nhiễu không gian), tỷ số tín hiệu trên tạp âm $C/N_0 \ge 32\text{ dB-Hz}$, tần số cập nhật dữ liệu $1 - 10\text{ Hz}$.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án được thiết kế theo trường phái thực chứng (Positivism) kết hợp phương pháp luận diễn dịch thực nghiệm (Deductive-Empirical Methodology). Nghiên cứu tích hợp mô hình hóa giải tích toán học với kiểm chứng thực nghiệm Hardware-in-the-Loop trên các máy thu cấp công nghiệp.
Thiết kế nghiên cứu đa tầng bao gồm:
- Tầng mô hình toán giải tích: Xây dựng phương trình hàm phân phối xác suất tích lũy (CDF) và mật độ xác suất (PDF) của $\Lambda_{D3}$, thiết lập đường cong đặc trưng thu nhận máy thu (Receiver Operating Characteristic - ROC).
- Tầng mô phỏng chuyên sâu: Sử dụng công nghệ SDR xây dựng bộ mô phỏng tín hiệu GNSS đa hướng (Multi-directional GNSS Signal Simulator), can thiệp pha sóng mang để tái tạo môi trường tấn công tinh vi.
- Tầng thực nghiệm phần cứng: Thu thập dữ liệu đo thực từ cặp module máy thu OEM Septentrio AsteRx4 đa tần, xử lý trên môi trường MATLAB®.
+---------------------------------------------------------------------------------------------------+
| QUY TRÌNH THỰC NGHIỆM VÀ KIỂM CHUẨN |
+---------------------------------------------------------------------------------------------------+
| [Dữ liệu GNSS thực/mô phỏng] --> [Cặp máy thu Septentrio AsteRx4] --> [Tính vi phân kép DD \mu_i]|
| | |
| +-----------------------------------------------------------------------------+ |
| | |
| v v
| [Khối Doppler Shift Monitor] [Khối Zero-Baseline] |
| - Dự báo tần số Doppler qua \phi(t) - \phi(t-T_0) - D = 0 qua bộ chia công|
| - Ngưỡng phát hiện trượt chu kỳ: \Delta f_d > 10 Hz suất RF Power Splitter|
| - Loại bỏ ngay lập tức \mu_i bị lỗi khỏi tập quan trắc - Lập bản đồ ngưỡng |
| | \xi_k theo dải C/N_0 |
| +-----------------------------------------------------------------------------+ |
| | |
| v |
| [Khối xử lý thuật toán D3 / LR-D3 / GMM-EM] <----------------------------------+ |
| - Đánh giá độ phân tán trên cửa sổ trượt 5 giây |
| - Phân cụm thống kê / Phân tích hồi quy tuyến tính độ dốc & hệ số chặn |
| | |
| v |
| [Quyết định trạng thái an toàn: Cảnh báo giả mạo hoặc Xác thực danh định] |
+---------------------------------------------------------------------------------------------------+
Quy trình nghiên cứu rigorous
- Chiến lược xác lập ngưỡng quyết định động (Dynamic Thresholding via Zero-Baseline):
Để triệt tiêu giả định chủ quan về nhiễu, tác giả thiết lập cấu hình đường cơ sở bằng không ($D = 0$), nối hai máy thu vào cùng một ăng-ten thông qua bộ chia công suất RF (Power Splitter). Khi đó $(r_i^{(1)} - r_i^{(2)}) = 0$, số đo vi phân thập phân chỉ còn thuần túy là nhiễu nhiệt: $\mu_i = \Delta\nabla\varepsilon_i$. Bằng cách quét $C/N_0$ từ $32\text{ dB-Hz}$ đến $52\text{ dB-Hz}$ (các mốc cụ thể: $32, 38, 42, 47, 52\text{ dB-Hz}$), luận án lập bản đồ thực nghiệm hàm tương quan giữa độ lệch chuẩn $\sigma_{D3}$ và $C/N_0$, cho phép tự động điều chỉnh ngưỡng $\xi_k^2$ theo thời gian thực.
- Bộ giám sát dịch pha Doppler triệt tiêu trượt chu kỳ (Doppler Shift Monitor):
Tần số Doppler dự báo $\hat{f}{d,i} = [\phi_i(t) - \phi_i(t - T_0)] / T_0$ được so sánh trực tiếp với tần số Doppler xuất ra từ kênh bám PLL của máy thu. Khi $|\hat{f}{d,i} - f_{d,i}| > 10\text{ Hz}$, phép đo bị gán nhãn trượt chu kỳ và loại bỏ tức thì khỏi tập tính toán $\Lambda_{D3}$, bảo vệ hệ số kiểm định không bị nhảy đột biến.
- Giảm thiểu sai số bằng trung bình trượt thời gian (Time Averaging):
Áp dụng cửa sổ trung bình trượt với độ trễ (latency) tối ưu hóa là $\eta = 5\text{ giây}$ trên tốc độ lấy mẫu $1\text{ Hz}$, tạo ra sự cân bằng hoàn hảo giữa tính nhạy cảnh báo và khả năng triệt tiêu nhiễu Gauss.
Data và phân tích
- Kịch bản thực nghiệm tĩnh (Static Tests): Kiểm thử trên tập dữ liệu GPS L1 C/A liên tục với cấu hình vệ tinh từ 8 đến 9 kênh, thiết lập tấn công giả mạo 3 vệ tinh hoặc 5 vệ tinh (vệ tinh tham chiếu thật hoặc bị giả mạo).
- Kịch bản thực nghiệm động (Dynamic Flight Trajectories):
- TRJ1: Mô phỏng máy bay bay thẳng gia tốc đều với vận tốc biến thiên, kiểm tra độ ổn định của vi phân kép dưới tác động của dịch chuyển Doppler động học lớn.
- TRJ2: Mô phỏng máy bay thực hiện động tác lượn vòng bán kính hẹp (High-G Turn), thay đổi liên tục góc hướng ăng-ten để kiểm tra tính vững của ma trận phương sai.
- Phương pháp phân tích nâng cao:
- Mô hình hồi quy tuyến tính LR-D3: Ước lượng độ dốc $\hat{a}_i$ và hệ số chặn $\hat{b}_i$ của chuỗi $\mu_i(t)$ qua xấp xỉ đoạn thẳng. So sánh đồng thời cả hai tham số $(\hat{a}_j - \hat{a}_k)^2 \le \xi_a^2$ và $(\hat{b}_j - \hat{b}_k)^2 \le \xi_b^2$.
- Mô hình hỗn hợp Gauss (GMM) 2D/3D: Tối ưu hóa hàm hợp lý Log-Likelihood qua 38 vòng lặp thuật toán EM (Expectation-Maximization), phân cụm chính xác các thành phần Gaussian tương ứng với các hướng phát giả mạo khác nhau.
Phát hiện đột phá và implications
Những phát hiện then chốt
+----------------------------------------------------------------------------------------------------+
| BẢNG TỔNG HỢP CÁC PHÁT HIỆN THEN CHỐT |
+-------------------+--------------------------------+-------------------------------+---------------+
| TÌNH HUỐNG / THỬ | HIỆN TRẠNG Y VĂN QUỐC TẾ | KẾT QUẢ THỰC NGHIỆM ĐỘT PHÁ | MỨC ĐỘ CẢI |
| NGHIỆM ĐỐI CHỨNG | (Thuật toán SoS gốc) | CỦA LUẬN ÁN (D3 & LR-D3) | THIỆN |
+-------------------+--------------------------------+-------------------------------+---------------+
| Bám hỗn hợp | Bỏ sót 100% cuộc tấn công | Bắt giữ chính xác 100% | Vượt trội |
| (Mixed Tracking) | (\Lambda_{SoS} vượt ngưỡng H1) | (Tập con \ge 3 vệ tinh bị gom)| hoàn toàn |
+-------------------+--------------------------------+-------------------------------+---------------+
| Trượt chu kỳ | Báo động giả liên tục | Triệt tiêu toàn bộ báo động | Khôi phục |
| (Cycle Slips) | khi có trượt chu kỳ pha | giả nhờ Doppler Monitor 10Hz | độ tin cậy |
+-------------------+--------------------------------+-------------------------------+---------------+
| Giao cắt quỹ đạo | Báo động giả cục bộ khi quỹ | Triệt tiêu hoàn toàn: | P_fa = 0% |
| vi phân pha 2 giờ | đạo DD thật cắt DD giả mạo | P_fa = 0.00%, P_md = 0.00% | P_md = 0% |
+-------------------+--------------------------------+-------------------------------+---------------+
| Tấn công tinh vi | Thất bại hoàn toàn trước | Phân cụm tự động qua GMM-EM | Tự động hóa |
| đa hướng | tín hiệu đa góc tới | sau 38 vòng lặp tối ưu | không cần tune|
+-------------------+--------------------------------+-------------------------------+---------------+
- Giải quyết triệt để bài toán bám hỗn hợp (Mixed Tracking):
Trong kịch bản tấn công trung cấp với 5 trên tổng số 8 vệ tinh bị giả mạo (hoặc 3 trên 9 vệ tinh), thuật toán SoS truyền thống hoàn toàn thất bại vì chỉ số tổng $\Lambda_{\text{SoS}}$ bị chi phối bởi các vệ tinh thật, dẫn đến việc kết luận sai lầm là hệ thống an toàn ($H_1$). Ngược lại, thuật toán D3 đã gom cụm chính xác toàn bộ 5 vệ tinh giả mạo vào vùng tương đồng $\Sigma_{cnt}^o$ với khoảng cách vi phân $|\mu_i - \mu_k|^2 \approx 0$, phát hiện cuộc tấn công với độ tin cậy tuyệt đối.
- Triệt tiêu ảnh hưởng của Cycle Slips bằng Doppler Monitor:
Khi xuất hiện hiện tượng trượt chu kỳ, chỉ số SoS lập tức tăng vọt hàng chục đơn vị gây báo động giả nghiêm trọng. Với bộ giám sát Doppler Shift tích hợp ngưỡng $10\text{ Hz}$, toàn bộ các bước nhảy pha dị thường bị phát hiện và loại bỏ trong đúng 1 epoch, khôi phục tính chính xác của chỉ số kiểm định.
- Loại bỏ hiện tượng giao cắt quỹ đạo bằng thuật toán LR-D3:
Trong các thử nghiệm mô phỏng kéo dài 2 giờ (Bảng 3.13), thuật toán D3 tiêu chuẩn đôi khi xuất hiện cảnh báo sai khi quỹ đạo vi phân pha của một vệ tinh thật (ví dụ PRN 23 hoặc PRN 25) tình cờ đi cắt ngang qua mức vi phân của các vệ tinh giả mạo. Thuật toán cải tiến LR-D3, thông qua việc đánh giá đồng thời độ dốc (Slope) và hệ số chặn (Intercept), đã chứng minh: Không có bất kỳ trường hợp báo động giả ($P_{fa} = 0%$) hay bỏ sót cảnh báo ($P_{md} = 0%$) nào xảy ra trong suốt 2 giờ thử nghiệm liên tục.
- Tự động hóa phát hiện giả mạo đa hướng bằng mô hình GMM:
Trong kịch bản tấn công tinh vi phát từ nhiều hướng khác nhau, mô hình GMM huấn luyện tự động qua thuật toán EM đã hội tụ chỉ sau 38 chu kỳ lặp, phân tách hoàn toàn các cụm Gauss của tín hiệu giả mạo mà không cần thiết lập ngưỡng tham số bằng tay.
Implications đa chiều
- Về mặt lý thuyết: Luận án hoàn thiện khung toán học về phân tích vi phân pha trên không gian đo lường phân tán, mở ra hướng tiếp cận mới trong việc kết hợp lý thuyết lọc không gian (Spatial Filtering) với học máy xác suất (Probabilistic Machine Learning) trong định vị an toàn.
- Về mặt phương pháp luận: Giải pháp sinh dữ liệu GNSS đa hướng trên nền tảng SDR phá vỡ sự phụ thuộc vào các hệ thống máy tạo tín hiệu phần cứng đa kênh cực kỳ đắt đỏ (như Spirent, IFEN, SkyDel), cung cấp quy trình nghiên cứu chi phí thấp nhưng có độ chính xác vật lý cao cho cộng đồng học thuật.
- Về mặt thực tiễn và công nghiệp: Thuật toán D3 và LR-D3 có độ phức tạp tính toán chỉ là $O(I^2)$ với $I$ là số lượng vệ tinh ($I \le 30$). Với tốc độ cập nhật $1-10\text{ Hz}$, thuật toán có thể nhúng trực tiếp vào firmware của các vi xử lý nhúng ARM Cortex-M hoặc DSP trên các module GPS thương mại hiện hữu (như U-blox F9P, Septentrio) mà không cần thay đổi phần cứng máy thu.
- Về mặt chính sách và an ninh quốc gia: Cung cấp cơ sở khoa học để các cơ quan quản lý hàng không (ICAO, FAA, EASA) và an toàn hàng hải xây dựng tiêu chuẩn kỹ thuật bắt buộc về khả năng tự phục hồi và cảnh báo giả mạo cho các phương tiện không người lái và cơ sở hạ tầng trọng yếu.
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 nghiên cứu:
- Ràng buộc về đường cơ sở vật lý (Baseline Constraint): Thuật toán yêu cầu khoảng cách giữa hai ăng-ten cố định ở mức xấp xỉ $2\text{ m}$. Cấu hình này rất phù hợp cho máy bay, tàu thủy, xe ô tô nhưng gặp khó khăn khi triển khai trên các thiết bị cầm tay hoặc drone siêu nhỏ (micro-UAV).
- Ảnh hưởng của môi trường đa đường (Multipath) đô thị khắc nghiệt: Trong các hẻm vực đô thị (Urban Canyons) với nhiều bề mặt phản xạ kim loại/kính cường lực, sai số đa đường $\Delta\nabla\varepsilon_i$ có thể vượt quá phân phối Gauss danh định, làm tăng độ giãn của vùng tương đồng $\Sigma_j^o$.
- Giới hạn đồng bộ trong mô phỏng đa hướng: Việc phát sóng giả mạo thực tế từ nhiều ăng-ten phân tán đòi hỏi hệ thống máy phát phải đồng bộ thời gian cấp pico-giây; nghiên cứu hiện tại mới kiểm chứng kịch bản này trên bộ mô phỏng tín hiệu phần mềm SDR.
Chương trình nghiên cứu tiếp nối trong giai đoạn 5 - 10 năm:
- Mở rộng 1 (Ăng-ten ảo trên phương tiện di động): Tận dụng chuyển động động học của phương tiện với một ăng-ten duy nhất để tạo ra đường cơ sở tổng hợp (Synthetic Baseline / Virtual Dual-Antenna).
- Mở rộng 2 (Học sâu thích nghi đa tần): Kết hợp mạng nơ-ron đồ thị (Graph Neural Networks - GNN) để mô hình hóa mối tương quan không gian - thời gian giữa các tần số L1, L2, L5 và nhiều chùm vệ tinh (GPS, Galileo, BeiDou).
- Mở rộng 3 (Tích hợp vệ tinh quỹ đạo thấp LEO PNT): Mở rộng thuật toán D3 để giám sát tính toàn vẹn cho các chùm vệ tinh LEO PNT thế hệ mới (như Starlink, Xona Space), nơi tốc độ dịch chuyển Doppler cao gấp hàng chục lần vệ tinh MEO.
Tác động và ảnh hưởng
+----------------------------------------------------------------------------------------------------+
| BẢN ĐỒ TÁC ĐỘNG VÀ ẢNH HƯỞNG |
+----------------------------------------------------------------------------------------------------+
| [HỌC THUẬT & KHOA HỌC] |
| - Công bố quốc tế: IEEE Transactions on Aerospace and Electronic Systems (2021), NAVITEC (2018). |
| - Giải quyết lỗ hổng lý thuyết 10 năm của bài toán bám hỗn hợp (Mixed Tracking). |
+----------------------------------------------------------------------------------------------------+
| [CÔNG NGHIỆP ĐỊNH VỊ & TỰ HÀNH] |
| - Tích hợp trực tiếp vào module OEM dân dụng (Septentrio, NovAtel, U-blox). |
| - Bảo vệ đội bay UAV giao hàng, taxi bay (eVTOL) và tàu container tự hành trước cướp quyền điều khiển|
+----------------------------------------------------------------------------------------------------+
| [HẠ TẦNG XÃ HỘI VÀ AN NINH QUỐC GIA] |
| - Đảm bảo an ninh đồng bộ thời gian cho mạng viễn thông 5G/6G và lưới điện thông minh. |
| - Ngăn chặn thiệt hại hàng tỷ USD do thao túng vị trí trong giao dịch tài chính tự động. |
+----------------------------------------------------------------------------------------------------+
- Tác động học thuật: Các bài báo khoa học xuất bản từ luận án trên tạp chí uy tín hàng đầu như IEEE Transactions on Aerospace and Electronic Systems (2021) và hội thảo ESA NAVITEC (Hà Lan) đã tạo lập một chuẩn mực mới trong việc đánh giá an ninh máy thu GNSS hai ăng-ten, mở ra hàng loạt trích dẫn và nghiên cứu tiếp nối trên toàn cầu.
- Tác động công nghiệp: Cung cấp cho các tập đoàn sản xuất thiết bị định vị (như Septentrio, Trimble, U-blox) một thuật toán phòng thủ nhẹ, có thể thương mại hóa ngay lập tức mà không làm tăng chi phí phần cứng (BOM cost).
- Tác động xã hội: Bảo vệ sự an toàn của tính mạng con người trong giao thông thông minh và vận tải hàng không dân dụng; ngăn chặn triệt để các hành vi phá hoại hệ thống định vị nhằm mục đích khủng bố hoặc trục lợi kinh tế.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật (Academic Researchers): Tiếp cận một khung phân tích toán học mẫu mực về xử lý tín hiệu vi phân sóng mang, lý thuyết ước lượng thống kê GLRT và mô hình học máy phân cụm GMM-EM ứng dụng trong bảo mật vô tuyến.
- Kỹ sư R&D trong ngành công nghiệp GNSS và Phương tiện tự hành: Sở hữu mã nguồn thuật toán tối ưu, dễ dàng lập trình trên các vi điều khiển nhúng thời gian thực để nâng cấp dòng sản phẩm máy thu an toàn cao cấp.
- Các nhà hoạch định chính sách và Quản lý Hàng không/Hàng hải: Có số liệu thực chứng cụ thể để ban hành các quy chuẩn kiểm định an toàn tín hiệu định vị cho hạ tầng giao thông thông minh quốc gia.
- Doanh nghiệp vận tải và Tài chính viễn thông: Giảm thiểu rủi ro gián đoạn dịch vụ, bảo vệ tài sản và tính liên tục của các giao dịch kinh tế số.
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 lý thuyết độc đáo nhất là việc mở rộng Lý thuyết vi phân pha sóng mang kép ngắn hạn (Short-baseline Carrier Phase Double Difference Theory) kết hợp với Kiểm định tỷ số hợp lý tổng quát (GLRT) để giải quyết bài toán phát hiện tập con giả mạo cục bộ ($S \subset {S \cup A}$). Luận án đã phá vỡ giả định toàn phần của mô hình Sum of Squares (Falletti & Nicola), chứng minh rằng hiệu vi phân bậc hai giữa các vệ tinh giả mạo luôn có kỳ vọng hình học đồng nhất ($M_{cnt}$), và phân phối của đại lượng kiểm định cặp $\Lambda_{D3}$ tuân theo quy luật Chi-square phi tâm với tham số phi tâm $\lambda = (m_j - m_k)^2 / (\sigma_j^2 + \sigma_k^2)$, cho phép tính toán giải tích chính xác xác suất phát hiện lỗi.
2. Điểm đột phá về phương pháp luận so với các nghiên cứu quốc tế trước đây là gì?
So với phương pháp SoS (Falletti & Nicola, 2018/2021) và phương pháp VSD (Psiaki et al., 2013):
- Không yêu cầu can thiệp vào tầng tương quan phần mềm phức tạp như VSD (vốn không thể chạy trên chip thương mại).
- Không yêu cầu đồng bộ xung nhịp giữa hai máy thu và không cần giải bất định số nguyên vòng pha ($N$).
- Bổ sung bộ giám sát Doppler Shift Monitor loại bỏ hoàn toàn hiện tượng cảnh báo sai do Cycle Slips.
- Nâng cấp lên mô hình hồi quy tuyến tính LR-D3, triệt tiêu $100%$ hiện tượng giao cắt quỹ đạo pha vốn là điểm yếu cố hữu của phương pháp D3 ban đầu.
3. Phát hiện thực nghiệm nào bất ngờ và có ý nghĩa thực tiễn lớn nhất?
Phát hiện bất ngờ nhất là hiện tượng giao cắt quỹ đạo vi phân pha thập phân (Fractional DD Crossing) giữa vệ tinh thật và chùm vệ tinh giả. Khi một vệ tinh thật chuyển động trên quỹ đạo và có góc tới thỏa mãn $\beta(\cos(\alpha_{authentic}) - \cos(\alpha_r)) \approx M_{cnt}$, giá trị $\mu_{authentic}$ sẽ tạm thời trùng với $\mu_{spoofed}$, khiến thuật toán D3 thông thường ngộ nhận vệ tinh thật này thuộc chùm giả mạo (gây cảnh báo sai). Bằng việc phát hiện ra rằng đạo hàm bậc nhất (độ dốc - slope) và hệ số chặn (intercept) của chúng hoàn toàn khác nhau do vận tốc góc vệ tinh khác biệt, tác giả đã thiết kế bộ dò LR-D3 đạt hiệu năng tuyệt đối: $P_{fa} = 0%$ và $P_{md} = 0%$ trong toàn bộ 2 giờ thử nghiệm.
4. Luận án có cung cấp quy trình nhân bản thực nghiệm (Replication Protocol) hoàn chỉnh không?
Hoàn toàn có. Luận án cung cấp chi tiết:
- Thông số phần cứng: Cặp máy thu Septentrio AsteRx4 OEM, khoảng cách baseline $D = 2\text{ m}$.
- Quy trình Zero-Baseline để đo chuẩn phương sai nhiễu $\sigma_{D3}$ tương ứng với các mức $C/N_0$ ($32, 38, 42, 47, 52\text{ dB-Hz}$).
- Ngưỡng lọc trượt chu kỳ Doppler Monitor cố định ở mức $10\text{ Hz}$.
- Cửa sổ lọc trung bình trượt thời gian $\eta = 5\text{ giây}$.
- Kiến trúc phần mềm mô phỏng SDR sinh tín hiệu GNSS đa hướng can thiệp pha sóng mang.
5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
Lộ trình 10 năm tập trung vào 3 mũi nhọn: (1) Phát triển thuật toán ăng-ten ảo (Synthetic Virtual Antenna) ứng dụng chuyển động học của xe/drone để loại bỏ ăng-ten thứ hai; (2) Tích hợp mạng nơ-ron học sâu thích nghi đa tần L1/L2/L5 và đa chùm vệ tinh (GPS/Galileo/BDS); (3) Mở rộng hệ thống bảo mật không gian cho các chùm vệ tinh định vị quỹ đạo thấp (LEO PNT Megaconstellations).
Kết luận
Luận án tiến sĩ của tác giả Nguyễn Văn Hiên đã giải quyết trọn vẹn và xuất sắc bài toán phát hiện, cô lập tín hiệu giả mạo GNSS, đóng góp cho khoa học và thực tiễn 5 giá trị cốt lõi:
- Phát triển thành công thuật toán D3 (Dispersion of Double Differences): Giải quyết dứt điểm điểm nghẽn 10 năm của bài toán bám hỗn hợp (Mixed Tracking), cho phép phát hiện chính xác tập con tín hiệu giả mạo từ một nguồn phát chung mà không cần đồng bộ xung nhịp máy thu.
- Thiết lập cơ chế triệt tiêu trượt chu kỳ bằng Doppler Shift Monitor: Đưa ra giải pháp loại bỏ xung đột pha sóng mang tức thời với ngưỡng $10\text{ Hz}$, đảm bảo hệ thống không bị kích hoạt báo động giả trong môi trường động học phức tạp.
- Sáng tạo thuật toán hồi quy tuyến tính LR-D3: Loại bỏ triệt để hiện tượng giao cắt quỹ đạo hình học giữa vệ tinh thật và vệ tinh giả, đạt hiệu suất lý tưởng $P_{fa} = 0%$ và $P_{md} = 0%$ trên các tập dữ liệu thử nghiệm động/tĩnh 2 giờ liên tục.
- Ứng dụng thành công mô hình học máy GMM và thuật toán EM: Tự động hóa quá trình học phân phối xác suất góc tới, giải quyết bài toán phát hiện tấn công giả mạo tinh vi đa hướng mà không cần tinh chỉnh tham số thủ công.
- Đề xuất phương pháp sinh dữ liệu GNSS đa hướng trên nền tảng SDR: Mở ra giải pháp kiểm chuẩn thực nghiệm kinh tế, độ chính xác cao, thay thế các hệ thống máy tạo tín hiệu phần cứng triệu USD.
Công trình tạo nên bước chuyển dịch quan trọng (Paradigm Shift) trong lĩnh vực an ninh định vị vô vệ tinh: Chuyển từ các hệ thống mảng ăng-ten cồng kềnh, đắt đỏ sang các giải pháp thuật toán thông minh, tinh gọn, có thể thương mại hóa diện rộng trên các thiết bị định vị dân dụng, góp phần bảo vệ an toàn cho các hạ tầng thông tin trọng yếu quốc gia và quốc tế.
