Mô phỏng các hệ thống thông tin vô tuyến sử dụng Matlab - Trần Xuân Nam

Hướng dẫn mô phỏng các hệ thống thông tin vô tuyến bằng MATLAB. Bài viết đi từ lý thuyết cơ bản đến ví dụ thực tế kèm code mẫu chi tiết.

• Tác giả: Trần Xuân Nam
• Trường học: Đại học Kỹ thuật Lê Quí Đôn
• Chuyên ngành: Vô tuyến điện tử
• Đề tài: Mô phỏng các hệ thống thông tin vô tuyến sử dụng matlab
• Loại tài liệu: Giáo trình
• Địa điểm: Hà Nội

Mục lục chi tiết

• 1. Chương 1: Giới thiệu Matlab
• 1.1. Matlab là gì?
• 1.2. Khởi động và Thoát khỏi MATLAB
• 1.3. Làm việc với MATLAB Desktop
• 1.4. Các lệnh MATLAB cơ bản
• 1.5. Các ký hiệu đặc biệt
• 1.6. Tài liệu tham khảo
• 2. Chương 2: Tính toán và Lập trình sử dụng Matlab
• 2.1. Các phép tính số học
• 2.2. Các toán tử so sánh
• 2.3. Các toán tử logic
• 2.4. Vector và Ma trận
• 2.4.1. Tạo vector và ma trận
• 2.4.2. Các phép toán đối với vector và ma trận
• 2.5. Lập trình với Matlab
• 2.5.1. Điều khiển luồng (flow control)
• 2.5.2. Tạo chương trình MATLAB bằng tệp
• 2.6. Sử dụng đồ hoạ trong MATLAB
• 2.7. Tài liệu tham khảo
• 3. Chương 3: Lý thuyết mô phỏng
• 3.1. Vai trò của mô phỏng
• 3.2. Mô phỏng vs.
• 3.2.1. Sơ đồ truyền dẫn số qua kênh AWGN
• 3.2.2. Sơ đồ truyền dẫn số qua kênh AWGN sử dụng các bộ lọc và KĐCS phi tuyến
• 3.2.3. Hệ thống truyền dẫn qua kênh thông tin vệ tinh
• 3.3. Xây dựng mô hình mô phỏng
• 3.4. Các phương pháp mô phỏng
• 3.5. BER vs Xác suất lỗi bit
• 3.6. Vai trò của mô phỏng
• 3.7. Tính toán quĩ tuyến và mô phỏng
• 3.8. Các tham số đánh giá phẩm chất hệ thống
• 3.9. Kiểm định mô hình Mô phỏng
• 3.10. Năng lượng và Công suất tín hiệu
• 3.11. Mô phỏng Monte-Carlo trong Truyền dẫn Số
• 4. Chương 4: Kênh thông tin vô tuyến
• 4.1. Kênh tạp âm AWGN
• 4.1.1. Tạp âm AWGN
• 4.1.2. Mô phỏng tạp âm AWGN
• 4.1.3. Mô phỏng truyền dẫn qua kênh AWGN
• 4.2. Kênh pha-đinh
• 4.2.1. Mô hình toán học của pha-đinh
• 4.2.2. Ảnh hưởng của chuyển động của MS
• 4.2.3. Hậu quả của truyền sóng pha-đinh đa đường
• 4.3. Kênh pha-đinh Rayleigh
• 4.4. Mô phỏng pha-đinh Rayleigh
• 4.4.1. Đặc tính thống kê
• 4.5. Tài liệu tham khảo
• 5. Chương 5: Điều chế số
• 5.1. Điều chế pha sóng mang
• 6. Chương 6: KỸ THUẬT THU PHÁT PHÂN TẬP KHÔNG GIAN-THỜI GIAN
• 6.1. Các phương pháp phân tập
• 6.1.1. Phân tập thời gian
• 6.1.2. Phân tập tần số
• 6.1.3. Phân tập phân cực
• 6.1.4. Phân tập không gian
• 6.2. Kỹ thuật kết hợp phân tập không gian thu
• 6.2.1. Mô hình tín hiệu
• 6.2.2. Kết hợp chọn lọc (Selection Combining)
• 6.2.3. Kết hợp tỷ lệ tối đa (Maximal Ratio Combining)
• 6.2.4. Kết hợp đồng độ lợi (Equal Gain Combining)
• 6.2.5. Kết hợp phân tập thu và tách sóng MLD
• 6.3. Kỹ thuật kết hợp phân tập không gian phát
• 6.3.1. Phân tập phát tỉ lệ tối đa (MRT)
• 6.3.2. Phân tập phát giữ chậm
• 6.3.3. Phân tập phát không gian-thời gian
• 6.4. Tài liệu tham khảo
• 7. Chương 7: CÁC HỆ THỐNG MIMO
• 7.1. Mô hình kênh MIMO
• 7.2. Dung lượng kênh truyền MIMO
• 7.2.1. Dung lượng kênh truyền cố định
• 7.2.2. Dung lượng kênh truyền Rayleigh pha-đinh
• 7.3. Các phương pháp truyền dẫn trên kênh truyền MIMO
• 7.4. Ghép kênh theo không gian
• 7.5. Các bộ tách tín hiệu tuyến tính
• 7.5.1. Bộ tách tín hiệu ZF
• 7.5.2. Bộ tách tín hiệu MMSE
• 7.5.3. Các tham số phẩm chất bộ tách tín hiệu tuyến tính
• 7.6. Các bộ tách tín hiệu phi tuyến
• 7.6.1. Bộ tách tín hiệu QRD
• 7.6.2. Bộ tách tín hiệu V-BLAST
• 7.6.3. Bộ tách tín hiệu có trợ giúp của phương pháp rút gọn cơ sở dàn
• 7.6.4. Bộ tách tín hiệu MLD
• 7.6.5. Bộ tách tín hiệu hình cầu (sphere detector)
• 7.7. Tài liệu tham khảo
• 8. Chương 8: MÃ KHÔNG GIAN-THỜI GIAN
• 8.2. Mã khối không gian-thời gian
• 8.2.1. Mã STBC cho tập tín hiệu thực
• 8.2.2. Mã STBC cho tập tín hiệu phức
• 8.3. Mã lưới không gian-thời gian
• 8.4. Mã không gian-thời gian cho các hệ thống đa người dùng
• 8.5. Tài liệu tham khảo

Tóm tắt

I. Vai trò của mô phỏng hệ thống thông tin vô tuyến là gì

Trong lĩnh vực viễn thông hiện đại, việc thiết kế và triển khai các hệ thống truyền dẫn ngày càng phức tạp. Các yêu cầu về tốc độ dữ liệu cao, băng tần và công suất hạn chế đòi hỏi việc áp dụng các kỹ thuật tiên tiến như điều chế bậc cao, mã hóa sửa sai và các thuật toán xử lý tín hiệu phức tạp. Đối với các hệ thống đơn giản, phương pháp giải tích toán học có thể đánh giá chính xác hiệu năng. Tuy nhiên, khi đối mặt với các kênh truyền thực tế như kênh di động tế bào, chịu ảnh hưởng của pha-đinh đa đường và nhiễu, việc phân tích giải tích trở nên bất khả thi. Đây chính là lúc mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB phát huy vai trò không thể thay thế. Mô phỏng, đặc biệt là mô phỏng Monte Carlo, cung cấp một 'phòng thí nghiệm ảo' cho phép các kỹ sư và nhà nghiên cứu kiểm tra, đánh giá và tối ưu hóa hệ thống trước khi chế tạo phần cứng. Công cụ mạnh mẽ như MATLAB và Simulink, đặc biệt là Communications Toolbox, cung cấp một môi trường tích hợp để xây dựng mô hình từ các khối cơ bản như nguồn phát, bộ điều chế, kênh truyền, đến khối thu và bộ giải điều chế. Thay vì giải các phương trình toán học phức tạp, mô phỏng cho phép quan sát trực quan các tham số quan trọng như dạng sóng, phổ tín hiệu, đồ thị chòm sao, và quan trọng nhất là tỷ lệ lỗi bit (BER). Việc này không chỉ tiết kiệm thời gian, chi phí mà còn cho phép thử nghiệm các kịch bản khác nhau, từ đó tìm ra cấu hình hệ thống tối ưu. Theo tài liệu của Trần Xuân Nam, 'động cơ quan trọng cho mô phỏng là do mô phỏng là công cụ quí giá cho phép tìm hiểu sâu về hoạt động của hệ thống'. Việc xây dựng một chương trình mô phỏng cho hệ thống cơ bản và kiểm nghiệm với kết quả lý thuyết là bước đệm vững chắc để mở rộng sang các hệ thống phức tạp hơn như MIMO hay OFDM.

1.1. Phân biệt giữa phương pháp giải tích và mô phỏng

Phương pháp giải tích dựa trên việc xây dựng các mô hình toán học và công thức để tính toán trực tiếp các chỉ số hiệu năng của hệ thống, ví dụ như xác suất lỗi bit. Phương pháp này rất hiệu quả cho các hệ thống đơn giản, tuyến tính, hoạt động trên kênh tạp âm AWGN. Ví dụ, xác suất lỗi của hệ thống BPSK trên kênh AWGN có thể được tính chính xác bằng hàm Q Gauss. Tuy nhiên, khi hệ thống trở nên phức tạp với các thành phần phi tuyến (bộ khuếch đại công suất) hoặc kênh truyền biến đổi theo thời gian (kênh Rayleigh), việc tìm ra một công thức giải tích là cực kỳ khó khăn hoặc không thể. Ngược lại, phương pháp mô phỏng không cố gắng giải một phương trình duy nhất. Thay vào đó, nó tái tạo hoạt động của hệ thống bằng cách cho một chuỗi dữ liệu lớn đi qua mô hình số hóa của từng khối chức năng. Bằng cách đếm số lỗi xảy ra sau khi truyền một lượng lớn bit (ví dụ, hàng triệu bit), mô phỏng ước tính được các chỉ số hiệu năng như BER. Phương pháp này có tính tổng quát cao, áp dụng được cho mọi loại hệ thống, dù phức tạp đến đâu.

1.2. Giới thiệu về MATLAB và Communications Toolbox

MATLAB (Matrix Laboratory) là một môi trường tính toán số và ngôn ngữ lập trình bậc cao, được tối ưu hóa cho các phép toán ma trận, vẽ đồ thị và thực thi thuật toán. Trong lĩnh vực viễn thông, MATLAB là công cụ không thể thiếu cho việc xử lý tín hiệu số và mô phỏng. Sức mạnh của nó được mở rộng đáng kể nhờ các 'toolbox' chuyên dụng. Communications Toolbox là một bộ sưu tập các hàm và ứng dụng được xây dựng sẵn, chuyên dùng để thiết kế và mô phỏng lớp vật lý của các hệ thống viễn thông. Nó cung cấp đầy đủ các khối chức năng cần thiết: từ việc tạo nguồn bit, các bộ điều chế số (BPSK, QPSK, 16-QAM), các mô hình kênh truyền (AWGN, Rayleigh, Rician), bộ mã hóa kênh, cho đến các thuật toán thu và đồng bộ. Việc sử dụng toolbox này giúp giảm đáng kể thời gian lập trình, cho phép người dùng tập trung vào việc phân tích và thiết kế thuật toán thay vì phải xây dựng lại các hàm cơ bản từ đầu. Đây là công cụ lý tưởng cho sinh viên làm đồ án viễn thông và các nhà nghiên cứu phát triển công nghệ mới.

II. Cách mô hình hóa các kênh thông tin vô tuyến phổ biến

Trái tim của việc mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB nằm ở việc xây dựng chính xác mô hình kênh truyền. Kênh truyền là môi trường mà tín hiệu đi qua từ máy phát đến máy thu, và nó là nguyên nhân chính gây ra suy hao và lỗi. Mỗi môi trường truyền dẫn lại có những đặc tính riêng, đòi hỏi các mô hình toán học và thuật toán mô phỏng khác nhau. Mô hình đơn giản nhất là kênh tạp âm trắng cộng Gauss (kênh tạp âm AWGN), đặc trưng cho tạp âm nhiệt không thể tránh khỏi trong các mạch điện tử. Trong MATLAB, kênh này thường được mô phỏng bằng hàm randn để tạo ra các số ngẫu nhiên theo phân bố chuẩn. Tuy nhiên, trong truyền thông di động, mô hình phức tạp hơn là kênh pha-đinh, gây ra bởi hiện tượng truyền sóng đa đường (multipath). Khi tín hiệu từ máy phát đến máy thu qua nhiều đường khác nhau (phản xạ từ nhà cửa, cây cối), các tín hiệu này có biên độ, pha và độ trễ khác nhau, gây ra hiện tượng giao thoa tăng cường hoặc triệt tiêu tại máy thu. Kênh pha-đinh Rayleigh là mô hình phổ biến khi không có đường truyền thẳng (Non-Line-of-Sight) giữa phát và thu. Ngược lại, kênh Rician được sử dụng khi có một đường truyền thẳng chiếm ưu thế. Việc mô phỏng chính xác các kênh này là yếu tố quyết định độ tin cậy của kết quả mô phỏng, giúp đánh giá đúng hiệu năng của các kỹ thuật như điều chế số hay các hệ thống MIMO và 5G.

2.1. Mô phỏng kênh tạp âm AWGN Additive White Gaussian Noise

Kênh tạp âm AWGN là mô hình kênh cơ bản nhất trong thông tin vô tuyến, mô tả tạp âm có mật độ phổ công suất không đổi trên toàn băng thông và có biên độ tuân theo phân bố Gauss (phân bố chuẩn). Tạp âm này có nguồn gốc chủ yếu từ chuyển động nhiệt của các electron trong các linh kiện điện tử. Trong mô phỏng, tín hiệu thu y được biểu diễn đơn giản bằng y = s + n, trong đó s là tín hiệu phát và n là mẫu tạp âm AWGN. Để mô phỏng tạp âm n trong MATLAB, người ta thường sử dụng hàm randn, hàm này tạo ra các số ngẫu nhiên với giá trị trung bình bằng 0 và phương sai bằng 1. Để có được tạp âm với phương sai mong muốn (liên quan đến tỷ số tín hiệu trên tạp âm - SNR), mẫu tạp âm từ randn sẽ được nhân với độ lệch chuẩn sigma. Ví dụ, để tạo tạp âm phức với công suất No, câu lệnh có thể là: n = sqrt(No/2) * (randn(1,N) + 1j*randn(1,N)). Việc mô phỏng trên kênh AWGN là bước đầu tiên và bắt buộc khi đánh giá bất kỳ hệ thống truyền thông nào.

2.2. Phân tích và mô phỏng kênh pha đinh Rayleigh

Kênh pha-đinh Rayleigh mô tả hiện tượng tín hiệu bị biến đổi ngẫu nhiên về biên độ và pha khi không có đường truyền thẳng (NLOS) giữa máy phát và máy thu. Tín hiệu thu là tổng hợp của rất nhiều tín hiệu phản xạ. Theo định lý giới hạn trung tâm, khi số lượng các đường truyền là lớn, thành phần đồng pha và vuông pha của đáp ứng kênh sẽ tuân theo phân bố Gauss độc lập, có giá trị trung bình bằng 0. Kết quả là biên độ của kênh sẽ tuân theo phân bố Rayleigh. Hiện tượng này gây ra sự suy giảm tín hiệu sâu và nhanh chóng, là một trong những thách thức lớn nhất của truyền thông di động. Để mô phỏng kênh Rayleigh, một phương pháp phổ biến là mô hình Jakes, dựa trên việc tổng hợp nhiều sóng hình sin với các tần số Doppler và pha ngẫu nhiên. Trong MATLAB, Communications Toolbox cung cấp các đối tượng kênh comm.RayleighChannel giúp việc mô phỏng trở nên đơn giản, người dùng chỉ cần định nghĩa các tham số như tần số Doppler tối đa và độ trễ của các đường truyền.

III. Hướng dẫn mô phỏng điều chế số cơ bản bằng MATLAB

Điều chế số là quá trình ánh xạ dữ liệu số (bit 0 và 1) thành các dạng sóng tương tự để truyền qua kênh vô tuyến. Việc mô phỏng các hệ thống thông tin vô tuyến không thể thiếu bước mô phỏng quá trình này và quá trình ngược lại là giải điều chế. MATLAB cung cấp môi trường lý tưởng để thực hiện và phân tích các kỹ thuật điều chế số từ cơ bản đến phức tạp. Các kỹ thuật phổ biến nhất bao gồm BPSK (Binary Phase Shift Keying), QPSK (Quadrature Phase Shift Keying), và 16-QAM (16-ary Quadrature Amplitude Modulation). Quá trình mô phỏng thường bắt đầu bằng việc tạo ra một chuỗi bit ngẫu nhiên. Sau đó, các bit này được nhóm lại và ánh xạ vào các ký hiệu phức (complex symbols) trên đồ thị chòm sao. Ví dụ, trong BPSK, bit '0' được ánh xạ thành '+1' và bit '1' thành '-1'. Trong QPSK, hai bit được nhóm lại và ánh xạ vào một trong bốn điểm {1+j, 1-j, -1+j, -1-j} (sau khi chuẩn hóa năng lượng). Các ký hiệu này sau đó được cho đi qua mô hình kênh (ví dụ kênh AWGN). Tại phía thu, quá trình giải điều chế bao gồm việc quyết định ký hiệu nào đã được phát dựa trên tín hiệu nhận được (đã bị nhiễu). Bằng cách so sánh chuỗi bit đã giải điều chế với chuỗi bit gốc, chúng ta có thể tính toán tỷ lệ lỗi bit (BER), một thước đo hiệu năng quan trọng. Toàn bộ quy trình này có thể được lập trình bằng các hàm trong MATLAB hoặc sử dụng các khối có sẵn trong Simulink.

3.1. Các bước tạo tín hiệu điều chế BPSK QPSK 16 QAM

Để tạo tín hiệu điều chế số trong MATLAB, quy trình chung như sau. Đầu tiên, tạo chuỗi dữ liệu nhị phân ngẫu nhiên, ví dụ: data = randi([0 1], N, 1). Đối với BPSK, phép ánh xạ đơn giản là s = 1 - 2*data, chuyển {0, 1} thành {+1, -1}. Đối với QPSK và 16-QAM, cần sử dụng các hàm điều chế có sẵn. Ví dụ, để điều chế QPSK, ta dùng hàm pskmod: s_qpsk = pskmod(data, 4, pi/4). Tham số M=4 chỉ ra đây là 4-PSK (tức QPSK) và pi/4 là độ lệch pha ban đầu để tạo chòm sao chuẩn. Tương tự, với 16-QAM, ta dùng hàm qammod: s_16qam = qammod(data, 16, 'UnitAveragePower', true). Tùy chọn 'UnitAveragePower' giúp chuẩn hóa công suất trung bình của các ký hiệu về 1, thuận tiện cho việc tính toán SNR và BER sau này. Các hàm này đều nằm trong Communications Toolbox.

3.2. Phương pháp tính Tỷ lệ lỗi bit BER qua mô phỏng

Tính Tỷ lệ lỗi bit (BER) là mục tiêu cuối cùng của hầu hết các chương trình mô phỏng hiệu năng. Phương pháp được sử dụng là mô phỏng Monte-Carlo. Ý tưởng rất đơn giản: truyền một số lượng lớn bit N qua mô hình hệ thống và đếm số bit bị lỗi Ne ở phía thu. Ước tính BER được tính bằng BER_est = Ne / N. Để kết quả có ý nghĩa thống kê, N phải đủ lớn. Một quy tắc kinh nghiệm là cần ghi nhận ít nhất 10-100 lỗi. Điều này có nghĩa là để đo được BER = 10^-5, cần truyền ít nhất 1-10 triệu bit. Trong MATLAB, sau khi có chuỗi bit gốc data và chuỗi bit thu được sau giải điều chế data_rx, ta có thể tính số lỗi bằng [number_of_errors, bit_error_rate] = biterr(data, data_rx). Vòng lặp for thường được sử dụng để mô phỏng hệ thống ở nhiều mức SNR khác nhau, từ đó vẽ được đường cong BER theo SNR, là đồ thị đặc trưng cho hiệu năng của mọi hệ thống viễn thông.

IV. Kỹ thuật mô phỏng hệ thống MIMO OFDM cho 5G và LTE

Các hệ thống truyền thông hiện đại như LTE và 5G đều dựa trên hai công nghệ nền tảng là MIMO (Multiple-Input Multiple-Output) và OFDM (Orthogonal Frequency Division Multiplexing). Việc mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB cho các công nghệ này đòi hỏi kiến thức sâu hơn về xử lý tín hiệu số. OFDM là một kỹ thuật điều chế đa sóng mang, chia một luồng dữ liệu tốc độ cao thành nhiều luồng tốc độ thấp hơn và truyền đồng thời trên các sóng mang con trực giao. Kỹ thuật này chống lại hiệu ứng pha-đinh chọn lọc tần số rất hiệu quả. Mô phỏng OFDM trong MATLAB bao gồm các bước chính: chuyển đổi chuỗi-song song, biến đổi Fourier ngược nhanh (IFFT), thêm tiền tố chu kỳ (Cyclic Prefix), đi qua kênh truyền, và thực hiện các quá trình ngược lại ở phía thu (loại bỏ CP, FFT). MIMO sử dụng nhiều anten ở cả phía phát và phía thu để cải thiện đáng kể hiệu năng hệ thống, hoặc bằng cách tăng độ tin cậy (phân tập không gian) hoặc tăng tốc độ dữ liệu (ghép kênh không gian). Mô phỏng một hệ thống MIMO yêu cầu mô hình hóa ma trận kênh H giữa các anten phát và thu. Tín hiệu thu được biểu diễn qua phương trình ma trận y = Hx + n. Việc mô phỏng kết hợp MIMO-OFDM là một chủ đề phức tạp nhưng vô cùng quan trọng, và MATLAB cung cấp các công cụ mạnh mẽ để thực hiện, ví dụ như các hàm để tạo ma trận kênh fading hoặc các thuật toán tách tín hiệu MIMO như Zero-Forcing (ZF) và MMSE.

4.1. Xây dựng mô hình hệ thống MIMO trong Simulink

Simulink, một môi trường đồ họa của MATLAB, rất phù hợp để xây dựng và mô phỏng các hệ thống MIMO phức tạp. Thay vì viết hàng trăm dòng code MATLAB, người dùng có thể kéo-thả các khối chức năng từ thư viện của Communications Toolbox và kết nối chúng lại với nhau. Một mô hình hệ thống MIMO 2x2 trong Simulink sẽ bao gồm các khối chính: khối 'Bernoulli Binary Generator' để tạo hai luồng dữ liệu độc lập; các khối 'QPSK Modulator' để điều chế; khối 'MIMO Fading Channel' để mô phỏng kênh truyền giữa 2 anten phát và 2 anten thu; khối 'AWGN Channel'; và một khối 'MIMO Detector' (ví dụ như ZF hoặc MMSE) để tách hai luồng tín hiệu tại phía thu. Simulink cho phép trực quan hóa tín hiệu tại bất kỳ điểm nào trong hệ thống bằng các khối 'Constellation Diagram' hay 'Error Rate Calculation'. Cách tiếp cận này giúp người dùng dễ dàng hiểu được luồng tín hiệu và tương tác giữa các thành phần, rất hữu ích cho mục đích giảng dạy và nghiên cứu ban đầu.

4.2. Mô phỏng hệ thống đa sóng mang OFDM cho LTE 5G

Mô phỏng một hệ thống OFDM trong MATLAB là một bài toán kinh điển trong xử lý tín hiệu số. Quá trình mô phỏng phía phát bao gồm: (1) Điều chế số (ví dụ QPSK) cho dữ liệu trên từng sóng mang con. (2) Áp dụng biến đổi IFFT (ifft trong MATLAB) lên vector các ký hiệu để chuyển từ miền tần số sang miền thời gian, tạo ra tín hiệu OFDM. (3) Thêm Cyclic Prefix (CP) bằng cách sao chép một phần cuối của ký hiệu OFDM và chèn vào đầu. CP giúp loại bỏ nhiễu giữa các ký hiệu (ISI) và nhiễu giữa các sóng mang con (ICI). Tại phía thu, quá trình ngược lại được thực hiện: (1) Loại bỏ CP. (2) Áp dụng FFT (fft trong MATLAB) để chuyển tín hiệu về miền tần số. (3) Cân bằng kênh để loại bỏ ảnh hưởng của kênh truyền. (4) Giải điều chế các ký hiệu để khôi phục dữ liệu. MATLAB cung cấp các đối tượng như comm.OFDMModulator và comm.OFDMDemodulator để tự động hóa các bước này, giúp việc mô phỏng hệ thống LTE hay 5G trở nên dễ dàng hơn.

V. Ứng dụng kết quả mô phỏng cho đồ án và luận văn

Kết quả từ việc mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB không chỉ là những con số, mà chúng là những bằng chứng khoa học quan trọng, đặc biệt trong môi trường học thuật như khi làm đồ án viễn thông hay luận văn mô phỏng. Một báo cáo hay luận văn thành công đòi hỏi việc trình bày kết quả một cách rõ ràng, trực quan và có tính thuyết phục. Các đồ thị là công cụ không thể thiếu. Đồ thị quan trọng nhất là đường cong BER theo SNR (hoặc Eb/No). Đồ thị này cho thấy hiệu năng của hệ thống thay đổi như thế nào khi chất lượng kênh truyền thay đổi. Việc so sánh đường cong BER của mô phỏng với đường cong lý thuyết (nếu có) là một bước quan trọng để kiểm chứng tính đúng đắn của chương trình mô phỏng. Một loại đồ thị hữu ích khác là đồ thị chòm sao (Constellation Diagram) tại phía thu. Đồ thị này cho thấy sự phân tán của các điểm ký hiệu nhận được xung quanh các vị trí lý tưởng do ảnh hưởng của tạp âm và pha-đinh. Độ phân tán càng lớn, hiệu năng hệ thống càng kém. Ngoài ra, các đồ thị về phổ tín hiệu, đáp ứng kênh, hay biểu đồ histogram cũng cung cấp những cái nhìn sâu sắc về hoạt động của hệ thống. Sử dụng thành thạo các hàm vẽ đồ thị của MATLAB như semilogy, plot, scatterplot là kỹ năng cần thiết để trình bày kết quả nghiên cứu một cách chuyên nghiệp. Việc phân tích và diễn giải ý nghĩa đằng sau các đồ thị này chính là phần cốt lõi, thể hiện sự hiểu biết sâu sắc của người thực hiện về hệ thống viễn thông.

5.1. Cách trình bày và phân tích đồ thị chòm sao

Đồ thị chòm sao là một công cụ trực quan hóa mạnh mẽ, biểu diễn các ký hiệu tín hiệu thu được trên mặt phẳng phức. Trong MATLAB, hàm scatterplot được sử dụng để vẽ đồ thị này. Một đồ thị chòm sao lý tưởng (không có nhiễu) sẽ chỉ bao gồm các điểm ký hiệu nằm chính xác tại vị trí được định nghĩa bởi kỹ thuật điều chế số (ví dụ, 4 điểm cho QPSK). Khi có tạp âm AWGN, các điểm nhận được sẽ tạo thành các cụm mờ xung quanh các điểm lý tưởng. Hình dạng và kích thước của các cụm này cho biết mức độ nhiễu. Nếu kênh truyền có các khiếm khuyết khác như méo pha, các cụm điểm sẽ bị xoay đi so với vị trí gốc. Nếu có méo biên độ, khoảng cách từ các cụm đến gốc tọa độ sẽ bị thay đổi. Bằng cách quan sát đồ thị chòm sao, người kỹ sư có thể nhanh chóng chẩn đoán các vấn đề trong hệ thống mà không cần phân tích các con số phức tạp.

5.2. Viết code MATLAB và tối ưu hóa thời gian mô phỏng

Viết code MATLAB hiệu quả là một kỹ năng quan trọng để giảm thời gian mô phỏng, đặc biệt khi cần tính BER ở các mức rất thấp. Một trong những nguyên tắc quan trọng nhất là 'vector hóa' (vectorization). Thay vì sử dụng vòng lặp for để xử lý từng mẫu tín hiệu, hãy tận dụng khả năng của MATLAB để thực hiện các phép toán trên toàn bộ vector hoặc ma trận trong một lệnh duy nhất. Ví dụ, thay vì lặp qua từng phần tử để cộng tạp âm, chỉ cần dùng lệnh y = s + n;. Kỹ thuật này nhanh hơn đáng kể. Ngoài ra, nên khởi tạo trước (pre-allocation) kích thước của các ma trận lớn trước khi điền dữ liệu vào chúng để tránh việc MATLAB phải thay đổi kích thước bộ nhớ liên tục. Đối với các mô phỏng Monte Carlo rất dài, có thể sử dụng Parallel Computing Toolbox để chạy các mô phỏng cho các mức SNR khác nhau trên nhiều lõi CPU song song, giúp giảm đáng kể tổng thời gian chờ đợi. Việc viết code sạch, có chú thích rõ ràng cũng giúp cho việc gỡ lỗi và phát triển luận văn mô phỏng sau này.

Nội dung chính

Báo cáo Phân tích và Mô tả Thư mục Học thuật: Giáo trình "Mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB"

Tổng quan về giáo trình

Giáo trình Mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB (tên tiếng Anh: Simulation of Radio Communication Systems using Matlab) do TS. Trần Xuân Nam thuộc Bộ môn Thông tin, Khoa Vô tuyến điện tử – Đại học Kỹ thuật Lê Quí Đôn biên soạn. Trong chương trình đào tạo kỹ sư và thạc sĩ ngành Kỹ thuật Điện tử – Viễn thông và Kỹ thuật Vô tuyến, tài liệu này đóng vai trò là giáo trình chuyên ngành cho các học phần về Mô phỏng hệ thống truyền thông, Xử lý tín hiệu trong vô tuyến, và Kỹ thuật thông tin di động.

Mục tiêu học tập của giáo trình là trang bị cho người học phương pháp luận mô hình hóa hệ thống truyền dẫn số, kỹ năng tính toán và lập trình mô phỏng trên môi trường MATLAB, cùng khả năng phân tích, đánh giá phẩm chất hệ thống thông tin vô tuyến (tỉ số tín hiệu trên tạp âm SNR, tỉ số lỗi bit BER, tỉ số lỗi ký hiệu SER, tỉ số lỗi khung FER). Người học được rèn luyện năng lực kiểm chứng giữa kết quả mô phỏng số thông qua phương pháp Monte-Carlo và các kết quả phân tích giải tích lý thuyết.

Cấu trúc giáo trình gồm 8 chương, được tổ chức theo quy trình từ việc làm chủ công cụ tính toán ma trận, nắm bắt lý thuyết mô phỏng thống kê, mô hình hóa các kênh truyền sóng vô tuyến (AWGN, pha-đinh Rayleigh, Rice), cho đến việc mô phỏng các sơ đồ điều chế số và các kỹ thuật truyền dẫn đa ăng-ten tiên tiến (phân tập không gian - thời gian, hệ thống MIMO, mã không gian - thời gian). Điểm đặc sắc của giáo trình là sự kết hợp chặt chẽ giữa cơ sở toán học giải tích và mã nguồn lập trình cụ thể, giúp người học chuyển đổi các mô hình toán học trừu tượng thành các thuật toán mô phỏng số có thể thực thi và kiểm chứng định lượng trên máy tính.

---

Nội dung kiến thức cốt lõi

Các chương và chủ đề chính

Nội dung học thuật của giáo trình được triển khai qua 8 chương có tính tiếp nối logic chặt chẽ:

• Chương 1: Giới thiệu MATLAB – Khái quát về 5 thành phần cốt lõi của hệ thống MATLAB (Môi trường phát triển, Thư viện hàm toán học, Ngôn ngữ MATLAB, Đồ họa khoa học kỹ thuật, và Giao diện chương trình ứng dụng API C/Fortran); hướng dẫn tương tác với Command Window, Command History, Current Directory, Workspace, các lệnh điều khiển cơ bản (ls, pwd, who, clear, clc) và các ký hiệu cú pháp đặc biệt.
• Chương 2: Tính toán và Lập trình sử dụng MATLAB – Trình bày các phép tính số học, toán tử so sánh, toán tử logic (&, |, ~, xor), các phép toán đại số ma trận và vector (nhân ma trận *, nhân từng phần tử .*, chuyển vị liên hợp phức ', chuyển vị thuần túy .', nghịch đảo ma trận inv, tính định thức det); cấu trúc điều khiển luồng (if-elseif-else, switch-case, for, while, break); kỹ thuật xây dựng hàm (function) và tệp lệnh (script) .m; kỹ thuật đồ họa kỹ thuật (plot, semilogy, định dạng đường, dấu và nhãn trục).
• Chương 3: Lý thuyết mô phỏng – Phân tích vai trò của mô phỏng trong việc giải quyết các bài toán truyền thông phi tuyến và kênh pha-đinh phức tạp; phân định sự khác biệt giữa mô phỏng và phân tích giải tích thông qua 3 sơ đồ hệ thống (kênh AWGN tuyến tính, kênh qua bộ khuếch đại công suất phi tuyến và bộ lọc gây nhiễu xuyên ký tự ISI, hệ thống vệ tinh hai chặng uplink/downlink); giới thiệu phương pháp mô phỏng bán giải tích (semi-analytical); định nghĩa tỉ số lỗi bit BER so với xác suất lỗi $P_E$; quy trình thiết kế "từ trên xuống dưới" (Top-down) và "từ dưới lên trên" (Bottom-up); phương pháp tính toán quỹ tuyến (link budget); nguyên lý mô phỏng Monte-Carlo và quy tắc kích thước mẫu $N > 10/P(s)$.
• Chương 4: Kênh thông tin vô tuyến – Thiết lập mô hình toán học cho tạp âm trắng Gauss cộng (AWGN) với mật độ phổ công suất $N_0/2$ và hàm tự tương quan $R_n(\tau) = (N_0/2)\delta(\tau)$; mô tả cơ chế truyền sóng pha-đinh đa đường giữa trạm gốc (BS) và trạm di động (MS); dẫn xuất hiện tượng trải Doppler (Doppler spread, tần số Doppler cực đại $f_D = \frac{v}{c}f_c$) và trải trễ (delay spread, mô hình kênh lọc FIR chọn lọc tần số); thiết lập hàm mật độ xác suất cho phân bố Rayleigh và phân bố Rice (với hệ số Rice $K$ và hàm Bessel sửa đổi bậc 0 loại 1 $I_0$); thuật toán mô phỏng pha-đinh Rayleigh và Rice dựa trên phương pháp tổng các hàm sin (Sum-of-Sinusoids - SOS) của Jakes và các hiệu chỉnh của Pop & Beaulieu.
• Chương 5: Mô phỏng Monte-Carlo một số hệ thống vô tuyến điển hình – Hiện thực hóa sơ đồ mô phỏng Monte-Carlo cho hệ thống điều chế BPSK qua kênh AWGN, so sánh đường cong BER thực nghiệm với hàm sai số bù $\frac{1}{2}\text{erfc}(\sqrt{E_b/N_0})$; xây dựng quy trình mô phỏng hệ thống điều chế $M$-PSK qua kênh pha-đinh Rayleigh sử dụng bộ tách sóng đồng bộ (coherent detection) và thuật toán tách tín hiệu hợp lẽ tối đa (MLD).
• Chương 6: Kỹ thuật thu phát phân tập không gian - thời gian – Phân loại các phương pháp phân tập (thời gian, tần số, phân cực, không gian); phân tích các kỹ thuật kết hợp phân tập thu gồm: Kết hợp chọn lọc (Selection Combining - SC), Kết hợp tỷ lệ tối đa (Maximal Ratio Combining - MRC), Kết hợp đồng độ lợi (Equal Gain Combining - EGC); trình bày các kỹ thuật phân tập phát gồm: Phân tập phát tỷ lệ tối đa (MRT), phân tập phát giữ chậm, và kỹ thuật mã khối không gian - thời gian (STBC) Alamouti $2 \times 1$ và $2 \times 2$.
• Chương 7: Các hệ thống MIMO – Mô hình hóa kênh truyền MIMO (SISO, MISO, SIMO, MIMO); tính toán dung lượng kênh truyền cố định và kênh pha-đinh Rayleigh; kỹ thuật ghép kênh theo không gian (MIMO-SDM, MIMO-SVD); các cấu trúc tách tín hiệu tuyến tính (Zero-Forcing - ZF, Minimum Mean Square Error - MMSE) và phi tuyến (tách tín hiệu QR Decomposition - QRD, V-BLAST, tách tín hiệu có trợ giúp của phương pháp rút gọn cơ sở dàn theo thuật toán LLL, MLD, và bộ tách tín hiệu hình cầu - Sphere Detector).
• Chương 8: Mã không gian - thời gian – Cấu trúc mã khối không gian - thời gian (STBC) cho tập tín hiệu thực và phức; mã lưới không gian - thời gian (STTC); mã không gian - thời gian cho hệ thống đa người dùng.

Kiến thức nền tảng được xây dựng

Giáo trình củng cố và hệ thống hóa các nền tảng lý thuyết:

• Lý thuyết xác suất và quá trình ngẫu nhiên: Các phân bố thống kê Gauss, Rayleigh, Rice, hàm $Q$-Gauss, hàm sai số bù $\text{erfc}$, tính dừng rộng (WSSUS), hàm mật độ phổ công suất (PSD) và hàm tự tương quan.
• Đại số tuyến tính ứng dụng: Không gian vector, ma trận phức, phép chuyển vị liên hợp Hermitian, phép nghịch đảo giả ma trận, phân tích giá trị kỳ dị (SVD), lý thuyết cơ sở dàn (lattice basis reduction).
• Lý thuyết truyền thông số: Năng lượng bit $E_b$, năng lượng ký hiệu $E_s$, mật độ tạp âm $N_0$, chòm sao tín hiệu (constellation), bộ lọc phối hợp (matched filter), nhiễu xuyên ký tự (ISI) và nguyên lý tách sóng tối ưu MLD.

Kỹ năng phát triển

• Kỹ năng kỹ thuật (Technical skills): Viết mã nguồn MATLAB chuẩn xác, tối ưu hóa các phép toán ma trận/mảng thay cho vòng lặp, sử dụng các hàm chuyên dụng (norm, sign, angle, floor, biterr), tạo chuỗi bit giả ngẫu nhiên, sinh biến ngẫu nhiên Gauss phức.
• Kỹ năng phân tích (Analytical skills): Phân tích mối tương quan giữa độ phức tạp mô hình, sai số tính toán và thời gian thực thi; phân tích suy hao đường truyền và tính toán quỹ tuyến (link budget); đánh giá dung lượng kênh và độ lợi phân tập.
• Kỹ năng thực hành (Practical competencies): Xây dựng chương trình mô phỏng đầu-cuối (end-to-end simulation), vẽ và phân tích đồ thị đường cong BER/SER trên thang logarith (semilogy), kiểm định mô hình mô phỏng thông qua đối sánh với các biểu thức giải tích chuẩn mực.

---

Phương pháp giảng dạy và học tập

Giáo trình áp dụng phương pháp sư phạm tiếp cận "từ trên xuống dưới" (Top-down approach). Mỗi chủ đề bắt đầu bằng việc thiết lập mô hình toán học giải tích trừu tượng của hệ thống hoặc hiện tượng vật lý, tiếp theo là thiết kế lưu đồ thuật toán xử lý số, và cuối cùng là hiện thực hóa thuật toán bằng các đoạn mã nguồn MATLAB minh họa hoàn chỉnh.

Các bài tập và ví dụ thực hành được thiết kế gắn liền với từng phần lý thuyết:

• Thực hành lập trình cơ bản: Xây dựng hàm giải phương trình bậc hai (discriminant.m, quadEqn.m), tạo chuỗi ngẫu nhiên nhị phân, tính kỳ vọng $E_x$, phương sai $\sigma_x^2$ và độ lệch chuẩn.
• Thực hành mô phỏng kênh truyền: Viết script tạo tạp âm AWGN phức; tạo kênh pha-đinh Rayleigh theo mô hình Jakes với số nhánh tán xạ $L = 35$, tần số Doppler $f_D = 100\text{ Hz}$ và vẽ phân bố biên độ bằng lệnh hist; mô phỏng kênh pha-đinh Rice với các hệ số $K = 5$ và $K = 10$.
• Thực hành đánh giá phẩm chất hệ thống: Viết chương trình mô phỏng điều chế BPSK qua kênh AWGN (BPSKAWGN.m), ánh xạ và giải điều chế QPSK (QPSKMap.m, QPSKMod.m), tính tỉ số BER của hệ thống BPSK qua kênh pha-đinh Rayleigh và so sánh với công thức lý thuyết $P_{E} = \frac{1}{2}\left(1 - \sqrt{\frac{E_b/N_0}{1 + E_b/N_0}}\right)$.

Phương pháp đánh giá học phần dựa trên:

1. Độ chính xác của lưu đồ thuật toán và cấu trúc mã nguồn mô phỏng.
2. Khả năng giải thích các hiện tượng vật lý thông qua biểu đồ (sự suy giảm BER theo $E_b/N_0$, hiện tượng sàn lỗi do pha-đinh hoặc méo phi tuyến).
3. Mức độ hội tụ giữa kết quả thực nghiệm mô phỏng Monte-Carlo và đường cong phân tích giải tích.

Hướng dẫn tự học: Sinh viên thực hiện kiểm thử từng đoạn lệnh trực tiếp trên cửa sổ Command Window để hiểu rõ cơ chế biến đổi dữ liệu, tiếp đó đóng gói thành các hàm chuyên biệt .m, và mở rộng kích thước mẫu thử nghiệm $N$ để quan sát sự cân bằng giữa độ chính xác thống kê và thời gian chạy máy tính.

---

Điểm nổi bật và cập nhật

Giáo trình giải quyết hạn chế của các tài liệu lý thuyết truyền thống vốn chỉ tiếp cận qua phương pháp giải tích trên kênh AWGN tuyến tính lý tưởng. Giáo trình tập trung vào các tình huống thực tế không thể giải bằng giải tích đơn thuần:

• Xử lý tính phi tuyến và nhiễu phức tạp: Mô hình hóa bộ khuếch đại công suất phi tuyến (KĐCS/HPA) gây méo hài, méo điều chế lẫn nhau (intermodulation) và mở rộng phổ tín hiệu; mô hình tuyến truyền dẫn vệ tinh hai chặng với hai nguồn tạp âm độc lập trên đường lên (uplink) và đường xuống (downlink).
• Mô hình kênh truyền di động thực tế: Mô hình hóa sự dịch chuyển của máy di động MS gây hiệu ứng Doppler spread và hiện tượng truyền sóng đa đường gây delay spread, làm rõ sự chuyển đổi giữa kênh pha-đinh phẳng (flat fading) và kênh chọn lọc tần số (frequency selective fading) thông qua mô hình bộ lọc FIR.
• Tích hợp kỹ thuật truyền thông đa ăng-ten hiện đại: Cập nhật đầy đủ các kỹ thuật không gian - thời gian then chốt gồm cấu trúc Alamouti STBC, hệ thống ghép kênh không gian MIMO-SDM, cùng các thuật toán tách tín hiệu hiện đại như V-BLAST, kỹ thuật rút gọn cơ sở dàn LLL (Lattice Reduction) và giải mã hình cầu (Sphere Decoding).
• Cơ sở tài liệu tham khảo học thuật chuẩn mực: Nội dung được xây dựng dựa trên các công trình nghiên cứu và tài liệu chuyên khảo quốc tế uy tín như W. Tranter et al. (2004), B. Sklar (1988), Jakes (1994), Cavers (2000), Stüber (2001), Rice (1948), Aulin (1981), Pop & Beaulieu (2001), Zheng (2003).

---

Đối tượng sử dụng giáo trình

• Sinh viên đại học và học viên cao học: Sinh viên năm thứ 3, năm thứ 4 và học viên cao học thuộc các chuyên ngành Kỹ thuật Điện tử – Viễn thông, Kỹ thuật Vô tuyến điện tử, Kỹ thuật Mạng và Xử lý Tín hiệu.
• Điều kiện tiên quyết (Prerequisites):
  • Toán học: Đại số tuyến tính (ma trận, định thức, không gian vector), Xác suất thống kê và Quá trình ngẫu nhiên.
  • Chuyên ngành: Tín hiệu và Hệ thống, Lý thuyết Truyền thông số (các phương thức điều chế BPSK, QPSK, M-PSK, cấu trúc máy thu lọc phối hợp).
  • Lập trình: Kiến thức cơ bản về tư duy thuật toán và lập trình máy tính.
• Giảng viên: Sử dụng làm tài liệu giảng dạy chính khóa cho môn học Mô phỏng hệ thống thông tin, hoặc làm giáo trình tham khảo cho các học phần Thông tin vô tuyến, Thông tin di động; khai thác hệ thống ví dụ và bài tập mẫu để xây dựng các bài thí nghiệm, đồ án môn học.
• Kỹ sư và nghiên cứu viên: Kỹ sư R&D trong lĩnh vực thiết kế hệ thống viễn thông sử dụng tài liệu làm tài liệu tham khảo để xây dựng công cụ mô phỏng thẩm định tham số, kiểm tra quỹ tuyến và đánh giá thuật toán trước khi hiện thực hóa trên phần cứng.

---

Câu hỏi thường gặp

1. Giáo trình này phù hợp với ai?

Giáo trình được biên soạn phục vụ sinh viên đại học giai đoạn chuyên ngành, học viên cao học, giảng viên kỹ thuật viễn thông và các kỹ sư nghiên cứu cần công cụ tính toán, mô phỏng số để khảo sát các hệ thống truyền dẫn vô tuyến.

2. Cần chuẩn bị những kiến thức nền tảng nào trước khi học?

Người học cần có kiến thức về Đại số ma trận, Xác suất và Quá trình ngẫu nhiên, các nguyên lý cơ bản của Kỹ thuật Thông tin số (điều chế, giải điều chế, tạp âm AWGN) và kỹ năng sử dụng máy tính căn bản.

3. Điểm khác biệt của giáo trình này so với các giáo trình lý thuyết thông tin thuần túy là gì?

Khác với các tài liệu chỉ cung cấp công thức giải tích trên các kênh truyền lý tưởng, giáo trình này hướng dẫn chi tiết cách mô hình hóa thuật toán và lập trình mô phỏng trên MATLAB cho cả các hệ thống phức tạp (phi tuyến, pha-đinh đa đường, MIMO), đồng thời cung cấp phương pháp định lượng sai số mô phỏng Monte-Carlo so với lý thuyết.

4. Làm sao để tự học và nghiên cứu giáo trình hiệu quả?

Người học nên kết hợp song song giữa việc đọc hiểu lý thuyết kênh truyền/điều chế và trực tiếp gõ, thực thi các đoạn mã nguồn mẫu trong MATLAB. Sau đó, tiến hành thay đổi các tham số đầu vào (tần số Doppler $f_D$, số nhánh tán xạ $L$, hệ số Rice $K$, tỉ số $E_b/N_0$, số lượng ăng-ten thu/phát) và tự giải các bài tập lập trình mở rộng ở cuối mỗi chương.

5. Giáo trình có tài liệu bổ trợ nào kèm theo không?

Nội dung giáo trình đã tích hợp đầy đủ hệ thống mã nguồn mẫu minh họa (m-files) cho từng chủ đề, bảng tham số định dạng đồ thị, cấu trúc thuật toán chi tiết và hệ thống danh mục tài liệu tham khảo chuyên khảo quốc tế đi kèm cuối mỗi chương.

---

Kết luận

Giáo trình Mô phỏng các hệ thống thông tin vô tuyến sử dụng MATLAB của TS. Trần Xuân Nam cung cấp một hệ thống kiến thức hoàn chỉnh và chuẩn mực về phương pháp luận mô hình hóa và kỹ thuật mô phỏng số trong lĩnh vực truyền thông vô tuyến.

Lộ trình học tập đề xuất:

1. Nắm vững cú pháp tính toán ma trận và kỹ thuật lập trình tệp .m trên MATLAB (Chương 1, 2).
2. Nghiên cứu nguyên lý mô phỏng Monte-Carlo, phương pháp bán giải tích và quy tắc tính kích thước mẫu (Chương 3).
3. Xây dựng các bộ sinh kênh truyền vô tuyến AWGN, Rayleigh và Rice (Chương 4).
4. Thực hiện mô phỏng đánh giá BER cho các sơ đồ điều chế số cơ bản (Chương 5).
5. Nghiên cứu và mô phỏng các cấu trúc đa ăng-ten nâng cao gồm phân tập không gian - thời gian, hệ thống MIMO và mã hóa không gian - thời gian (Chương 6, 7, 8).

Tài liệu bổ trợ nghiên cứu mở rộng gồm các công trình chuyên khảo của W. Tranter et al. (Principles of Communication Systems Simulation with Wireless Applications), B. Sklar (Digital Communications), W. C. Jakes (Microwave Mobile Communications) và tài liệu hướng dẫn kỹ thuật từ The MathWorks Inc.

Trích đoạn nội dung tài liệu

Chương 1. Giới thiệu Matlab Window), Lịch sử Lệnh (Command History), Chương trình Soạn thảo và Gỡ rối (Editor and Debugger), và một Trình duyệt (Browser) để xem trợ giúp, Không gian Làm việc (Workspace), các Tệp, và Đường dẫn Tìm kiếm (Search Path). • Thư viện Hàm Toán học (Mathematical Function Library). Đây là một tập hợp các thuật toán tính toán trải rộng từ các hàm cơ cở như cộng, trừ, sin, cos, các phép tính số học phức, tới các hàm phức tạp hơn như đảo ma trận, tính giá trị riêng (eigenvalue) của ma trận, các hàm Bessel, và các phép biến đổi nhanh.

• Ngôn ngữ MATLAB (MATLAB Language). Đây là ngôn ngữ ma trận/mảng bậc cao với các khai báo luồng điều khiển, các hàm số, các cấu trúc dữ liệu, vào/ra, các đặc điểm lập trình hướng đối tượng. Nó cho phép viết cả các chương trình gọn nhẹ hay các chương trình ứng dụng phức tạp.) MATLAB có nhiều phương tiện hiển thị vector và ma trận ở dạng đồ thị, cũng như sửa đổi và in các đồ thị này. Nó bao gồm các hàm bậc cao để trực giác hóa các dữ liệu hai và ba chiều, xử lý ảnh, hoạt hình, và biểu diễn đồ họa.

Nó cũng bao gồm cả các hàm bậc thấp cho phép tùy biến hóa đồ họa cũng như xây dựng các giao diện đồ họa hoàn chỉnh cho các ứng dụng MATLAB của người sử dụng. • Giao diện Chương trình Ứng dụng MATLAB (MATLAB Ap- plication Program Interface [API]). Đây là một thư viện cho phép viết các chương trình C và Fortran tương tác với MATLAB. Nó có cả các phương tiện để gọi các thường trình (routine) từ MATLAB, dùng MATLAB như là động cơ tính toán, và dùng để đọc và viết MAT-files.

MATLAB cung cấp một họ các giải pháp theo từng ứng dụng, được gọi la hộp công cụ (toolbox). Hộp công cụ MATLAB bao gồm một tập hợp đầy đủ các hàm MATLAB ở dạng tệp "m" (m-file) dùng để mở rộng môi trường MATLAB cho việc giải quyết các loại vấn đề cụ thể. Các ví dụ về phạm vi ứng dụng của các hộp công cụ MATLAB là xử lý tín hiệu, hệ thống điều khiển, mạng nơ-ron, fuzzy logic, wavelet, mô phỏng, và nhiều ứng dụng khác.2 Khởi động và Thoát khỏi MATLAB Để khởi động MATLAB từ Windows, nhắp đúp (double-click) vào biểu tượng MATLAB trên màn hình desktop của windows. Sau khi khởi động xong màn hình sẽ hiện ra cửa sổ Môi trường Làm việc của MATLAB gồm 3 phần chính là: Thư mục Hiện thời (Current Directory), cửa sổ Lịch sử Câu lệnh (Command History) và Cửa sổ Câu lệnh (Command Windows) như ở Hình 1.

Làm việc với MATLAB Desktop 3 Hình 1.1: Môi trường làm việc của MATLAB Để kết thúc MATLAB có thể thực hiện bằng cách nhắp phím chuột trái vào ô đóng cửa sổ hình dấu sao ( ) phía trên và bên tay trái cửa sổ MAT- LAB. Ngoài ra cũng có thể kết thúc MATLAB bằng cách nhập vào câu lệnh » quit ở cửa sổ Command Windows rồi bấm Enter.3 Làm việc với MATLAB Desktop MATLAB Desktop bao gồm một Thanh Công cụ (Tool Bar) với các menu File, Edit, Debug, Desktop, Windows và Help. Bên cạnh Thanh Công cụ là một menu kéo xuống (Pull-down Menu) cho phép xem và thay đổi thư mục làm việc hiện thời. Nội dung của thư mục làm việc hiện thời được hiển thị ở cửa sổ Current Menu.

Phía dưới của sổ Current Menu là cửa sổ Command History hiển thị các câu lệnh MATLAB đã được nhập trước đó. Tiếp theo cửa sổ Command History xuống phía dưới có phím , cho phép truy nhập nhanh đến các thư viện của MATLAB, Simulink và các cài đặt Desktop Tools hay là các lựa chọn Preferences. Cửa sổ to nhất trong MATLAB Desktop là cửa sổ câu lệnh Command Window dùng để nhập các câu lệnh MATLAB hay chạy các chương trình cho trước. Giới thiệu Matlab 1.4 Các lệnh MATLAB cơ bản Các câu lệnh của MATLAB gần giống với các câu lệnh Unix.

Một số câu lệnh cơ bản của MATLAB được tóm tắt lại dưới đây: » ls Liệt kê nội dung của thư mục làm việc hiện thời. Cũng có thể dùng câu lệnh dir thay cho câu lệnh ls. temp liệt kê thư mục con temp bên trong thư mục làm việc hiện thời của MATLAB. » pwd Hiển thị đường dẫn của thư mục hiện tại.

Ví dụ: » pwd ans = C:\MATLAB701\work chỉ ra đường dẫn của thư mục làm việc hiện tại là C:\MATLAB701\work » who chỉ ra các biến đang được lưu ở bộ nhớ. Ví dụ: » a=1 a = 1 » b=2 b = 2 » who Your variables are: a b liệt kê hai biến a và b đang được lưu trữ ở bộ nhớ chương trình. Các ký hiệu đặc biệt 5 » clear [tên biến] xóa biến có tên được khai báo khỏi bộ nhớ. Ví dụ: » clear a » who Your variables are: b Để xóa hết tất cả các biến đang được lưu tại bộ nhớ, sử dụng lệnh »clear all »clc lệnh xóa toàn bộ thông tin trên Command Windows và đưa con trỏ trở về vị trí ban đầu.5 Các ký hiệu đặc biệt ( ) dấu ngoặc tròn được sử dụng để chỉ ra thứ tự ưu tiên trong các biểu thức số học hoặc bao quanh đối số của một hàm số.

Dấu ngoặc đơn cũng được dùng để bao quanh chỉ số phần tử trong một vector hay ma trận. Ngoài ra, dấu ngoặc đơn này còn được sử dụng để bao quanh các chỉ số dưới (subscript) logic. Ví dụ: A(2) chỉ ra phần tử thứ 2 của A. A([1 2 3]) liệt kê các phần tử thứ nhất, hai và ba của A.5) liệt kê các phần tử của A lớn hơn 0.

Giới thiệu Matlab A = 2 6 3 1 2 3 định nghĩa một ma trận với sáu phần tử cho trước. { } dấu ngoặc móc được sử dụng để tạo ra các mảng tế bào (cell array). Bộ dấu ngoặc móc này tương tự như bộ ngoặc vuông ngoại trừ các cấp độ nesting được bảo toàn. ′ biểu diễn phép toán chuyển vị liên hợp phức của một ma trận.

Ví dụ, A′ là ma trận chuyển vị liên hợp phức của A còn A.′ là ma trận chuyển vị của A. dấu chấm biểu diễn phân cách giữa phần nguyên và phần thập phân của một số thập phân. ; dấu chấm phảy dùng để ngăn cách các hàng khi khai báo ma trận, hoặc ngăn không hiển thị kết quả một phép toán trên màn hình. % dấu phần trăm dùng để tạo chú thích.

Tất cả các câu lệnh viết sau dấu phần trăm này đều bị bỏ qua. dấu 3 chấm dùng để nối hai phần của một câu lênh trên 2 dòng với nhau. Một câu lệnh dài có thể viết trên 2 dòng cho tiên theo dõi. Khi đó, dấu 3 chấm được sử dụng để nối 2 dòng với nhau.

Tài liệu tham khảo [1] Getting started with Matlab. The Mathworks Inc. Breiner, Matlab for Engineers. Chương 2 Tính toán và Lập trình sử dụng Matlab 2.1 Các phép tính số học Bốn phép tính số học cơ bản gồm cộng, trừ, nhân, chia được thể hiện tương ứng bằng các ký hiệu +, −, ∗, /.

Ví dụ » 2 + 1 ans= 3 » 3 − 1 ans= 2 » 2 ∗ 3 ans= 6 » 6/3 ans= 2 Với các phép tính phức tạp hơn có dấu ngoặc thì dấu ngoặc đơn (gồm cả mở và đóng) được sử dụng để phân cách thứ tự ưu tiên. Ví dụ, phép tính [(2 + 3) − (15 − 3)][7 + 5 − 4] (2.1) 2 được biểu diễn trong Matlab như sau 7 8 Chương 2. Tính toán và Lập trình sử dụng Matlab » ((2+3)-(15-3))*(7+5-4)/2 trong đó dấu ngoặc đơn được sử dụng thay cho dấu ngoặc vuông đã được mặc định dùng cho vector và ma trận trong Matlab. Trong Matlab phép tính lấy mũ được biểu diễn bởi ký hiệu bởi dấu mũ ˆ như: » 5ˆ 2 ans= 25 2.2 Các toán tử so sánh Trong Matlab các toán tử so sánh được biểu diễn như sau: nhỏ hơn (<), lớn hơn (>), nhỏ hơn hoặc bằng (<=), lớn hơn hoặc bằng (>=), bằng (trùng) nhau (==), khác nhau (∼=).

Khi hai mảng có cùng kích thước được so sánh với nhau thì toán tử so sánh sẽ thực hiện việc so sánh từng phần tử với nhau. Các toán tử <, >, <=, and >= chỉ so sánh phần thực của các toán hạng với nhau. Các toán tử == and ∼= thực hiện so sánh cả phần thực và phần ảo của hai toán hạng. Kết quả của phép toán so sánh cho ta 1 nếu phép so sánh là TRUE và ngược lại 0 nếu FALSE.

Một số ví dụ về toán tử so sánh được trình bày ở dưới đây » 1==2 ans = 0 » 3 > 1 ans = 1 » 4 <= 10 ans = 1 » 3 ∼= 7 ans = 1 2. Các toán tử so sánh 9 Trong trường hợp so sánh 2 vector hay 2 ma trận với nhau thì toán sử so sánh thực hiện và cho ta kết quả so sánh của từng phần tử tương ứng ở 2 vector hay ma trận với nhau. Ví dụ » A=[1 3 4; 2 8 7; 6 9 5] A = 1 3 4 2 8 7 6 9 5 » B=[3 1 4; 7 8 2; 6 5 9] B = 3 1 4 7 8 2 6 5 9 » A==B ans = 0 0 1 0 1 0 1 0 0 » A∼=B ans = 1 1 0 1 0 1 0 1 1 » A >= B ans = 0 1 1 0 1 1 10 Chương 2. Tính toán và Lập trình sử dụng Matlab 1 1 0 » A < B ans = 1 0 0 1 0 0 0 0 1 2.3 Các toán tử logic Các ký hiệu &, |, và ∼ được sử dụng để biểu diễn toán tử logic AND, OR, và NOT.

Các toán tử này làm việc với từng phần tử của mảng, với 0 biểu diễn FALSE còn 1 hay bất kỳ phần tử khác 0 nào biểu diễn TRUE. Các toán tử logic trả lại một mảng logic với các phần tử 0 (FALSE) hoặc 1 (TRUE). Các ký hiệu trên cũng có thể thay thế bằng cách sử dụng các hàm Matlab ở dạng and(A, B), or(A, B), hay not(A). Hàm HOẶC tuyện đối được biểu diễn như sau xor(A, B).

Trình tự ưu tiên của các toán tử logic là NOT, OR và AND. Phép toán sử dụng toán tử AND (&) cho kết quả TRUE nếu cả hai toán hạng đều TRUE về mặt logic. Nói theo thuật ngữ số, thì phép toán AND cho ta kết quả TRUE nếu cả hai toán hạng đều khác 0. Ví dụ: » a=[3 5 0 4 0] a = 3 5 0 4 0 » b=[3 1 0 0 2] b = 3 1 0 0 2 » a & b ans = 1 1 0 0 0 Các số 1 chỉ ra các phần tử tương ứng khác không của cả a và b .