tailieunhanh - Thực hành xử lý tín hiệu với Matlab part 7

Tham khảo tài liệu 'thực hành xử lý tín hiệu với matlab part 7', công nghệ thông tin, kỹ thuật lập trình phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả | Hình . Mật độ phổ công suất của tín hiệu theo phương pháp Burg. toán MUSIC Thuật toán MUSIC Multiple signal Classification là thuật toán dùng để xác định phổ của các tín hiệu ngẫu nhiên dựa vào sự phân tích phổ các véctơ riêng của ma trận tương quan. Thuật toán này phân tích không gian các trị riêng của ma trận tương quan của tín hiệu để xác định các thành phần tần số của tín hiệu. Phương pháp rất hữu hiệu cho các tín hiệu sin nhiễm tạp nhiễu trắng dạng Gauss. Trong phương pháp này hàm mật độ phổ công suất của tín hiệu được xác định từ biểu thức 1 1 N e 0 ZXv e f v p l k p l Trong đó N là chiều dài của các véctơ riêng Vk của ma trận tương quan . Véctơ e f 1 ej2 r F . ei2nr N I F T Do đó tích nội VkHe f đại diện cho biến đổi Fourier để tính mật độ phổ công suất. 259 Trong MATLAB có hai hàm xác định phổ dựa theo thuật toán này Đó là hàm pmusic và . Mỗi hàm có rất nhiều phiên bản ở đây ta chỉ sử dụng một hàm thông dụng nhất là hàm pmusic. Hàm này có các cú pháp như sau S w pmusic x p S w pmusic . nfft S f pmusic x p nfft fs S f pmusic . corr S f pmusic x p nfft fs nwin noverlap Hàm íS wJ pmusic x p thực thi thuật toán MUSIC và cho ra s là PSD của tín hiệu X và véctơ tần số đã được chuẩn hóa w trong rad sample . Phổ được tính bằng cách xác định các véctơ riêng của ma trận tương quan của dữ liệu lối vào X. Ví dụ randn state O n 0 199 X cos pi n sin pi n randn size n pmusic x 4 Kết quả cho PSD trên hình . Cũng cùng tín hiệu nhưng bây giờ dử dụng pmusic để tạo nên ma trận dữ liệu 100x7 bằng cách sử dụng cửa sổ và FFT có chiều dài 512 điểm. Tức là ta dùng pmusic x 2 512 7 0 . Ta sẽ được kết quả như trên hình . Hình . Mật độ phổ công suất của cùng tín hiệu được xác định bởi thuật toán MUSIC dùng pmusic x 4 . 260 Hình . Mật độ phổ công suất của cùng tín hiệu được xác định bởi thuật toán MUSIC dùng pmusic x 2 512 7 0 . Để so sánh các phương pháp xác định phổ công suất dùng phương pháp cổ điển phương