Nghiên cứu, ứng dụng phân tích thành phần độc lập vào hệ thống DS-CDMA

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (2.55 MB, 15 trang )

<span class="text_page_counter">Trang 1</span><div class="page_container" data-page="1">

(PT

<small>—— GÌ GGEEZẸFĐ</small>

NGƠ VĂN THÁI

NGHIÊN CỨU, UNG DUNG PHAN TÍCH THÀNH PHAN ĐỘC

</div><span class="text_page_counter">Trang 2</span><div class="page_container" data-page="2">

<small>Luận văn được hồn thành tại:</small>

HỌC VIỆN CƠNG NGHỆ BƯU CHÍNH VIÊN THƠNG

Người hướng dẫn khoa học: TS. Đặng Hồi Bắc

Phản biện 1: TS. Đinh Chí Hiếu

Phản biện 2: TS. Nguyễn Quốc Khương

Luận văn sẽ được bảo vệ trước Hội đồng chấm luận văn thạc sĩ tại Học viện Cơng nghệ Bưu chính Viễn thơng

<small>Vào lúc: II giờ 45 ngày 20 tháng 9 năm 2015</small>

</div><span class="text_page_counter">Trang 3</span><div class="page_container" data-page="3">

MO DAU

Bài toán phân tách nguồn mù trong nhiều năm trở lại đây đã được nghiên cứu rat nhiều và áp dụng cho nhiều lĩnh vưc: viễn thông, y tế, quân sự... Đây là một phương pháp được sử dụng phổ biến cho mục đích đánh giá các nguồn tín hiệu ban đầut ừ tập hợp các tín hiệu trộn của chúng, mà khơng có thơng tin(hay có rat ít thơng tin) về q trình trộn cũng như là tín hiệu gốc. Phương pháp được biết đến nhiều nhất và có tính ưu việt là phương pháp phân tích thành phần

<small>độc lập (Independent Component Analysis).</small>

Sau mot qua trinh tim hiéu va nghiên cứu em quyết định lựa chọn luận văn “Nghiên cứu,

ứng dụng phân tích thành phần độc lập vào hệ thống DS-CDMA”. Luận văn tập trung nghiên

<small>cứu bộ thu đa người sử dụng sử dụng phương pháp ICA phục vụ cho quá trình đánh giá kênh</small>

truyền hướng xuống của hệ thống DS-CDMA. Nghiên cứu bộ thu đa người sử dụng ICA, đề xuất bộ thu ICA-SUD để có được kết quả tốt nhất.

Bồ cục luận văn bao gồm ba chương:

Chương 1: Nghiên cứu chung về phân tích thành phần độc lập (ICA)

Nghiên cứu các cơ sở lý thuyết của phân bố Gauss và phi Gauss, khái niệm độc lập,

phương pháp phân tích thành phần độc lập (ICA)

Chương 2: Nghiên cứu về thuật toán FASTICA

Nghiên cứu hàm tối ưu Negentropy, phương pháp ICA bằng cực đại hóa Negentropy, thuật tốn FastICA, kết quả mơ phỏng tách nguồn mù, đánh giá chất lượng phân tách

Chương 3: Ứng dụng phương pháp ICA cho quá trình đánh giá kênh truyền hướng xuống của hệ thống DS-CDMA

Trình bày tổng quan về các công nghệ đa truy nhập, mô hình tín hiệu của hệ thống trải phơ chuỗi trực tiếp (DSSS), mô phỏng bộ thu hướng xuống DS-CDMA nền tảng ICA

</div><span class="text_page_counter">Trang 4</span><div class="page_container" data-page="4">

CHƯƠNG 1: NGHIÊN CỨU CHUNG VỀ PHÂN TÍCH THÀNH PHẢN

ĐỘC LAP (ICA)

Phương pháp ICA giả định rằng các ngn tín hiệu là độc lập với nhau và có tính phi

Gauss. Khi mà gia thiết các thành phân nguồn độc lập được thỏa mãn, quá trình phân tách

nguồn mù sử dụng phương pháp ICA sẽ cho kết quả tốt. Một điều kiện quan trọng là nếu có N nguồn tín hiệu, thì phải có ít nhất N tín hiệu trộn đề có thể đánh giá được các nguồn tín hiệu

ban đâu.

1.1 Phân bố Gauss và phi Gauss

Phân bố Gauss hay còn gọi là phân bố chuẩn là một phân bố xác suất cực ky quan trọng trong nhiều lĩnh vực. Nó là họ phân bố có dạng tông quát giống nhau , chỉ khác tham số kỳ vọng p1 và phương sai ø?. Một biến ngẫu nhiên X có kỳ vọng p và phương sai ổ2 được gọi là có

<small>phân bơ Gauss nêu hàm mật độ xác suât của nó dạng:</small>

exp (--) (1.1)

[@) =

<small>Supper Gauss</small>

Ham phân bố dang Supper Gauss là hàm phân bố mật độ xác suất đặc biệt: hàm khá là lớn khi các biến có giá trị tương đối gần khơng hoặc là có giá trị lớn, tại các giá trị trung bình thì hàm mật độ xác suất khá là nhỏ Một ví dụ điển hình của Supper Gauss là phân bố Laplace. Phân bố Laplace có cơng thức như sau:

p(x| ,b) =z-exp(——) (1.2)

<small>Sub Gauss</small>

Ham phân bố dang Sub Gauss có phân bố xác suất det, hàm phân bồ rất ít ở các gid trị lớn và giá trị gần 0, phân bố nhiều ở các giá trị trung bình. Một ví dụ điển hình là phân bố uniform,

<small>thường được dùng trong giả lập máy tính:</small>

f@ =| (13)

</div><span class="text_page_counter">Trang 5</span><div class="page_container" data-page="5">

Khái niệm độc lập được định nghĩa bởi hàm mật độ xác suất (probability density function

- pdf). Các biến ngẫu nhiên s¿, s¿...., s„ là độc lập khi và chi khi:

ĐẤS, Sạ, ..., Sn) = Ð1(S1)P2 (S2) ... Pn (Sn) (1.4)

với p(S¡, 82, ..., s„) ham mật độ xác suất của $1, Sp, ...,S, và ps; là hàm mật độ xác suất

<small>riêng của Sj.</small>

<small>1.3 Phương pháp ICA</small>

Nguyên lý thực hiện ICA v ới giả thiết là các thành phần độc lập với nhau, dựa trên

định lý giới hạn trung tâm , đó là “Ham phân bé của tông nhiều biến ngẫu nhiên độc lập luôn

<small>hội tụ tới phân bô Gaussian”.</small>

<small>1.3.1 Định lý giới hạn trung tâm</small>

L4 N \: = , ' , w

<small>Tiến tới Gauss phi Gauss phi Gauss phi Gauss</small>

<small>Hình 1.1- Minh họa định lý giới hạn trung tâm.</small>

Dinh lý gi ới hạn trung tâm chỉ ra rằng : do x;,i = 1, N là tổ hợp tuyến tinh của các tin

<small>và ngược lại s; sẽ có tính phi</small>

<small>hiệu nguồn %,ƒj = 1N, nên x; sẽ có tính Gaussian hon Sjs</small>

<small>Gaussian hơn x;.</small>

1.3.2 Các giả thiết của phương pháp ICA

e Các nguồn nguyên thuỷ ban đầu phải độc lập thống kê với nhau .

e Khơng có thành phan độc lập nào (nguồn) có phân bố Gaussian , hoặc tối đa chỉ có 1 nguồn có phân bố Gaussian.

<small>e Ma trận trộn A là ma trận vng (sơ lượng tín hiệu nguồn và tín hiệu trộn băng nhau).</small>

</div><span class="text_page_counter">Trang 6</span><div class="page_container" data-page="6">

Với các giả thiết trên , mục tiêu của bài toán phân tách nguồn mù là tìm ra các vector w

thoả mãn tơ hợp tuyến tính y = w”x đạt tính phi Gaussian cực đại , khi đó y chính là một

nguồn tín hiệu được mong đợi là tín hiệu ban đầu. Tập hợp các vector w chính là ma trận tách

1.3.3 Một số bat định trong mơ hình ICA tuyễn tính

e Khơng thê xác định lại được chính xác năng lượng ban đầu của của các nguồn tín hiệu nguyên thuỷ do cả s và 4 đều không biết nên:

x=As= (2) (ks) = (aA) (:) =e (1.5)

Do đó trong mơ hình ICA người ta ln giả thiết mọi nguồn tín hiệu nguyên thuy —s; đều

có năng lượng (phương sai) xác định, thoả mãn: E(sể) = 1 hay E(ss’) = I với I là ma trận

e Không thê xác định được thứ tự ban đầu các thành phan độc lập khi phân tach do ca s và

A đều không biết nên khi đổi vị trí các hàng trong s và A m6 hình ICA khơng thay đơi.

e Khơng xác định được dấu của các thành phần độc lập ban đầu.

1.5 Kết luận chương

Một trong các phương pháp thông dụng nhất dé giải quyết bài toán phân tách nguồn tin hiệu mù chính là phương pháp phân tích các thành phần độc lập (Independent component Analys-ICA). Phương pháp ICA giả định rang các nguồn tín hiệu là độc lập với nhau và có tinh phi Gauss. Khi mà giả thiết các thành phần nguồn độc lập được thỏa mãn, quá trình phân tách nguồn mù sử dụng phương pháp ICA sẽ cho kết quả tốt. Ứng dụng đơn giản chính là van dé tách âm thanh tại một bữa tiệc, ta cần phải có giả thiết khi dùng ICA để giải quyết bài tốn

<small>này, đó là khơng có thời gian trễ cho các tín hiệu âm thanh cũng như khơng có độ vọng. Và</small>

một điều kiện quan trọng nữa là nếu có N nguồn tín hiệu, thì phải có ít nhất N tín hiệu trộn dé

</div><span class="text_page_counter">Trang 7</span><div class="page_container" data-page="7">

CHƯƠNG 2: NGHIÊN CỨU VE THUẬT TOÁN FASTICA

Chương 2 của luận văn tơi sẽ trình bày về thuật tốn FastICA, một thuật toán dựa vào hàm toi ưu Negentropy để ước lượng tính phi Gauss, và dùng thuật tốn tối wu Newton dé giúp

<small>thuật toán hội tụ nhanh.</small>

<small>2.1 Negentropy</small>

Một kết quả cơ sở của lý thuyết thông tin là một biến với phân bố Gauss có entropy lớn

nhất trong tất cả các biến ngẫu nhiên của biển đổi cân bang. Diéu nay co nghia 1a entropy vi

phân có thé được sử dụng như một phép do tính phi Gauss. Cách do này có thé đưa ra sự khác nhau giữa entropy của một biến ngẫu nhiên Gauss y„„„;; và một biến phi Gauss y . Phép do này

được gọi là negentropy, được biểu diễn là J, và được định nghĩa là :

JO) = Hsaus;) — HO) (2.1)

Trong đó H là ham entropy, yGauss là một vector ngẫu nhiên N chiều có phân bố Gaussian,

<small>và có chung ma trận hiệp phương sai với vector y.</small>

Vậy negentropy là một phương pháp đề đo tính phi Gauss của tín hiệu.

2.2 Phương pháp ICA bằng cực đại hóa Negentropy

Phương pháp này gồm 3 bước : Tiền xử ly dit liệu, xấp xỉ hóa negentropy, tối ưu hóa hàm

xap xi negentropy.

2.2.1 Q trình tiền xử lý

Lua chọn những thuật tốn cho ICA dùng qui tâm, trắng hóa và giảm chiều là những bước tiền xử lý để làm đơn giản và giảm tính phức tạp của vấn đề lặp thuật tốn.

Đề khơng làm mat tính tổng quát, chúng ta có thể giả thiết rằng cả các biến trộn và các thành phần độc lập đều có trung bình bang 0. Giả thiết này làm đơn giản hóa các lý thuyết và

thuật tốn rất nhiều.Nếu các tín hiệu chưa có giá trị trung bình bằng 0, chúng ta có thể thực

hiện q trình tiền xử lý, gọi là phép quy tâm tức trừ phân bố của các biến được khảo sát với

<small>các giá trị trung bình của chúng:</small>

</div><span class="text_page_counter">Trang 8</span><div class="page_container" data-page="8">

Trong đó x là vecto ngẫu nhiên chưa có trung bình là khơng. Sau khi đã ước lượng ma

trận A và các thành phan s ta có thé thêm trở lại các trị trung bình của chúng:

Khi vecto ngẫu nhiên x (hoặc s) có trị trung bình bằng khơng thi hiệp phương sai va tương quan của nó giống nhau.

° Trắng hóa

Sự trắng hố với mục dich làm cho dir liệu trở nên bat tương quan . Nếu ma trận vector có các vector thành phan x; đều thoả mãn điều kiện phương sai bằng 1 hay ma trận hiệp phương sai đồng nhất với ma trận đơn vị I: E{xxT = 1} thì ma trận vector x được xem là “trắng”.

Q trình trắng hố thự c chat là một phép biến đối tuyến tính : z =V.x . Trong đó x là

dữ liệu cần làm trắng, V là ma trận trắng hoá, z là dit liệu đã trang hố.

® Vai trị của q trình trắng hod

Với một tín hiệu đã trắng hố _, nhiệm vụ cịn lại của ICA là tìm ra một vector w sao cho y¡ = wTz dat giá trị phi Gaussian cực đại dưới điều kiện ràng buộc E(y?) = 1. Khi đó y,

<small>sẽ tương ứng với một ngn tín hiệu ban đâu.</small>

2.2.2 Xap xi hóa Negentropy

Đối với một biến ngẫu nhiên phân bố Gaussian, negentropy luôn bằng 0 và với tat cả các

loại biến cịn lại (phi Gaussian) negentropy ln có giá trị dương . Tuy nhiên vấn đề ở chỗ

chúng ta khơng thể tính negentropy một cách trực tiếp, mà phải đánh giá negentropy thơng qua

<small>2.3 Thuật toán FastICA</small>

Thuật toán FastICA theo phương pháp tách nguồn tuần tự đối voi dir liệu đã trang hoá

1. Chọn m là số nguồn tín hiệu cần đánh giá, gan p © 1.

2. Khởi tạo ngẫu nhiên giá trị vector đơn vị w, thoả man: \|w, lš =1

3.Gán w, — E{zg(w7z)} — E{g'(wjz)}w,

4. Thực hiện trực giao hoá theo phương pháp trực giao tuần tự (deflation) .

</div><span class="text_page_counter">Trang 9</span><div class="page_container" data-page="9">

6. Nếu thuật toán chưa hội tu:

|lewz*1,wz)| — 1| < €, trở lại bước 3.

7. Gắn p — p + 1 nếu p <m trở lại bước 2.

Sau đây là lưu đồ thuật toán cho q trình mơ phỏng thuật tốn FastICA.

</div><span class="text_page_counter">Trang 10</span><div class="page_container" data-page="10">

2.4 Kết quả mô phỏng

2.4.1 Mô phỏng tách nguồn mù âm thanh

<small>2.4.2 Mơ phóng tách mù hình ảnh</small>

2.4.3 Mơ phóng loại bó nhiễu trong ảnh 2.5 Đánh giá chất lượng phân tách

2.6 Kết luận chương

Phân tích thành phần độc lập (ICA) là một kỹ thuật thống kê tổng quát, ở đó đữ liệu độc lập được xem xét dé tách bang cach dua cac thanh phan độc lập lớn nhất, đã ước lượng được từ những thành phần khác nhau dưới giả thiết là khơng có hoặc có tối đa một nguồn có phân bố Gauss. ICA được thực hiện với nhiều thuật toán khác nhau, nhưng nguyên tắc cơ bản của nó đều dựa vào sự bất tương quan phi tuyến bằng cực đại hoặc cực tiểu negentropy hoặc cực đại

thơng tin tương hỗ.

<small>Thuật tốn FastlCA là một thuật toán đa nhiệm của ICA. Đây là thuật toán dựa trên</small>

phương pháp điểm cố định và có những ưu điểm sau:

<small>e Hộitụnhanh</small>

<small>e Thich hợp cho cả Super-Gauss và Sub-Gauss.</small>

ICA đã chứng minh được rằng nó hồn tồn có thé giải quyết van dé bài toán Cocktail và

một số van dé BSS khác. ICA có thé ứng dụng trong nhiều trường hợp đa dạng, ví dụ như phan tách mù trong điện não đồ (EEG) và trường não đồ (MEG), phân tích chuỗi thời gian kinh tế,

xử lý ảnh...Mới đây, ICA đã được ứng dụng rộng dai trong viễn thông dé phân tách tín hiệu

của người sử từ nhiễu của những người sử dụng khác trong DS-CDMA. Sử dụng ICA như là

một công cụ xử lý cho máy thu nối tiếp thơng thường, ICA có khả năng làm giảm sóng

jamming liên tục trong hệ thong DS-CDMA. ICA có tiềm năng sử dung rất lớn. Vì vậy nghiên

cứu phân tách những tín hiệu trong tac động của tạp âm có thé được nghiên cứu.

</div><span class="text_page_counter">Trang 11</span><div class="page_container" data-page="11">

CHƯƠNG 3: ỨNG DỤNG PHƯƠNG PHÁP ICA CHO QUÁ TRÌNH ĐÁNH GIÁ KENH TRUYEN HUONG XUONG CUA HE THONG

Trong chương này, em xin được trình bày tổng quan về các cơng nghệ da truy nhập, mơ hình tín hiệu của hệ thống DSSS và mô phỏng bộ thu SUD, bộ thu ICA và bộ thu kết hợp ICA-SUD. Với kết quả thu được, chúng ta sẽ đánh giá dé tìm ra một bộ thu toi wu cho hệ thống.

Đa truy nhập có thể hiểu là sự chia sẻ tài nguyên vô tuyến hoặc là sử dụng chung môi trường truyền dẫn chung (thời gian, tần số) giữa nhiều người sử dụng. Các công nghệ đa truy nhập được phát triển với mục đích cải thiện hiệu suất sử dụng tài ngun vơ tuyến (vì tài

ngun vơ tuyến là giới hạn)

3.2 Mơ hình tín hiệu của một hệ thống DSSS

Phan này sẽ giới thiệu mơ hình tín hiệu của hệ thống DSSS sử dụng điều chế BPSK.

Trong hệ thống này dữ liệu b(t) được trải phổ bang tín hiệu PN c(t) nhờ việc nhân hai tin

<small>hiệu nay với nhau. Tín hiệu nhận được b(t)c(t) là tín hiệu nhị phân lưỡng cực biên độ +1 sau đó</small>

sẽ dùng dé điều chế cho sóng mang sử dung BPSK, cho ta tín hiệu DSSS-BPSK s(t). Ta có thé thấy một bit bản tin bằng một chu kỳ của tín hiệu PN (trong tín hiệu này P=7).

(BPS V

Hình 3.1- Sơ đồ khối máy phát DSSS-BPSK

</div><span class="text_page_counter">Trang 12</span><div class="page_container" data-page="12">

v(ĐofÐ Ị [| » t

Hình 3.2 —Sơ đồ tín hiệu bên phát 3.2.1 Điều chế tín hiệu

3.2.2 Giải điều chế tín hiệu

3.3 Mơ phỏng bộ thu hướng xuống DS-CDMA nền tang ICA 3.3.1 Mô phỏng bộ thu DS-CDMA cổ điển

Chúng ta thực hiện mơ phỏng giải điều chế DSSS có sự xuất hiện của nhiễu AWGN, và mơ

hình mơ phỏng đó được thê hiện ở hình 3.3

</div><span class="text_page_counter">Trang 13</span><div class="page_container" data-page="13">

<small>3.3.2 Mơ phóng bộ thu cơ sở ICA</small>

Trong phan này, chúng ta mô phỏng máy thu hướng xuống DS-CDMA nên tảng ICA (máy thu ICA). Đầu tiên, chúng ta thực hiện mô phỏng một máy thu cơ sở ICA, sau đó nghiên cứu kich thước vector được yêu cầu bởi thuật toán ICA, xem xét một số yếu tố khác như là giá trị khởi tạo của quá trình lặp, tiêu chuẩn hội tụ. Cuối cùng mô phỏng máy thu liên kết

ICA-SUD để so sánh xác suất lỗi với máy thu ICA hoặc máy thu ICA-SUD thơng thường.

<small>3.3.2.2 Mô phỏng bộ thu ICA-SUD</small>

</div><span class="text_page_counter">Trang 14</span><div class="page_container" data-page="14">

3.4 Kết luận chương

Chương này nghiên cứu những khái niệm cơ bản của đa truy nhập và DSSS. Cuối cùng là mô phỏng bộ thu hướng xuống DS-CDMA nền tang ICA, kết qua chỉ ra rằng bộ thu ICA có thê giải quyết van đề ước lượng kí tự với mã trải phổ được u cau, thơng qua mã trải phơ có thé dùng dé nhận ra mỗi người sử dung. Do đó bộ thu ICA-SUD đã được đưa ra và tốc độ lỗi kí tự của nó thấp hơn bộ thu SUD qui chuẩn, điều này được kết luận từ số những thực nghiệm. Thậm chí nếu cơng suất của tín hiệu là giống nhau, ICA cịn có thể làm giảm nhỏ ảnh hưởng của nhiễu đa truy nhập, do đó cải thiện được hiệu suất của SUD.

</div><span class="text_page_counter">Trang 15</span><div class="page_container" data-page="15">

KET LUẬN

Đặc điểm nổi bật của ICA, phân tách nguồn bán mù tự nhiên, được sử dụng hữu hiệu trong giải điều chế CDMA. Khi bên thu có nhiều hoặc ít hơn thông tin chỉ dẫn trong hệ thống truyền thông, điển hình khi mã trải phố của người sử dụng được biết ít nhất, thì ICA có thé bỏ qua quá trình nhận dạng người dùng và tập trung giải mã. Thơng tin ưu tiên này có thể được truyền theo một hướng thích hợp với cơng nghệ ICA mù nhằm đạt kết quả tối ưu nhất. Đặc điểm thiết kế quan trọng mà những thuật toán khác trong thực tế khơng thé thực hiện được, đó là khả năng tính tốn vượt q u cau, dé có thé thực hiện trong thời gian thực, đây là cái ma

tương lai đang cần.

Đầu tiên chúng ta bắt đầu với những khái niệm cơ bản của ICA. FastlCA được đưa vào

xem xét kĩ lưỡng. FastICA là một thuật toán thực thi đa mục đích rất tốt của ICA. Đây là thuật tốn dựa trên phương pháp lặp điểm cố định, nó có một số những ưu điểm là hội tụ nhanh,

thích hợp cho cả những thành phần super-Gauss và sub-Gauss. Cuối cùng là mô phỏng bộ thu

hướng xuống DS-CDMA nền tảng ICA, kết quả chi ra rang bộ thu ICA có thể giải quyết van dé

ước lượng kí tự với mã trải phổ được yêu cầu, thông qua mã trải phổ có thé dùng dé nhận ra

mỗi người sử dụng. Do đó bộ thu ICA-SUD đã được đưa ra và tốc độ lỗi kí tự của nó thấp hơn bộ thu SUD qui chuẩn, điều này được kết luận từ số những thực nghiệm. Thậm chí nếu cơng

suất của tín hiệu là giống nhau, ICA cịn có thể làm giảm nhỏ ảnh hưởng của nhiễu đa truy

nhập, do đó cải thiện được hiệu suất của SUD.

Do thời gian có hạn và những hạn chế không tránh khỏi của việc hiểu biết các vấn đề dựa trên lý thuyết là chính nên báo cáo luận văn của tôi chắc chăn không tránh khỏi những thiếu sót.Tơi rat mong có được những ý kiến đánh giá, góp ý của các thầy và các bạn dé luận

<small>văn thêm hoàn thiện.</small>

</div>