Các quy định về chuyển đổi Hilbert Huang


7

Sau khi nghiên cứu các phương pháp phân tích tín hiệu dựa trên phân rã chế độ thực nghiệm (EMD), tôi thấy rằng các phát triển gần đây chủ yếu liên quan đến biến đổi Hilbert Huang (HHT) và phương pháp Phân tích trung bình cục bộ (LMD).

Tôi đã đọc một vài bài viết về chủ đề này và tôi muốn có ý kiến ​​của bạn về HHT.

EMD dẫn đến các chức năng chế độ nội tại đặc biệt phù hợp với biến đổi Hilbert. HHT dường như thực sự được sử dụng rộng rãi cho nhiều loại ứng dụng công nghiệp hoặc học thuật khác nhau.

Tôi có đúng không khi nghĩ rằng HHT thực chất là biến đổi EMD + Hilbert? Bạn có xem xét tính đặc hiệu của HHT nằm trong quy trình sàng lọc dựa trên phép nội suy spline (ngược lại với LMD bằng thuật toán trung bình di động chẳng hạn)?

Câu trả lời:


3

Sau khi xem qua các tài liệu liên quan đến HHT và EMD, tôi thấy rằng phần "Huang" của HHT xuất phát từ thực tế rằng anh ấy là người đề xuất EMD ngay từ đầu. Điều đó giải thích tên của phương thức ...

Để biết thêm về sự phát triển liên quan đến EMD và HHT, tôi đề xuất các bài báo của Rilling et al. "Về phân rã chế độ thực nghiệm và các thuật toán của nó". Đối với những người may mắn nói tiếng Pháp, Tiến sĩ của Rilling. luận án về EMD dường như không còn có sẵn trực tuyến; tài liệu này có vẻ rất đầy đủ và chứa một phân tích toán học rất chi tiết về EMD. Các bài viết liên quan cũng có sẵn bằng tiếng Anh tại đây trên Google Scholar.

Người ta có thể tóm tắt HHT theo cách này:

  1. EMD : phân tách tín hiệu ban đầu dưới dạng danh sách các Hàm Chế độ Nội tại (IMF);
  2. Hilbert tranform tính toán các tần số tức thời liên quan đến IMF (chính xác rất phù hợp để chuyển đổi như vậy)
  3. Phổ Hilbert , có nghĩa là biểu diễn biên độ của IMF trong miền tần số / thời gian sử dụng tần số tức thời.

Dưới đây là một ví dụ đơn giản dựa trên tín hiệu thời gian: cho s. Tín hiệu này có sự tăng trưởng bậc hai của biên độ của nó kết hợp với sự tăng trưởng tuyến tính của tần số. Đây là những gì một HHT của tín hiệu dẫn đến:y(t)= =t2STôin(t2)t[0; 16]

tín hiệu quan tâm

nhập mô tả hình ảnh ở đây

IMF thu được từ EMD

IMF đầu tiên, rất giống với <span class =y(t) ">

IMF thứ hai, gần như 0 trong ví dụ này

tần số tức thời

tần số tức thời liên quan đến IMF 1

Phổ Hilbert (trắng đến đen từ 0 đến biên độ cực đại)

nhập mô tả hình ảnh ở đây

Một số peaople xem HHT là một Biến đổi Fourier tổng quát theo nghĩa là sự phân rã của tín hiệu quan tâm bởi HHT dẫn đến cả tín hiệu thời gian thay đổi biên độ và tần số.

Một nhược điểm đáng chú ý của HHT nằm ở độ nhạy của nó đối với các hiệu ứng cạnh (điều xảy ra đối với tín hiệu gần với ranh giới bên trái và bên phải của nó trong 1D). Một số kỹ thuật tồn tại để giảm thiểu những hiệu ứng này. Các bài viết của Rilling đã đề cập ở trên đi với các chiến lược phản chiếu trong khi các chiến lược định hướng kỹ thuật khác có thể liên quan đến các kỹ thuật khớp dạng sóng.


Chào mừng đến với DSP.se! Thật tuyệt khi bạn quay lại và trả lời câu hỏi của chính mình khi bạn tìm ra câu trả lời (thậm chí bạn có thể chấp nhận nó). Tuy nhiên, sẽ còn tuyệt vời hơn (và có lợi cho người khác) nếu câu trả lời của bạn có thể trả lời rõ ràng câu hỏi bạn đã hỏi: ví dụ: viết một so sánh 1 hoặc 2 câu về HHT và EMD, và trên thực tế, tính cụ thể của HHT nằm ở đâu trong những gì bạn giả định. Hiện tại bạn đã cung cấp một câu trả lời OK , tôi nghĩ rằng thêm một vài dòng sẽ làm cho nó trở thành một câu trả lời tuyệt vời .
Penelope

Cảm ơn thông tin tốt của bạn; Tôi cần mã Matlab để áp dụng tín hiệu EMD. Bạn vui lòng cho tôi một số lời khuyên?

1
Chà, điều này phụ thuộc vào thời gian bạn dành cho nó ... Nếu bạn muốn tự viết mã EMD, tôi khuyên bạn nên đọc bài viết ( ncbi.nlm.nih.gov/pmc/articles/PMC1618495 ) của Smith. Mặt
Alain

@ Alain - cả hai liên kết của bạn đã chết hoặc không thể truy cập được.
DuckDucking

@SpaceDog Liên kết đầu tiên vẫn hoạt động với tôi. Dường như người thứ hai đã chết tuy nhiên ...
Alain
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.