(Luận văn) đánh giá hiệu năng hệ thống truyền hình số mặt đất dvb t2 sử dụng kỹ thuật mimo ofdm

HỌC VIỆN CƠNG NGHỆ BƯU CHÍNH VIỄN THƠNG - lu an va n Nguyễn Mạnh Tiến p ie gh tn to ĐÁNH GIÁ HIỆU NĂNG HỆ THỐNG TRUYỀN HÌNH nl w SỐ MẶT ĐẤT DVB – T2 SỬ DỤNG KỸ THUẬT d oa MIMO - OFDM u nf va an lu ll LUẬN VĂN THẠC SĨ KỸ THUẬT m oi (Theo định hướng ứng dụng) z at nh z m co l gm @ an Lu HÀ NỘI – 2019 n va ac th si HỌC VIỆN CÔNG NGHỆ BƯU CHÍNH VIỄN THƠNG - lu an n va Nguyễn Mạnh Tiến SỐ MẶT ĐẤT DVB – T2 SỬ DỤNG KỸ THUẬT p ie gh tn to ĐÁNH GIÁ HIỆU NĂNG HỆ THỐNG TRUYỀN HÌNH d oa nl w MIMO - OFDM Chuyên ngành: Kỹ Thuật Viễn Thông an lu ll u nf va Mã Số: 8.52.02.08 oi m LUẬN VĂN THẠC SĨ KỸ THUẬT z at nh NGƯỜI HƯỚNG DẪN KHOA HỌC: z PGS.TS VŨ VĂN SAN m co l gm @ HÀ NỘI – 2019 an Lu n va ac th si i 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 lu Nguyễn Mạnh Tiến an n va p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si ii LỜI CẢM ƠN Đầu tiên, học viên xin gửi lời cảm ơn chân thành đến tất thầy cô Khoa Đào tạo Sau Đại học - Học viện Cơng nghệ Bưu Viễn thơng ln nhiệt tình hướng dẫn, truyền đạt kiến thức suốt thời gian học tập Học viện, tảng giúp học viên thực luận văn tốt nghiệp Học viên xin chân thành cảm ơn Thầy giáo, PGS.TS Vũ Văn San – Giám đốc Học viện tận tình hướng dẫn học viên hồn thành luận văn Học viên xin chân thành cảm ơn bạn bè sát cánh giúp học viên có lu an kết ngày hơm va n Đề tài nghiên cứu luận văn có nội dung bao phủ rộng Tuy nhiên, thời gian tn to nghiên cứu cịn hạn hẹp Vì vậy, luận văn có thiếu sót Học viên ie gh mong nhận đóng góp ý kiến thầy cô bạn p Xin chân thành cảm ơn! w d oa nl Tác giả luận văn u nf va an lu Nguyễn Mạnh Tiến ll oi m z at nh z m co l gm @ an Lu n va ac th si iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT v DANH MỤC CÁC BẢNG viii DANH MỤC CÁC HÌNH ix lu MỞ ĐẦU .1 an n va CHƯƠNG 1-TỔNG QUAN VỀ TRUYỀN HÌNH SỐ MẶT ĐẤT DVB – T Tổng quan hệ thống truyền hình số Các đặc điểm chung hệ thống truyền hình số 1.1.1 ie gh tn to 1.1 Các tiêu chuẩn truyền hình số tiêu biểu p 1.1.2 Xử lý truyền dẫn tín hiệu truyền hình số 1.2 oa nl w 1.1.3 Truyền hình số mặt đất theo tiêu chuẩn DVB – T .11 d lu Đặc tính kỹ thuật DVB – T .12 1.2.2 Đặc điểm DVB – T 14 u nf va Truyền hình số mặt đất theo tiêu chuẩn DVB – T2 15 ll 1.3 an 1.2.1 oi m Đặc tính kỹ thuật DVB – T2 15 1.3.2 Đặc điểm DVB – T2 .19 Kết luận chương .21 z @ 1.4 z at nh 1.3.1 l gm CHƯƠNG – KỸ THUẬT MIMO – OFDM TRONG TRUYỀN HÌNH SỐ MẶT ĐẤT DVB – T2 23 m co Tổng quan kỹ thuật MIMO .23 2.2 Mơ hình kênh MIMO dung lượng kênh MIMO .24 an Lu 2.1 n va ac th si iv 2.3 Kỹ thuật MIMO – OFDM 26 2.4 Kỹ thuật ghép kênh không gian 32 2.5 Hoạt động hệ thống truyền hình số mặt đất DVB – T2 sử dụng kỹ thuật MIMO – OFDM 35 Kết luận chương .37 2.6 CHƯƠNG – PHÂN TÍCH, ĐÁNH GIÁ HIỆU NĂNG HỆ THỐNG DVB – T2 SỬ DỤNG KỸ THUẬT MIMO – OFDM 38 lu an n va Mơ hình hệ thống DVB – T2 sử dụng kỹ thuật MIMO – OFDM .38 3.2 Phân tích, đánh giá kết 42 3.3 Khuyến nghị, đề xuất 49 3.4 Kết luận chương .49 gh tn to 3.1 p ie KẾT LUẬN .50 w DANH MỤC TÀI LIỆU THAM KHẢO 51 d oa nl PHỤ LỤC 52 ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si v DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT Viết tắt Tiếng Anh ACE Active Constellation Extension Mở rộng chịm tích cực Advanced Television System Ủy ban hệ thống truyền Committee hình tiên tiến ATSC Tiếng Việt Additive White Gaussian Noise Tạp âm Gausse trắng cộng AWGN British Broadcasting BBC Hiệp hội phát Anh quốc Corporation Bit Error Rate Tỷ lệ lỗi bit BPSK Binary Phase Shift Keying Điều chế pha nhị phân CP Cyclic Prefix Tiền tố tuần hồn C/N Carrier/Noise Sóng mang/tạp âm Coding Othogonality Fequency Ghép kênh phân chia theo tần số Dvision Mltiplexing trực giao có mã Digital Terrestrial Television Truyền hình số mặt đất Digital Video Broadcasting Truyền hình số lu BER an n va gh tn to p ie COFDM DTT Digital Video Broadcasting - d Truyền hình kỹ thuật số qua cáp Cable va an lu Digital Video Broadcasting - Truyền hình số mặt đất ll Terrestrial u nf DVB-T oa DVB-C nl w DVB m Digital Video Broadcasting - z at nh Second Generation Terrestrial NGH Handheld Ủy ban truyền thông liên bang-Mỹ an Lu Commission m co Federal Communications Truyền hình kỹ thuật số - vệ tinh l Satellite gm FCC Digital Video Broadcasting - Thiết bị cầm tay hệ @ DVB - Next Generation z DVB- DVB-S Truyền hình số mặt đất hệ thứ oi DVB-T2 n va ac th si vi FEC Forward Error Correction Sửa lỗi trước (thuận) FFT Fast Fourrier Transform Biến đổi Fourrier nhanh High Definitiom Television Truyền hình độ phân giải cao HDTV Independent and Identically IID Phân bố độc lập đồng dạng Distributed Inverse Fast Fourrier IFFT Biến đổi nhanh Fourrier ngược Transform Internet Protocol Television Truyền hình IP IRD Integrated Decoder Bộ giải mã tích hợp Integrated Service Didital Truyền hình số tích hợp dịch vụ Broadcasting – Terrestrial mặt đất Inter-Symbol Interference Nhiễu liên ký tự Low Density Parity Check Mã kiểm tra chắn lẻ mật độ (Codes) thấp Low Priority Độ ưu tiên thấp Multi Input Multi Output Nhiều đầu vào nhiều đầu lu IPTV an ISDB-T n va LDPC p ie gh tn to ISI LP Giao diện điều chế Modulator Interface d oa MI nl w MIMO Nhiều đầu vào đầu Multi Input Single Output MPEG Moving Pictures Experts Group u nf va an lu MISO Hội đồng hệ thống truyền hình Committee quốc gia Mỹ oi m z at nh Othogonality Fequency OFDM ảnh động National Television System ll NTSC Nhóm chuyên gia nghiên cứu Dvision Ghép kênh phân chia theo tần số z trực giao @ OSI Open Systems Interconnection Mơ hình tham chiếu PAL Phase Alternative Line Đảo pha theo dịng m co Tỷ số cơng suất đỉnh/cơng suất trung bình an Lu Peak to Average Power Ratio l PAPR gm Mltiplexing n va ac th si vii PLP Ống lớp vật lý Physical Layer Pipes Quadrature Amplitude QAM Điều chế biên độ vng góc Modulation QPSK RF SDTV Quadrature Phase Shift Keying Khóa dịch pha vng góc Radio Frequence Tần số vơ tuyến Truyền hình độ phân giải tiêu Standard Definition Television chuẩn an Single Frequence Network Mạng đơn tần SNR Signal to Noise Ratio Tỷ số tín hiệu tạp âm SISO Single Input Single Output Một đầu vào đầu Standard-Definition Độ nét tiêu chuẩn STBC Space–time block coding Mã hóa không gian – thời gian STTC Space Time Trellis Code Mã hóa lưới khơng gian - thời gian SM Spatial Multiplexing Ghép kênh không gian p lu SFN TR Tone Reservation Hạn chế âm sắc TV Television Truyền hình n va SD ie gh tn to nl w Ultra-high-definition d Ultra High Frequency VHF Very High Frequency Tần số siêu cao (Siêu cao tần) an lu UHF Tần số cao va Kỹ thuật tách tín hiệu phía thu, u nf Vertical-Bell Laboratories phát triển phịng thí nghiệm ll Layered Space-Time m V_BLAST Siêu nét oa UHD oi Bell Vietnam Television Điều biên cụt Tổng Công ty Truyền thông đa z VTC Vestigal Side Band z at nh VSB phương tiện Việt Nam m co l gm @ Corporation an Lu n va ac th si viii DANH MỤC CÁC BẢNG Bảng 1.1 Đặc điểm ATSC .5 Bảng 1.2 Các thông số truyền dẫn cho ISDB-T với độ rộng kênh truyền MHz .7 Bảng 3.1 Các thông số mô hệ thống .41 Bảng 3.2 Các thông số đầu vào hệ thống 44 Bảng 3.3 Thống kê kết mô trường hợp .45 Bảng 3.4 Thống kê kết mô trường hợp .47 lu Bảng 3.5 Thống kê kết mô trường hợp .48 an n va p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si 54 close all%xoa tat cac hinh clc%xoa man hinh Dlen=1024; % Tat ca du lieu symbol duoc truyen Flen=64; % 64 du lieu symbol tren khung duoc bien doi IFFT IFFTlen=2*Flen; % chieu dai 128 symbol cho IFFT F=Dlen/Flen; % so luong frame G=16; % chieu dai chuoi bao ve x1=zeros(Dlen/2,1);%tao chuoi symbol sau duoc ma hoa qpsk y=zeros(IFFTlen,1); lu an y1=zeros(IFFTlen,1); n va z=zeros(F*IFFTlen,1);%chuoi symbol sau duoc chuyen doi song song noi tiep % tao nguon tin hieu gh tn to ofdm_sg=zeros(F*IFFTlen+G,1);%chuoi symbol sau duoc ghep chuoi bao ve p ie x=randi(Dlen,1,2);%tao chuoi tin hieu va bat ki % chuyen sang mien phuc d oa for i=1:F nl w %x1=qpsk_dc(x);%ham chuyen doi qpsk va for j=1:Flen/2 an lu % chuyen doi noi tiep song song ll % chuyen doi song song noi tiep z at nh y1=ifft(y,length(y));%bien doi ifft oi % dua ve cac song truc giao m end u nf y(j+IFFTlen/2-Flen/4)=x1((i-1)*Flen/2+j); z an Lu %chen chuoi bao ve m co end l end gm z((i-1)*IFFTlen+j)=y1(j); @ for j=1:IFFTlen n va ac th si 55 for i=1:G ofdm_sg(i)=z(F*IFFTlen-G+i); end for i=1:F*IFFTlen ofdm_sg(G+i)=z(i); end % Do thi mien thoi gian %figure(1); %subplot(2,1,1); lu an %stem(real(z),'r.'); n va %xlabel('chi so thoi gian');ylabel('chi so bien do'); %subplot(2,1,2); gh tn to %title('phan thuc'); p ie %stem(imag(z),'r.'); %xlabel('chi so thoi gian');ylabel('chi so bien do'); %figure(2); d oa nl w %title('phan ao'); an lu %f = linspace(-Flen,Flen,length(z)); u nf mien thoi gian va %plot(f,abs(z));%bieu dien gia tri tuyet doi cua tin hieu sau dieu che ofdm theo ll %title('gia tri tuyet doi cua tin hieu mien thoi gian'); m oi y2 = fft(z);%chuyen tin hieu sang mien tan so z at nh %xlabel('chi so thoi gian'); %ylabel('chi so tuyet doi cua tin hieu'); z an Lu end m co %y2(j)=0.01; l if abs(y2(j)) < 0.01 gm for j=1:F*IFFTlen; @ % neu Y2 la be hon 0.01 thi Y2=0.001 n va ac th si 56 end %y2 = 10*log10(abs(y2));%chuyen y2 sang dB % Do thi mien tan so %figure(3); %f = linspace(-Flen,Flen,length(y2)); %plot(f,y2); %axis([-IFFTlen/2 IFFTlen/2 -20 20]); %xlabel('truc tan so'); %ylabel('chi so tuyet doi cua tin hieu theo dB'); lu an %title('pho cua tin hieu mien tan so'); n va %%%%%%%%%%%%%%%%%%%%%%%% %function y=dc64qam(x) gh tn to %Dieu che 64-QAM%%%%% p ie for i=1:4:length(x) if x(i)==0 & x(i+1)==0 & x(i+2)==0 & x(i+3)==0 nl w y((i+3)/4)=exp(j*pi/4); d oa elseif x(i)==1 & x(i+1)==0 & x(i+2)==0 & x(i+3)==0 an lu y((i+3)/4)=exp(j*3*pi/4); y((i+3)/4)=exp(j*5*pi/4); u nf va elseif x(i)==1 & x(i+1)==1 & x(i+2)==0 & x(i+3)==0 ll elseif x(i)==0 & x(i+1)==1 & x(i+2)==0 & x(i+3)==0 oi m y((i+3)/4)=exp(j*7*pi/4); z at nh elseif x(i)==0 & x(i+1)==0 & x(i+2)==1 & x(i+3)==0 y((i+3)/4)=2*exp(j*pi/6); z y((i+3)/4)=2*exp(j*pi/4); y((i+3)/4)=2*exp(j*pi/3); an Lu elseif x(i)==1 & x(i+1)==0 & x(i+2)==0 & x(i+3)==1 m co elseif x(i)==0 & x(i+1)==0 & x(i+2)==0 & x(i+3)==1 l gm @ elseif x(i)==0 & x(i+1)==0 & x(i+2)==1 & x(i+3)==1 n va ac th si 57 y((i+3)/4)=2*exp(j*2*pi/3); elseif x(i)==1 & x(i+1)==0 & x(i+2)==1 & x(i+3)==1 y((i+3)/4)=2*exp(j*3*pi/4); elseif x(i)==1 & x(i+1)==0 & x(i+2)==1 & x(i+3)==0 y((i+3)/4)=2*exp(j*5*pi/6); elseif x(i)==1 & x(i+1)==1 & x(i+2)==1 & x(i+3)==0 y((i+3)/4)=2*exp(j*7*pi/6); elseif x(i)==1 & x(i+1)==1 & x(i+2)==1 & x(i+3)==1 y((i+3)/4)=2*exp(j*5*pi/4); lu an elseif x(i)==1 & x(i+1)==1 & x(i+2)==0 & x(i+3)==1 n va y((i+3)/4)=2*exp(j*4*pi/3); y((i+3)/4)=2*exp(j*5*pi/3); gh tn to elseif x(i)==0 & x(i+1)==1 & x(i+2)==0 & x(i+3)==1 p ie elseif x(i)==0 & x(i+1)==1 & x(i+2)==1 & x(i+3)==1 y((i+3)/4)=2*exp(j*7*pi/4); nl w elseif x(i)==0 & x(i+1)==1 & x(i+2)==1 & x(i+3)==0 z at nh n = 2e5; oi k = log2(M); m M = 32; ll clc; u nf va clear all; an end lu end d oa y((i+3)/4)=2*exp(j*11*pi/6); z an Lu hDemod = modem.qamdemod(hMod); m co hMod.SymbolOrder = 'Gray'; l hMod.InputType = 'Bit'; gm hMod = modem.qammod(M); @ nSyms = n/k; n va ac th si 58 x = randi([0 1],n,1); tx = modulate(hMod,x); EbNo = 0:30; SNR = EbNo +5- 10*log10(k); rx = zeros(nSyms,length(SNR)); bit_error_rate = zeros(length(SNR),1); for i=1:length(SNR) rx(:,i) = awgn(tx,SNR(i),'measured'); end lu an rx_demod = demodulate(hDemod,rx); n va for i=1:length(SNR) end gh tn to [~,bit_error_rate(i)] = biterr(x,rx_demod(:,i)); p ie EbNo1 = 0:30; % In dB SNR1 = EbNo1 +8- 10*log10(k); nl w rx = zeros(nSyms,length(SNR1)); d oa bit_error_rate1 = zeros(length(SNR1),1); an lu for i=1:length(SNR1) u nf end va rx(:,i) = awgn(tx,SNR1(i),'measured'); ll rx_demod = demodulate(hDemod,rx); oi m for i=1:length(SNR1) z at nh [~,bit_error_rate1(i)] = biterr(x,rx_demod(:,i)); end z an Lu for i=1:length(SNR2) m co bit_error_rate2 = zeros(length(SNR2),1); l rx = zeros(nSyms,length(SNR2)); gm SNR2 = EbNo2 +10- 10*log10(k); @ EbNo2 = 0:29; % In dB n va ac th si 59 rx(:,i) = awgn(tx,SNR2(i),'measured'); end rx_demod = demodulate(hDemod,rx); for i=1:length(SNR2) [~,bit_error_rate2(i)] = biterr(x,rx_demod(:,i)); end EbNo3 = 0:30; % In dB SNR3 = EbNo3 +12- 10*log10(k); rx = zeros(nSyms,length(SNR3)); lu an bit_error_rate3 = zeros(length(SNR3),1); n va for i=1:length(SNR3) end gh tn to rx(:,i) = awgn(tx,SNR3(i),'measured'); p ie rx_demod = demodulate(hDemod,rx); for i=1:length(SNR3) an lu figure; d oa end nl w [~,bit_error_rate3(i)] = biterr(x,rx_demod(:,i)); va semilogy(EbNo, bit_error_rate,'^-',EbNo1, bit_error_rate1, 'o-',EbNo2, ll m grid on; u nf bit_error_rate2, '*-',EbNo3, bit_error_rate3, '+-'); oi legend('Kenh TU6','Kenh P1','Kenh F1','Kenh Gaussian'); ylabel('Bit Error Rate'); z at nh xlabel('Eb/No (dB)'); z clear all an Lu n = 2e5; m co k = log2(M3); l M3 = 32; gm @ title('BER CHO CAC KENH MO PHONG KHAC NHAU'); n va ac th si 60 nSyms = n/k; hMod = modem.qammod(M3); hMod.InputType = 'Bit'; hMod.SymbolOrder = 'Gray'; hDemod = modem.qamdemod(hMod); x = randi([0 1],n,1); tx = modulate(hMod,x); EbNo7 = 0:30; SNR7 = EbNo7 +6.3- 10*log10(k); lu an rx = zeros(nSyms,length(SNR7)); n va bit_error_rate7 = zeros(length(SNR7),1); rx(:,i) = awgn(tx,SNR7(i),'measured'); gh tn to for i=1:length(SNR7) p ie end rx_demod = demodulate(hDemod,rx); nl w for i=1:length(SNR7) EbNo8 = 0:30; % In dB va an lu end d oa [~,bit_error_rate7(i)] = biterr(x,rx_demod(:,i)); u nf SNR8 = EbNo8 +11.4- 10*log10(k); ll rx = zeros(nSyms,length(SNR8)); m oi bit_error_rate8 = zeros(length(SNR8),1); z at nh for i=1:length(SNR8) rx(:,i) = awgn(tx,SNR8(i),'measured'); z rx_demod = demodulate(hDemod,rx); an Lu end m co [~,bit_error_rate8(i)] = biterr(x,rx_demod(:,i)); l for i=1:length(SNR8) gm @ end n va ac th si 61 EbNo9 = 0:29; % In dB SNR9 = EbNo9 +11.6- 10*log10(k); rx = zeros(nSyms,length(SNR9)); bit_error_rate9 = zeros(length(SNR9),1); for i=1:length(SNR9) rx(:,i) = awgn(tx,SNR9(i),'measured'); end rx_demod = demodulate(hDemod,rx); for i=1:length(SNR9) lu an [~,bit_error_rate9(i)] = biterr(x,rx_demod(:,i)); n va end semilogy(EbNo7, bit_error_rate7,'^-',EbNo8, bit_error_rate8, 'o-',EbNo9, gh tn to figure; p ie bit_error_rate9, '*-'); grid on; nl w legend('SISO','MISO','MIMO'); d oa xlabel('Eb/No (dB)'); an lu ylabel('Bit Error Rate'); ll M2 = 32; z at nh clear; oi clc; m clear all u nf Eb/No'); va title('BER TREN KENH MIMO VA MISO TOT HON KENH SISO TREN CUNG z an Lu hMod.InputType = 'Bit'; m co hMod = modem.qammod(M2); l nSyms = n/k; gm n = 2e5; @ k = log2(M2); n va ac th si 62 hMod.SymbolOrder = 'Gray'; hDemod = modem.qamdemod(hMod); x = randi([0 1],n,1); tx = modulate(hMod,x); EbNo4 = 0:30; SNR4 = EbNo4 +6.3- 10*log10(k); rx = zeros(nSyms,length(SNR4)); bit_error_rate4 = zeros(length(SNR4),1); for i=1:length(SNR4) lu an rx(:,i) = awgn(tx,SNR4(i),'measured'); n va end for i=1:length(SNR4) gh tn to rx_demod = demodulate(hDemod,rx); end p ie [~,bit_error_rate4(i)] = biterr(x,rx_demod(:,i)); nl w EbNo5 = 0:30; % In dB d oa SNR5 = EbNo5 +7- 10*log10(k); an lu rx = zeros(nSyms,length(SNR5)); u nf for i=1:length(SNR5) va bit_error_rate5 = zeros(length(SNR5),1); ll rx(:,i) = awgn(tx,SNR5(i),'measured'); oi m end z at nh rx_demod = demodulate(hDemod,rx); for i=1:length(SNR5) z end an Lu rx = zeros(nSyms,length(SNR6)); m co SNR6 = EbNo6 +12.5- 10*log10(k); l EbNo6 = 0:29; % In dB gm @ [~,bit_error_rate5(i)] = biterr(x,rx_demod(:,i)); n va ac th si 63 bit_error_rate6 = zeros(length(SNR6),1); for i=1:length(SNR6) rx(:,i) = awgn(tx,SNR6(i),'measured'); end rx_demod = demodulate(hDemod,rx); for i=1:length(SNR6) [~,bit_error_rate6(i)] = biterr(x,rx_demod(:,i)); end figure; lu an semilogy(EbNo4, bit_error_rate4,'^-',EbNo5, bit_error_rate5, 'o-',EbNo6, n va bit_error_rate6, '*-'); legend('SISO','MISO','MIMO'); gh tn to grid on; p ie xlabel('Eb/No (dB)'); ylabel('Bit Error Rate'); ll hMod = modem.qammod(M4); z at nh nSyms = n/k; oi n = 2e5; m k = log2(M4); u nf M4 = 32; va clc; an lu clear all d oa Eb/No'); nl w title('BER TREN KENH MIMO TOT HON KENH MISO VA SISO TREN CUNG z an Lu tx = modulate(hMod,x); m co x = randi([0 1],n,1); l hDemod = modem.qamdemod(hMod); gm hMod.SymbolOrder = 'Gray'; @ hMod.InputType = 'Bit'; n va ac th si 64 EbNo10 = 0:30; SNR10 = EbNo10 +12- 10*log10(k); rx = zeros(nSyms,length(SNR10)); bit_error_rate = zeros(length(SNR10),1); for i=1:length(SNR10) rx(:,i) = awgn(tx,SNR10(i),'measured'); end rx_demod = demodulate(hDemod,rx); for i=1:length(SNR10) lu an [~,bit_error_rate10(i)] = biterr(x,rx_demod(:,i)); n va end SNR11 = EbNo11 +11.3- 10*log10(k); gh tn to EbNo11 = 0:30; % In dB p ie rx = zeros(nSyms,length(SNR11)); bit_error_rate11 = zeros(length(SNR11),1); nl w for i=1:length(SNR11) an lu end d oa rx(:,i) = awgn(tx,SNR11(i),'measured'); u nf for i=1:length(SNR11) va rx_demod = demodulate(hDemod,rx); ll [~,bit_error_rate11(i)] = biterr(x,rx_demod(:,i)); oi m end z at nh EbNo12 = 0:29; % In dB SNR12 = EbNo12 +11.7- 10*log10(k); z an Lu end m co rx(:,i) = awgn(tx,SNR12(i),'measured'); l for i=1:length(SNR12) gm bit_error_rate12 = zeros(length(SNR12),1); @ rx = zeros(nSyms,length(SNR12)); n va ac th si 65 rx_demod = demodulate(hDemod,rx); for i=1:length(SNR12) [~,bit_error_rate12(i)] = biterr(x,rx_demod(:,i)); end figure; semilogy(EbNo10, bit_error_rate10,'^-',EbNo11, bit_error_rate11, 'o-',EbNo12, bit_error_rate12, '*-'); grid on; legend('SISO','MISO','MIMO'); lu an xlabel('Eb/No (dB)'); n va ylabel('Bit Error Rate'); Eb/No'); gh tn to title('BER TREN KENH MIMO, MISO VA SISO NHU NHAU TREN CUNG M5 = 32; nl w clc; p ie clear all va nSyms = n/k; an lu n = 2e5; d oa k = log2(M5); oi m hMod.SymbolOrder = 'Gray'; ll hMod.InputType = 'Bit'; u nf hMod = modem.qammod(M5); z at nh hDemod = modem.qamdemod(hMod); x = randi([0 1],n,1); z an Lu bit_error_rate = zeros(length(SNR13),1); m co rx = zeros(nSyms,length(SNR13)); l SNR13 = EbNo13 +6.3- 10*log10(k); gm EbNo13 = 0:30; @ tx = modulate(hMod,x); n va ac th si 66 for i=1:length(SNR13) rx(:,i) = awgn(tx,SNR13(i),'measured'); end rx_demod = demodulate(hDemod,rx); for i=1:length(SNR13) [~,bit_error_rate13(i)] = biterr(x,rx_demod(:,i)); end EbNo14 = 0:30; % In dB SNR14 = EbNo14 +11.3- 10*log10(k); lu an rx = zeros(nSyms,length(SNR14)); n va bit_error_rate14 = zeros(length(SNR14),1); rx(:,i) = awgn(tx,SNR14(i),'measured'); gh tn to for i=1:length(SNR14) p ie end rx_demod = demodulate(hDemod,rx); nl w for i=1:length(SNR14) an lu end d oa [~,bit_error_rate14(i)] = biterr(x,rx_demod(:,i)); va EbNo15 = 0:29; % In dB u nf SNR15 = EbNo15 +12- 10*log10(k); ll rx = zeros(nSyms,length(SNR15)); m oi bit_error_rate15 = zeros(length(SNR15),1); z at nh for i=1:length(SNR15) rx(:,i) = awgn(tx,SNR15(i),'measured'); z rx_demod = demodulate(hDemod,rx); an Lu end m co [~,bit_error_rate15(i)] = biterr(x,rx_demod(:,i)); l for i=1:length(SNR15) gm @ end n va ac th si 67 figure; semilogy(EbNo13, bit_error_rate13,'^-',EbNo14, bit_error_rate14, 'o-',EbNo15, bit_error_rate15, '*-'); grid on; legend('SISO','MISO','MIMO'); xlabel('Eb/No (dB)'); ylabel('Bit Error Rate'); title('BER TREN KENH MISO TIEM CAN KENH MIMO'); clear all lu an Dlen=1024; n va Flen=64; F=Dlen/Flen; gh tn to IFFTlen=2*Flen; p ie G=16; x1=zeros(Dlen/2,1); nl w y=zeros(IFFTlen,1); d oa y1=zeros(IFFTlen,1); an lu z=zeros(F*IFFTlen,1); oi m for j=1:Flen/2 ll for i=1:F u nf x=randi(Dlen,1,2); va ofdm_sg=zeros(F*IFFTlen+G,1); z at nh y(j+IFFTlen/2-Flen/4)=x1((i-1)*Flen/2+j); end z an Lu end m co end l z((i-1)*IFFTlen+j)=y1(j); gm for j=1:IFFTlen @ y1=ifft(y,length(y)); n va ac th si 68 for i=1:G ofdm_sg(i)=z(F*IFFTlen-G+i); end for i=1:F*IFFTlen ofdm_sg(G+i)=z(i); end lu an n va p ie gh tn to d oa nl w ll u nf va an lu oi m z at nh z m co l gm @ an Lu n va ac th si