Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu giải pháp nâng cao chất lượng định vị nguồn âm sử dụng nguyên lý TDOA" do nghiên cứu sinh Trần Công Thìn thực hiện dưới sự hướng dẫn khoa học của TS Bùi Ngọc Mỹ và TS Nguyễn Huy Hoàng, bảo vệ tại Viện Khoa học và Công nghệ quân sự năm 2022, thuộc chuyên ngành Kỹ thuật điện tử (mã số: 9 52 02 03).

Tính cấp thiết và khoảng trống nghiên cứu

Kỹ thuật định vị nguồn âm thanh dựa trên nguyên lý khác biệt thời gian đến (Time Difference of Arrival - TDOA) sở hữu nhiều ưu điểm như yêu cầu số lượng cảm biến tương đối nhỏ, độ chính xác cao và thời gian triển khai nhanh, rất phù hợp cho các bài toán giám sát ngoài trời cự ly xa (như trinh sát hỏa lực, cảnh giới xâm nhập, định vị điểm nổ hoặc cảnh báo máy bay không người lái). Tuy nhiên, khi hoạt động trong môi trường tự nhiên ngoài trời với khoảng cách xa, kỹ thuật TDOA bộc lộ các hạn chế kỹ thuật:

  • Tỷ số tín trên tạp (Signal to Noise Ratio - SNR) tại đầu vào cảm biến suy giảm mạnh theo khoảng cách, gây khó khăn cho việc phát hiện sự kiện âm thanh mục tiêu và làm tăng sai số ước tính trễ thời gian.
  • Tín hiệu âm thanh bị ảnh hưởng bởi tạp âm môi trường và các nguồn âm thanh khác đan xen.
  • Vận tốc lan truyền của âm thanh trong không khí không phải là một hằng số cố định mà biến đổi liên tục theo nhiệt độ, độ ẩm và gió, dẫn đến sai lệch lớn trong quá trình giải hệ phương trình định vị hình học nếu giả định vận tốc âm là hằng số.
  • Các công trình và giải pháp kỹ thuật chi tiết phục vụ mục đích quân sự thường được bảo mật, hạn chế công bố rộng rãi.

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

Luận án tập trung xây dựng các thuật toán xử lý tín hiệu nhằm nâng cao chất lượng định vị nguồn âm ngoài trời cự ly xa theo nguyên lý TDOA, cụ thể gồm 4 mục tiêu:

  1. Nghiên cứu, xây dựng thuật toán phát hiện sự kiện âm thanh trong bài toán định vị nguồn âm theo nguyên lý TDOA.
  2. Nghiên cứu, xây dựng thuật toán nâng cao độ chính xác ước lượng khác biệt thời gian tới.
  3. Nghiên cứu, xây dựng thuật toán giải phương trình định vị nguồn âm trong điều kiện vận tốc âm thanh là một biến số.
  4. Nghiên cứu, xây dựng mô hình hệ thống định vị nguồn âm thanh ngoài trời, cự ly định vị xa, theo nguyên lý TDOA.

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

  • Đối tượng nghiên cứu: Các giải pháp kỹ thuật và thuật toán định vị nguồn âm thanh dựa trên nguyên lý TDOA.
  • Phạm vi và giới hạn nghiên cứu:
    • Môi trường định vị ngoài trời, chịu tác động của điều kiện môi trường tự nhiên (nhiệt độ, độ ẩm, gió).
    • Cự ly định vị thực nghiệm và mô phỏng được giới hạn trong phạm vi dưới $500\text{ m}$.
    • Nguồn âm thanh mục tiêu được xác định trước, dữ liệu mô phỏng và thực nghiệm sử dụng tín hiệu tiếng nổ đầu nòng của súng AK47.
    • Tín hiệu truyền từ nguồn âm đến cảm biến chịu ảnh hưởng của tạp âm Gauss.

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

Các hướng nghiên cứu tiền đề

Luận án tổng thuật lịch sử phát triển của các kỹ thuật định vị âm thanh từ thế kỷ 19 đến nay:

  • Cơ sở định vị thính giác: Nghiên cứu của Lord Rayleigh (1876, 1907) xác định tỷ lệ cường độ âm giữa hai tai là cơ sở định hướng, đồng thời chứng minh âm thanh tần số thấp dưới $256\text{ Hz}$ phụ thuộc chủ yếu vào sự chênh lệch pha tương đối. Silvanus Thompson (1878) kiểm chứng thực nghiệm sự lệch pha hai tai. Các công trình nửa đầu thế kỷ 20 của Jeffress, Newman, Rosenzweig, Stevens, Von Hornbostel và Wertheimer (1920) đặt nền móng cho việc xác định sai lệch thời gian đến từ đặc tính pha, tạo tiền đề cho kỹ thuật TDOA hiện đại.
  • Phân loại phương pháp định vị nguồn âm:
    • Dựa trên đặc tính thời gian (Time Based Localization - TBL): Phương pháp thời gian đến (Time of Arrival - ToA), phương pháp khác biệt thời gian đến (TDOA), và phương pháp định hướng phản hồi năng lượng (Steered Response Power - SRP / Global Coherence Field - GCF).
    • Dựa trên đặc tính năng lượng (Energy Based Localization - EBL): Phương pháp tỷ lệ năng lượng (Energy Ratio - ER) và phương pháp xác định hướng đến của tín hiệu (Direction of Arrival - DoA).
  • Các kỹ thuật phát hiện sự kiện âm thanh (SKAT): Ban đầu dựa vào ngưỡng biên độ tín hiệu (phù hợp định vị tiếng ồn đô thị hoặc hỏng hóc máy móc); tiếp đến là kỹ thuật tương quan tín hiệu; gần đây là các mô hình học máy như mô hình hỗn hợp Gauss kết hợp Markov ẩn (GMM-HMM), mạng nơ-ron tích chập (CNN), mạng nơ-ron hồi tiếp (RNN).
  • Kỹ thuật ước tính trễ thời gian (TDOA / TDE):
    • Miền thời gian: Tìm cực trị hàm tương quan chéo kinh điển (CCC), cấu trúc đơn giản nhưng suy giảm độ chính xác khi SNR thấp.
    • Miền tần số: Nhóm nghiên cứu của G. Clifford Carter tại Trung tâm Tác chiến Dưới nước Hải quân Hoa Kỳ (NUWC) đề xuất phương pháp tương quan chéo tổng quát (Generalized Cross Correlation - GCC) với các hàm trọng số như PHAT (Phase Transform), SCOT (Smoothed Coherence Transform).
  • Kỹ thuật giải phương trình vị trí nguồn âm: Phương pháp dự đoán hợp lý cực đại (Maximum Likelihood Estimator - MLE), phương pháp bình phương trung bình tối thiểu (Least Squares Estimator - LSE), phương pháp tuyến tính hóa kết hợp ràng buộc hình học hoặc phối hợp với bản đồ xác suất SRP.

Khoảng trống luận án lựa chọn giải quyết

Phần lớn các nghiên cứu hiện có tập trung vào môi trường trong nhà (âm học phòng, hội nghị truyền hình, tương tác người - máy) hoặc chưa tối ưu hóa cho điều kiện ngoài trời cự ly xa, nơi tạp âm lớn và vận tốc âm thanh thay đổi liên tục theo thời tiết. Luận án lựa chọn giải quyết đồng thời chuỗi bài toán: phát hiện tín hiệu mục tiêu xác định trước khi bị nhiễu che lấp, nâng cao độ chính xác trễ TDOA ở cự ly xa, và giải bài toán vị trí với giả định vận tốc âm thanh là biến số cần ước tính.

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

Khung lý thuyết và mô hình toán học

  • Mô hình kênh truyền: Hệ thống định vị nguồn âm được mô hình hóa dưới dạng hệ SIMO (Single-Input Multiple-Output). Trong không gian tự do ngoài trời, tín hiệu thu tại cảm biến thứ $i$ được biểu diễn: $$x_i(t) = h * s(t) + b_i(t) = h^T s(t) + b_i(t)$$ Với mô hình lan truyền trực tiếp (suy hao đa đường nhỏ), phương trình suy biến thành: $$x_i(t) = \alpha_i s(t - \tau_i) + b_i(t), \quad i = 1, 2, \dots, K$$ Trong đó $\alpha_i \in [0, 1]$ là hệ số suy giảm công suất, $\tau_i$ là thời gian lan truyền từ nguồn âm $x_s = [x_s, y_s, z_s]^T$ tới cảm biến thứ $i$ tại vị trí $m_i = [x_i, y_i, z_i]^T$, và $b_i(t)$ là tạp âm Gauss.
  • Quan hệ hình học TDOA: Khác biệt thời gian đến giữa hai cảm biến $i$ và $j$ là $\tau_{ij} = \tau_i - \tau_j$, quan hệ với khoảng cách hình học qua phương trình hypeboloit: $$\tau_{ij} = \frac{|x_s - m_i| - |x_s - m_j|}{v}$$ Để giải tọa độ 3 chiều $x_s$, cần tối thiểu 3 mặt hypeboloit độc lập, tương ứng với ít nhất 4 cảm biến âm thanh: $$\begin{cases} v\tau_{12} = |x_s - m_1| - |x_s - m_2| \ v\tau_{13} = |x_s - m_1| - |x_s - m_3| \ v\tau_{14} = |x_s - m_1| - |x_s - m_4| \end{cases}$$
                +------------------------------+
                |     Tín hiệu âm thanh vào    |
                |   (Cụm 08 cảm biến INMP401)  |
                +--------------+---------------+
                               |
                               v
                +------------------------------+
                |     Tiền xử lý tín hiệu      |
                | (Lọc BPF + Phân tích ICA)   |
                +--------------+---------------+
                               |
                               v
                +------------------------------+
                | Phát hiện sự kiện âm thanh   |
                | (Lọc tương quan mẫu AK47)    |
                +--------------+---------------+
                               |
                               v
                +------------------------------+
                | Trích xuất cửa sổ tín hiệu   |
                |    đồng bộ trên các kênh     |
                +--------------+---------------+
                               |
                               v
                +------------------------------+
                |    Ước tính trễ TDOA (tau)   |
                | (Thuật toán GCC-PHAT-beta-TN)|
                +--------------+---------------+
                               |
                               v
                +------------------------------+
                |   Tính toán vị trí nguồn âm  |
                | (Vận tốc âm v là biến số /   |
                |    kết hợp đa cụm cảm biến)  |
                +------------------------------+

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

  • Phương pháp lý thuyết: Vận dụng lý thuyết lan truyền sóng âm, xác suất thống kê, đại số tuyến tính, giải tích và xử lý tín hiệu số để thiết lập mô hình toán.
  • Phương pháp mô phỏng: Thực hiện mô phỏng Monte-Carlo trên môi trường MATLAB để kiểm nghiệm thuật toán với các mức cự ly từ $100\text{ m}$ đến $500\text{ m}$.
  • Phương pháp thực nghiệm:
    • Thiết kế cấu trúc phần cứng thu thập dữ liệu gồm 08 cảm biến âm thanh vi cơ (MEMS) INMP401 của hãng InvenSense (độ nhạy $-42\text{ dBV}$, dải đáp ứng tần số $60\text{ Hz} - 15\text{ kHz}$, dòng tiêu thụ $< 250\ \mu\text{A}$).
    • Cụm 8 cảm biến được bố trí trên hai hình tứ diện đều lồng nhau có độ dài cạnh $1\text{ m}$, đảm bảo tính đồng đều búp sóng thu theo mọi hướng không gian.
    • Tín hiệu thử nghiệm: Xung âm thanh tạo giả tiếng nổ đầu nòng của súng AK47, phát qua loa nén công suất $30\text{ W}$ tại các cự ly $100\text{ m}$, $200\text{ m}$, $300\text{ m}$, $400\text{ m}$, $500\text{ m}$ trong môi trường thực địa.
    • Dữ liệu âm thanh được số hóa và ghi nhận đồng bộ tại tần số lấy mẫu $f_s = 20\text{ kHz}$.

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

Chương 1: Tổng quan về định vị nguồn âm sử dụng nguyên lý TDOA

Chương 1 hệ thống hóa lịch sử phát triển, khả năng ứng dụng trong dân dụng (hệ thống SoundCompass với 52 cảm biến MEMS định vị ô nhiễm tiếng ồn, SoundCam của hãng CAE chẩn đoán hỏng hóc cơ khí xe hơi, hệ thống hội nghị truyền hình) và quân sự (hệ thống ASLS của Rheinmetall trên xe bọc thép, thiết bị mảng Bionic XS-56 của CAE cảnh giới drone, khí tài định vị điểm nổ pháo binh và trinh sát bắn tỉa).

Chương này phân tích chi tiết sơ đồ khối chức năng của hệ thống định vị nguồn âm gồm 3 bước tuần tự:

  1. Tiền xử lý và phát hiện sự kiện âm thanh: Tách tín hiệu cần định vị khỏi tạp âm nền, kích hoạt trích xuất cửa sổ dữ liệu.
  2. Ước tính khác biệt thời gian đến ($\tau_{ij}$): So khớp tương quan giữa các cặp cảm biến để xác định giá trị trễ.
  3. Tính toán vị trí nguồn âm: Giải hệ phương trình phi tuyến dựa trên các giá trị $\tau_{ij}$, tọa độ cảm biến $m_i$ và vận tốc âm $v$.

Chương 1 cũng xác lập mô hình thực nghiệm với cụm 08 cảm biến INMP401 xếp dạng tứ diện kép lồng nhau, ghi nhận dữ liệu thực tế tại $f_s = 20\text{ kHz}$ và chỉ ra rằng ở cự ly $> 300\text{ m}$, phổ năng lượng tín hiệu nổ AK47 bị chìm hoàn toàn trong tạp âm nền, đặt ra yêu cầu phải cải tiến thuật toán.

Chương 2: Đề xuất giải pháp phát hiện sự kiện âm thanh và nâng cao độ chính xác ước lượng khác biệt thời gian đến

Chương 2 tập trung giải quyết hai mắt xích xử lý tín hiệu đầu vào:

  • Phát hiện sự kiện âm thanh (SKAT):
    • Đánh giá phương pháp phát hiện theo ngưỡng biên độ dựa trên bổ đề Neyman-Pearson: Hiệu quả ở cự ly gần ($100\text{ m}$), nhưng thất bại ở cự ly xa ($300\text{ m} - 500\text{ m}$) do biên độ tín hiệu giảm sâu dưới mức đỉnh của tạp âm.
    • Đề xuất giải pháp kết hợp bộ lọc thông dải (Band Pass Filter - BPF) và phân tích thành phần độc lập (Independent Component Analysis - ICA) sử dụng thuật toán FastICA dựa trên hàm độ đo phi Gauss Negentropy $J(z)$: $$J(z) \propto \left[ E{G(z)} - E{G(\nu)} \right]^2$$ Giải pháp này giúp tách tín hiệu mục tiêu khỏi nguồn nhiễu đan xen, làm tăng rõ rệt hệ số tương quan cực đại $r_{MAX}$ và nâng cao xác suất phát hiện đúng tín hiệu.
  • Nâng cao độ chính xác ước tính trễ TDOA:
    • Phân tích thuật toán tương quan chéo tổng quát biến đổi pha (GCC-PHAT) kinh điển: $$R_{x_i x_j}^{PHAT}(\tau) = \int_{-\infty}^{\infty} \frac{S_{x_i x_j}(f)}{|S_{x_i x_j}(f)|} e^{j 2 \pi f \tau} df$$
    • Đề xuất giải thuật GCC-PHAT-$\beta$-TN, đưa vào hàm trọng số thích nghi với hệ số loại bỏ biên độ $\beta \in [0, 1]$: $$\Phi_{\beta}(f) = \frac{1}{|S_{x_i x_j}(f)|^\beta}$$ Khi $\beta = 0$, hàm trọng số trở về dạng tương quan kinh điển (CCC); khi $\beta = 1$, trở thành GCC-PHAT tiêu chuẩn. Bằng cách chọn $\beta$ thích nghi theo điều kiện SNR thực tế, hàm tương quan làm sắc nét đỉnh cực trị, triệt tiêu các đỉnh giả do nhiễu ở cự ly xa ($500\text{ m}$).
Cự ly định vị Phương pháp phát hiện / Xử lý tín hiệu Hệ số tương quan trung bình Khả năng phát hiện sự kiện
$100\text{ m}$ Ngưỡng biên độ thuần túy Tín hiệu phân tách rõ Tốt
$100\text{ m}$ Lọc thông dải (BPF) + Tương quan Cải thiện biên độ đỉnh Rất tốt
$100\text{ m}$ Tiền xử lý ICA + GCC-PHAT-$\beta$-TN Đỉnh tương quan sắc nhọn Tối ưu
$300\text{ m}$ Ngưỡng biên độ thuần túy Bị chìm trong tạp âm Kém / Bỏ sót mục tiêu
$500\text{ m}$ Lọc thông dải (BPF) Tạp âm còn ảnh hưởng lớn Trung bình
$500\text{ m}$ Tiền xử lý ICA + GCC-PHAT-$\beta$-TN Loại bỏ triệt để đỉnh giả Ổn định

Chương 3: Xây dựng giải pháp định vị nguồn âm với vận tốc âm thanh là biến số

Chương 3 giải quyết bài toán tính toán tọa độ không gian:

  • Phân tích tác động môi trường tới vận tốc âm: Vận tốc âm thanh trong không khí $v$ phụ thuộc chặt chẽ vào nhiệt độ tuyệt đối $T$ theo hệ thức $v = \sqrt{\frac{\gamma R T}{M}}$ (trong đó $\gamma$ là hằng số đoạn nhiệt, $R$ là hằng số khí, $M$ là khối lượng phân tử). Ngoài ra, độ ẩm, áp suất và vận tốc/hướng gió làm biến dạng phương truyền sóng. Sai số giả định vận tốc âm cố định gây ra độ lệch vị trí nguồn âm tăng dần theo cự ly.
  • Phương pháp định vị coi vận tốc âm là biến số:
    • Chia nhỏ không gian định vị thành các phân vùng hình học để giới hạn vùng tìm kiếm.
    • Thiết lập bài toán tối ưu hóa giải đồng thời 4 ẩn số bao gồm tọa độ nguồn âm $[x_s, y_s, z_s]^T$ và vận tốc truyền âm thực tế $v$.
  • Mô hình định vị đa cụm cảm biến:
    • Khảo sát sự suy giảm sai số định vị khi mở rộng số lượng cụm cảm biến từ 1 cụm, 2 cụm đến 4 cụm cảm biến phân tán.
    • Việc kết hợp dữ liệu TDOA từ 4 cụm cảm biến giúp giảm thiểu hiện tượng suy thoái hình học (Geometric Dilution of Precision - GDOP), duy trì sai số định vị ở mức thấp ngay cả ở cự ly tiếp cận $500\text{ m}$.
  • Thiết kế tổng thể hệ thống: Xây dựng mô hình hệ thống gồm các cụm cảm biến ngoại vi (chứa mảng MEMS, khối tiền khuếch đại, ADC, xử lý sơ bộ) truyền dữ liệu về cụm xử lý trung tâm (chạy thuật toán GCC-PHAT-$\beta$-TN và giải hệ phương trình tối ưu hóa đa biến).

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

Đóng góp mới về mặt lý luận và khoa học

  1. Xây dựng giải pháp phát hiện sự kiện âm thanh kết hợp tiền xử lý BPF và phân tích thành phần độc lập (FastICA): Cung cấp giải pháp phân tách tín hiệu mục tiêu khỏi nguồn nhiễu môi trường, nâng cao xác suất phát hiện sự kiện âm thanh dạng xung nổ ở khoảng cách xa với SNR thấp.
  2. Đề xuất giải thuật ước tính trễ thời gian GCC-PHAT-$\beta$-TN: Đưa hệ số điều chỉnh biên độ $\beta$ thích nghi vào hàm trọng số miền tần số, giúp làm mịn phổ chéo, làm sắc nét đỉnh tương quan chính và triệt tiêu các đỉnh tương quan giả do tạp âm ở cự ly lên tới $500\text{ m}$.
  3. Đề xuất giải pháp giải bài toán định vị nguồn âm coi vận tốc truyền âm là biến số: Loại bỏ sai số hệ thống phát sinh do giả định vận tốc âm thanh không đổi dưới tác động của biến thiên thời tiết ngoài trời; kết hợp phương pháp phân vùng không gian và mở rộng đa cụm cảm biến để nâng cao độ chính xác hội tụ tọa độ.

Đóng góp thực tiễn và kiến nghị ứng dụng

  • Thiết kế cụm 08 cảm biến MEMS dạng tứ diện đều kép lồng nhau: Cung cấp mô hình bố trí cảm biến vật lý tối ưu về kích thước ($1\text{ m}$), đảm bảo tính đối xứng không gian và thu nhận tín hiệu đồng đều từ mọi hướng.
  • Khả năng chuyển giao ứng dụng quân sự và dân sự:
    • Quân sự: Ứng dụng trong các tổ hợp trinh sát âm thanh phát hiện vị trí hỏa lực bắn tỉa, định vị điểm nổ pháo binh phục vụ hiệu chỉnh phần tử bắn, cảnh giới phát hiện phương tiện bay không người lái (UAV/Drone) tầm thấp.
    • Dân sự: Ứng dụng trong các mạng cảm biến giám sát tiếng ồn đô thị, hệ thống cứu hộ cứu nạn tìm kiếm nạn nhân qua âm thanh, và giám sát an ninh khu vực trọng yếu.

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

Hạn chế của luận án

  • Cự ly thử nghiệm và khảo sát hiệu năng mới giới hạn trong phạm vi $\le 500\text{ m}$.
  • Thực nghiệm tập trung chủ yếu vào một dạng tín hiệu chuẩn đặc thù là tiếng nổ đầu nòng của súng bộ binh (AK47); chưa mở rộng tập dữ liệu thực nghiệm cho các nguồn âm phức tạp khác.
  • Các điều kiện môi trường khí tượng (nhiệt độ, gió, độ ẩm) trong quá trình đo đạc thực nghiệm được coi là hằng số cục bộ tức thời, chưa mô hình hóa đầy đủ trường hợp gió giật hoặc gradient nhiệt độ phân tầng mạnh.
  • Giả thiết tạp âm tác động chủ yếu tuân theo phân bố Gauss.

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

  • Nghiên cứu mở rộng cự ly định vị vượt trên $500\text{ m}$ và đánh giá thuật toán với nhiều loại nguồn âm mục tiêu khác nhau (tiếng cánh quạt drone, động cơ xe cơ giới, đạn pháo cỡ lớn).
  • Tích hợp thêm các cảm biến đo thông số môi trường tức thời (nhiệt độ, tốc độ gió) tại từng trạm cảm biến để bổ trợ cho thuật toán ước lượng biến số vận tốc âm.
  • Nghiên cứu hiện thực hóa toàn bộ thuật toán trên các phần cứng nhúng xử lý tín hiệu số thời gian thực tốc độ cao (như chip DSP, FPGA) và tích hợp đồng bộ với hệ thống camera quang điện tử / bản đồ số.

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cao đối với các nhóm đối tượ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 Kỹ thuật điện tử, Xử lý tín hiệu, Kỹ thuật viễn thông và Kỹ thuật điều khiển. Có thể tham khảo phương pháp xây dựng mô hình toán kênh truyền SIMO, thuật toán FastICA tách nguồn âm và biến thể hàm trọng số thích nghi $\text{GCC-PHAT-}\beta\text{-TN}$.
  • Kỹ sư nghiên cứu và phát triển (R&D) trong lĩnh vực quốc phòng: Cung cấp cấu trúc thiết kế hình học cụm mảng cảm biến MEMS (tứ diện đều kép), giải pháp đồng bộ thời gian và quy trình xử lý tín hiệu phát hiện hỏa lực bắn tỉa, định vị điểm nổ.
  • Kỹ sư và nhà quản lý môi trường đô thị / công nghiệp: Tham khảo nguyên lý định vị nguồn âm TDOA để xây dựng các trạm quan trắc tự động phát hiện nguồn ô nhiễm tiếng ồn hoặc giám sát sự cố rung động âm thanh trong nhà máy.

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

1. Cấu trúc hình học của cụm cảm biến thực nghiệm trong luận án được thiết kế như thế nào?

Cụm thu nhận thực nghiệm sử dụng 08 cảm biến âm thanh MEMS INMP401, được bố trí tại các đỉnh của hai hình tứ diện đều lồng nhau với độ dài cạnh là $1\text{ m}$. Thiết kế tứ diện đều giúp đảm bảo tính đối xứng hình học, cho phép thu nhận tín hiệu âm thanh đồng đều từ mọi góc hướng tới trong không gian 3 chiều mà không bị suy giảm độ nhạy theo hướng.

2. Thuật toán GCC-PHAT-$\beta$-TN cải tiến yếu tố gì so với GCC-PHAT kinh điển?

GCC-PHAT kinh điển sử dụng hàm trọng số $\Phi(f) = \frac{1}{|S_{x_i x_j}(f)|}$ để làm trắng hoàn toàn biên độ, chỉ giữ lại thông tin pha. Tuy nhiên, khi tỷ số SNR thấp ở cự ly xa ($500\text{ m}$), việc làm trắng biên độ khiến tạp âm nền bị khuếch đại, gây ra nhiều đỉnh tương quan giả. Thuật toán GCC-PHAT-$\beta$-TN đưa vào hệ số điều chỉnh $\beta \in [0, 1]$ với hàm trọng số $\Phi_{\beta}(f) = \frac{1}{|S_{x_i x_j}(f)|^\beta}$, cho phép cân bằng thích nghi giữa việc làm sắc nét đỉnh trễ pha và kìm hãm nhiễu nền, giúp xác định chính xác vị trí cực trị trễ thời gian $\tau_{ij}$.

3. Vì sao luận án cần giải bài toán định vị với vận tốc âm thanh là biến số?

Trong môi trường ngoài trời, vận tốc âm thanh biến thiên liên tục theo nhiệt độ môi trường ($v \propto \sqrt{T}$), độ ẩm và vận tốc gió. Ở cự ly xa ($100\text{ m} - 500\text{ m}$), nếu sử dụng một hằng số vận tốc cố định (như $331\text{ m/s}$ hoặc $340\text{ m/s}$) trong phương trình TDOA, sai số vận tốc sẽ nhân với khoảng cách và thời gian trễ, dẫn đến sai số định vị tọa độ rất lớn hoặc làm hệ phương trình vô nghiệm. Việc coi vận tốc âm $v$ là biến số cần giải đồng thời với tọa độ $[x_s, y_s, z_s]^T$ giúp loại bỏ sai số hệ thống do điều kiện thời tiết gây ra.

4. Hiệu quả của việc tăng số lượng cụm cảm biến từ 1 cụm lên 4 cụm là gì?

Khi chỉ sử dụng 01 cụm cảm biến đơn lẻ ở cự ly xa, độ mở góc và cơ sở đo hình học bị hạn chế, dẫn đến sai số định vị tăng nhanh theo khoảng cách (hiện tượng suy thoái độ chính xác hình học GDOP). Khi kết hợp mạng lưới 04 cụm cảm biến phân tán trong không gian, hệ thống có thêm nhiều cặp trễ TDOA độc lập từ các hướng khác nhau, tạo ra các mặt hypeboloit giao cắt ở góc tối ưu, giúp giảm đáng kể sai số định vị trung bình và duy trì độ ổn định tại cự ly $500\text{ m}$.

Kết luận

Luận án tiến sĩ của tác giả Trần Công Thìn đã giải quyết một cách có hệ thống các vấn đề kỹ thuật cốt lõi trong bài toán định vị nguồn âm ngoài trời cự ly xa theo nguyên lý TDOA. Thông qua việc đề xuất bộ lọc BPF kết hợp FastICA để phát hiện sự kiện âm thanh, thuật toán tương quan thích nghi GCC-PHAT-$\beta$-TN nâng cao độ chính xác ước tính trễ thời gian, và mô hình giải hệ phương trình định vị coi vận tốc âm là biến số trên mạng mảng cảm biến đa cụm, công trình đã hoàn thiện cơ sở lý thuyết và giải pháp kỹ thuật có độ tin cậy cao. Các kết quả nghiên cứu đã được công bố trên 05 công trình khoa học chuyên ngành, đóng góp thiết thực cho sự phát triển của kỹ thuật điện tử - xử lý tín hiệu âm học phục vụ quốc phòng và dân sinh tại Việt Nam.