Nghiên cứu công nghệ ofdm và ứng dụng trong truyền hình số mặt đất dvb t2

LỜI CAM ĐOAN Tơi cam đoan cơng trình nghiên cứu riêng Các số liệu, kết nêu luận văn trung thực chưa cơng bố cơng trình khác Tác giả luận văn Vũ Văn Chiến LỜI CẢM ƠN Để hồn thành luận văn này, tơi xin gửi lời cảm ơn chân thành tới PGS TS Nguyễn Hữu Trung, người trực tiếp hướng dẫn, bảo, định hướng giúp đỡ tận tình suốt thời gian qua Tơi muốn nói lời cảm ơn tới thầy giáo Viện đại học mở Hà Nội tận tình bảo tạo điều kiện giúp đỡ, động viên tơi suốt q trình học tập Cuối xin gửi lời cảm ơn đến bố mẹ, anh chị em tất người thân bạn bè, người giành tốt đẹp suốt trình học tập để có đến ngày hơm Do có hạn chế thời gian nên luận văn không tránh khỏi thiếu sót Tơi mong góp ý thầy để kiến thức nghiên cứu công nghệ OFDM truyền hình số mặt đất DVB- T2 tơi phong phú MỤC LỤC MỞ ĐẦU CHƯƠNG TỔNG QUAN VỀ TRUYỀN HÌNH SỐ MẶT ĐẤT THEO TIÊU CHUẨN DVB-T 1.1 Giới thiệu truyền hình số mặt đất DVB – T 1.1.1 Sơ đồ khối 1.1.2 Đặc điểm 1.2 Điều chế truyền hình số 1.2.1 Tổng quan kỹ thuật điều chế COFDM 1.2.2 Đặc tính COFDM 10 1.2.3 Nguyên tắc COFDM 11 1.2.4 Điều chế số 13 1.2.5 Đồng kênh 18 1.2.6 Khoảng bảo vệ 19 1.2.7 Mạng đơn tần SFN 21 1.3 Mã hóa kênh điều chế DVB – T 22 1.3.1 Phân tán lượng 22 1.3.2 Mã tráo 23 1.3.3 Mã 24 1.3.4 Tráo 25 1.3.5 Hàm tín hiệu COFDM chuẩn DVB – T 26 1.3.6 Máy thu DVB – T thực tế 28 1.4 Các thông số đo kiểm tra 30 1.4.1 Độ xác tần số RF 30 1.4.2 Độ chọn lọc 31 1.4.3 Phạm vi điều khiển tự động tần số 31 1.4.4 Công suất RF/IF 31 1.4.5 Công suất tạp nhiễu 32 1.4.6 Độ nhạy máy thu/ dải động kênh gaussian 32 1.4.7 Hiệu suất công suất 33 1.4.8 Can nhiễu liên kết 33 1.4.9 Quan hệ ber tỉ số c/n thay đổi công suất máy phát 33 CHƯƠNG 2: TRUYỀN HÌNH SỐ MẶT ĐẤT THEO TIÊU CHUẨN DVB-T2 35 2.1 Giới thiệu chung tiêu chuẩn truyền hình số mặt đất DVB-T2 35 2.2 Cấu trúc DVB-T2 35 2.3 Lớp vật lý DVB - T2 37 2.4 Những giải pháp kỹ thuật 37 2.4.1 Ống lớp vật lý (Physical Layer Pipes - PLPs) 37 2.4.2 Băng tần phụ (1,7 Mhz 10 Mhz) 37 2.4.3 Mode sóng mang mở rộng (đối với 8K, 16K, 32K) 38 2.4.4 MISO dựa Alamouti (trên trục tần số) 38 2.4.5 Symbol khỏi đầu (P1 P2) 39 2.4.6 Mẫu hình tín hiệu Pilot (Pilot Pattern) 39 2.4.7 Phương thức điều chế 256-QAM 39 2.4.8 Chòm xoay (Rotated Constellation) 39 2.4.9 16K, 32K FFT tỷ lệ khoảng bảo vệ 1/128 40 2.4.10 Mã sửa sai LDPC/BCH 40 2.4.11 Kỹ thuật giảm thiểu tỷ số cơng suất đỉnh/cơng suất trung bình (Peak - to Average Power Ratio - PAPR) 40 2.4.12 Tráo bit, ánh xạ bit lên đồ thị chòm 40 2.4.13 Ánh xạ lên đồ thị chòm I/Q 41 2.4.14 Tráo tế bào, tráo thòi gian 41 2.4.15 Điều chế mã sửa sai liệu lớp 41 2.4.16 Cấu trúc khung tín hiệu DVB-T2 41 2.4.17 Những tiêu chuẩn DVB-T2 42 CHƯƠNG ĐẶC ĐIỂM VÀ TÍNH NĂNG KĨ THUẬT HỆ THỐNG OFDM 45 A Lịch sử phát triển công nghệ OFDM 45 B Đặc điểm kỹ thuật điều chế OFDM 45 3.1 Đa sóng mang (MULTICARRIER) 47 3.1.1 Sự trực giao (ORTHOGONAL) 48 3.1.2 Mơ tả tốn học OFDM 52 3.1.3 Trực giao miền tần số 55 3.2 Tạo thu OFDM 55 3.2.1 Nối tiếp - song song 57 3.2.2 Điều chế tải phụ 57 3.2.3 Điều chế RF 58 3.2.4 Thuật ngữ nhiễu pha vấn đề liên quan đế 58 3.3 Khoảng bảo vệ (GUARA PERIOD) 60 3.3.1 Bảo vệ chống lại OFFSET thời gian 61 3.3.2 Bảo vệ chống lại ISI 62 3.3.3 Giới hạn băng thông OFDM cửa số 63 3.3.4 Lọc băng thông 64 3.3.5 Độ phức tạp tính tốn lọc băng thơng FIR 66 3.3.6 Ảnh hưởng lọc băng thông tới tiêu kỹ thuật OFDM 67 3.4 Khoảng bảo vệ COSIN tăng RC (RAISED COSINE GUARD PERIOD) 69 3.5 Ảnh hưởng nhiễu GAUSS trắng cộng AWGN (ADDITIVE WHITE GAUSIAN NOISE) đến OFEM 69 3.5.1 Các sơ đồ điều chế 70 3.5.2 So sánh truyền OFDM với truyền sóng mang đơn 71 3.5.3 Các giới hạn điều chế hệ thống 72 3.5.4 Mã Gray 73 3.5.5 Điều chế kết hợp 75 3.5.6 Điều chế pha vi sai 76 3.5.7 Visai QAM 77 3.6 Ảnh hưởng méo tới OFDM 81 3.6.1 Mơ hình hố méo 82 3.6.2 Mở rộng phổ cắt méo 84 3.6.3 SNR hiệu dụng từ cắt méo 85 3.7 Ảnh hưởng lỗi đồng thời gian 86 3.8 Ảnh hưởng lỗi đồng tần số 87 CHƯƠNG MƠ PHỎNG Q TRÌNH ĐIỀU CHẾ VÀ ỨNG DỤNG OFDM TRONG TRUYỀN HÌNH SỐ DVB-T2 89 4.1 Giới thiệu chương 89 4.2 Mô hệ thống OFDM simulink 89 4.3 Một số lưu đồ thuật tốn chương trình 92 4.3.1 Lưu đồ mô kênh truyền 92 4.3.2 Lưu đồ mơ thu phát tín hiệu OFDM 93 4.3.3 Lưu đồ mơ thu phát tín hiệu QAM 94 4.3.4 Lưu đồ mơ thuật tốn tính BER 96 4.4 Kết chương trình mơ 97 4.4.1 So sánh tín hiệu QAM OFDM 97 4.4.2 So sánh tín hiệu âm điều chế QAM OFDM 98 4.5 Kết đo thực tế tín hiệu Việt Nam 98 4.5.1 Vùng phủ sóng theo chuẩn DVB-T2 98 4.5.2 So sánh kết mạng đơn tần mô đo dạc thực tế 99 KẾT LUẬN 101 TÀI LIỆU THAM KHẢO 102 PHỤ LỤC 103 DANH MỤC CHỮ VIẾT TẮT 2G 3G 4G AC ADSL Alpha AM APR ASK AWGN b/s/Hz BER Bps BPSK BS CD CDMA CF DAB DB DBc dBm DC DDS DFT DMT DPSK DSBSC DS-CDMA DSP DSSS DVB DVB-C DVB-S Second Generation mobile phone system (GSM, IS-95) Third Generation mobile phone system Fourth Generation mobile phone system Alternating Current ( Frequency above Hz) Asymmetric Digital Subscriber Line Path loss exponent ( rate of path loss change with distance) Amplitude Modulation Access Point Repeater Amplitude Shift Keying Additive White Gaussian Noise Bits per second per hertz ( unit of spectral efficiency) Bits Error Rate Bits per second Binary Phase Shift Keying Base Station Compact Station Code Division Multiple Access Crest Factor (peak to average power ratio of the RF envelope) Digital Audio Broadcasting Decibel ( ration in log scale) Decibel relative to main signal power Decibel relative to milliwatt Direct Current (0 Hz) Direct Digital Synthesiser Discrete Fourier Transform Discrete Multi – Tone Differential Phase Shift Keying Double Side Band Suppressed Carrier Direct Sequence Code Division Multiole Access Digital Signal Processing Direct Sequence Spread Spectrum Digital Video Broadcasting Digital Video Broadcasting – Cable Digital Video Broadcasting – Satellite DVB-T EBNR FDM FEC FFT FIR FM Fs FSK GA GHz GMSK GSM HDTV Digital Video Broadcasting – Terrestrial Energyper Bit to Noise Ratio Frequency Division Mulyiplexing Forward Error Correction Fast Fourier Transform Finite Impulse Response ( digital filter) Frequency Modulation Sample Frequency Frequency Shift Keying Genetic Algorithm Gigahertz – 109 Hz Gaussian Minimum Shift Keying Global System for Mobile Communications High Definition Televison High Performance Radio Local Area Network, WLAN satandard ( Europe) HiperLAN2 Based on OFDM, with a maximum data rate of 54 Mbps Similar to IEEE802.11a Hz Hertz (cycles per second) ICI Inter- Carrier Interference IEEE802.11a WLAN standard( U.S) base on OFMD, with a maximum data rate of 54 MBps.Similar to HiperLAN2 IEEE802.11b WLAN standard( U.S) base on DSSS, with a maximum data rate of 11 MBps IF Intermediate Frequency IFFT Inverse Fast Fourier Transform IMD Inter-Modulation Distortion IQ Inphase Quadreature ISI Inter-symbol Interference ISM Industrial Scientific Medical IS-95 Mobile phone standard using CDMA transmission metho JCU Kelvin Kbps Kilo bits per secon ( 103 bps) Km Kilometer( 103m) λ Lambda- Rfwavelength LO Local Oscillator LOS Line Of Sight M Metre Mbps Mega bits per second( 106 bps) Mhz Megahert – 106 Hz MPEG NF OBO OFDM PAPR PC PDA π PLL PM PRC PRS PSK QAM QOS QPSK RAKE RC RF RMS SHARC SNR SSB SSPA TDM TDMA TWTA µm UMTS µs VSB W W-CDMA WLAN Wll Moving Picture Experts Grop ( video compression standard) Receiver Noise Figure Output power BackOff Ortheogonal Frequency Division Multiplexing Peak to Average Power Ratio Personal Computer Personal Digital Assistant Pi(3,14…) Phase Locked Loop Phase Modulation Peak Reduction Carriers Pseudo Random Sequence Phase Shift Keying Quandreature Amplitude Modulation Quanlity Of Sevice Quandrature Plase Shift Keying Multiple tap multipath euqalisation for CDMA Raised Cosine (Guard Period) Radio Frequency Root Mean Squared Super Harvard ARChitecture, Digital Signal Processor by Analog Devices) Signal to Noise Ratio Single Side band SolidState Power Amplifier Time Division Multiplexing Time Division Multiple Access Travelling Wave Tube Amplifier Micrometre( 10-6m) Universal Mobile Telecommunication System Microsecond( 10-6s) Vestigal Side band Watt( energy per unit time, one joule per second) Wide-band Code Dvision Multiple Access Wireless Local Area Network Wireless Local Loop DANH MỤC BẢNG BIỂU Bảng 1.1 Mô tả thông số mode làm việc DVB – T Bảng 1.2 Tín hiệu điều chế QPSK 14 Bảng 1.3 Các hoán vị mode 2k 26 Bảng 1.4 Các hoán vị mode 8k 26 Bảng 1.5 Các giá trị khoảng phòng vệ 27 Bảng 4.1 Kết đo Cường độ trường vùng phủ sóng 99 10 Bảng 4.1 Kết đo Cường độ trường vùng phủ sóng Kết đo cường độ trường thông số tham khảo kết đo thực tế tín hiệu Việt Nam - Tơ Thị Thu Trang (trường ĐH Quốc Gia Hà Nội) 4.5.2 So sánh kết mạng đơn tần mô đo dạc thực tế 99 Kết luận: Trong chương cuối mô hệ thống OFDM simulink Matlab, với scope để thị tín hiệu giúp cho việc phân tính đánh giá tác động kênh truyền đến tín hiệu, tác dụng ước lượng bù kênh Tuy nhiên, simulink dừng lại mức độ đơn giản, tức mô hệ thống OFDM băng gốc với phương thức điều chế QPSK Trong chương so sánh tín hiệu OFDM tín hiệu QAM, file âm chúng để thấy rõ ưu điểm OFDM Từ thực tế đo được, ta nhận thấy kết đo kiểm thực tế điểm đo thỏa mãn yêu cầu đặt ra, Cùng điểm đo với cường độ trường theo chuẩn DVB- T2 lớn khoảng 50% so với DVB-T Kết đo kiểm thực tế mạng đơn tần theo chuẩn DVB-T2 có kết gần tiệm cận với kết theo tính tốn Với tính trội tiêu chuẩn DVB-T2, việc công ty AVG thực thành cơng phát sóng truyền hình số mặt đất theo chuẩn DVB-T2 tạo hội việc cung cấp dịch vụ truyền hình độ nét cao (HDTV), 3DTV, dịch vụ truyền hình di động, dịch vụ khác tương lai Đặc biệt, đóng góp đáng kể việc thực mạng SFN nhằm tiết kiệm tài nguyên tần số quốc gia 100 KẾT LUẬN OFDM công nghệ cho phép ghép kênh tiết kiệm phổ Chúng ta biết tín hiệu truyền kênh truyền mặt đất bị ảnh hưởng nhiều hiệu ứng như: nhiễu xung, fading, nhiễu kênh kề, nhiễu kênh chung trễ lan truyền đa đường… vấn đề quan tâm ảnh hưởng tín hiệu trễ lên tín hiệu gốc gây giao thoa kí hiệu ISI Trễ tạo lan truyền theo nhiều đường tín hiệu mặt đất, ảnh hưởng trễ đến chất lượng tín hiệu thu mạnh trễ biến đổi nhanh Hiện công nghệ OFDM nghiên cứu ứng dụng lĩnh vực thông tin vô tuyến, công nghệ kết hợp phương pháp điều chế cổ điển phương pháp đa truy cập vô tuyến Việc triển khai ứng dụng truyền hình số mặt đất DVB- T2 được diễn trung tâm truyền hình lớn nước ta Xin chân thành cảm ơn Thầy giáo hướng dẫn PGS.TS Nguyễn Hữu Trung, Thầy hội đồng bảo vệ, thầy, cô Viện Đại học Mở Hà Nội Khoa đào tạo sau đại học, bạn bè, đồng nghiệp giúp đỡ em hoàn thành luận văn 101 TÀI LIỆU THAM KHẢO Richard Van Nee & Ramjee Prasad.OFDM For Wireless multimedia communications (Artech House, 2000) ETSI EN 300744 (2004/06) Digital Video Broadcasting (DVB); Frame structure, channel coding and modulation for digital terrestrial television Dambacher, P (1998) Digital Video Broadcasting The international Standard for Digital Television Digital Broadcasting Seminar, Hanoi 21 july 2003 (ROHDER & SCHWARZ) DVB’S Comments to the Submission for the ANATEL CONSULTA PUBLICA N0.291 - BRASILIAN DIGITAL TERRESTRIAL TELEVISION (DTV) SYSTEM The How and why of COFDM, J.H Stott, BBC Research and Development ( 1995) http://wikipedia.com, http://DVB.org… http://www.citd.edu.vn Phạm Đắc Bi, Lê Trọng Bằng, Đỗ Anh Tú, ”Các đặc điểm máy phát số DVB-T”, Tạp chí Bưu Chính Viễn Thơng & Cơng Nghệ Thơng Tin, (8/2004) 10 Ngơ Thái Trị - Truyền hình số Nhà xuất Đại học Quốc gia Hà Nội 2004 102 PHỤ LỤC • Chuẩn hóa liệu trước sử dụng If channel_on=1 Disp(‘Simulating Channel’) Norm factor = max(abs(recv)); Recv = (1/norm_factor)*recv; ch_clipping ch_multipat ch_noise recv = norm_factor*recv; end • Mơ hiệu ứng xén tín hiệu For I = 1:length(recb) if recv(i) > clip_level recv(i) = clip_level; end if recv(i) < -clip_level recv(i) = -clip_level; end end • Mơ hiệu ứng đa đường Copy1=zeros(size(recv)); For I = 1+d1:length(recv) Copy1(i)=a1*recv(i-d1); End Copy2=zeros(size(recv)); For i=1+d2:length(recv) Copy2(i)=a2*recv(i-d2); End 103 Recv=recv+copy1+copy2; • Tính tốn nhiễu( thực phía thu) If already_made_noise==0 % only generrate once and use for both QAM and OFDM Noise = (rand(1,length(recv))-0.5)*2*noise_level; Already_made_noise = 1; End Recv = recv+noise; % khơi phục biên độ liệu • Phát symbol OFDM Disp(‘Transmiting’) Read % Đọc liệu vào Data_in_pol = bin2pol(data_in) % Chuyen doi du lieu nhi phan du lieu phan cuc Tx_chunk • Thực IFFT để tạo dạng sóng miền thời gian biểu diễn liệu Td_sets = zeros(num_chunks,fft_size); For I = 1:num_chunks Td_sets(I,1:fft_size) = real(ifft(spaced_chunks(I,1:iff_size))); End Tx_dechunk • Đổi liệu nhị phân(0,1) thành liệu cực(-1,1) Y = ones(1,length(x)); For i= 1:length(x) If x(i) == Y(i)= -1 End End 104 • Đổi liệu nhị phân thành hexa Y=0 K=0 For I = 1:8 Y = y+ x(8-k)*2^k; K = k+1 End • **************************** % Simulation of digital M-PSK modulation schemes over an AWGN channel % November 2004 Robert Morelos-Zaragoza San José State University % Needs the Communications toolbox clear Fd= 1; Fs = 1; N = 100000 ; method-psk'; set(l) = 2; set(2) = 4; set(3) = 8; for j=l:l:3 M = set(j); i=i; for esno=0:l:18 sigma = sqrt(10A(-esno/10)/2); x = floor(M*rand(N,l)); y = modmap(x,Fd,Fs,method,M); ynoisy = y + sigma*randn(N*Fs,2); z = demodmap(ynoisy,Fd,Fs,method,M); s = symerr(x,z); ber(j,i) = (s/N)/log2(M); snr(i) - esno; i=i+l; end 105 j end semilogy(snr,ber(l,:),'-bA',snr,ber(2,:),'-bo',snr,ber(3,:),'-bs'); grid on; ylabel('BER'); xlabel('E/N_0 (dB)'); legend('BPSK', 'QPSK', ’8PSK',1); title('Simulated error performance of M-PSK SJSU - Fall 2004.'); • ******************************************** % Simulation of digital M-QAM modulation schemes over an AWGN channel % November 2004 Robert Morelos-Zaragoza San Jose State University % Needs the Communications toolbox clear Fd= 1; Fs = 1; N= 100000 ; method-qask'; set(l) = 4; set(2) =16; set(3) = 64; for j=l:l:3 M = set(j); 1=1:1 :M; aux = sum(abs(modmap(l-l,Fd,Fs,method,M)).A2)/M; energy(j) = aux(l)+aux(2); i=l; for esno=0:2:26; sigma = sqrt(10A(-esno/10)/2)*sqrt(energy(j)); x = floor(M*rand(N,l)); y = modmap(x,Fd,Fs,method,M); ynoisy = y + sigma*randn(N*Fs,2); z = demodmap(ynoisy,Fd,Fs,method,M); s = 106 symerr(x,z); ber(j,i) = (s/N)/log2(M); snr(i) esno; i=i+l; end j end semilogy(snr,ber(l,:),'-bA',snr,ber(2,:),'-bo',snr,ber(3,:),'-bs'); grid on; ylabel('BER'); xlabel('E/N_0 (dB)'); legend('QPSK', '16-QAM', ’64-QAM',l); title('Cac kieu dieu che M-QAM khac kenh truyen AWGN'); • *********************************************** function y = eight2bin(x) % eight2bin % % Converts eight bit data (0-255 decimal) to a binary form for processing y = zeros(l,8); k = 0; while x > y(8-k) = rem(x,2); k k+1; x = floor(x/2); end • *********************************************** % Run OFDM simulation tic % Start stopwatch to calculate how long QAM simulation takes disp(' '),disp(' ') disp('OFDM Simulation') tx eh rx % Stop stopwatch to calculate how long QAM 107 simulation takes OFDM_simulation_time - toe; if OFDM_simulation_time > 60 disp(strcat('Time for OFDM simulation-, num2str(OFDM_simulation_time/60),' minutes.')); else disp(strcat('Time for OFDM simulation-, num2str(OFDM_simulation_time),' seconds.')); end function y = pol2bin(x) % pol2bin % % Chuyen doi cac so phan cue (-1,1) cac so nhi phan (0,1) % Chap nhan mot mang -D cua cac so phan cue % Loai bo cac zeros, chung khong hop le % % Loai bo cac zeros Khong can cho giai ma y = ones(l,length(x)); for i = :length(x) if x(i) —— -1 y(i) = 0; end end • ******************************************* % QAM.m So sanh OFDM (Da song mang) voi QAM da muc (Don song mang) % Khi chung ta phat cung mot so luong bit giong ngau tren mot chu ky thoi % gian 108 Read % Doc du lieu cho QAM - Khong anh huong den OFDM data_in_pol = bin2pol(data_in); % Chuyên doi du lieu nhi phan du lieu phan cue % Kiem tra so song mang co phai la luy thua cua is_pow_2 = num_carriers; tempdoQAM = 0; if is_pow_2 ~= while temp_do_QAM — tempdoQAM = rem(is_pow_2,2); is_pow_2 = is_pow_2/2; if is_pow_2 — tempdoQAM = -99; end end else temp_do_QAM = -99; % la luy thua cua end if tempdoQAM ~= -99 do_QAM = 0; % Khong the thuc hien disp(' '),disp('ERROR: Cannot run QAM because num_carriers is not valid.') dispC Please see "setup.m" for details.') end if doQAM = tic % Bat dau de tinh toan thoi gian mo phong thuc hien QAM disp(‘ ‘), disp(' ‘) disp(‘QAM simulation'), disp('Transmitting') • Them cac muc zeros de du lieu duoc chia cac phan bang data_length = length(data_in_pol); r = rem(data_length,num_carriers); if r ~= 109 for i = :num_carriers-r data_in_pol(data_length+i) = 0; %Them dau vao voi cac zeros vao tap hop du lieu end %Toc CO the duoc cai thien end data_length = length(data_in_pol); %Cap nhat sau them num_OFDM_symbols = ceil(data_length / (2*num_carriers)); % So ky hieu QAM duoc bieu dien bang so luong cua du lieu tren mot ky % hieu OFDM num_QAM_symbols = num_carriers / 2; % So mau tren ky hieu QAM num_symbol_samples - fft_size / num_QAM_symbols; • Chuyen doi du lieu phan cue [-1, 1] du lieu muc [-3, -1, 1, 3] data_in_4 = zeros(l,data_length/2); for i = 1:2:data_length data_in_4(i - (i-1 )/2) = data_in_pol(i)*2 + data_in_pol(i+l); end % Dinh ro diem lay mau giua va 2*pi ts = linspace(0, 2*pi*QAM_periods, num_symbol_samples+l); % Phat du lieu 16-QAM % Tong dai cua truyen dan 16-QAM tx_length - num_OFDM_symbols * num_QAM_symbols * num_symbol_samples; QAM_tx_data = zeros(l,tx_length); for i = l:2:data_length/2 for k = :num_symbol_samples QAM_tx_data(k+((i-1)/2)*num_symbol_samples) data_in_4(i)*cos(ts(k)) + data_in_4(i+l)*sin(ts(k)); end 110 end % Do channel simulation on QAM data xmit = QAM _tx _data; % ch dung du lieu ‘xmit’ va tra ve ‘ recv’ch QAM_ rx _data - recv; % Luu du lieu QAM sau mo phong kenh clear recv % Loai bo ‘recv’ cho no khong nhieu voi OFDM clear xmit % Loai bo ‘xmit’ cho no khong nhieu voi ODFM disp('Receiving') % Khoi phuc du lieu nhi phan (Giai ma QAM) cos_temp - zeros( ,num_symbol_samples); % sin_temp - cos_temp; %xxx = zeros(l,data_length/4); % Khoi tao muc khong cho toe yyy - xxx; % QAM_data_out_4 = zeros( ,data_length/2); % for i = l:2:data_length/2 % "cheating" for k = :num_symbol_samples • Yang so song mang de tao tan so cao va du lieu goc cos_temp(k) - QAM_rx_data(k+((i- l)/2)*num_symbol_samples) * cos(ts(k)); sin_temp(k) - QAM_rx_data(k+((il)/2)*num_symbol_samples) * sin(ts(k)); end % LPF va xac dinh - chung ta se rat don gian LPF bang phep trung % binh xxx(l+(i-l)/2) = mean(cos_temp); yyy(l+(i-l)/2) = mean(sin_temp); 111 • Khoi phuc du lieu dang noi tiep QAM_data_out_4(i) = xxx(l+(i-l)/2); QAM_data_out_4(i+l) = yyy(l+(i-l)/2); end %Tinh toan giua cac muc khong zeros_between = ((fft_size/2) - (num_carriers + num_zeros))/(num_carriers + num_zeros); spaced_chunks = zeros(num_chunks,fft_size); %Them vao giua cac muc khong i= l; for k = zeros_between +1 :zeros_between +1 :fft_size/2 spaced_chunks(l:num_chunks,k) = padded_chunks(l:num_chunks,i); i = i+1; end % Gap du lieu de tao mot ham le cho dau vao IFFT for i = :num_chunks % Chu y: chi muc = that la tan so mot chieu de ifft -> no khong % tao ban len true y thi spaced_chunks(i,fft_size:-l :fft_size/2+2) = conj (spaced_chunks(i,2: fft_size/2)); end • Thuc hien xac dinh giua cac muc [-3, -1, 1, 3] for i = :data_length/2 if QAM_data_out_4(i) >= 1, QAM_data_out_4(i) - 3; elseif QAM_data_out_4(i) >= 0, QAM_data_out_4(i) - 1; elseif QAM_data_out_4(i) >= -1, QAM_data_out_4(i) = -1; else QAM_data_out_4(i) - -3; end 112 • Chuyen doi du lieu muc [-3, -1, 1, 3] ve du lieu phan cuc [-1, i] QAM_data_out_pol = zeros(l,data_length); % "cheating" for i = 1:2:data_length switch QAM_data_out_4(l + (i-l)/2) case -3 QAM_data_out_pol(i) = -1; QAM_data_out_pol(i+l) = -1; case -1 QAM_data_out_pol(i) = -1; QAM_data_out_pol(i+l) - 1; case QAM_data_out_pol(i) - 1; QAM_data_out_pol(i+l) = -1; case QAM_data_out_pol(i) - 1; QAM_data_out_pol(i+l) - 1; otherwise disp('Error detected in switch statment - This should not be happening.'); end end QAM_data_out = pol2bin(QAM_data_out_pol); % Chuyen doi ve du lieu nhi phan 113