intTypePromotion=1
zunia.vn Tuyển sinh 2024 dành cho Gen-Z zunia.vn zunia.vn
ADSENSE
 » 
 » 

Báo cáo bài tập lớn: Xử lý thống kê tín hiệu ngẫu nhiên

Chia sẻ: Du Hành Kẻ | Ngày: | Loại File: PDF | Số trang:25

Thêm vào BST
Báo xấu
55
lượt xem 11
download
 
  Download Vui lòng tải xuống để xem tài liệu đầy đủ

Nội dung của tiểu luận trình bày ước lượng phổ tín hiệu sử dụng phương pháp Periodogram và MUSIC; dùng phương pháp Monte Carlo để đánh giá độ chính xác của phương pháp ước lượng phổ Periodogram và MUSIC; dùng phương pháp Monte Carlo để đánh giá chất lượng của phương pháp ước lượng kênh truyền dùng bộ lọc thích nghi với giải thuật RLS so với dùng bộ lọc Wiener...

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Báo cáo bài tập lớn: Xử lý thống kê tín hiệu ngẫu nhiên

  1. ĐẠI HỌC QUỐC GIA TP.HCM - TRƯỜNG ĐẠI HỌC BÁCH KHOA Chương trình kỹ sư chất lượng cao Việt Pháp BÁO CÁO BÀI TẬP LỚN XỬ LÝ THỐNG KÊ TÍN HIỆU NGẪU NHIÊN SPECTRAL ESTIMATION AND ADAPTIVE FILTERS GVHD: PGS. TS Đỗ Hồng Tuấn Nhóm 1 - Sinh viên thực hiện: Võ Nguyễn Gia Luật – 1611944 Lâm Minh Quân – 16 TP. Hồ Chí Minh, ngày 25 tháng 12 năm 2020
  2. MỤC LỤC I. PROBLEM 1: SPECTRAL ESTIMATION ............................................................... 2 1. Yêu cầu .................................................................................................................... 2 2. Thực hiện ................................................................................................................. 2 2.1. Ước lượng phổ tín hiệu sử dụng phương pháp Periodogram và MUSIC ......... 2 2.2. Dùng phương pháp Monte Carlo để đánh giá độ chính xác của phương pháp ước lượng phổ Periodogram và MUSIC .................................................................. 9 II. PROBLEM 2: ADAPTIVE FILTERS ..................................................................... 12 1. Yêu cầu .................................................................................................................. 12 2. Thực hiện ................................................................................................................ 12 2.1. Ước lượng thông số của kênh truyền sử dụng bộ lọc thích nghi dùng giải thuật RLS ......................................................................................................................... 13 2.2. Ước lượng số lượng mẫu cần thiết để đạt được sự hội tụ ............................... 19 2.3. So sánh kết quả ước lượng thông số kênh truyền sử dụng bộ lọc thích nghi dùng giải thuật RLS so với bộ lọc Wiener ............................................................. 21 2.4. Dùng phương pháp Monte Carlo để đánh giá chất lượng của phương pháp ước lượng kênh truyền dùng bộ lọc thích nghi với giải thuật RLS so với dùng bộ lọc Wiener. ................................................................................................................... 22 1
  3. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn I. PROBLEM 1: SPECTRAL ESTIMATION 1. Yêu cầu Consider a sum of sinusoids in white Gaussian noise. The sinusoids are with frequencies of Fs/2, Fs/2 and Fs/2 (𝜶 = 𝟎. 𝟏, 𝜷 = 𝟎. 𝟒,  = 𝟎. 𝟑𝟕 ) and with amplitudes of a, b and c (𝒂 = 𝟐, 𝒃 = 𝟓, 𝒄 = 𝟓), respectively. The Fs is sampling frequency. Number of signal samples is d (𝒅 = 𝟏𝟎𝟎), DFT size is e (𝒆 = 𝟐𝟎𝟒𝟖). The white Gaussian noise is with variance of 1 and zero mean. a) Plot the estimated spectrum of sum of sinusoids in white Gaussian noise using Periodogram and MUSIC. b) Evaluate the accuracy of frequency estimation for different values of noise variance from 0.5 to 5 (using Monte Carlo simulations). 2. Thực hiện 2.1. Ước lượng phổ tín hiệu sử dụng phương pháp Periodogram và MUSIC - Sử dụng lý thuyết về ước lượng phổ dựa trên phương pháp Periodogram và MUSIC để thực hiện ước lượng phổ của tín hiệu gồm 3 thành phần hình sin có dạng: 𝑥 [𝑛] = 2cos(0.1. Fs/2) + 5cos(0.4. Fs/2) + 5cos(0.37. Fs/2) Nhiễu cộng vào là AWGN với 𝑚𝑒 = 0, 𝜎𝑛2 = 1, số mẫu lấy vào là 100 mẫu, số điểm DFT là 2048 điểm. Thực hiện mô phỏng với Matlab tìm được phổ ước lượng của tín hiệu: 2
  4. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn  Nhận xét: Phổ ước lượng khi sử dụng phương pháp Periodogram có thành phần sidelobe nhiều nên nếu tín hiệu có biên độ nhỏ sẽ rất khó phân biệt do sẽ lẫn vào trong các sidelobe. Tuy nhiên, thành phần mainlobe có độ rộng hẹp nên độ phân giải khá tốt, nhất là khi hai tín hiệu có tần số gần nhau như 0.37 Fs/2 và 0.4 Fs/2. Đối với phương pháp MUSIC, phổ nhận được là phổ giả nên thành phần sidelobe rất bé nên hầu như không ảnh hưởng đến việc nhận biết các tín hiệu có biên độ bé, tuy nhiên do số lượng mẫu lấy vào chỉ là 100 mẫu nên phổ tại tần số mong muốn không tiến về vạch phổ được do đó việc nhận dạng các tín hiệu gần nhau không tốt như đối với phương pháp Periodogram. - Khảo sát lại phổ với a = 0.1 để thấy được sự giảm biên độ ảnh hưởng đến hai phương pháp như thế nào: 3
  5. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Dựa vào kết quả mô phỏng ta thấy khi biên độ tín hiệu giảm, nếu dựa vào phương pháp Periodogram thì kết quả nhận dạng không còn chính xác tín tần số của hiệu ứng với biên độ a = 1 nữa do thành phần sidelobe có biên độ lớn nên rất khó để phân biệt nó so với tín hiệu thực. Phương pháp MUSIC mặc dù hơi khó nhìn nhưng vẫn có thể tìm ra phổ chính xác ứng với biên độ lớn nhất, vì ảnh hưởng của sidelobe là không nhiều. - Giả sử, số mẫu lấy vào tăng lên là d = 2000 mẫu, khảo sát lại phổ trên với a = 1: 4
  6. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Rõ ràng, khi số mẫu tăng lên rất nhiều thì cả 2 phương pháp đều cho kết quả rất chính xác dù tín hiệu có tần số khá gần nhau hay biên độ có giảm đáng kể. Trở lại với phổ của tín hiệu với các thông số ban đầu. Dựa vào phổ tín hiệu, nếu nhận biết bằng mắt ta sẽ tìm ra các tần số của tín hiệu vào tại các vị trí có đỉnh phổ cao nhất. Tuy nhiên do hai tín hiệu ở tần số 0.4 Fs/2 và 0.37 Fs/2 gần nhau nên khi sử dụng phương pháp MUSIC giữa chúng có thể xuất hiện một đỉnh nhỏ nhô cao hơn so với 2 đỉnh ứng với 2 tần số chính xác: - Nếu dùng mắt để quan sát, ta có thể phỏng đoán được đỉnh cao nhất này không ứng với tần số thực do phổ hẹp hơn nhiều so với hai đỉnh xung quanh. Vì vậy, vẫn có thể phỏng đoán được 3 tần số đầu vào từ việc quan sát phổ tín hiệu thu được. - Để có thực hiện việc nhận biết phổ dựa trên code Matlab thì cần thêm một đoạn chương trình để tìm các đỉnh của phổ. Từ các đỉnh này, lựa ra 3 đỉnh có biên độ lớn nhất: 5
  7. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Nhận thấy rằng, nếu dựa vào đỉnh cao nhất để xác định sẽ có trường hợp chọn nhằm điểm nhô lên giữa hai tần số 0.4 Fs/2 và 0.37 Fs/2 khi sử dụng phương pháp MUSIC: Để có loại bỏ đỉnh nhô lên này, cần thêm điều kiện ràng buộc trong chương trình về độ phân giải nhỏ nhất của tần số như sau : for i = 2:K-1 if((abs(BTmp(i)-BTmp(i-1))
  8. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Đoạn chương trình này chỉ lấy những đỉnh mà khoảng cách giữa nó với hai đỉnh lân cận lớn hơn 0.028 Fs/2, đồng nghĩa với nếu trong một đoạn tần số rộng 0.056 Fs/2 mà có chứa 3 đỉnh thì sẽ bỏ qua đỉnh ở giữa. Kết quả đạt được là xác định được khá chính xác tần số của tín hiệu ban đầu: - Thực hiện ước lượng phổ tín hiệu sử dụng Periodogram và MUSIC 2000 lần và lấy giá trị DFT trung bình tại mỗi điểm, ta được kết quả: 7
  9. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Phổ trung bình của 2000 lần mô phỏng độc lập này rõ ràng hơn và dễ quan sát hơn so với việc thực hiện riêng lẻ từng lần một. 8
  10. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 2.2. Dùng phương pháp Monte Carlo để đánh giá độ chính xác của phương pháp ước lượng phổ Periodogram và MUSIC - Sau mỗi lần ước lượng phổ, xác định được 3 tần số cần ước lượng: + 𝛼̂𝑃 , β̂𝑃 ,  ̂𝑃 khi dùng Periodogram. + 𝛼̂𝑀𝑈 , β̂𝑀𝑈 ,  ̂𝑀𝑈 khi dùng MUSIC. - Để đánh giá chất lượng của phương pháp ước lượng không thể chỉ dựa vào một lần thực nghiệm mà phải lặp đi lặp lại thí nghiệm nhiều lần để xác định sai số RMS (Root mean square error) của từng thông số ước lượng so với giá trị thực tế. Đây chính là phương pháp Monte Carlo: 𝜀αP/𝑀𝑈 = √𝐸(|𝛼̂𝑃 − 𝛼𝑃 |2 ) 2 𝜀βP/𝑀𝑈 = √𝐸 (|β̂𝑃 − β𝑃 | ) 2 ̂𝑃 − 𝑃 | ) 𝜀P/𝑀𝑈 = √𝐸 (| - Nếu số lần thực hiện phép lặp đủ lớn, có thể xấp xỉ trung bình thống kê bởi trung bình như sau: 𝑀 1 𝜀αP/𝑀𝑈 = √ ∑|𝛼̂𝑃 − 𝛼𝑃 |2 𝑀 𝑖=1 𝑀 1 2 𝜀βP/𝑀𝑈 = √ ∑|β̂𝑃 − β𝑃 | 𝑀 𝑖=1 𝑀 1 2 𝜀P/𝑀𝑈 ̂𝑃 − 𝑃 | = √ ∑| 𝑀 𝑖=1 - Dựa vào công thức trên, thực hiện lặp lại thí nghiệm 2000 lần với phương sai của nhiễu 𝜎𝑛2 = 1, ta được kết quả: 𝜀αP = 0.98% 9
  11. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 𝜀βP = 0.36% 𝜀P = 0.46% 𝜀αMU = 0.98% 𝜀βMU = 0.41% 𝜀MU = 1.85%  Nhận xét: Kết quả trên cho thấy khi phương sai nhiễu là 1, phương pháp Periodogram và MUSIC đều ước lượng phổ với độ chính xác cao. Đối với tần số của tín hiệu có biên độ thấp, phương pháp MUSIC cho sai số thấp hơn. Đối với các tín hiệu có tần số gần nhau nhìn chung phương pháp Periodogram cho kết quả tốt hơn với số lượng mẫu hữu hạn 100. - Thay đổi phương sai trong khoảng 0.1 đến 5, bước 0.1 khảo sát kết quả đạt được: 10
  12. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 11
  13. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn  Nhận xét: Khi phương sai của nhiễu nhỏ hơn 1, sai số RMS của cả hai phương pháp rất nhỏ cho thấy độ chính xác của cả 2 phương pháp rất tốt. Tuy nhiên khi phương sai của nhiễu tăng lên, sai số của cả 2 phương pháp cũng tăng lên nhanh. Đặc biệt khi phương sai của nhiễu lớn hơn 2, sai số trong việc ước lượng phổ trở nên lớn hơn, dẫn đến chất lượng của việc ước lượng giảm. Nhìn chung thì đối với tín hiệu có biên độ nhỏ, phương pháp MUSIC cho sai số nhỏ hơn tức là khả năng nhận biết tốt hơn do đặc tính chứa sidelobe nhỏ của MUSIC so với Periodogram. Với các tín hiệu có tần số gần nhau thì Periodogram cho kết quả nhận dạng tốt hơn với một lượng mẫu quan sát hữu hạn. II. PROBLEM 2: ADAPTIVE FILTERS 1. Yêu cầu Consider the use of Adaptive FIR filter (AF) using Recursive Least Squares (RLS) algorithm for Channel estimation. The signal and noise are stochastic processes with Gaussian distribution and zero mean. The signal and noise processes are with variances of 𝝈𝟐𝒔 = 𝟏. 𝟐 and 𝝈𝟐𝒏 = 𝟎. 𝟑𝟓 , respectively. The number of samples can be selected arbitrarily. a) Estimate the weight vector w of Adaptive FIR filter. b) Estimate the number of samples for the convergence. c) Compare the estimated weight vector w with the optimum weight vector using Wiener FIR filter. d) Evaluate the root mean squared error (RMSE) between the estimated weight vector w and the optimum weight vector for different values of 𝜎𝑛2 using Monte Carlo simulations. 12
  14. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 2. Thực hiện 2.1. Ước lượng thông số của kênh truyền sử dụng bộ lọc thích nghi dùng giải thuật RLS. - Sử dụng giải thuật RLS (Recursive Least Squares) để cập nhật lại thông số bộ lọc AF sau mỗi vòng lặp theo công thức: Khởi tạo các giá trị trọng số ban đầu của bộ lọc cũng như nghịch đảo của ma trận tự tương quan của giá trị đầu vào lúc ban đầu, thông thượng chọn 𝛿 nhỏ để các giá trị x ban đầu hầu như không tương quan nhau: 𝑤0 = 0 𝑃(0) = 𝛿 −1 𝐼 Thực hiện vòng lặp đệ quy với 𝑛 = 1,2, … , 𝑁 để cập nhật trọng số của bộ lọc dựa trên các giá trị trọng số và sai số trước đó. 𝑧(𝑛) = 𝑃(𝑛 − 1)𝑥 ∗ (𝑛) 1 𝑔 (𝑛 ) = 𝑧(𝑛) λ + 𝑥 𝑇 (𝑛 )𝑧 (𝑛 ) 𝑇 𝛼(𝑛) = 𝑑 (𝑛) − 𝑤𝑛−1 𝑥(𝑛) 𝑤𝑛 = 𝑤𝑛−1 + 𝛼(𝑛)𝑔(𝑛) 1 [𝑃(𝑛 − 1) − 𝑔(𝑛)𝑧 𝐻 (𝑛)] 𝑃 (𝑛 ) = λ Trong đó hệ số λ được gọi là hệ số quên (0 ≤ λ ≤ 1). Giá trị λ thể hiện sự ảnh hưởng của mẫu trước đó đến giá trị của ma trận hiệp phương sai. Giá trị λ gây ảnh hưởng lên khả năng bám theo của thuật toán thích ứng. Khi λ càng bé thì 13
  15. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn cập nhật các thông số càng nhanh, tuy nhiên thuật toán lại nhạy cảm với nhiễu. Trong thực tế λ thường được chọn trong khoảng từ 0,98 đến 1. - Giải thuật RLS trong Matlab: function [W, alpha] = rls(x,d,nord,lambda) delta = 0.001; X=convm(x,nord); [M,N] = size(X); if nargin < 4 lambda = 1.0; end P = eye(N)/delta; W(1,:)= zeros(1,N); alpha(1)= d(1) - X(1,:)*W(1,:).' for k=2:M-nord+1; z = P*X(k,:)'; g = z/(lambda+X(k,:)*z); alpha(k) = d(k) - X(k,:)*W(k-1,:).'; W(k,:) = W(k-1,:)+alpha(k)*g.'; P = (P-g*z.')/lambda; end end - Sử dụng giải thuật RLS đã nêu ở trên để thực hiện chương trình ước lượng kênh truyền với số mẫu là 1000, hệ số quên chọn là 1, giả sử kênh truyền FIR thực tế có dạng: 𝐻 (𝑧) = 1 + 2𝑧 −1 + 5𝑧 −2 14
  16. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn So sánh với việc thực hiện giải thuật LMS: Nhận xét: Phương pháp sử dụng bộ lọc thích nghi với giải thuật RLS cho kết quả ước lượng kênh truyền tốt, kết quả đạt được sau 1000 mẫu là [1.0233 1.9965 5.0126] hầu như đúng với giá trị thực. So với sử dụng giải thuật LMS, giải thuật RLS cho kết quả chính xác hơn, tốc độ hội tụ cao hơn nhiều mặc dù thuật toán có phức tạp hơn sơ với LMS. - Quan sát sai số giữa tín hiệu đầu ra khi đi qua kênh truyền thực tế và đầu ra của bộ lọc thích nghi. Thực hiện phép thử 200 lần để tìm trung bình bình phương của sai số, kết quả cho thấy sai số trung bình bình phương của tín hiệu đầu ra giữa bộ lọc thích nghi và kênh truyền thực tế là hằng số và bằng 0.59822 (gần về 0) cho thấy phương pháp ước lượng dựa trên AF RMS này là phù hợp để ước lượng kênh truyền. 15
  17. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn - Thực hiện lại bộ lọc AF RLS với kênh truyền FIR thực tế có dạng: 𝐻 (𝑧) = 1 + 1.4𝑧 −1 + 5.6𝑧 −2 + 2.5𝑧 −3 + 4.7𝑧 −4 16
  18. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn - Thực hiện lại với kênh truyền FIR thực tế có dạng: 𝐻 (𝑧) = 1 + 3𝑧 −1 + 2. 8𝑧 −2 + 1.5𝑧 −3 + 2.3𝑧 −4 + 5.6𝑧 −5 17
  19. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Nhận xét: Sau khi khảo sát nhận dạng các kênh truyền có bậc khác nhau cho thấy phương pháp sử dụng bộ lọc thích nghi với giải thuật RLS cho kết quả ước lượng kênh truyền có độ chính xác cao và tốc độ hội tụ nhanh. Tín hiệu ngõ ra của bộ lọc AF và tín hiệu ngõ ra thực tế của kênh truyền chưa biết khi có nhiễu là giống nhau, cho thấy bộ lọc đã ước lượng được một cách chính xác kênh truyền chưa biết. Đây chính là những ưu điểm vượt trội của RLS so với các giải thuật đơn giản như LMS. - Thực hiện lại bộ lọc AF RLS có phương sai nhiễu tăng lên 𝜎𝑛2 = 0.5 với kênh truyền FIR thực tế có dạng: 𝐻 (𝑧) = 1 + 1.4𝑧 −1 + 5.6𝑧 −2 + 2.5𝑧 −3 + 4.7𝑧 −4 18
  20. Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Nhận xét: Kết quả đạt được cho thấy sai số MMSE của bộ lọc AF RLS là 0.70589 lớn hơn so với trường hợp khảo sát với 𝜎𝑛2 = 0.35 ở trên (MMSE đạt 0.62701). Điều này cho thấy khi phương sai nhiễu tăng lên, khả năng ước lượng của bộ lọc AF cho chất lượng ước lượng giảm xuống với cùng một lượng mẫu hữu hạn đã chọn. - Thực hiện lại bộ lọc AF RLS có hệ số quên λ = 0.98 với kênh truyền FIR thực tế có dạng: 𝐻 (𝑧) = 1 + 1.4𝑧 −1 + 5.6𝑧 −2 + 2.5𝑧 −3 + 4.7𝑧 −4 Nhận xét: Hệ số quên có ảnh hưởng đến quá trình ước lượng RLS. Khi hệ số quên càng giảm, sai số của phép ước lượng càng lớn, do bộ lọc AF lúc này nhạy cảm hơn với nhiễu của kênh truyền. 2.2. Ước lượng số lượng mẫu cần thiết để đạt được sự hội tụ. - Ở câu a đã tiến hành thực hiện ước lượng kênh truyền với số mẫu cố định là 1000. Tuy nhiên, số mẫu thực tế tùy thuộc vào nhiễu kênh truyền , số tham số ước lượng và sai số mục tiêu đặt ra để đạt được sự hội tụ. Giả sử, sai số lớn nhất ̂𝑘 ước lượng và các giá trị thực tế của hàm truyền yêu có thể giữa các giá trị 𝑤 cầu cho sự hội tụ là 10−3 : ̂𝑘 − 𝑤𝑘 )/𝑤𝑘 } ≤ 10−3 𝑚𝑎𝑥 {𝑎𝑏𝑠(𝑤 19
ADSENSE

CÓ THỂ BẠN MUỐN DOWNLOAD

LV.09: Bộ Luận Văn Tốt Nghiệp Chuyên Ngành Quản Trị Kinh Doanh
81 tài liệu
1643 lượt tải
THÔNG TIN
TRỢ GIÚP
HỖ TRỢ KHÁCH HÀNG
Theo dõi chúng tôi

Chịu trách nhiệm nội dung:

Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA

LIÊN HỆ

Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM

Hotline: 093 303 0098

Email: support@tailieu.vn

Giấy phép Mạng Xã Hội số: 670/GP-BTTTT cấp ngày 30/11/2015 Copyright © 2022-2032 TaiLieu.VN. All rights reserved.

 

Đồng bộ tài khoản
2=>2