thực hiện thuật toán music với các cấu trúc anten mảng ula

Khóa luận tốt nghiệp 2011 i TÓM TẮT NỘI DUNG KHÓA LUẬN Mục tiêu chính của khóa luận là thực thi thuật toán MUSIC trên kit DSP TMS320C6713 của Texas Instrument, với các cấu trúc anten mảng khác nhau nhƣ ULA, UCA Vì vậy trọng tâm của khóa luận là mô phỏng thuật toán MUSIC với các cấu trúc mảng khác nhau, từ đó tìm các giải thuật để triển khai thuật toán trên kit DSP. Chi tiết khóa luận bao gồm 5 chƣơng: Chƣơng 1: Thực hiện thuật toán MUSIC với các cấu trúc anten mảng ULA, UCA: thực hiện việc mô phỏng thuật toán sử dụng chƣơng trình Matlab, từ đó đánh giá ƣu, nhƣợc điểm của thuật toán với các cấu trúc anten khác nhau. Chƣơng 2: Giới thiệu cơ bản về kit DSP TMS320C6713. Chƣơng này phân tích các đặc điểm của kit, từ đó nêu lý do sử dụng kit để thực thi thuật toán MUSIC. Tiếp đó giới thiệu môi trƣờng giao tiếp với kit là chƣơng trình Code Composer Studio, và nghiên cứu khả năng tích hợp giữa Code Composer Studio và Matlab qua tính năng “Real Time Data Exchange”. Chƣơng 3: Thực thi thuật toán MUSIC trên kit TMS320C6713: nêu lên các công việc thực hiện đƣợc và các kết quả đạt đƣợc trong thực tế. Cuối cùng là phần kết luận và hƣớng phát triển tiếp theo cho khóa luận. Khóa luận tốt nghiệp 2011 ii LỜI CAM ĐOAN Tôi xin cam đoan khóa luận tốt nghiệp: “Thực thi thuật toán MUSIC trên kit DSPTMS320C6713” là công trình nghiên cứu của bản thân. Những phần sử dụng tài liệu tham khảo trong khóa luận đã đƣợc nêu rõ trong phần tài liệu tham khảo. Các số liệu, kết quả trình bày trong khóa luận là hoàn toàn trung thực, nếu sai tôi xin chịu hoàn toàn trách nhiệm và chịu mọi kỷ luật của khoa và nhà trƣờng đề ra. Tác giả khóa luận Nguyễn Đức Anh Khóa luận tốt nghiệp 2011 iii LỜI CẢM ƠN Để có thể hoàn thành tốt khóa luận này, em xin chân thành gửi lời cảm ơn tới Ths.Trần Thị Thúy Quỳnh, ngƣời đã hƣớng dẫn tận tình và giúp đỡ em rất nhiều trong quá trình thực hiện bài khóa luận của mình. Ngoài ra, trong quá trình thực hiện khóa luận em còn nhận đƣợc rất nhiều sự động viên và giúp đỡ từ phía gia đình, ngƣời thân và tập thể các bạn trong lớp. Do đó kết quả cũng nhƣ tính khả dụng của bài luận văn này trong thực tế là lời cảm ơn sâu sắc nhất của em gửi tới mọi ngƣời và là nguồn động lực để em có thể tự tin vào các kiến thức mình đã thu đƣợc sau khi tốt nghiệp. Khóa luận tốt nghiệp 2011 iv MỤC LỤC CHƢƠNG 1 LÝ THUYẾT VÀ MÔ HÌNH THUẬT TOÁN MUSIC 2 1.1 Giới thiệu về thuật toán MUSIC 2 1.2 Thuật toán MUSIC đối với dàn anten ULA 3 1.2.1 Mô hình dàn anten ULA 3 1.2.2 Tín hiệu thu đƣợc sau khi qua dàn anten ULA 3 1.2.3 Thuật toán MUSIC với dàn anten ULA 5 1.2.4 Mô phỏng thuật toán MUSIC với dàn anten ULA trên Matlab 8 1.3 Thuật toán MUSIC đối với dàn anten UCA 10 1.3.1 Mô hình dàn anten UCA 10 1.3.2 Tín hiệu thu đƣợc sau khi qua dàn anten ULA 10 1.3.3 Thuật toán MUSIC với dàn anten UCA 11 1.3.4 Mô phỏng thuật toán MUSIC với dàn anten UCA 11 CHƢƠNG 2 TỔNG QUAN VỀ KIT DSP TMS320C6713 13 2.1 Lựa chọn phần cứng 13 2.2 Giới thiệu chung về kit 13 2.3 Bảng mạch DSP 14 2.4 Code Composer Studio 16 2.5 Tích hợp CCS với Matlab thông qua tính năng Real Time Data Exchange . 18 CHƢƠNG 3 THỰC THI THUẬT TOÁN MUSIC TRÊN KIT TMS320C6713 21 3.1 Thiết kế thuật toán 21 3.2 Lập trình thuật toán 21 3.3 Thực thi thuật toán trên kit 26 3.3.1 Mô tả quá trình 27 3.3.2 Phân vùng bộ nhớ của kit cho dữ liệu và chƣơng trình 28 3.3.3 Kết quả thực thi thuật toán: 30 Khóa luận tốt nghiệp 2011 v DANH MỤC CÁC CHỮ VIẾT TẮT Thuật ngữ Tiếng anh Tiếng Việt DOA Direction of Arrival Hƣớng sóng tới MUSIC MUltiple Signal Classification Phƣơng pháp phân lớp đa tín hiệu ESPRIT Estimation of Signal Parameters via Rotational Invariance Techniques Phƣơng pháp đánh giá các thông số tín hiệu thông qua kỹ thuật quay bất biến SNR Signal to Noise Ratio Tỷ số tín hiệu trên nhiễu ULA Uniform Linear Array Hệ anten sắp xếp theo đƣờng thẳng UCA Uniform Circular Array Hệ anten sắp xếp theo đƣờng tròn DSP Digital Signal Proccessor Bộ xử lý tín hiệu số FPGA Field-programmable gate array Vi mạch cấu trúc mảng phần tử khả trình RTDX Real Time Data Exchange Trao đổi dữ liệu thời gian thực API Application programming interface Giao diện lập trình ứng dụng JTAG Joint Test Action Group CCS Code Composer Studio Khóa luận tốt nghiệp 2011 vi DANH MỤC BẢNG BIỂU Bảng 1: Lịch sử phát triển của thuật toán tìm hƣớng sóng đến 2 Bảng 3: Thông số mô phỏng dàn anten 12 Bảng 4: Các dạng file sử dụng trong CCS 17 Bảng 5: Các lệnh cơ bản trong thƣ viện RTDX.h để liên kết CCS→Matlab 19 Bảng 6: Các lệnh cơ bản của toolbox để liên kết Matlab → CCS 20 Bảng 7: So sánh thuật toán MUSIC trên ngôn ngữ Matlab và ngôn ngữ C 22 Bảng 8: Mô hình hệ thống thực thi thuật toán MUSIC với anten ULA 30 Khóa luận tốt nghiệp 2011 vii DANH MỤC HÌNH VẼ Hình 1: Mô hình cấu trúc anten ULA 3 Hình 2: Lƣu đồ thuật toán MUSIC 5 Hình 3: Dạng phân bố các giá trị riêng của ma trận tƣơng quan 7 Hình 5: Hệ anten ULA trong trƣờng hợp hai tín hiệu đối xứng 9 Hình 6: Mô hình cấu trúc anten UCA 10 Hình 7:Kết quả mô phỏng hệ anten UCA 12 Hình 8: Kit TMS320C6713 14 Hình 9: Sơ đồ khối kit TMS320C6713 15 Hình 10: Sơ đồ và địa chỉ vùng nhớ L2 của kit 16 Hình 11: Giao diện chƣơng trình Code Composer Studio 17 Hình 12: Mô hình giao tiếp giữa kit và PC thông qua RTDX 18 Hình 13: Quá trình liên kết từ kit đến máy tính 19 Hình 14: Quá trình liên kết từ máy tính đến kit 20 Hình 15: Mô hình xây dựng thuật toán 21 Hình 16: Lƣu đồ thuật toán Jacobi 26 Hình 17: Quá trình thực thi thuật toán MUSIC trên kit DSP 27 Hình 18: Thực thi hệ thống DSP với mảng anten ULA 30 Hình 19: So sánh kết quả thực thi thuật toán và mô phỏng trên dàn anten ULA 31 Hình 20: Thực thi hệ thống DSP với mảng anten UCA 32 Hình 21: So sánh kết quả thực thi thuật toán và mô phỏngtrên dàn anten UCA 32 Hình 22: Số nhịp CPU của chip chạy đến khi kết thúc nhận dữ liệu 33 Hình 23: Số nhịp CPU của chip chạy đến khi bắt đầu gửi dữ liệu 33 Khóa luận tốt nghiệp 2011 1 MỞ ĐẦU Anten thông minh là một công nghệ mới đƣợc quan tâm nhiều trong thời gian gần đây với các ƣu điểm: cải thiện vùng phủ sóng, giảm nhiễu, tăng dung lƣợng, mở rộng phạm vi của hệ thống [1]. Nói một cách đơn giản, anten thông minh gồm nhiều phần tử anten đơn giản kết hợp với bộ xử lý tín hiệu để có thể thay đổi búp sóng một cách tự động. Các cấu trúc anten mảng phổ biến là cấu trúc dạng đƣờng thẳng cách đều (ULA), cấu trúc dạng chữ nhật cách đều (URA) hay cấu trúc dạng đƣờng tròn cách đều (UCA). Trong các cấu trúc anten thông minh, việc xử lý tín hiệu đƣợc thực thi bằng cách thay đổi các trọng số tại mỗi phần tử anten theo một tham số tối ƣu xác định. Việc xác định sơ bộ hƣớng sóng đến (DOA) là một trong số các trọng số này. MUSIC là một trong những thuật toán xác định DOA với nhiều ƣu điểm nhƣ độ chính xác cao, khả năng áp dụng với nhiều cấu trúc anten mảng. Vì vậy thuật toán MUSIC đang ngày càng phát triển và đƣợc ứng dụng rộng rãi trong các hệ thống quân sự, điều khiển, thông tin liên lạc, truyền thông… Để phục vụ cho việc đo lƣờng, lƣu trữ và xử lý các tín hiệu thu thập từ anten trong thực tế, thuật toán MUSIC đƣợc thực thi trên các bộ xử lý số DSP. Các bộ xử lý số DSP đƣợc lựa chọn bởi khả năng xử lý tín hiệu số rất mạnh cùng khả năng tái lập trình đơn giản. Khóa luận tốt nghiệp thực hiện việc thƣc thi thuật toán MUSIC trên kit DSP TMS320C6713 nhằm mục đích tiến tới thực tiễn hóa hệ thống xác định hƣớng sóng tới. Khóa luận tốt nghiệp Chƣơng 3 2 CHƢƠNG 1 LÝ THUYẾT VÀ MÔ HÌNH THUẬT TOÁN MUSIC 1.1 Giới thiệu về thuật toán MUSIC MUSIC là một thuật toán xác định hƣớng sóng tới (DOA-Direction Of Arrival) dựa vào những tín hiệu thu thập đƣợc từ mảng anten. Trên thế giới việc nghiên cứu phƣơng pháp tìm hƣớng sóng bắt đầu từ những năm 70 của thế kỷ trƣớc, và rất nhiều thuật toán đƣợc tìm ra và tiếp tục đƣợc nghiên cứu phát triển đến tận bây giờ. Lịch sử phát triển các thuật toán có thể đƣợc trình bày tóm tắt qua bảng sau [5]: Bảng 1: Lịch sử phát triển của thuật toán tìm hướng sóng đến Năm tìm ra Tên thuật toán Tác giả thuật toán 1967 Maximum Entropy Method Burg 1969 Maximum Likehood Method Capon 1973 Covariance Method Pisarenco 1977 Subspace Method –MUSIC Schmidt 1979 Subspace Method –MUSIC Bienvenu 1989 Subspace Method –ESPRIT Richchard Roy,Thomas Kailath Từ bảng ta có thể thấy: các phƣơng pháp tìm hƣớng sóng đến hiện nay cơ bản dựa trên hai phƣơng pháp không gian con: phƣơng pháp phân lớp đa tín hiệu (MUltiple Signal Classification – MUSIC) và phƣơng pháp đánh giá các thông số tín hiệu thông qua kỹ thuật quay bất biến (Estimation Of Signal Parameters Via Rotational Invariance Techniques – ESPRIT). Phƣơng pháp ESPRIT ra đời sau nên tốc độ thực hiện nhanh hơn MUSIC, tuy nhiên MUSIC lại có ƣu điểm là đơn giản, linh hoạt hơn trong việc áp dụng cho các cấu trúc anten khác nhau nên đƣợc nghiên cứu và ứng dụng thực tế nhiều hơn. Đƣợc tìm ra bởi những nghiên cứu độc lập của Schmidt (1977) và Bienvenu (1979), thuật toán này có thể xác định số lƣợng tín hiệu đến anten cũng nhƣ hƣớng của các tín hiệu đó, chỉ dựa trên vectơ tín hiệu đầu vào. Khóa luận tốt nghiệp Chƣơng 3 3 1.2 Thuật toán MUSIC đối với dàn anten ULA 1.2.1 Mô hình dàn anten ULA Anten ULA hay dàn chấn tử đồng pha ULA là loại anten mảng đơn giản nhất. Dạng hình học của dàn anten này đƣợc biểu diễn qua hình vẽ sau: … 0 Y X d 1 2 M-1 θ Hình 1: Mô hình cấu trúc anten ULA Cấu trúc của dàn anten ULA bao gồm M chấn tử đặt song song trên cùng một trục thẳng với cùng khoảng cách d. Mỗi phần tử trong hệ anten đóng vai trò là một nguồn đẳng hƣớng. Các chấn tử này hoạt động cùng pha với nhau tạo nên một hƣớng bức xạ duy nhất và để cho tín hiệu tới bộ xử lí giữ nguyên đƣợc pha và biên độ so với tín hiệu tới hệ anten. 1.2.2 Tín hiệu thu được sau khi qua dàn anten ULA Nhƣ vậy, với một nguồn sóng đến, tín hiệu nhận đƣợc tại phần tử 0 (phần tử tham chiếu) đƣợc biểu diễn dƣới dạng: 0 ( ) ( ) os(2 ( ) ) c u t a t c f t t       (1.1) với: a(t): biên độ của tín hiệu f c : tần số sóng mang   t  : thành phần mang tin  : pha của tín hiệu. Giả sử rằng khoảng cách từ anten tới nguồn tín hiệu r >> d. Do vậy tia sóng ở [...]... (1.35) 1.3.3 Thuật toán MUSIC với dàn anten UCA Thuật toán MUSIC trong dàn anten UCA cũng được áp dụng giống như trong dàn anten ULA 1.3.4 Mô phỏng thuật toán MUSIC với dàn anten UCA Nhằm so sánh, đánh giá hiệu quả thuật toán MUSIC với các cấu trúc dàn anten ULA và UCA, việc mô phỏng đƣợc thực hiện với cùng thông số về môi trƣờng và nguồn sóng với dàn anten ULA - Mô hình hệ mô phỏng 11 Khóa luận tốt nghiệp... thuật toán MUSIC với dàn anten ULA trên Matlab - Mô hình hệ mô phỏng 8 Khóa luận tốt nghiệp Chƣơng 3 Thuật toán đƣợc mô phỏng với dàn anten ULA với các thông số của nguồn cũng nhƣ của hệ anten đƣợc cho trong bảng: Nhận xét: Kết quả mô phỏng thuật toán MUSIC với hệ anten ULA xác định đƣợc 8 nguồn sóng tới tại các góc lần lƣợt là [20 0, 400, 600, 800, 2800, 3000, 3200, 3400] So với thông số các nguồn... dẫn đến việc, khi xây dựng một hệ ULA cần phải thiết kế vị trí đặt hệ thích hợp để có thể che khuất nửa mặt phẳng mà hệ không phân biệt đƣợc, hoặc phải kết hợp các hệ ULA lại sao cho có thể tìm đƣợc vị trí tín hiệu cần nhận biết 1.3 Thuật toán MUSIC đối với dàn anten UCA 1.3.1 Mô hình dàn anten UCA Khác với cấu trúc anten ULA, cấu trúc anten sắp xếp theo đƣờng tròn cách đều UCA đƣợc nghiên cứu từ lâu,... anten 1    j d cos  e  A        j  M 1d cos  e    (1.11) si (t ) là các tín hiệu tới dàn anten với i  1  D n(t ) là tín hiệu nhiễu M chiều Công thức thu gọn của (1.10): U  A S  N (1.12) Đây cũng là công thức thể hiện tín hiệu thu nhận đƣợc sau khi qua dàn anten 1.2.3 Thuật toán MUSIC với dàn anten ULA Thuật toán MUSIC thực hiện việc tìm hƣớng sóng tới DOA dựa trên các. .. kit DSP Các bƣớc xây dựng thuật toán trên nền tảng DSP bao gồm 4 bƣớc đƣợc mô tả bởi lƣu đồ dƣới đây: Thiết kế thuật toán Lập trình thuật toán Thực thi thuật toán trên Kit Phân tích, sửa lỗi, tối ƣu Hình 14: Mô hình xây dựng thuật toán 3.1 Thiết kế thuật toán Cụ thể về thuật toán MUSIC đã đƣơc nghiên cứu ở chƣơng đầu Tuy nhiên khi nhúng vào DSP, ta phải làm sao cho chƣơng trình không quá nặng nề với chip... 1.3.2 Tín hiệu thu được sau khi qua dàn anten ULA So với hệ anten mảng thẳng ULA, hệ anten mảng tròn chỉ khác cách bố trí anten trong không gian Do vậy dạng tín hiệu nhận đƣợc tại mỗi anten riêng lẻ cũng nhƣ cấu trúc tín hiệu ở đầu vảo hệ anten là giống nhau, chỉ khác biệt ở thành phần 10 Khóa luận tốt nghiệp Chƣơng 3 vectơ lái: Vectơ lái của dàn anten ULA: A i ULA 1    j d cosi  e    ... ngôn ngữ C STT Các bƣớc thực hiện 1 Code Matlab Khó khăn Code C Thu nhân tín hiệu Có sẵn trƣờng Không có vào số phức Phải xây dựng trƣờng số phức với các phép tính thông thƣờng và các phép biến đổi ma trận 2 Tính ma trận tự Thuận tiện Sử dụng cấu Các biến xử lý đều là tƣơng quan trong việc tính trúc mảng với các ma trận kích cỡ toán ma trận, có thể thực hiện bƣớc này trong 1 dòng lệnh các giá trị lớn,... tông các phần tử ngoài đường chéo của Ruu O=O*Ruu Đúng END Giá trị riêng là phần tử đường chéo của Ruu Vector riêng là cột của O Sai Tổng >độ chính xác Hình 15: Lưu đồ thuật toán Jacobi 3.3 Thực thi thuật toán trên kit Với mục đích tập trung vào thuật toán và khả năng thực thi của thuật toán trên kit DSP, khóa luận không đi sâu vào thực hiện các phần cứng khác của hệ thống, ví dụ nhƣ hệ thống anten. .. phỏng thuật toán cho ta giá trị chính xác 4 góc tới là [200, 400, 600, 800] Ngoài ra trong phổ không gian còn xuất hiện thêm 4 góc tới khác ở các hƣớng [280 0, 3000, 3200, 3400] Nguyên nhân của hiện tƣợng này là do phƣơng trình biểu diễn vectơ lái hay do cấu trúc hình học của hệ anten Từ công thức tính vectơ lái của mảng anten (1.8) ta có: A(i )  e j d cos(i ) (1.30) Trong cấu trúc dàn anten ULA, ... Lại do cấu trúc dàn anten ULA có các phần tử sắp xếp cách đều theo đƣờng thẳng nên vectơ lái của mảng với các góc θ và (–θ) là nhƣ nhau, dẫn đến kết quả phổ không gian đánh giá DOA ở góc θ cũng giống với góc (–θ) Vấn đề này đƣợc thể hiện rõ ràng hơn ở hình sau: S1 Y θ 1 2 -θ X S2 Hình 4: Hệ anten ULA trong trường hợp hai tín hiệu đối xứng 9 Khóa luận tốt nghiệp Chƣơng 3 Ví dụ nhƣ ở hình trên, hệ anten . cấu trúc mảng khác nhau, từ đó tìm các giải thuật để triển khai thuật toán trên kit DSP. Chi tiết khóa luận bao gồm 5 chƣơng: Chƣơng 1: Thực hiện thuật toán MUSIC với các cấu trúc anten mảng. ULA 3 1.2.3 Thuật toán MUSIC với dàn anten ULA 5 1.2.4 Mô phỏng thuật toán MUSIC với dàn anten ULA trên Matlab 8 1.3 Thuật toán MUSIC đối với dàn anten UCA 10 1.3.1 Mô hình dàn anten UCA 10. dàn anten UCA cũng được áp dụng giống như trong dàn anten ULA. 1.3.4 Mô phỏng thuật toán MUSIC với dàn anten UCA Nhằm so sánh, đánh giá hiệu quả thuật toán MUSIC với các cấu trúc dàn anten ULA