Tổng quan về luận án

Nghiên cứu của tác giả Ngô Duy Tân (2018) tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, dưới sự hướng dẫn của PGS. TS Thái Quang Vinh và TS. Bùi Trọng Tuyên, giải quyết bài toán cốt lõi trong công nghệ vũ trụ: "Ước lượng tư thế vệ tinh nhỏ quan sát Trái đất bằng việc hợp nhất hóa dữ liệu của cảm biến tốc độ góc và cảm biến sao". Trong bối cảnh ngành công nghệ vũ trụ Việt Nam tiếp cận và làm chủ công nghệ vệ tinh sau thành công của dự án VNREDSat-1 (Vietnam Natural Resources, Environment and Disaster Monitoring Satellite), phân hệ xác định và điều khiển tư thế (Attitude Determination and Control Subsystem - ADCS) được định vị là phân hệ phức tạp, tinh vi nhất, quyết định trực tiếp tới chất lượng quang học và độ chính xác hình học của ảnh chụp bề mặt Trái đất.

Khoảng trống nghiên cứu (research gap) trọng tâm xuất phát từ sự mâu thuẫn giữa yêu cầu định hướng chính xác cao ở chế độ chụp ảnh linh hoạt và các ràng buộc vật lý khắc nghiệt trên dòng vệ tinh nhỏ (khối lượng dưới 500 kg, đặc biệt là phân khúc khoảng 100 kg hoạt động ở quỹ đạo thấp LEO dưới 1000 km, chu kỳ quay quỹ đạo $T_o \approx 90$ phút). Các bộ xử lý trên vệ tinh như dòng máy tính OBC750 của hãng Surrey Satellite Technology (SST) chỉ sở hữu năng lực tính toán khiêm tốn: "Bộ vi xử lý : IBM PPC750FL, 2 MIPS; Hệ điều hành : RTEM (hệ điều hành thời gian thực); Bộ nhớ trong : 256MB; Công suất : 10W". Hơn nữa, môi trường không gian gây ra hiện tượng đảo bit dữ liệu (Single Event Upset - SEU) do bức xạ vũ trụ và biến đổi nhiệt độ cực đoan. Khi cảm biến quang học gặp điểm mù (bị lóa sáng do Mặt Trời, Mặt Trăng hoặc Trái Đất đi vào trường nhìn) hoặc cảm biến quán tính bị trôi dạt dữ liệu, các thuật toán lọc Kalman phi tuyến truyền thống dễ bị phân kỳ hoặc mất ổn định.

Luận án thiết lập hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) chặt chẽ:

  • RQ1: Làm thế nào để hợp nhất tối ưu dữ liệu từ cảm biến tốc độ góc (gyroscope) và cảm biến sao (star tracker - SST) khi có sự suy giảm chất lượng đo lường và trôi dạt độ trượt (bias drift)?
  • RQ2: Cơ chế thích nghi nào có thể tự động tái cấu hình hệ số lọc theo thời gian thực mà vẫn tối ưu hóa tài nguyên tính toán trên chip FPGA/SoC?
  • RQ3: Phương pháp biểu diễn hình học mới Pivot vector có thể ứng dụng ra sao trong mô hình hóa động thái quay phi suy biến của vệ tinh?
  • H1: Việc bổ sung biến trạng thái độ trượt cảm biến tốc độ góc $\beta$ vào bộ lọc Kalman mở rộng (Extended Kalman Filter - EKF) sẽ triệt tiêu sai số tích lũy dạng bước đi ngẫu nhiên (random walk) khi cảm biến sao mất tín hiệu.
  • H2: Việc tích hợp bộ suy luận logic mờ (Fuzzy Logic) để giám sát và hiệu chỉnh ma trận hiệp biến nhiễu theo thời gian thực sẽ tạo ra khả năng kháng lỗi (fault-tolerant) vượt trội trước các đột biến đo lường.
  • H3: Cấu trúc biểu diễn Pivot vector của Russell P. Patera cho phép biểu diễn các phép quay không gian liên tiếp mà không gặp điểm kỳ dị toán học (gimbal lock) như góc Euler.

Khung lý thuyết của công trình tích hợp lý thuyết động lực học vật rắn Euler-Newton, đại số Quaternion 4 thành phần $q = [v^T, q_4]^T$, bài toán tối ưu hóa Wahba, lý thuyết lọc ngẫu nhiên Kalman/EKF, và lý thuyết điều khiển mờ Mamdani/Takagi-Sugeno. Phạm vi nghiên cứu bao quát chu trình khép kín của phân hệ ADCS với vòng điều khiển phản hồi PID và cơ cấu chấp hành cụm 4 bánh xe động lượng (reaction wheels), mang lại đóng góp đột phá về khả năng tự chủ công nghệ điều khiển bay cho vệ tinh viễn thám thế hệ mới.

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự tiến hóa mạnh mẽ của các phương pháp xác định tư thế vệ tinh qua nhiều thập kỷ. Trường phái xác định tư thế tiền định (deterministic methods) bắt đầu từ bài toán tối ưu hóa hàm mục tiêu bình phương tối thiểu $J(A) = \frac{1}{2} \sum a_i |b_i - A r_i|^2$ do Grace Wahba đề xuất năm 1965. Lời giải cho bài toán này được hiện thực hóa qua thuật toán TRIAD (Black, 1964) và bước nhảy vọt thuật toán ước lượng Quaternion - QUEST (Quaternion Estimator) của Shuster và Oh (1981). Dù QUEST tối ưu hóa tốc độ tính toán cho các phép đo vector tức thời, phương pháp này không có khả năng lọc nhiễu động học và không tận dụng được lịch sử truyền trạng thái từ cảm biến quán tính.

Trường phái lọc ngẫu nhiên thống kê (stochastic filtering) do Kalman (1960) khởi xướng đã được Lefferts, Markley và Shuster (1982) chuẩn hóa cho bài toán không gian qua cấu trúc Bộ lọc Kalman mở rộng dạng nhân (Multiplicative EKF - MEKF). MEKF sử dụng quaternion 4 thành phần cho việc cập nhật động học toàn cục nhưng biểu diễn sai số góc dưới dạng vector 3 chiều vi phân nhằm tránh suy biến ma trận hiệp biến $P$. Sau đó, Julier và Uhlmann (1997) cùng Crassidis và Markley (2003) phát triển Bộ lọc Kalman không mùi (Unscented Kalman Filter - UKF) để xử lý các hàm phi tuyến cao bậc mà không cần tính ma trận Jacobi.

Tuy nhiên, các tranh luận học thuật sâu sắc nảy sinh quanh sự đánh đổi giữa độ chính xác và tính khả thi phần cứng:

  • Tranh luận 1 (Độ phức tạp tính toán vs. Tài nguyên hạn chế): Nghiên cứu của Si Mohammed và cộng sự (2013) khi đánh giá các thuật toán trên máy tính nhúng vệ tinh OBC750 chỉ ra rằng UKF đòi hỏi năng lực tính toán quá lớn và tiêu tốn bộ nhớ gấp nhiều lần EKF, khiến UKF khó đáp ứng tính thời gian thực trong chu kỳ điều khiển khắt khe của vệ tinh chụp ảnh.
  • Tranh luận 2 (Tính bền vững trước lỗi vs. Tính tối ưu tĩnh): Các bộ lọc EKF/MEKF cổ điển phụ thuộc hoàn toàn vào ma trận hiệp biến nhiễu quá trình $Q$ và nhiễu đo $R$ cố định. Khi cảm biến sao bị lóa sáng do nguồn quang học ngoại lai hoặc cảm biến con quay bị dịch chuyển hệ số tỉ lệ (scale factor error), EKF truyền thống sẽ suy giảm chất lượng nghiêm trọng hoặc phân kỳ hoàn toàn.

Công trình của Mokhtar Aboelaze và cộng sự (2007) về triển khai điều khiển vệ tinh dùng logic mờ trên chip FPGA Xilinx và các công bố kiến trúc FPGA của Georg Grillmayer (2013) cùng Muhammad Yasir (2011) đã mở ra hướng đi mới: ứng dụng trí tuệ tính toán để giám sát thích nghi. Luận án của Ngô Duy Tân định vị chính xác tại điểm giao thoa này: kết hợp tính tinh gọn, chuẩn xác của cấu trúc MEKF với sự linh hoạt của bộ suy luận logic mờ, tạo nên một hệ thống ước lượng kháng lỗi tự động, trực tiếp giải quyết sự suy giảm đo lường thực tế mà không làm quá tải bộ xử lý trên khoang.

Đó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 điều khiển và ước lượng trạng thái cho các hệ thống hàng không vũ trụ phi tuyến:

  1. Mở rộng mô hình không gian trạng thái bù độ trượt: Tác giả tích hợp phương trình vi phân trôi dạt cảm biến con quay vào vector trạng thái động học, biến đổi phương trình liên tục $\omega = \omega_{true} + \beta + \eta_v$ và $\dot{\beta} = \eta_u$ sang mô hình sai số vi phân 6 chiều $\Delta x = [\delta \theta^T, \Delta \beta^T]^T$. Mô hình này cho phép thuật toán ước lượng đồng thời cả sai lệch góc quay lẫn tốc độ trượt bias ngẫu nhiên theo thời gian thực.
  2. Hình thức hóa phương pháp biểu diễn hình học Pivot vector: Kế thừa cơ sở toán học của Russell P. Patera (2005, 2014), luận án thiết lập ánh xạ hình học biểu diễn phép quay không gian như một cung trên đường tròn lớn của mặt cầu đơn vị thông qua cặp vector trực chuẩn $p_1, p_2$ trên mặt phẳng quay vuông góc với trục Euler $e$. Cấu trúc này mở rộng khả năng tổng hợp các chuỗi quay liên tiếp bằng phép đại số giao tuyến mặt phẳng, loại bỏ hoàn toàn các điểm kỳ dị mà không đòi hỏi chuẩn hóa liên tục như quaternion.
  3. Mô hình hóa toán học các mệnh đề điều khiển:
    • Mệnh đề 1 (Tính hội tụ ước lượng độ trượt): Khi hệ thống quan sát được trạng thái thông qua các bản đo vector của cảm biến sao ($q_{Star}$), sai số ước lượng độ trượt con quay $\tilde{\beta}_k = \beta_k - \hat{\beta}_k$ hội tụ tiệm cận về lân cận gốc tọa độ với ma trận hiệp biến bị chặn.
    • Mệnh đề 2 (Tính ổn định kháng lỗi thích nghi): Bộ hiệu chỉnh logic mờ thay đổi trọng số và ma trận $R_k$ theo hàm độ lệch dư (residual innovation), đảm bảo tính ổn định Lyapunov cho sai số bám tư thế ngay cả khi cảm biến sao ngắt quãng dữ liệu trong khoảng thời gian $\Delta t > 10$ giây.
                  +--------------------------------------------------+
                  |              DỮ LIỆU CẢM BIẾN                    |
                  |  +---------------------+  +-------------------+  |
                  |  | Cảm biến tốc độ góc |  |  Cảm biến sao     |  |
                  |  | (Gyroscope - FOG)   |  |  (Star Tracker)   |  |
                  |  +----------+----------+  +---------+---------+  |
                  +-------------|-----------------------|------------+
                                | (\omega)              | (q_Star)
                                v                       v
                  +--------------------------------------------------+
                  |      KHUNG PHÂN TÍCH ƯỚC LƯỢNG VÀ THÍCH NGHI     |
                  |                                                  |
                  |   +------------------------------------------+   |
                  |   |    Cơ chế Giám sát & Kháng lỗi Fuzzy    |   |
                  |   |    (Đánh giá độ lệch dư Innovation)      |   |
                  |   +--------------------+---------------------+   |
                  |                        |                         |
                  |                        v (\alpha, R_k)           |
                  |   +------------------------------------------+   |
                  |   |  Bộ lọc Kalman Mở rộng Bù trượt (MEKF)   |   |
                  |   |  Trạng thái: \Delta x = [\delta \theta, \Delta \beta]^T |   |
                  |   +--------------------+---------------------+   |
                  +------------------------|-------------------------+
                                           | q_Est, \omega_Est
                                           v
                  +--------------------------------------------------+
                  |         BỘ ĐIỀU KHIỂN TƯ THẾ PHẢN HỒI PID        |
                  |    u = -K_p * q_e - K_d * \omega_e               |
                  +------------------------+-------------------------+
                                           | Moment điều khiển (N_c)
                                           v
                  +--------------------------------------------------+
                  |         CƠ CẤU CHẤP HÀNH & ĐỘNG LỰC HỌC          |
                  |    Cụm 4 Bánh xe động lượng -> Thân Vệ tinh      |
                  +--------------------------------------------------+

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

Khung phân tích của luận án tích hợp đa tầng giữa động học, động lực học và các phép biến đổi không gian:

  • Tầng quy chiếu không gian: Thiết lập mối quan hệ hình học chuyển đổi giữa Hệ tọa độ quán tính quy tâm Trái Đất (ECI - $i$), Hệ tọa độ quỹ đạo (Orbit frame - $o$, với vận tốc góc quay quỹ đạo $\omega_{oi}^o = [0, -\omega_o, 0]^T$), và Hệ tọa độ liên kết thân vệ tinh (Body frame - $b$).
  • Tầng động lực học phi tuyến khép kín: Phương trình động lực học tổng quát tính đến mô men động lượng của cụm bánh xe $h_w$ và các nguồn mô men nhiễu môi trường: $$\dot{\omega} = -I_s^{-1} S(\omega) I_s \omega - I_s^{-1} S(\omega) h_w + I_s^{-1} N_c + I_s^{-1} N_e$$ Trong đó $S(\omega)$ là ma trận phản đối xứng của vector vận tốc góc, $I_s$ là ma trận quán tính $3 \times 3$ của vệ tinh, $N_c$ là mô men điều khiển và $N_e$ là tổng các ngoại lực nhiễu (nhiễu gradient trọng trường $\tau_g$, áp lực bức xạ Mặt Trời, khí quyển tầng cao).

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết lý thực chứng suy diễn (positivism/deductive engineering approach). Hệ thống phân hệ ADCS được thiết kế theo cấu trúc đa tầng mô phỏng số học chi tiết nhằm phản ánh chính xác các điều kiện hoạt động khắc nghiệt trên quỹ đạo đồng bộ Mặt Trời (Sun Synchronous Orbit - SSO).

Thiết kế nghiên cứu bao gồm 3 chế độ vận hành đặc thù:

  1. Chế độ duy trì hướng Trái Đất (Nadir-pointing mode): Trục $z_b$ bám sát tâm Trái Đất với tốc độ quay danh định $\omega_0 = 2\pi / T_o \approx 0.00116$ rad/s quanh trục pitch ($y_b$).
  2. Chế độ chụp ảnh bề mặt (Imaging scan mode): Vệ tinh thực hiện cơ động góc lớn, thay đổi góc nghiêng (roll) và góc chúc ngóc (pitch) theo quỹ đạo quét 3 giai đoạn: tăng gia tốc góc, duy trì vận tốc góc chụp ảnh ổn định, và giảm tốc hội tụ.
  3. Chế độ sự cố cảm biến (Sensor failure/occlusion mode): Cắt tín hiệu hoặc bơm nhiễu xung vào cảm biến sao trong khoảng thời gian xác định ($5-10$ giây) để kiểm thử độ bền vững thuật toán.

Quy trình nghiên cứu rigorous

Mô hình hóa toán học các cảm biến đầu vào được thực hiện theo tiêu chuẩn hàng không vũ trụ quốc tế:

Loại Cảm biến Bản chất vật lý Đặc tính toán học & Mô hình sai số Ưu điểm & Ràng buộc vận hành
Cảm biến Tốc độ góc (Gyroscope / FOG) Cảm biến quán tính (quang sợi / cơ học) $\omega_k = \omega_{true, k} + \beta_k + \eta_{v, k}$
$\beta_{k+1} = \beta_k + \eta_{u, k}$
$E[\eta_v \eta_v^T] = \sigma_v^2 I_{3\times 3}, E[\eta_u \eta_u^T] = \sigma_u^2 I_{3\times 3}$		 Ưu điểm: Đo liên tục, tần số cao, chống nhiễu tốt.
Ràng buộc: Tồn tại độ trôi dạt (bias) và nhiễu ngẫu nhiên random walk cộng dồn theo thời gian.
Cảm biến Sao (Star Tracker - SST) Cảm biến quang học nhận dạng bản đồ sao $q_s = q \otimes \delta q$
$E[\delta q] = 0, E[\delta q \delta q^T] = \sigma_s^2 I_{4\times 4}$		 Ưu điểm: Độ chính xác cực cao (vài arcsecond), không trôi dạt sai số.
Ràng buộc: Tần số thấp, dễ mù quang học khi bị nguồn sáng ngoại lai chiếu thẳng.

Khái niệm hợp nhất dữ liệu được định nghĩa chuẩn xác: "Hợp nhất dữ liệu đa cảm biến là quá trình kết hợp dữ liệu từ nhiều cảm biến khác nhau với mục đích cho ra kết quả cuối cùng tốt hơn khi sử dụng từng cảm biến độc lập". Trong đó, bản chất của sai lệch con quay được xác định rõ: "Bias là giá trị đo được khi mà chuyển động quay quán tính thực tế là 0o/s. Sai lệch này có một thành phần tương đối ổn định, chỉ thay đổi trong thời gian dài hoặc trong các quá trình tắt/bật".

Quy trình hợp nhất tuyến tính cơ bản dựa trên trọng số được mô hình hóa: $$q_{Est} = \alpha q_{Gyro} + (1 - \alpha) q_{Star}, \quad \alpha \in [0, 1]$$ Khi chất lượng đo lường biến động, hệ số $\alpha$ được điều chỉnh thích nghi thông qua khối suy luận mờ nhằm tránh suy thoái nghiệm.

Data và phân tích

Toàn bộ thuật toán được mô phỏng và kiểm chứng trên phần mềm MATLAB/Simulink với bộ thông số vệ tinh viễn thám thực tế (tương thích cấu hình vệ tinh VNREDSat-1):

  • Ma trận mô men quán tính vệ tinh $I_s$: $$I_s = \begin{bmatrix} 15.0 & 0.2 & 0.1 \ 0.2 & 17.5 & 0.15 \ 0.1 & 0.15 & 20.0 \end{bmatrix} \text{ kg}\cdot\text{m}^2$$
  • Cơ cấu chấp hành bánh xe động lượng: Cụm 4 bánh xe bố trí dạng hình chóp (pyramid configuration) với ma trận phân phối lực $A_w$ kích thước $3 \times 4$.
  • Luật điều khiển phản hồi: Bộ điều khiển vi tích phân tỉ lệ dạng quaternion: $$u = -K_p q_e - K_d \omega_e$$ Với $q_e = q_s^{-1} \otimes q_r$ là quaternion sai lệch giữa tư thế ước lượng $q_s$ và tư thế đặt mong muốn $q_r$.
  • Bộ hiệu chỉnh Fuzzy Logic: Sử dụng 2 biến đầu vào là độ lệch dư trung bình (innovation residual) và đạo hàm độ lệch dư, ánh xạ qua các hàm thuộc tính tam giác/hình thang để sinh ra hệ số hiệu chỉnh ma trận hiệp biến đo $R_k$.

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

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

  1. Sự sụp đổ độ chính xác của phương pháp trọng số tĩnh: Mô phỏng thực nghiệm chứng minh khi cảm biến sao bị mất tín hiệu trong khoảng thời gian từ $5$ đến $10$ giây, phương pháp trọng số tĩnh tích phân trực tiếp từ cảm biến con quay làm sai số trỏ hướng bùng nổ theo hàm phi tuyến do tích lũy độ trượt bias. Sai số góc lệch vượt ngưỡng cho phép của thiết bị chụp ảnh viễn thám ($> 0.5^\circ$).
  2. Khả năng triệt tiêu độ trôi của bộ lọc EKF bù trượt: Khi triển khai bộ lọc EKF với vector trạng thái mở rộng chứa biến $\beta$, sai số ước lượng độ trượt nhanh chóng hội tụ về giá trị thực sau chưa đầy $15$ giây từ khi có tín hiệu cảm biến sao. Trong giai đoạn cảm biến sao bị mất dấu $5-10$s, bộ lọc duy trì sai số trỏ hướng dưới $0.02^\circ$, cải thiện độ chính xác hơn $25$ lần so với phương pháp trọng số.
  3. Hiệu năng vượt trội của bộ lọc Fuzzy-EKF kháng lỗi: Trong kịch bản cảm biến sao bị nhiễu đột biến cường độ cao (giả lập hiện tượng phản xạ ánh sáng Mặt Trời), bộ giám sát mờ đã tự động tăng ma trận hiệp biến $R_k$, giảm độ lợi Kalman $K_k$, chuyển quyền ưu tiên sang mô hình lan truyền quán tính. Nhờ đó, hệ thống loại bỏ hoàn toàn các đỉnh gai sai số (error spikes), giữ sai số xác định tư thế luôn nằm trong dải ổn định $0.005^\circ - 0.01^\circ$.
  4. Tính khả thi của mô hình Pivot vector: Chứng minh thành công rằng cấu trúc biểu diễn 4 tham số độc lập của cặp vector Pivot $p_1, p_2$ thực hiện các phép quay liên tiếp hoàn toàn tương đương với tích quaternion nhưng giảm bớt các thao tác lượng giác phức tạp, mở ra khả năng tiết kiệm chu kỳ lệnh trên chip xử lý chuyên dụng.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp một phương pháp luận hoàn chỉnh về thiết kế bộ lọc ước lượng phi tuyến thích nghi kết hợp giữa giải thuật ngẫu nhiên và trí tuệ nhân tạo mờ, đặt nền tảng cho việc nghiên cứu các hệ thống tự hành không gian sâu.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn để kiểm thử, mô phỏng và đánh giá phân hệ ADCS khép kín từ mô hình cảm biến, cơ cấu chấp hành đến động lực học vật rắn, có thể tái sử dụng cho nhiều lớp vệ tinh khác nhau.
  • Về mặt thực tiễn công nghiệp: Đề xuất giải pháp thuật toán nhẹ, tối ưu tài nguyên tính toán, cho phép nạp trực tiếp vào các dòng chip FPGA hoặc máy tính nhúng công suất thấp ($10\text{W}$) như OBC750 mà không cần nâng cấp phần cứng đắt đỏ.
  • Về mặt chính sách: Đóng góp trực tiếp vào mục tiêu tự chủ công nghệ vũ trụ quốc gia, giảm thiểu sự phụ thuộc vào các gói phần mềm điều khiển thương mại nước ngoài trong các chương trình vệ tinh tiếp nối của Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, nghiên cứu thẳng thắn chỉ ra các giới hạn kỹ thuật:

  1. Giả thiết vật rắn tuyệt đối: Luận án giả định thân vệ tinh là khối cứng lý tưởng, bỏ qua dao động đàn hồi của các tấm pin Mặt Trời khi cơ động nhanh và sự dịch chuyển trọng tâm do chất lỏng tiêu hao trong phân hệ đẩy (thruster sloshing).
  2. Tuyến tính hóa quanh điểm làm việc danh định: Mô hình trạng thái phục vụ phân tích tuyến tính được xấp xỉ quanh tốc độ quay hướng Trái Đất $[0, -\omega_o, 0]^T$, có thể làm giảm nhẹ độ chính xác trong các pha cơ động phi tuyến cực lớn.
  3. Môi trường thử nghiệm dừng ở mức mô phỏng số: Các thuật toán mới được kiểm chứng thông qua mô phỏng phần mềm MATLAB/Simulink, chưa thực hiện đo kiểm thực tế trên bàn quay 3 bậc tự do mô phỏng không trọng lực (Hardware-in-the-Loop Simulation - HILS).

Chương trình nghiên cứu tương lai định hướng các nội dung then chốt:

  • Mở rộng mô hình sang hệ vật rắn biến dạng mềm (flexible multibody dynamics) tích hợp hiệu ứng dao động tấm pin quang điện.
  • Ứng dụng cấu trúc lọc Unscented Kalman Filter thích nghi trên phần cứng chip đa nhân System on Chip (SoC) thế hệ mới chống bức xạ.
  • Thực nghiệm kiểm chứng mã nguồn thuật toán nhúng trên bàn thử nghiệm HIL với cảm biến sao quang học thực và hệ thống mô phỏng vòm trời quang học.

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

Nghiên cứu tạo ra tác động sâu rộng trên nhiều bình diện:

  • Học thuật: Định hình hướng nghiên cứu chuyên sâu về điều khiển tự động và ước lượng tối ưu trong ngành kỹ thuật không gian tại Việt Nam; cung cấp tài liệu tham khảo chuẩn mực cho các chương trình đào tạo tiến sĩ, thạc sĩ ngành Tự động hóa và Kỹ thuật Vũ trụ.
  • Công nghiệp hàng không vũ trụ: Cung cấp giải pháp công nghệ trực tiếp cho Trung tâm Vũ trụ Việt Nam (VNSC) và Viện Công nghệ Vũ trụ (STI) thuộc Viện Hàn lâm KHCNVN trong việc thiết kế phân hệ ADCS cho các vệ tinh viễn thám thế hệ tiếp theo lớp Micro/Mini-satellite.
  • Kinh tế - Xã hội: Nâng cao độ chính xác định vị hình ảnh viễn thám, hỗ trợ đắc lực công tác giám sát tài nguyên rừng, theo dõi biến đổi khí hậu, dự báo thiên tai lũ lụt và bảo đảm an ninh quốc phòng toàn vẹn lãnh thổ.
  • Hội nhập quốc tế: Khẳng định năng lực nghiên cứu của các nhà khoa học Việt Nam trong việc làm chủ các bài toán công nghệ cao tiệm cận tiêu chuẩn của các cơ quan vũ trụ lớn như ESA (Châu Âu) và JAXA (Nhật Bản).

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận khung phương pháp luận mẫu mực về mô hình hóa quaternion, xử lý tín hiệu ngẫu nhiên, và thiết kế bộ lọc Kalman mờ thích nghi cho các hệ thống động lực học phức tạp.
  • Kỹ sư thiết kế hệ thống ADCS: Ứng dụng trực tiếp các thuật toán bù trượt con quay và cơ chế kháng lỗi vào phần mềm bay (flight software) của các dự án CubeSat và vệ tinh nhỏ quan sát Trái Đất.
  • Các viện nghiên cứu & Doanh nghiệp công nghệ cao: Sở hữu tài sản trí tuệ và thuật toán lõi để thương mại hóa các giải pháp tích hợp hệ thống nhúng tự hành trên máy bay không người lái (UAV), tàu lặn tự hành (AUV) và tên lửa đẩy.
  • Cơ quan quản lý nhà nước: Có cơ sở khoa học vững chắc để xây dựng lộ trình đầu tư phát triển công nghệ vệ tinh quan sát Trái Đất quốc gia đến năm 2030 và tầm nhìn 2045.

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ì?

Đóng góp độc đáo nhất là việc xây dựng thành công kiến trúc ước lượng trạng thái tích hợp đa cấp: mở rộng lý thuyết MEKF kinh điển của Lefferts-Markley-Shuster bằng cách gắn kèm vector trôi dạt con quay 3 trục $\beta$, đồng thời bổ sung lớp giám sát mờ thích nghi (Fuzzy supervisor). Cấu trúc này phá vỡ hạn chế phụ thuộc vào tham số nhiễu tiên nghiệm cố định của bộ lọc Kalman truyền thống.

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

So với nghiên cứu thuần EKF của Si Mohammed et al. (2013) trên máy tính OBC750 và nghiên cứu điều khiển mờ thuần túy của Mokhtar Aboelaze et al. (2007), luận án đã kết hợp xuất sắc cả hai phương pháp: duy trì sự tinh gọn toán học của EKF để đảm bảo tính thời gian thực trên chip FPGA, đồng thời tận dụng bộ luật mờ phi tham số để tự động điều chỉnh độ lợi lọc khi cảm biến sao gặp sự cố mù quang học, đạt độ bền vững tối đa mà không gây quá tải tài nguyên nhúng.

3. Phát hiện bất ngờ nhất trong quá trình mô phỏng thực nghiệm là gì?

Phát hiện ấn tượng nhất là trong điều kiện cảm biến sao bị gián đoạn hoàn toàn tín hiệu tới $10$ giây, bộ lọc EKF có bù trượt vẫn duy trì sai số xác định tư thế ở mức dưới $0.02^\circ$, trong khi phương pháp hợp nhất trọng số cổ điển lập tức phân kỳ vượt quá $0.5^\circ$. Điều này chứng minh rằng việc ước lượng chính xác động thái trôi dạt $\beta$ có vai trò quyết định hơn nhiều so với việc chỉ tăng tần số lấy mẫu cảm biến con quay.

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

Có. Luận án trình bày chi tiết toàn bộ hệ phương trình vi phân trạng thái, ma trận quán tính $Is$, ma trận phân phối cơ cấu chấp hành $Aw$, mô hình toán học nhiễu Gauss rời rạc của cảm biến, bảng luật hợp suy mờ Mamdani và thông số các bộ điều khiển PID. Mọi nhà nghiên cứu độc lập đều có thể tái lập hoàn toàn kết quả mô phỏng trên các nền tảng tính toán khoa học như MATLAB/Simulink hoặc Python.

5. Lộ trình nghiên cứu phát triển 10 năm được vạch ra như thế nào?

Lộ trình 10 năm bao gồm: (1) Tối ưu hóa mã nguồn nhúng VHDL/Verilog cho chip FPGA kháng bức xạ; (2) Tích hợp thuật toán với hệ thống cảm biến quang học đa đầu đo (Multi-Star Tracker); (3) Tiến hành thử nghiệm Hardware-in-the-Loop trên bàn quay không gian 3 bậc tự do; (4) Thử nghiệm bay thực tế (In-orbit demonstration) trên vệ tinh viễn thám thế hệ mới của Việt Nam.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán ước lượng tư thế vệ tinh nhỏ quan sát Trái Đất trong điều kiện ràng buộc tài nguyên tính toán phần cứng khắt khe bằng phương pháp hợp nhất dữ liệu cảm biến đa nguồn tiên tiến.
  2. Thiết lập thành công mô hình bộ lọc Kalman mở rộng bù độ trượt con quay (MEKF-Bias), giúp triệt tiêu sai số tích lũy dạng bước đi ngẫu nhiên và nâng cao độ chính xác trỏ hướng gấp nhiều lần phương pháp trọng số cổ điển.
  3. Đề xuất và chứng minh tính hiệu quả vượt trội của cơ chế kháng lỗi thích nghi sử dụng logic mờ, đảm bảo hệ thống ADCS vận hành bền vững, ổn định khi cảm biến sao bị mất tín hiệu hoặc nhiễu lóa quang học.
  4. Mở ra hướng tiếp cận hình học mới mẻ trong biểu diễn không gian thông qua phương pháp Pivot vector của Russell P. Patera, tạo tiền đề cho các thuật toán quay không kỳ dị tối ưu phần cứng.
  5. Xây dựng hoàn chỉnh mô hình mô phỏng phân hệ ADCS khép kín với luật điều khiển PID và cụm bánh xe động lượng, tương thích hoàn toàn với các thông số kỹ thuật thực tế của vệ tinh viễn thám Việt Nam.
  6. Mở ra 3 nhánh nghiên cứu chuyên sâu mới: lọc phi tuyến thích nghi trên chip SoC/FPGA, điều khiển tư thế hệ vệ tinh linh hoạt có kết cấu biến dạng, và tích hợp thuật toán định hướng tự hành cho các chòm vệ tinh nhỏ bay theo đội hình (satellite constellation formation flying).