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

ĐẠ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 MỤC LỤC I. PROBLEM 1: SPECTRAL ESTIMATION ............................................................... 2 1. Yêu cầu .........

pdf25 trang | Chia sẻ: huong20 | Ngày: 12/01/2022 | Lượt xem: 284 | Lượt tải: 0download
Tóm tắt tài liệu 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, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
............................................................................................................ 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 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) 2 Nhiễu cộng vào là AWGN với 푚푒 = 0, 휎푛 = 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 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 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 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 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))<0.028)&&(abs(BTmp(i)-BTmp(i+1))<0.028)) BTmp2(i-x) = []; pks21temp(i-x) = []; x = x+1; end end 6 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 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 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: (| |2) 휀αP/푀푈 = √퐸 훼̂푃 − 훼푃 ̂ 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 휀 = √ ∑|훼̂ − 훼 |2 αP/푀푈 푀 푃 푃 푖=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 2 của nhiễu 휎푛 = 1, ta được kết quả: 휀αP = 0.98% 9 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 휀βP = 0.36% 휀 = 0.46% P 휀αMU = 0.98% 휀βMU = 0.41% 휀 = 1.85% MU  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 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 11 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 2 vector w and the optimum weight vector for different values of 휎푛 using Monte Carlo simulations. 12 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 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 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 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 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 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. 2 - Thực hiện lại bộ lọc AF RLS có phương sai nhiễu tăng lên 휎푛 = 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 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 2 lớn hơn so với trường hợp khảo sát với 휎푛 = 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ó thể giữa các giá trị 푤̂푘 ước lượng và các giá trị thực tế của hàm truyền yêu cầu cho sự hội tụ là 10−3: −3 푚푎푥{푎푏푠(푤̂푘 − 푤푘)/푤푘} ≤ 10 19 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Thực hiện vòng lặp while for đến khi đạt được sai số bé hơn sai số hội tụ đã chọn thì dừng. Giả sử 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: Khi chọn điều kiện hội tụ nhỏ, thời gian hội tụ sẽ lâu nhưng bù lại được độ chính xác cao hơn. Do đó tùy vào mục đích bài toán mà chọn điều kiện hội tụ phù hợp. 20 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn 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. - Thực hiện ước lượng kênh truyền dùng bộ lọc Wiener, với bộ lọc thực tế có dạng: 퐻(푧) = 1 + 1.4푧−1 + 5.6푧−2 + 2.5푧−3 + 4.7푧−4 - So sánh kết quả của bộ lọc thực tế với bộ lọc Winer, với giá trị mẫu N = 3005 tìm được ở câu b: 21 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn Nhận xét: Kết quả so sánh cho thấy phương pháp ước lượng sử dụng bộ lọc thích nghi RLS sau khi hội tụ cho kết quả bám rất tốt so với giá trị ước lượng từ bộ lọc Wiener. Có thể kết luận phương pháp ước lượng dùng bộ lọc thích nghi RLS vừa cho kết quả tiệm cận giá trị tối ưu của bộ lọc Wiener, vừa có thể cập nhật thông số một cách thích nghi tùy thuộc vào sự thay đổi đặc tính thống kê của tín hiệu. 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. - Giả sử hàm truyền thực tế có dạng: 퐻(푧) = 1 + 1.4푧−1 + 5.6푧−2 + 2.5푧−3 + 4.7푧−4 - Sau mỗi lần ước lượng kênh truyền đạt được giá trị hội tụ, xác định được các thông số cần ước lượng: + 푤̂1, 푤̂2, 푤̂3, 푤̂4, 푤̂5 khi dùng bộ lọc thích nghi với giải thuật RLS. + 푤̂1표푝푡, 푤̂2표푝푡, 푤̂3표푝푡, 푤̂4표푝푡, 푤̂5표푝푡 khi dùng bộ lọc Wiener. - Để đá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ị chuẩn là giá trị có được từ bộ lọc Winer. Đây chính là phương pháp Monte Carlo: 2 휀W1 = √퐸 (|푤̂1 − 푤̂1표푝푡| ) 2 휀W2 = √퐸 (|푤̂2 − 푤̂2표푝푡| ) 2 휀W3 = √퐸 (|푤̂3 − 푤̂3표푝푡| ) 2 휀W4 = √퐸 (|푤̂4 − 푤̂4표푝푡| ) 2 휀W5 = √퐸 (|푤̂5 − 푤̂5표푝푡| ) 22 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn - 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 2 휀 = √ ∑|푤̂ − 푤̂ | W1 푀 1 1표푝푡 푖=1 푀 1 2 휀 = √ ∑|푤̂ − 푤̂ | W2 푀 2 2표푝푡 푖=1 푀 1 2 휀 = √ ∑|푤̂ − 푤̂ | W3 푀 3 3표푝푡 푖=1 푀 1 2 휀 = √ ∑|푤̂ − 푤̂ | W4 푀 4 4표푝푡 푖=1 푀 1 2 휀 = √ ∑|푤̂ − 푤̂ | W5 푀 5 5표푝푡 푖=1 - Dựa vào công thức trên, thực hiện lặp lại thí nghiệm 500 lần với phương sai của 2 nhiễu 휎푛 = 0.35, ta được kết quả: 휀W1 = 0.9905% 휀W2 = 1.0122% 휀W3 = 0.9624% 휀W4 = 1.0446% 휀W5 = 1.0066% 23 Stochastic Signal Processing PGS.TS Đỗ Hồng Tuấn  Nhận xét: Khi phương sai của nhiễu là 0.35, sai số RMS của phép ước lượng so với giá trị tối ưu của bộ lọc Wiener là rất nhỏ, cho thấy phương pháp sử dụng bộ lọc thích nghi RLS ước lượng trọng số kênh truyền với độ chính xác cao. - Thay đổi phương sai trong khoảng 0.05 đến 5, bước 0.05 khảo sát kết quả đạt được:  Nhận xét: Khi nhiễu AWGN tăng phương sai, sai số của bộ lọc thích nghi tăng lên nhưng giá trị sai số vẫn nhỏ trong tầm phương sai nhiễu từ 0.05 đến 5. Điều đó cho thấy khả năng hoạt động tốt của bộ lọc thích nghi khi điều kiện môi trường thay đổi. 24

Các file đính kèm theo tài liệu này:

  • pdfbao_cao_bai_tap_lon_xu_ly_thong_ke_tin_hieu_ngau_nhien_spect.pdf