Đề tài Tách sóng đa truy nhập dùng mạng Hopfield trong hệ thống CDMA

Tách sóng đa truy nhập dùng mạng Hopfield trong hệ thống CDMA Multiuser Detection using a Hopfield network for CDMA system ThS. Phạm Hùng Kim Khánh Tóm tắt Bài này giới thiệu một loại máy thu đa truy nhập sử dụng mạng Hopfield bằng cách kết hợp thuật toán cận tối ưu với khả năng hội tụ nhanh của mạng neural. Abstract This paper introduces a novel multi-user receiver, Hopfield network receiver, which combines fast convergence of neural network with the asymptotically optimum algorithm. I. GIỚI THIỆU Trong thời gian gần đây, hệ thống thông tin trải phổ (Spread Spectrum Communication System) đã được xem xét rộng rãi do quá trình thực hiện đã dễ dàng hơn. Hệ thống thông tin SS có khả năng chống fading đa đường và chống nhiễu cao. CDMA (Code Division Multiple Access) là một ứng dụng của hệ thống trải phổ. Trong hệ thống CDMA, tách sóng đa truy nhập (MUD - multi-user detection) là một trong những bài toán cơ sở của thông tin di động. Do tính thực thi kém của máy thu đơn kênh cổ điển nên MUD trở thành vấn đề chính trong hệ thống CDMA. Verdú đã chứng minh rằng lời giải tối ưu tương đương với quá trình cực tiểu hóa bài toán bậc hai nhưng quá trình tìm nghiệm tối ưu sẽ rất phức tạp khi số lượng user trong hệ thống tăng cao. Do đó, thông thường các thuật toán chỉ thực hiện tìm các lời giải cận tối ưu. Mạng Hopfield là một trong những phương pháp phổ biến sử dụng cho mục đích này. II. MÔ HÌNH HỆ THỐNG Xét kênh truyền có K user với nhiễu cộng Gaussian: y(t) = K k k k k 1 A b s (t) n(t) = + σ∑ , t ∈ [0,T] (1) trong đó T: chu kỳ symbol sk(t): dạng sóng nhận dạng của cho user thứ k, chuẩn hoá với năng lượng đơn vị. T 2 k k 0 s s (t)dt 1= =∫ (2) Giả sử tín hiệu nhận dạng bằng 0 ở ngoài khoảng [0,T] nên không tồn tại giao thoa liên ký tự (intersymbol interference). Ak: biên độ thu được của user k. Giá trị 2kA là năng lượng của user k. bk∈{-1,+1}: chuỗi bit truyền của user thứ k n(t): nhiễu Gaussian với mật độ phổ đơn vị. Nó thường dùng cho mô hình nhiễu nhiệt và nguồn nhiễu không liên quan đến quá trình truyền. Công suất nhiễu trên băng thông B là σ2B. Mô hình máy thu dùng bộ lọc thích hợp, còn gọi là máy thu RAKE cổ điển (Conventional RAKE detector): Hình 1 - Mô hình ngõ ra của bộ lọc thích hợp MF user 1 MF user 2 MF user K r(t) y1 y2 yK Ngõ ra bộ lọc thích hợp của user thứ k: T k k 0 y y(t)s (t)dt= ∫ k k j j jk j k A b A b n(k) ≠ = + ρ +∑ (3) Ta có thể biểu diễn phương trình (3) ở dạng vector như sau: y = RAb + n (4) trong đó: y = {y1,..,yK}T b = {b1,..,bK}T A = 1 2 K A 0 0 0 A 0 0 0 A ⎡ ⎤⎢ ⎥⎢ ⎥⎢ ⎥⎢ ⎥⎣ ⎦ K K M M K M K n là vector ngẫu nhiên Gaussian có trung bình bằng 0. Trong trường hợp kênh truyền không trực giao thì sẽ tồn tại ảnh hưởng của các user còn lại (MAI – Multiple Access Interference). Theo Verdú, xác suất lỗi của hệ thống ứng với trường hợp 2 user và K user là: 1 2 1 2A A A A1 1Q Q 2 2 − ρ + ρ⎛ ⎞ ⎛ ⎞+⎜ ⎟ ⎜ ⎟σ σ⎝ ⎠ ⎝ ⎠ (4) 1 j K jk j jkK 1 e { 1} e { 1} e { 1} j k j k AA1 ... ... Q e 2 − ∈ ± ∈ ± ∈ ± ≠ ≠ ⎛ ⎞⎜ ⎟+ ρ⎜ ⎟σ σ⎜ ⎟⎝ ⎠ ∑ ∑ ∑ ∑ (5) III. TÁCH SÓNG ĐA TRUY NHẬP DÙNG MẠNG HOPFIELD A. Mạng Hopfield Mạng Hopfield là mạng neural đơn lớp theo dạng lan truyền ngược. Mỗi neural nhận tổng các hoạt động từ các neural khác trong mạng và cập nhật theo quy luật sau: Vi = g(Ui) = i ij j j i g J T V ≠ ⎛ ⎞⎜ ⎟+⎜ ⎟⎜ ⎟⎝ ⎠ ∑ (6) Trong đó Tij là trọng số liên kết giữa neural thứ j và thứ i, Ji là trạng thái hiện tại của neural thứ i. Nếu kết nối giữa các neural là đối xứng (nghĩa là Tij = Tji) thì phương trình cập nhật của hệ thống sẽ hội tụ về trạng thái ổn định (giá trị tại ngõ ra sẽ là hằng số). Còn trong trường hợp các phần tử trên đường chéo Tii = 0 thì trạng thái ổn định của mạng gồm N neural sẽ hội tụ về giá trị cực tiểu địa phương (gọi là hàm năng lượng): E = N N N ij i j i i i 1 j 1 i 1 1 T V V V J 2 = = = − −∑∑ ∑ (7) Phương trình cập nhật của neural thứ i có thể biểu diễn như sau: i iij j i i i j dU UE T V J dt V ≠ ∂= − = + −∂ τ∑ (8) B. Máy thu HNN (Hopfield neural network) Máy thu cổ điển (CD – Conventional Detector) sử dụng một bank các bộ lọc thích hợp (matched filter) để nhận dạng từng user và ước lượng các bit thông tin chỉ dựa vào ngõ ra của các bộ lọc này: k k (i 1)T (i) k k k iT y r(t)s (t iT )dt + −τ −τ = − − τ∫ (9) (i) (i)CD kbˆ sign(y )= (10) Một phương pháp khác là dùng máy thu đa truy nhập tối ưu (OMD – Optimum Multiuser Detector). Phương pháp này thực hiện ước lượng các bit thông tin bằng cách cực đại hoá hàm khả năng (likelihood function). Trong trường hợp đồng bộ: { }TK(i) (i) TOMD b { 1, 1}b arg max 2y b b Rb∈ + −= − (11) Tuy nhiên khi sử dụng phương pháp này thì độ phức tạp của quá trình tính toán sẽ thay đổi theo hàm mũ của số user. Như vậy, khi số lượng user lớn thì quá trình thực hiện là không khả thi. Do đó, ta chỉ dùng các sơ đồ tối ưu phụ, đó là máy thu đa tầng (MSD – Multistage Detector). MSD gồm có một tập hợp các tầng, mỗi tầng dùng để ước lượng các bit thông tin như sau: (i)MSDb (m+1) = sign(y (i) – (R – I) (i)MSDb (m)) (12) Ngõ vào của tầng thứ nhất chính là ngõ ra của một bộ tách sóng cổ điển. MSD có số tầng vô hạn và có thể hội tụ tới cực tiểu địa phương của hàm đối tượng OMD. Đối với trường hợp bất đồng bộ thì bài toán máy thu tối ưu giải quyết bằng cách biểu diễn giống như phương trình (11) nhưng ma trận tương quan chéo R lúc này có dạng: R(0) R( 1) 0 0 R(1) R(0) R( 1) R 0 R(1) R(0) 0 R( 1) 0 0 R(1) R(0) −⎡ ⎤⎢ ⎥−⎢ ⎥⎢ ⎥= ⎢ ⎥−⎢ ⎥⎢ ⎥⎣ ⎦ K K M % O M O O K (13) Các bit thông tin nhận được lúc này dùng ước lượng có dạng như sau: { }T( 2 M 1) K(i) (i) TOMD b { 1, 1}b arg max 2y b b Rb+∈ + −= −%% % % %%% (14) Như vậy, độ phức tạp trong trường hợp này sẽ lớn hơn rất nhiều so với trường hợp đồng bộ. Từ (11), ta thấy rằng hàm đối tượng OMD tương tự như hàm năng lượng của HNN. Mà (11) có thể viết ở dạng như sau: { } { } { } T K T K T K (i) (i) T1 OMD 2b { 1, 1} (i) T T1 1 2 2b { 1, 1} (i) T1 2b { 1, 1} b arg min y b b Rb arg min y b b (R I)b b Ib arg min y b b (R I)b ∈ + − ∈ + − ∈ + − = − + = = − + − + = − + − (15) (do bTIb luôn là số dương). Như vậy, ta có thể chuyển trực tiếp thành hàm năng lượng của mạng Hopfield với ma trận trọng số T = –(R – I) và trạng thái ban đầu của mạng J = y(i). Xét trường hợp rời rạc trên miền thời gian: i i ij j j i V (m 1) sign y V (m) ≠ ⎛ ⎞⎜ ⎟+ = − ρ⎜ ⎟⎜ ⎟⎝ ⎠ ∑ (16) hay có thể viết ở dạng ma trận như sau: V(m+1) = sign(y – (R – I)V(m)) (17) ™ Mạng Hopfield dừng: Để tránh trường hợp cực tiểu địa phương, chúng ta có thể thay thế HNN bằng mạng Hopfield dừng (SHN – Stochastic Hopfield Network). Phương trình (17) của mạng Hopfield có thể thay đổi bằng phương trình: i i ij j j i V (m 1) sign y V (m) (m) ≠ ⎛ ⎞⎜ ⎟+ = − ρ + ν⎜ ⎟⎜ ⎟⎝ ⎠ ∑ (18) trong đó ν(m) là biến ngẫu nhiên độc lập với trung bình bằng 0 và hàm phân phối F(x,m). Levendovzky đã xác định được một dạng hàm F cho biến ngẫu nhiên ν như sau: F(x) = x 1 1 e−α+ (19) Đối với dạng hàm F như trên, phương sai của biến ngẫu nhiên phụ thuộc vào giá trị của α (giá trị của α càng lớn thì phương sai càng nhỏ). IV. KẾT QUẢ MÔ PHỎNG Kết quả mô phỏng thực hiện trên ngôn ngữ MatLab ứng với chuỗi giả ngẫu nhiên có chiều dài cực đại. Chương trình mô phỏng thực hiện so sánh một số loại máy thu tiêu biểu ứng với chuỗi giả ngẫu nhiên là chuỗi nhị phân có chiều dài cực đại (Maximal Length Binary Sequence). Do chất lượng của hệ thống thông tin được đánh giá dựa cơ sở trên tỉ số lỗi bit (Bit Error Rate) nên quá trình mô phỏng sẽ thực hiện tính toán thông số này. Trong tất cả các hình vẽ, trục tung biểu thị BER theo dB và trục hoành biểu diễn giá trị tỉ số tín hiệu trên nhiễu (SNR – Signal-to-Noise Ratio). Hình 2 – So sánh BER của các phương pháp ứng với số user K = 60, chiều dài chuỗi giả ngẫu nhiên M = 1024 bit Đối với máy thu RAKE cổ điển, có thể tính BER theo công thức (5), nhưng công thức này khá phức tạp, đòi hỏi quá trình tính toán khá lớn nên tác giả đã giả lập một kênh truyền tồn tại nhiễu trắng Gaussian (AWGN – Addition White Gaussian Noise). Đối với máy thu HNN, để tránh tình trạng cực tiểu địa phương, tác giả đã hạn chế số vòng lặp tối đa là 100 (theo chương trình thực hiện thì thông thường mạng sẽ hội tụ dưới 10 lần lặp). Hình 3 – So sánh BER của các phương pháp ứng với số user K = 30, chiều dài chuỗi giả ngẫu nhiên M = 1024 bit Dựa vào các hình vẽ, ta thấy rằng khi tỉ số SNR nhỏ thì BER của các phương pháp tách sóng không sai khác nhiều, nhưng khi thực hiện với SNR lớn hơn thì máy thu HNN sẽ cho BER tốt hơn các phương pháp khác. Hình 4 - So sánh BER của các phương pháp ứng với số user K = 60, chiều dài chuỗi giả ngẫu nhiên M = 4096 bit Hình 5 - So sánh BER của các phương pháp ứng với số user K = 30, chiều dài chuỗi giả ngẫu nhiên M = 4096 bit V. TÀI LIỆU THAM KHẢO [1] Sergio Verdú, "Multiuser Detection", ISBN 0521593735, Cambridge University Press, 1998 [2] Matti Latva-aho, "Advanced Receivers for WCDMA system", ISBN 9514250397, University of Oulu, Oulu Finland, 1998 [3] Habib Fathallah, "Narrowband Interference Suppression in Spread Spectrum CDMA Communications via Blind Equalization", 1997