(Đồ án hcmute) áp dụng phương pháp ước lượng kênh dùng tín hiệu pilot trên hệ thống ofdm trong môi trường fading

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT TP.HCM KHOA ĐÀO TẠO CHẤT LƯỢNG CAO ĐỒ ÁN TỐT NGHIỆP NGÀNH CÔNG NGHỆ KỸ THUẬT ĐIỆN TỬ, TRUYỀN THÔNG ÁP DỤNG PHƯƠNG PHÁP UỚC LUỢNG KÊNH DÙNG TÍN HIỆU PILOT TRÊN HỆ THỐNG OFDM TRONG MƠI TRUỜNG FADING GVHD: PGS.TS PHẠM HỒNG LIÊN SVTH: HUỲNH HOÀNG NAM MSSV: 12141342 SKL 0 Tp Hồ Chí Minh, tháng 08 năm 2016 an BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH KHOA ĐÀO TẠO CHẤT LƯỢNG CAO ĐỒ ÁN TỐT NGHIỆP ÁP DỤNG PHƯƠNG PHÁP ƯỚC LƯỢNG KÊNH DÙNG TÍN HIỆU PILOT TRÊN HỆ THỐNG OFDM TRONG MƠI TRƯỜNG FADING SVTH: HUỲNH HỒNG NAM MSSV: 12141342 Khoá: 2012 Ngành: CNKT Điện Tử, Truyền Thơng GVHD: PGS.TS PHẠM HỒNG LIÊN Tp Hồ Chí Minh, tháng năm 2016 an CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc *** -Tp Hồ Chí Minh, ngày……tháng… năm…… NHIỆM VỤ ĐỒ ÁN TỐT NGHIỆP Họ tên sinh viên: Huỳnh Hoàng Nam MSSV:12141342 Ngành: CNKT Điện Tử, Truyền Thông Lớp: 12141CLVT Giảng viên hướng dẫn: PGS.TS Phạm Hồng Liên ĐT: 093 979 21 24 Ngày nhận đề tài: 03/2016 Ngày nộp đề tài: 23/07/2016 Tên đề tài: Áp dụng phương pháp ước lượng kênh dùng tín hiệu pilot hệ thống OFDM môi trường Fading Nội dung thực đề tài: Mô hệ thống OFDM sử dụng phương pháp ước lượng kênh dùng tín hiệu pilot mơi trường trung Fading Sản phẩm: Mô hệ thống OFDM sử dụng phương pháp ước lượng kênh dựa vào tín hiệu pilot môi trường kênh truyền bị fading phần mềm Matlab TRƯỞNG NGÀNH GIẢNG VIÊN HƯỚNG DẪN i an CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh Phúc ******* PHIẾU NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN Họ tên sinh viên: Huỳnh Hồng Nam MSSV: 12141342 Ngành: CNKT Điện Tử, Truyền Thơng Tên đề tài: Áp dụng phương pháp ước lượng kênh dùng tín hiệu pilot hệ thống OFDM mơi trường Fading Họ tên giáo viên hướng dẫn: PGS.TS Phạm Hồng Liên NHẬN XÉT Về nội dung đề tài & khối lượng thực hiện: Ưu điểm: Khuyết điểm: Đề nghị cho bảo vệ hay không? Đánh giá loại: Điểm:……………….(Bằng chữ: ) Tp.Hồ Chí Minh, ngày……tháng… năm… Giáo viên hướng dẫn (Ký & ghi rõ họ tên) ii an CỘNG HOÀ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh Phúc ******* PHIẾU NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN Họ tên Sinh viên: Huỳnh Hoàng Nam MSSV: 12141342 Ngành: CNKT Điện Tử, Truyền Thông Tên đề tài: Áp dụng phương pháp ước lượng kênh dùng tín hiệu pilot hệ thống OFDM môi trường Fading Họ tên giáo viên phản biện: NHẬN XÉT Về nội dung đề tài & khối lượng thực hiện: Ưu điểm: Khuyết điểm: Đề nghị cho bảo vệ hay không? Đánh giá loại: Điểm:……………….(Bằng chữ: ) Tp Hồ Chí Minh, ngày…….tháng…….năm…… Giáo viên phản biện (Ký & ghi rõ họ tên) iii an LỜI CẢM ƠN Để hồn thành chương trình học ngành Điện Tử Viễn Thông thuộc khoa Đào Tạo Chất Lượng Cao trường Đại Học Sư Phạm Kỹ Thuật Thành Phố Hồ Chí Minh, em thực đề tài “ÁP DỤNG PHƯƠNG PHÁP ƯỚC LƯỢNG KÊNH DỰA VÀO TÍN HIỆU PILOT TRÊN HỆ THỒNG OFDM TRONG MƠI TRƯỜNG FADING” Với lịng kính trọng biết ơn sâu sắc em xin gửi lời cảm ơn chân thành đến: Em xin chân thành gửi lời cảm ơn đến Cô PGS.TS PHẠM HỒNG LIÊN, người trực tiếp hướng dẫn em thực đề tài Trong trình thực đề tài, Cơ tận tình bảo, theo dõi hỗ trợ em để em hoàn thành tốt đề tài Em xin gửi lời cảm ơn đến Thầy/Cô khoa Đào Tạo Chất Lượng Cao trường hỗ trợ, tận tình dạy em suốt thời gian em theo học trường, cung cấp cho em kiến thức để em thực đề tài Bên cạnh đó, em xin gửi lời cảm ơn đến gia đình, nhà trường bạn bè tạo điều kiện, hỗ trợ, động viên em suốt trình học làm đề tài tốt nghiệp Tuy nhiên, kiến thức em hạn chế nên chắn thiếu sót đề tài Rất mong nhận góp ý q Thầy/Cơ, bạn bè để em hoàn thiện Xin trân trọng cảm ơn Người thực đề tài Huỳnh Hồng Nam iv an TĨM TẮT Ngày yêu cầu truyền thông không dây ngày phát triển mạnh, địi hỏi kỹ tḥt điều chế khắt khe để cung cấp tốc độ truyền liệu cao Đặc biệt ứng dụng địi hỏi tốc độ cao truyền thơng đa phương tiện không dây, truy cập Internet không dây hay hệ thiết bị tương lai… Vì kỹ thuật OFDM đời kỹ thuật hứa hẹn giải yêu cầu Lợi ích OFDM thực thi đơn giản hạn chế ảnh hưởng kênh truyền chọn lọc tần số Ngồi vấn đề phải đối mặt hệ thống truyền thông không dây ảnh hưởng kênh truyền vô tuyến Để giải vấn đề thông thường người ta sử dụng ước lượng kênh truyền Mà phương pháp truyền thống cho ước lượng kênh truyền sử dụng chuỗi liệu biết trước (tín hiệu pilot) để xác định hàm truyền kênh truyền Kỹ thuật hiệu để loại bỏ ảnh hưởng kênh truyền Đồ án thực ứng dụng phương pháp ước lượng kênh dựa vào tín hiệu pilot với giải thuật LS hệ thống OFDM môi trường kênh truyền fading v an ABSTRACT Recently, increasing interest has been concentrated on modulation techniques that provide high data rates over broadband wireless channels for applications, including wireless multimedia, wireless Internet access, and future-generation mobile communication systems Orthogonal frequency division multiplexing (OFDM) is a promising digital modulation scheme to simplify the equalization in frequency selective channels The main benefit is that it simplifies implementation, and it is against the frequency-selective fading channels Thus, OFDM systems can provided high-performance transmission Moreover, one of the main problems faced in wireless communications is the influence of radio channels To resolve common problems people will use the estimated channels That traditional methods for estimating channels using sequence data are known in advance (pilot signal) to determine the channel transfer function That technique is very effective to eliminate the influence of channels This thesis is focused on working out a equalizer based on LS algoritm for OFDM system in fading channel Key words: OFDM, channel estimation with pilot signal, LS algorithm,Fading vi an MỤC LỤC NHIỆM VỤ ĐỒ ÁN TỐT NGHIỆP i PHIẾU NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN ii PHIẾU NHẬN XÉT CỦA GIÁO VIÊN PHẢN BIỆN iii TÓM TẮT v ABSTRACT vi MỤC LỤC vii DANH MỤC CÁC CHỮ VIẾT TẮT .x DANH MỤC CÁC BẢNG BIỂU xi DANH MỤC CÁC BIỂU ĐỒ VÀ HÌNH ẢNH xi CHƯƠNG 1: GIỚI THIỆU .1 1.1 Vai trò đề tài 1.2 Mục đích đề tài .1 1.3 Phạm vi nghiên cứu giới hạn đề tài .1 1.4 Bố cục đề tài CHƯƠNG 2: LÝ THUYẾT LIÊN QUAN .3 2.1 Kênh truyền vô tuyến 2.1.1 2.1.1.1 Suy hao theo khoảng cách truyền 2.1.1.2 Ảnh hưởng phản xạ 2.1.1.3 Suy hao vật cản 2.1.2 2.2 Mô hình suy giảm diện rộng Mơ hình fading diện hẹp hiệu ứng đa đường 2.1.2.1 Hiệu ứng đa đường 2.1.2.2 Hiệu ứng Doppler 2.1.2.3 Các thông số kênh truyền đa đường 2.1.2.4 Phân loại kênh truyền Small scale fading .10 Kỹ thuật điều chế OFDM 14 2.2.1 Nguyên lí OFDM 14 2.2.2 Sơ đồ khối hệ thống OFDM .17 2.2.2.1 Ánh xạ điều chế .18 2.2.2.2 Bộ chuyển đổi nối tiếp – song song 22 vii an 2.2.2.3 Chuyển đổi miền tần số sang miền thời gian 23 2.2.2.4 Chèn khoảng bảo vệ 24 2.2.2.5 Điều chế RF 25 2.2.2.6 Máy thu OFDM .25 2.2.3 Ưu nhược điểm kỹ thuật điều chế OFDM 28 2.2.3.1 Ưu điểm kỹ thuật điều chế OFDM 28 2.2.3.2 Nhược điểm kỹ thuật điều chế OFDM .29 CHƯƠNG 3: PHƯƠNG PHÁP ƯỚC LƯỢNG KÊNH TRUYỀN DÙNG TÍN HIỆU PILOT .30 3.1 Cấu trúc tín hiệu pilot 30 3.1.1 Pilot dạng khối 30 3.1.2 Pilot dạng lược 31 3.1.3 Pilot dạng Lattice .31 3.2 Ước lượng kênh truyền sử dụng tín hiệu pilot 32 3.2.1 Ước lượng kênh truyền sử dụng pilot dạng khối 32 3.2.1.1 Ước lượng LS 33 3.2.1.2 Ước lượng MMSE 33 3.2.1.3 Ước lượng MMSE cải tiến 34 3.2.2 3.2.2.1 Ước lượng kênh truyền sử dụng pilot dạng lược 36 Ước lượng LS nội suy chiều 37 Nội suy tuyến tính (Linear Interpolation) 37 Nội suy bậc (Second Order) .37 Nội suy lowpass 38 Nội suy spline cubic 38 3.2.2.2 Ước lượng ML 38 3.2.2.3 Ước lượng PCMB 39 CHƯƠNG 4: KẾT QUẢ MÔ PHỎNG 40 4.1 Mô 40 4.1.1 pilot Sơ đồ khối hệ thống OFDM áp dụng phương pháp ước lượng kênh dùng tín hiệu 40 4.1.2 Mơ hình kênh truyền ITU sử dụng cho WiMax di động 41 4.2 Kết mô 43 viii an 190 -19.2 800 -4.9 410 -22.8 1200 -8 - - 2300 -7.8 - - 3700 -23.9 Bảng 4.3 Mơ hình kênh truyền vehicular Kênh truyền A Kênh truyền B Tap Trễ (s) Công suất (dB) Trễ (s) Công suất (dB) 0 0 0.4 -1 0.8 -1 0.8 -9 1.6 -9 1.2 -10 2.2 -10 1.8 -15 3.6 -15 2.6 -20 5.2 -20 Ngoài chịu tác động đa đường mơ kênh truyền cịn chịu tác động nhiễu Gauss Do việc ước lượng xác hàm truyền kênh truyền vơ tuyến phức tạp kênh truyền chịu tác động nhiều tác nhân thường xuyên thay đổi nên việc nội suy kênh truyền khó ta khơng biết quy luật kênh truyền 4.2 Kết mô  So sánh BER ước lượng LS sử dụng điều chế 16QAM 16PSK: 43 an Bảng 4.4 Kết BER ước lượng LS sử dụng điều chế 16QAM 16PSK SNR(dB) 10 12 14 16 18 20 0.3611 0.3207 0.2760 0.2244 0.1693 0.1150 0.0692 0.0329 0.0128 0.0043 0.0012 0.4069 0.3729 0.3357 0.2958 0.2550 0.2082 0.1564 0.1044 0.0591 0.0248 0.0074 Map 16QAM 16PSK SO SANH BER CUA DIEU CHE QAM VA PSK 10 16QAM 16PSK -1 BER 10 -2 10 -3 10 -4 10 10 12 SNR(dB) 14 16 18 20 Hình 4.2 So sánh BER 16QAM 16PSK Hình 4.2 kết mô so sánh tỉ lệ lỗi bit hai phương pháp ánh xạ điều chế khác điều chế biên độ pha QAM điều chế pha PSK, mức điều chế 16 Mô với vận tốc thiết bị v = 5km/h điều kiện LOS Kết mô cho thấy tỉ lệ lỗi bit phương pháp điều chế QAM thấp so với điều chế PSK, điều với lý thuyết Vì điều chế QAM điều chế pha biên độ nên khác biệt nhiều so với điều chế PSK điều chế pha tín hiệu 44 an  So sánh BER ước lượng LS sử dụng loại điều chế khác nhau: Bảng 4.5 Kết BER ước lượng LS sử dụng loại điều chế khác SNR(dB) 10 12 14 16 18 20 0.2027 0.1360 0.0761 0.0341 0.0341 0.0102 0.0017 0.0002 0 0.3340 0.2838 0.2287 0.1638 0.1075 0.0595 0.0279 0.0096 0.0027 0.0005 0.0001 0.3589 0.3184 0.2728 0.2271 0.1710 0.1154 0.0687 0.0340 0.0139 0.0043 0.0012 0.2728 0.1917 0.1135 0.0513 0.0149 0.0028 0.0002 0.0009 0 0.3712 0.3239 0.2647 0.2029 0.1358 0.0764 0.0332 0.0096 0.0016 0.0001 0.4060 0.3735 0.3371 0.2966 0.2547 0.2063 0.1565 0.1051 0.0581 0.0254 0.0071 Map 4QAM 8QAM 16QAM 4PSK 8PSK 16PSK SO SANH BER CUA CAC KIEU VA CAC MUC DIEU CHE 10 -1 10 -2 BER 10 -3 10 4QMA 4PSK 8QAM 8PSK 16QAM 16PSK -4 10 -5 10 -6 10 10 12 SNR(dB) 14 16 18 20 Hình 4.3 So sánh BER kiểu mức điều chế khác Hình 4.3 kết mơ so sánh tỉ lệ lỗi bit hai phương pháp ánh xạ điều chế khác điều chế biên độ pha QAM điều chế pha PSK, với ba mức điều chế khác 4, 16 Mô với vận tốc thiết bị v = 5km/h điều kiện LOS Kết mô cho thấy tỉ lệ lỗi bit phương pháp điều chế QAM 45 an mức điều chế thấp so với điều chế PSK với mức tương ứng Ở phương pháp điều chế tỉ lệ lỗi bit tỷ lệ thuận với mức điều chế - mức điều chế tăng tỉ lệ lỗi bit cao  So sánh BER ước lượng LS sử dụng hàm nội suy khác nhau: SO SANH BER CUA CAC PHUONG PHAP NOI SUY 10 -1 BER 10 LinearQAM LinearPSK Second OrderQAM Second OrderPSK LowpassQAM LowpassPSK SplineQAM SplinePSK -2 10 -3 10 -4 10 10 12 SNR(dB) 14 16 18 20 Hình 4.4 So sánh BER với hàm nội suy khác Hình 4.4 kết mơ so sánh tỉ lệ lỗi bit hai phương pháp ánh xạ điều chế khác điều chế biên độ pha QAM điều chế pha PSK, với mức điều chế 16 sử dụng hàm nội suy khác việc ước lượng kênh truyền Mô với vận tốc thiết bị v = 5km/h điều kiện LOS Kết mô cho thấy tỉ lệ lỗi bit sử dụng hàm nội suy tuyến tính thấp nhất, hàm nội suy bậc cao hàm nội suy spline cubic Tỷ lệ lỗi bit chịu ảnh hưởng nhiều phương pháp điều chế, ta thấy dùng hàm nội suy tuyến tính với phương pháp điều chế PSK tỉ lệ lỗi bit cao so với dùng hàm nội suy spline cubic sử dụng điều chế QAM 46 an  So sánh BER ước lượng LS môi trường Single Path (LoS) với vận tốc khác nhau: Bảng 4.6 Kết BER ước lượng LS môi trường Single Path (LoS) với vận tốc khác SNR(dB) V(km/h) 20 60 100 10 15 20 0.2034 0.2046 0.2033 0.2038 0.2192 0.0520 0.0530 0.0569 0.0774 0.0819 0.0015 0.0016 0.0017 0.0033 0.0249 0 0.00001 0.0001 0 0 SO sanh BER moi truong Single Path voi cac van toc khac 10 1Km/h 5Km/h 20Km/h 60Km/h 100Km/h -1 10 -2 BER 10 -3 10 -4 10 -5 10 -6 10 10 12 SNR(dB) 14 16 18 20 Hình 4.5 So sánh BER với vận tốc khác điều kiện LOS Hình 4.5 kết mơ so sánh tỉ lệ lỗi bit với trường hợp vận tốc thiết bị di chuyển khác Mô sử dụng điều chế 4QAM với hàm nội suy tuyến tính điều kiện LOS Kết mô cho thấy tỉ lệ lỗi bit tăng 47 an vận tốc di chuyển thiết bị tăng Điều giải thích hiệu ứng Doppler  So sánh BER ước lượng LS môi trường khác nhau: Bảng 4.7 Kết BER ước lượng LS môi trường khác SNR(dB) Môi trường Single Path Indoor Pedestrian Vehicular 10 15 20 0.2022 0.2619 0.3473 0.4401 0.0524 0.1338 0.2704 0.3711 0.0016 0.0573 0.2316 0.3849 0.0440 0.2368 0.3880 0.0329 0.2033 0.3779 SO SANH BER TRONG CAC MOI TRUONG KHAC NHAU 10 -1 BER 10 -2 10 SinglePath4QAM Indoor4QAM Pedestrian4QAM Vehicular4QAM -3 10 -4 10 10 12 SNR(dB) 14 16 18 20 Hình 4.6 So sánh BER với mơi trường khác Hình 4.6 kết mô so sánh tỉ lệ lỗi bit môi trường kênh truyền khác sử phương pháp ánh xạ điều chế điều chế biên độ pha QAM, với mức điều chế Mô với vận tốc thiết bị v = 1km/h sử dụng hàm nội suy tuyến tính Kết mơ cho thấy tỉ lệ lỗi bit truyền với đường truyền LOS cho 48 an BER thấp nhất, BER cao truyền môi trường Indoor cao truyền môi trường Vehicular Điều ảnh hưởng mặt độ trễ độ lợi công suất suất môi trường khác nhau, mơi trường Indoor bị ảnh hưởng thông số so với môi trường Pedestrian nhiều môi trường Vehicular 49 an CHƯƠNG 5: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ĐỀ TÀI Chương kết luận luận văn đề hướng phát triển đề tài 5.1 Kết luận Với luận văn “Áp dụng phương pháp ước lượng kênh dùng tín hiệu pilot hệ thống OFDM môi trường Fading” em thực số vấn đề như: Tìm hiểu loại mơ hình kênh truyền vơ tuyến ảnh hưởng đến kênh truyền Tìm hiểu kỹ thuật điều chế OFDM Tìm hiểu phương pháp ước lượng kênh dùng pilot Áp dụng phương pháp ước lượng tìm hiểu vào mô Matlab với điều kiện môi trường khác từ thấy ảnh hưởng môi trường kênh truyền vô tuyến Tuy thực vấn đề hạn chế kiến thức thời gian thực nên chắn đề tài cịn nhiều thiếu sót chưa đánh giá cách cụ thể chất lượng phương pháp,… 5.2 Hướng phát triển đề tài Đề tài thực số phương pháp ước lượng kênh truyền đơn giản nên chưa cho thấy nhìn tổng quát phương pháp ước lượng kênh Tuy nhiên, với tìm hiểu đề tài tạo tảng giúp em thực nghiên cứu sâu giúp hiểu ước lượng kênh truyền Em thực nghiên cứu theo hướng như: - Ước lượng kênh với giải thuật tốt Tiến hành mã hoá kênh để làm giảm tỷ lệ lỗi bit Áp dụng mô phương pháp ước lượng kênh truyền dùng tín hiệu pilot khác MMSE, ML,… Áp dụng phương pháp ước lượng kênh truyền vào hệ thống MIMO OFDM 50 an TÀI LIỆU THAM KHẢO Tiếng Việt [1] Đặng Ngọc Khoa , Phạm Hồng Liên, Trần Thanh Phương (2006), MATLAB ứng dụng viễn thông, NXB Đại học Quốc Gia TPHCM Tiếng Anh [2] Coleri, S., Ergen, M., Puri, A and Bahai (2002), Channel estimation techniques based on pilot arrangement in OFDM systems, IEEE Transactions on Broadcasting.Volume 48, No 03, [Pg 769-770] [3] Mohammad Ali, Rehan Elahi, Ali Saljuk, Sheheryar Bukhari (2008), Channel Estimation of MIMO OFDM Systems [4] Yong Soo Cho, Jaekwon Kim, Won Young Yang, Chung-Gu Kang (2010), MIMO – OFDM Wireless Communications with MALAB.[Pg 187-190] [5] Raj Jain (2007), Channel Models A Tutorial [6] ITU-R Recommendation M.1225 (1997), Guidelines for evaluation of radio transmission technologies for IMT-2000.[Pg 28] [7] Yushi Shen, Ed Martinez (2006), Channel Estimation in OFDM Systems,[Pg.7-9] [8] Ramjee Prasad (2004), OFDM for wireless communications system, Chapter 1, Universal personal communications [9] R van Nee and R Prasad (2000), OFDM for Wireless Multimedia Communications, Artech House 51 an PHỤ LỤC  Hàm transmit function [base_ofdm,data_Map,numMap,data_pilot_comb] = transmitLS(bits,Tye_map,M_map,numcarr,Pilot_Distance,NFFT,G) [data_Map,numMap] = Map_Data(bits,Tye_map,M_map); len=length(data_Map); Np = ceil(numcarr/Pilot_Distance); numcarr_total = numcarr; numcarr = numcarr - Np; r = rem(len,numcarr); if r ~= for i = 1:numcarr-r data_Map(len+i) = 0; end end len=length(data_Map); num_OFDM=len/numcarr; data_pattern=[]; for i=0:num_OFDM-1 Map_tem=[]; for n=1:numcarr Map_tem=[Map_tem,data_Map(i*numcarr+n)]; end data_pattern=[data_pattern;Map_tem]; clear QAM_tem; end  Hàm receiver function [bits_output] = receiveLS(data_channel,numcarr,NFFT,G,numQAM,data_pilot_comb,Pilot_Distance ,algorithm,Tye_map,M_map) num_OFDM = length(data_channel)/(NFFT+G); base_ofdm = reshape(data_channel,NFFT+G,num_OFDM); base_ofdm = conj(base_ofdm'); receiver_data = []; for i=1:num_OFDM Demodulated_signal = OFDM_Demodulator(base_ofdm(i,:),NFFT,G,numcarr); [data_MAP,pilot_receiver] = Extractor_pilot_comb(Demodulated_signal,numcarr,Pilot_Distance); 52 an H_pilot = pilot_receiver./data_pilot_comb if(algorithm==1) %Ham noi suy tuyen tinh H_estimated = interpolation_linear(H_pilot,numcarr,Pilot_Distance); elseif(algorithm==2) % Ham noi suy bac H_estimated = interpolation_secondorder(H_pilot,numcarr,Pilot_Distance); elseif(algorithm==3) % Ham noi suy lowpass H_estimated = interpolation_lowpass(H_pilot,numcarr,Pilot_Distance); elseif(algorithm==4) % Ham noi suy spline cubic H_estimated = interpolation_spline(H_pilot,numcarr,Pilot_Distance); end receiver_OFDM = Demodulated_signal./H_estimated; [receiver_OFDM,pilot_receiver] = Extractor_pilot_comb(receiver_OFDM,numcarr,Pilot_Distance); receiver_data = [receiver_data;receiver_OFDM]; end y=conj(receiver_data'); y=reshape(y,size(y,1)*size(y,2),1); receiver_QAM=y(1:numQAM); bits_output = Demap_Data(receiver_QAM,Tye_map,M_map); clear y receiver_QAM;  Hàm Mod_OFDM function [y]=OFDM_Modulator(data,NFFT,G,numcarr); x=zeros(1,NFFT); for i=1:numcarr x(1,i+NFFT/2-numcarr/2)=data(1,i); end a=ifft(x); y=[a(NFFT-G+1:NFFT),a]; clear data x a  Hàm Demod_OFDM function [y]=OFDM_Demodulator(data,NFFT,G,numcarr); x_remove_guard_interval=[data(G+1:NFFT+G)]; x=fft(x_remove_guard_interval); for i=1:numcarr y(i)=x(1,i+NFFT/2-numcarr/2); end clear data x_remove_guard_interval x;  Hàm Insert Pilot 53 an function [IPS,data_pilot_comb] = Insert_Pilot_comb(QAM_symbol,numcarr,Np,num_ofdm,Tye_map,M_map,Pilot_Distan ce) L = log2(M_map); source_pilot = randi([0 1],1,Np*L); [Map_pilot,~] = Map_Data(source_pilot,Tye_map,M_map); data_pilot_comb = Map_pilot.'; IPS = []; for i = 1:num_ofdm TS2_BeforeIFFT = []; m=1; n=1; for k = 0:numcarr-1 if (mod(k, Pilot_Distance)==0) TS2_BeforeIFFT = [TS2_BeforeIFFT, data_pilot_comb(m)]; m=m+1; else TS2_BeforeIFFT = [TS2_BeforeIFFT,QAM_symbol(i,n)]; n=n+1; end; end; IPS = [IPS;TS2_BeforeIFFT]; end clear QAM_symbol source_pilot QAM_pilot num_ofdm TS2_BeforeIFFT;  Hàm Extract Pilot function [data_QAM,pilot_estimated] = Extractor_pilot_comb(data_afterFFT,numcarr,Pilot_Distance) pilot_estimated = []; data_QAM = []; for i = 0:numcarr-1 if (mod(i, Pilot_Distance)==0) pilot_estimated = [pilot_estimated,data_afterFFT(i+1)]; else data_QAM = [data_QAM,data_afterFFT(i+1)]; end end clear data_afterFFT;  Hàm channel function [outsignal] = channel(rs,snr,v,B,fc,delay_vector,gain_vector) if(v>=0) %kenh truyen co da duong c=3e8; v=v*10^3/3600; %doi km/h m/s f_dmax=v*fc/c; %tan so dich Doppler max chan = rayleighchan(1/B,f_dmax,delay_vector,gain_vector); outsignal = filter(chan,rs); else 54 an outsignal = rs; end outsignal = awgn(outsignal,snr,'measured','dB'); clear rs; clear delay_vector gain_vector; clear c v f_dmax;  Hàm Map function [Data_map,Num_map] = Map_Data(bits,tye_map,level) L = log2(level); len = length(bits); n = rem(len,L); if n ~= re = zeros(1,L-n); bits = [bits,re]; end len = length(bits); k = len/L; Num_map = k; data = reshape(bits,L,k); data = data'; symbols = bi2de(data,'left-msb'); if tye_map == Tx = modem.qammod(level); elseif tye_map == Tx = modem.pskmod(level); end Data_map = modulate(Tx,symbols); clear bits data symbols; end  Hàm Demap function [bits_out] = Demap_Data(code_Map,tye_map,level) if tye_map == Rx = modem.qamdemod(level); elseif tye_map == Rx = modem.pskdemod(level); end symbols = demodulate(Rx,code_Map); data = de2bi(symbols,log2(level),'left-msb'); data = data'; s = size(data); m = s(1)*s(2); data = reshape(data,1,m); 55 an len = length(data); n = rem(len,log2(level)); bits = zeros(1,len-n); for k=1:(len-n) bits(k)=data(k); end bits_out = bits; clear bits data symbols; end 56 an an