Bạn nên lưu ý rằng việc ước tính phổ công suất bằng cách sử dụng biểu đồ không được khuyến nghị và trên thực tế đã không được thực hiện kể từ ~ 1896. Đây là một công cụ ước tính không nhất quán cho bất cứ điều gì ít hơn hàng triệu mẫu dữ liệu (và thậm chí sau đó ...), và nói chung thiên vị. Điều tương tự chính xác áp dụng cho việc sử dụng các ước tính tiêu chuẩn của tự động tương quan (ví dụ Bartlett), vì chúng là các cặp biến đổi Fourier. Miễn là bạn đang sử dụng một công cụ ước tính nhất quán, có một số tùy chọn có sẵn cho bạn.
Điều tốt nhất trong số này là ước tính nhiều cửa sổ (hoặc độ côn) của phổ công suất. Trong trường hợp này, bằng cách sử dụng các hệ số của mỗi cửa sổ với tần suất quan tâm, bạn có thể tính toán Thống kê Harmonic F dựa trên giả thuyết khống về tiếng ồn trắng. Đây là một công cụ tuyệt vời để phát hiện các thành phần dòng trong tiếng ồn, và rất được khuyến khích. Đây là lựa chọn mặc định trong cộng đồng xử lý tín hiệu để phát hiện các tạp âm trong tạp âm theo giả định về sự ổn định.
Bạn có thể truy cập cả phương pháp ước lượng đa phổ và kiểm tra F liên quan thông qua multitapergói trong R (có sẵn thông qua CRAN). Các tài liệu đi kèm với gói phải đủ để đưa bạn đi; kiểm tra F là một tùy chọn đơn giản trong lệnh gọi hàm spec.mtm.
Tài liệu tham khảo ban đầu xác định cả hai kỹ thuật này và đưa ra các thuật toán cho chúng là Ước tính phổ và Phân tích hài hòa , DJ Thomson, Kỷ yếu của IEEE, tập. 70, pg. 1055-1096, 1982.
Dưới đây là một ví dụ sử dụng tập dữ liệu đi kèm với multitapergói.
require(multitaper);
data(willamette);
resSpec <- spec.mtm(willamette, k=10, nw=5.0, nFFT = "default",
centreWithSlepians = TRUE, Ftest = TRUE,
jackknife = FALSE, maxAdaptiveIterations = 100,
plot = TRUE, na.action = na.fail)
Các tham số bạn cần lưu ý là k và nw : đây là số lượng cửa sổ (được đặt thành 10 ở trên) và sản phẩm băng thông thời gian (5.0 ở trên). Bạn có thể dễ dàng để những thứ này ở các giá trị gần như mặc định này cho hầu hết các ứng dụng. Lệnh centreWithSlepians loại bỏ một ước tính mạnh mẽ về giá trị trung bình của chuỗi thời gian bằng cách sử dụng phép chiếu lên các cửa sổ Slepian - điều này cũng được khuyến nghị, vì để lại giá trị trung bình tạo ra nhiều năng lượng ở tần số thấp.
Tôi cũng sẽ khuyên bạn nên vẽ sơ đồ phổ từ 'spec.mtm' trên thang ghi nhật ký, vì nó làm sạch mọi thứ đáng kể. Nếu bạn cần thêm thông tin, chỉ cần đăng và tôi rất vui được cung cấp.
bootspecdenscó thể hữu ích.