.24 Lưu đồ thuật tốn ước tính khác biệt thời gian đến

âm thanh các bước tính tốn tiếp theo được thực hiện, cụ thể các bước tính tốn nhằm xây dựng hàm tương quan rx1x2 (τ ) được thực hiện bao gồm

các bước:

S Có sự kiện âm thanh

Đ

Đọc tín hiệu thu được trên cảm biến thứ i xi (t) và thứ j xj (t) BẮT ĐẦU Nhập dữ liệu đầu vào Đọc tín hiệu chuẩn XS

Biến đổi Fourier tín hiệu thu được từ cảm biến

Xi ( f )  xi (t)

X j ( f )  xj (t)

KẾT THÚC

i j

ˆij  arg max rx x ( )

Ước tính khác biệt thời gian đến:

 xi xj  x1x2 r (Tính hàm tương quan: ) = 1( f )S ( f ) xi xj 1 . | S ( f ) | Tính tốn hàm trọng số ( f )  i1

Tính tốn hệ số chuẩn hóa cơng suất, với N=2:

ris  N N ( f ) j X   i i j x x S ( f )  X ( f )

Tính tốn phổ tương quan phổ năng lượng

*

j

i* và X ( f ).*

X ( f )

Σ

(τ ) = F

Biến đổi Fourier tín hiệu thu được trên các cảm biến; xác định liên hiệp phức của các tín hiệu sau biến đổi Fourier; tính tốn phổ tương quan năng lượng theo cơng thức (2.19); tính tốn hệ số chuẩn hóa cơng suất β theo cơng thức (2.32); tính tốn hàm trọng số Φ(f ). Cuối cùng rx1x2 (τ ) được

tính tốn thơng qua cơng thức (2.17) sử dụng phép biến đổi Fourier ngược. Từ hàm tương quan rx1x2 (τ ) thu được, giá trị khác biệt thời gian

đến τˆij được tính tốn thơng qua cơng thức (2.16).

Trên hình 2.24 là lưu đồ thuật tốn ước tính khác biệt thời gian đến của một cặp cảm biến, đối với một hệ thống định vị nguồn âm hồn chỉnh, cần thiết phải có q trình tính tốn, lọc giá trị phù hợp theo các giải pháp đã trình bày trong mục 2.2.1.3.

2.2.3. Mơ phỏng đánh giá giải pháp nâng cao chất lượng ước tính khácbiệt thời gian đến biệt thời gian đến

2.2.3.1. Phương pháp mô phỏng đánh giá

Để đánh giá hiệu quả của phương pháp được đề xuất, một chương trình mơ phỏng trên Matlab được xây dựng. Trên cơ sở triển khai công thức (2.25) với tín hiệu đã được số hóa thu được trên cảm biến thứ i và cảm biến thứ j, đã được phân tách thành các cửa sổ sau khi phát hiện sự kiện âm thanh:

Trong đó:

rx1x2 −1 Xi(f ) [Xj(f )]∗

|Xi(f ) [Xj(f )]∗| Σ

(2.33)

Xi(f ) = Fxi(t) là biến đổi Fourier của tín hiệu thu được trên cảm

biến thứ i;

Xj(f ) = Fxj(t) là biến đổi Fourier của tín hiệu thu được trên cảm

biến thứ j;

Xi∗(f ) và Xj∗(f ) là liên hiệp phức của Xi(f ) và Xj(f );

Dựa trên cơng thức (2.33), mơ hình mơ phỏng được xây dựng như trên hình 2.25. Trước hết, tín hiệu có tồn tại sự kiện âm thanh trên một cặp cảm biến thứ

Nạp tín hiệu trên các cảm biết, thiết lập thời gian trễ FFT FFT IFFT Tính tốn hàm trọng số Tính tốn tương quan phổ năng lượng Tính tốn sai số Xác định trễ thời gian đến Sai số e