1
E.3. Tính toán của phổ tần số bằng phương pháp Blackman-Tukey
Phổ tần số thường được ước tính theo phương pháp Blackman-Tukey,
mà là dựa trên các định lý Wiener-Khinchine. Thủ tục Blackman- Tukey
đòi hỏi tính toán trong các bước sau đây (Massel, 1996):
1. Trừ một giá trị trung bình từ dữ liệu kỹ thuật số
ü
^
n , xu hướng loại bỏ và lọc
(nếu cần thiết).
2. Tính năng tự tương quan:
N r
n
n r
^
n
^
N r
K rt
1
1 , r = 0,1,2, ..., m, (E.9)
nơi r được gọi là số tụt hậu và m là số lag tối đa (mgiá trị m cung cấp các ước tính tối ưu cho hàm tự tương quan. Một giá trị hữu hạn
của m hàm ý rằng độ cao bề mặt ) t (, vào những thời điểm t mt là không tương quan.
3. Sự ngăn chặn sự rò rỉ phổ sử dụng một cửa sổ cho các tương quan
chức năng. Cửa sổ như vậy thon chức năng tự tương để loại bỏ sự gián đoạn
tại kết thúc chức năng K. có rất nhiều cửa sổ như trong sử dụng Một cửa sổ điển hình là.
cửa sổ Hanning:
m
r
u cos giờ
1
2
. 1 (E.10)
chức năng tự tương biến đổi trở thành :
giờ
~
K rt K rtu (E.11).
4. Tính toán mật độ tần số quang phổ bằng tích phân số của
hàm tự tương Krt ~
:
ü
K cos mt (k)
m
rk KK r () t cos
t S
m ~
r
~ ~
k
1
1
0 2 (E.12)
cho các tần số
m t
K kk
, k = 0,1,2, ..., m.
2
các S ước tính mô tả thời gian trung bình của về
2
) t (^
về của nó
thành phần tần số đặt bên trong các băng tần:
2 2
e Be, B
, chia
bởi độ rộng băng Be (rad / s). Phương trình (E.12) cho m
2
1 ước tính độc lập
của quang phổ. Các ước tính được phân cách bởi gia số tần số nhỏ hơn
m
f c 2
được
tương quan
Đối với một băng thông cho Bé, các yêu cầu số lag tối đa m là:
B t
m
e
2 . (E.13)
Sai số chuẩn dự toán phổ thường được trình bày như là một chức năng của
số bậc tự do n:
N
m
n 2
, (E.14)
nơi
m
N
n
2 (E. 15)
cần lưu ý rằng, để tồn tại nhất định, khi số lượng lag tối đa m nhỏ, thì
lỗi cũng nhỏ. Tối thiểu tổng chiều dài kỷ lục t, cần thiết để đạt được một mong muốn
lỗi về các thông số, được tính bằng:
2 2
2
mt
B
N
t
e
. (E.16)
đang được dịch, vui lòng đợi..
