Nâng cao chất lượng mạng tính toán di động thế hệ sau

ĐẠI HỌC QUỐC GIA H À NỘI TRƯỜNG ĐẠI HỌC CÔNG NGHỆ LƯƠNG DUY HIẾU NÂNG CAO CHẤT LƯỢNG MẠNG TÍNH TỐN DI ĐỘNG THẾ HỆ SAU LUẬN VĂN THẠC SĨ HÀ NỘI - 2010 ĐẠI HỌC QUỐC GIA H À NỘI TRƯỜNG ĐẠI HỌC CÔNG NGHỆ LƯƠNG DUY HIẾU NÂNG CAO CHẤT LƯỢNG MẠNG TÍNH TỐN DI ĐỘNG THẾ HỆ SAU NGÀNH CÔNG NGHỆ THÔNG TIN CHUYÊN NGÀNH TRUYỀN DỮ LIỆU VÀ MẠNG MÁY TÍNH MÃ SỐ 60 48 15 LUẬN VĂN THẠC SĨ GIÁO VIÊN HƯỚNG DẪN: PGS.TS TRẦN HỒNG QUÂN HÀ NỘI - 2010 LỜI CAM ĐOAN Tơi xin cam đoan luận văn hồn thành thân tự nghiên cứu không chép từ luận văn Nội dung luận văn có trích dẫn, tham khảo nêu rõ phần “Tài liệu tham khảo” Tôi xin cam đoan lời khai đúng, thông tin sai lệch xin hoàn toàn chịu trách nhiệm trước Hội đồng Hà Nội, ngày 18 tháng năm 2010 Người cam đoan Lương Duy Hiếu MỤC LỤC LỜI CẢM ƠN BẢNG CÁC TỪ VIẾT TẮT BẢNG THỐNG KÊ CÁC HÌNH VẼ CHƯƠNG LỊCH SỬ HỆ THỐNG THÔNG TIN DI ĐỘNG V À MẠNG 4G .3 1.1 LỊCH SỬ MẠNG DI ĐỘNG V À XU HƯỚNG PHÁT TRIỂN 1.2 GIỚI THIỆU MẠNG 4G 1.3 CÁC ỨNG DỤNG VÀ DỊCH VỤ TRONG MẠNG 4G CHƯƠNG CÁC ĐẶC TÍNH CỦA KỸ THUẬT OFDM TRONG HỆ THỐNG THÔNG TIN DI ĐỘNG 4G 13 2.1 MỘT SỐ KHÁI NIỆM 13 2.2 LỊCH SỬ PHÁT TRIỂN OFDM 14 2.3 NGUYÊN LÝ CƠ BẢN CỦA OFDM 15 2.4 MƠ HÌNH TỐN HỌC CỦA OFDM 19 2.5 CÁC KỸ THUẬT ĐIỀU CHẾ TRONG OFDM 19 2.6 ƯU, NHƯỢC ĐIỂM CỦA OFDM 22 CHƯƠNG NÂNG CAO HIỆU NĂNG CỦA HỆ THỐNG THÔNG TIN DI ĐỘNG BẰNG GIẢI PHÁP KẾT HỢP OFDM -MIMO 24 3.1 GIỚI THIỆU 24 3.2 MƠ HÌNH MÁY PHÁT DI D ỘNG OFDM 26 3.3 MƠ HÌNH KÊNH 26 3.4 MƠ HÌNH MÁY THU DI D ỘNG OFDM 28 3.5 CÁC THUẬT TOÁN ĐIỀU KHIỂN ANTEN 29 3.6 KẾT QUẢ MÔ PHỎNG 31 3.7 MƠ HÌNH HỆ THỐNG MIMO 33 3.7.1 Khái niệm 33 3.7.2 Mơ hình tốn học 33 3.7.3 Dung lượng hệ thống MIMO 34 3.8 MƠ HÌNH HỆ THỐNG OFDM - MIMO KẾT HỢP 35 3.8.1 Đề xuất mơ hình hệ thống MIMO-OFDM 36 3.9 TÍNH HIỆU NĂNG CỦA HỆ THỐNG OFDM -MIMO 38 3.10 XÁC SUẤT LỖI ĐỐI VỚI HỆ THỐNG VÔ TUYẾN ĐIỀU CHẾ MPSK 40 3.11 XÁC SUẤT LỖI ĐỐI VỚI HỆ THỐNG ĐIỀU CHẾ MQAM 41 3.12 TÍNH XÁC SUẤT LỖI ĐỐI V ỚI HỆ THỐNG ĐIỀU CHẾ QPSK 42 CHƯƠNG MƠ PHỎNG, PHÂN TÍCH ĐÁNH GIÁ HIỆU NĂNG GIẢI PHÁP KẾT HỢP OFDM-MIMO 44 4.1 CÔNG CỤ MÔ PHỎNG 44 4.2 MÔ PHỎNG BỘ PHÁT VÀ BỘ THU OFDM 44 4.3 MÔ PHỎNG GIẢI PHÁP KẾT HỢP OFDM-MIMO 47 4.3.1 Mơ hình giải pháp OFDM-MIMO 47 4.3.2 Lưu đồ phát ký tự mơ hình OFDM-MIMO 47 4.3.3 Lưu đồ thu ký tự mơ hình OFDM-MIMO 48 4.3.4 Lưu đồ mô thuật tốn tính BER 49 4.3.5 Số liệu mô 49 4.3.6 Xác suất lỗi hệ thống vô tuyến điều chế MPSK 50 4.3.7 Xác suất lỗi hệ thống vô tuyến điều chế MQAM v MPAM 51 4.3.8 Xác suất lỗi hệ thống vô tuyến điều chế QPSK 53 4.4 ĐÁNH GIÁ 54 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ĐỀ T ÀI 56 TÀI LIỆU THAM KHẢO 56 PHẦN PHỤ LỤC 59 LỜI CẢM ƠN Tôi xin chân thành cảm ơn Thầy giáo PGS.TS Trần Hồng Quân, Viện khoa học kỹ thuật Bưu điện tận tình hướng dẫn tơi suốt thời gian nghiên cứu luận văn Xin cảm ơn Thầy giáo, Cô giáo Khoa Công nghệ thông tin - Đại học Công nghệ - Đại học Quốc gia Hà Nội tận tình giảng dạy, trang bị vốn kiến thức kinh nghiệm quý báu để tơi có kết tốt học tập Tơi xin bày tỏ lòng biết ơn sâu sắc đến gia đình, bạn bè giúp đỡ, động viên tơi hồn thành tốt chương trình học đề tài nghiên cứu Hà Nội, ngày 18 tháng năm 2010 BẢNG CÁC TỪ VIẾT TẮT A AM Amplitude Modulation AWGN Additive White Gaussian Noise B BER Bit Error Rate BPSK Binary Phase Shift Keying BS Base Station C CDMA Code Division Multiple Access CP Cyclic Prefix D DC Direct Current (0 Hz) DFT Discrete Fourier Transform DPLL Digital Phase Look Loop DS-CDMA Direct Sequence CDMA DSP Digital Signal Processor DVB Digital Video Broadcasting F FDM Frequency Division Multiplexing FEC Forward Error Correcting FFT Fast Fourier Transform FIR Finite Impulse Response (digital filter) FM Frequency Modulation FOE Frequency Offset Estimation FSC Frame Synchronization Code FSK Frequency Shift Keying G GI Guard Interval I ICI InterChannel Interference ICI InterCarrier Interference ISI InterSymbol Interference IDFT Inverse Discrete Fourier Transform IEEE Institute of Electrical and Electronic Engneers IFFT Inverse FFT IMD Inter-Modulation Distortion ISI InterSymbol Interference O OFDM Orthogonal Frequency Division Multiplexing P P/S Parallel to Serial PAPR Peak to Average Power Ratio PM Phase Modulation PN Pseudo Noise PSK Phase-Shift Keying Q QAM Quadrature Amplitude Modulation QPSK Quadrature Phase-Shift Keying S S/P Serial to Parallel SC Single Carrier SNR Signal to Noise Ratio W Wimax Worldwide Interoperability for Microwave Access BẢNG THỐNG KÊ CÁC HÌNH VẼ HÌNH SỐ NỘI DUNG HÌNH VẼ Hình 1.1 Trục phát triển hệ thống thơng tin di động Hình 1.2 Các hệ di động Hình 1.3 Bảng so sánh tham số cơng nghệ c Hình 1.4 Lịch sử hệ di động theo cơng nghệ phát triển Hình 1.5 Mơ hình dịch vụ mạng 4G Hình 2.1 Truyền dẫn sóng mang đơn Hình 2.2 Cấu trúc hệ thống truyền dẫn đa sóng mang Hình 2.3 Các sóng mang trực giao Hình 2.4 Sơ đồ hệ thống OFDM Hình 2.5 OFDM với sóng mang Hình 3.1 Mơ hình máy phát thu hệ thống OFDM Hình 3.2 Mơ hình véc tơ kênh khơng gian cho ULA kênh Hình 3.3 AWGN kênh Rayleigh Hình 3.4 Phân tích giả lập kết BER cho kênh AWGN kênh Rayleigh Hình 3.5 Mơ hình kết hợp MIMO OFDM Hình 3.6 N kênh truyền Gauss trắng song song HÌNH SỐ NỘI DUNG HÌNH VẼ Hình 3.7 Các chuẩn thơng tin khơng dây IEEE Hình 3.8 Mơ hình hệ thống MIMO Hình 3.9 Mơ hình hệ thống OFDM-MIMO với phát thu cải tiến Hình 3.10 Mơ hình thực biến đổi Hình 4.1 Mơ hình hệ thống OFDM-MIMO Hình 4.2 Lưu đồ mơ phát ký tự mơ h ình OFDM-MIMO Hình 4.3 Phổ tín hiệu OFDM truyền Hình 4.4 Phổ tín hiệu OFDM nhận Hình 4.5 Dạng sóng tín hiệu OFDM truyền Hình 4.6 Dạng sóng tín hiệu OFDM nhận Hình 4.7 Mơ hình hệ thống OFDM-MIMO Hình 4.8 Lưu đồ mô phát ký tự mô h ình OFDM-MIMO Hình 4.9 Lưu đồ mô thu ký tự mơ hình OFDM-MIMO Hình 4.10 Lưu đồ mơ thuật tốn BER Hình 4.11 Đồ thị mơ tả xác suất lỗi điều chế MPSK Hình 4.12 Đồ thị mô tả xác suất lỗi điều chế MQAM Hình 4.13 Đồ thị mơ tả xác suất lỗi điều chế MPAM Hình 4.14 Đồ thị mô tả xác suất lỗi điều chế QPSK 58 [11] L.HANZO,W.WEBB,and T.KELLER,"Single-and Multi-Carrier Quadrature Amplititude Modulation" New York: IEEE Press/ Wiley, Apr.2000 [12] Richard van Nee, Ramjee Prasad, " OFDM for wireless multimedia communications", Artech House, 2000 [13] Ahmad R.S Bahai, Burton R Saltzberg, “Multicarier Digital Communications Theory and Applications of OFDM”, Kluwer Academic Publishers, 2002 [14] C.Williams, M.A.Beach, and S.McLaughlin, BRobustOFDM timing synchronisation [15] C.-N.Chuah, D.N.C.Tse, J.M.Kahn, and R.A.Valenzuela, Bcapacity scaling in MIMO wireless systems under correlated fading, [ IEEETrans.Inf.Theory] 59 PHẦN PHỤ LỤC MÔ PHỎNG BỘ PHÁT VÀ THU OFDM % -% Simulation program OFDM transmission/receive system % File name: OFDM.m % Development: Lương Duy Hiếu % -% Preparation part para=128; % So kenh truyen fftlen=128; %FFF length noc=128; %So song mang nd=6; %So ky hieu OFDM ml=2; %Muc dieu che: QPSK sr=250000; %Ty le br=sr.*ml; %Ty le bit tren song mang gilen=32; %Chieu dai khoang bao ve ebn0=3; %Eb/N0 % Main part nloop=100; %So vong gia lap noe = 0; % So ty le loi nod = 0; %So du lieu truyen eop = 0; % So goi loi nop = 0; % So goi truyen for iii=1:nloop %Transmitter %Data generation seldata=rand(1,para*nd*ml) >0.5; %Convert to parallel paradata=reshape(seldata,para,nd*ml); % QPSK modulelation [ich,qch]=qpskmod(paradata,para,nd,ml); kmod=1/sqrt(2); ich1=ich.*kmod; qch1=qch.*kmod; %IFFT x=ich1+qch1.*i; y=ifft(x); ich2=real(y); 60 qch2=imag(y); %Guard [ich3,qch3]=giins(ich2,qch2,fftlen,gilen,nd); fftlen2=fftlen+gilen; %Attenuation spow=sum(ich3.^2+qch3.^2)/nd./para; attn=0.5*spow*sr/br*10.^( -ebn0/10); attn=sqrt(attn); %Receiver %AWGN add [ich4,qch4]=comb(ich3,qch3,attn); %Guard [ich5,qch5]=girem(ich4,qch4,fftlen2,gilen,nd); %FFT rx=ich5+qch5.*i; ry=fft(rx); ich6=real(ry); qch6=imag(ry); %Demodule ich7=ich6./kmod; qch7=qch6./kmod; [demodata]=qpskdemod(ich7,qch7,para,nd,ml); %Parallel demodata1=reshape(demodata,1,para*nd*ml); %Bit Error noe2=sum(abs(demodata1 -seldata)); nod2=length(seldata); noe=noe+noe2; nod=nod+nod2; %Calculating PER if noe2~=0 eop=eop+1; else eop=eop end eop; nop=nop+1; fprintf('%d\t%e\t%d\n',iii,noe2/nod2,eop); end %Output per=eop/nop; 61 ber=noe/nod; fprintf('%f\t%e\t%e\t%d\t\n',ebn0,ber,per,nloop); fid=fopen('BERofdm.dat','a'); fprintf(fid,'%f\t%e\t%e\t%d\t\n',ebn0,ber,per,nloop); fclose(fid); % -% Function qpskmod % File name: qpskmod.m % Development: Lương Duy Hiếu % -function [iout,qout]=qpskmod(paradata,para,nd,m1) %paradata: input data (para -by-nd matrix) %iout: output Ich data %qout: output Qch data %para: Number of parallel channels %nd: Number of data %ml: Number of modulation levels %(QPSK-2 16QAM-4) m2=m1./2; paradata2=paradata.*2 -1; count2=0; for jj=1:nd isi=zeros(para,1); isq=zeros(para,1); for ii=1:m2 isi=isi+2.^(m2-ii) *paradata2((1:para),ii+count2); isq=isq+2.^(m2-ii) *paradata2((1:para),m2+ii+count2); end iout((1:para),jj)=isi; qout((1:para),jj)=isq; count2=count2+m1; end % -% Function comb % File name: comb.m % Development: Lương Duy Hiếu % -function [iout,qout] = co mb (idata,qdata,attn) 62 %****************** variables ************************* % idata : input Ich data % qdata : input Qch data % iout output Ich data % qout output Qch data % attn : attenuation level caused by Eb/No or C/N %****************** ************************************ iout = randn(1,length(idata)).*attn; qout = randn(1,length(qdata)).*attn; iout = iout+idata(1:length(idata)); qout = qout+qdata(1:length(qdata)); % -% Function qpskdemod % File name: qpskdemod.m % Development: Lương Duy Hiếu % -function [demodata]=qpskdemod(idata,qdata,para,nd,ml) %****************** variables ************************* % idata :input Ich d ata % qdata :input Qch data % demodata: demodulated data (para -by-nd matrix) % para : Number of paralell channels % nd : Number of data % ml : Number of modulation levels % (QPSK ->2 16QAM -> 4) % ************************************************* **** demodata=zeros(para,ml*nd); demodata((1:para),(1:ml:ml*nd 1))=idata((1:para),(1:nd))>=0; demodata((1:para),(2:ml:ml*nd))=qdata((1:para),(1:nd))>=0; % -% Function giins % File name: giins.m % Development: LD Hiếu % -function [iout,qout]=giins(idata,qdata,fftlen,gilen,nd); idata1=reshape(idata,fftlen,nd); qdata1=reshape(qdata,fftlen,nd); idata2=[idata1(fftlen -gilen+1:fftlen,:);idata1]; 63 qdata2=[qdata1(fftlen -gilen+1:fftlen,:);idata1]; iout=reshape(idata2,1,(fftlen+gilen)*nd); qout=reshape(qdata2,1,(fftlen+gilen)*nd); % -% Function girem % File name: girem.m % Development: LD Hi ếu % - %girem.m function [iout,qout]=girem(idata,qdata,fftlen2,gilen,nd); idata2=reshape(idata,fftlen2,nd); qdata2=reshape(qdata,fftlen2,nd); iout=idata2(gilen+1:fftlen2,:); qout=qdata2(gilen+1:fftlen2,:); MÔ PHỎNG GIẢI PHÁP OFDM-MIMO % -% Channel estimation for a MIMO -OFDM % Development: Lương Duy Hiếu % -clear all; %close all; Rayleigh=1; AWGN=0; % for AWGN channel MMSE=0; % estimation technique Nsc=64; % Number of subcarriers Ng=16; % Cyclic prefix length SNR_dB=[0 10 15 20 25 30 35 40]; % Signal to noise ratio Mt=2; % Number of Tx antennas Mr=2; % Number of Rx antennas pilots=[1:Nsc/Ng:Nsc]; % pilot subcarriers DS=15; % Delay spread of channel iteration_max=200; %%%%%%%%%%%%%%%%%%%%% %%%%%%% % Channel impulse response % %%%%%%%%%%%%%%%%%%%%%%%%%%%% if (Rayleigh) N=50; fm=100; %??? B=20e3; fd=(rand(1,N)-0.5)*2*fm; %?? theta=randn(1,N)*2*pi; 64 c=randn(1,N); c=c/sum(c.^2); t=0:fm/B:10000*fm/B; Tc=zeros(size(t)); Ts=zeros(size(t)); for k=1:N Tc=c(k)*cos(2*pi*fd(k)*t+theta(k))+Tc; Ts=c(k)*sin(2*pi*fd(k)*t+theta(k))+Ts; end r=ones(Mt*Mr,1)*(Tc.^2+Ts.^2).^0.5; index=floor(rand(Mt*Mr,DS)*5000+1); end MEE1=zeros(1,length(SNR_dB )); MEE2=zeros(1,length(SNR_dB)); for snrl=1:length(SNR_dB) % snrl estimation_error1=zeros(Mt*Mr,Nsc); estimation_error2=zeros(Mt*Mr,Nsc); R1=besselj(0,2*pi*fm*(Nsc+Ng)/B); sigma2=10^(-SNR_dB(snrl)/10); aa=(1-R1^2)/(1-R1^2+sigma2); bb=sigma2*R1/(1-R1^2+sigma2); for iteration=1:iteration_max %iteration if AWGN==1 h=ones(Mt*Mr,1); else phi=rand*2*pi; h=r(index+iteration)*exp(i*phi); %h=rand(Mt*Mr,DS); h=h.*(ones(Mt*Mr,1)*(exp( -0.5).^[1:DS])); h=h./(sqrt(sum(abs(h).^2,2))*ones(1,DS)); end CL=size(h,2); % %channel length data_time=zeros(Mt,Nsc+Ng); data_qam=zeros(Mt,Nsc); data_out=zeros(Mr,Nsc); output=zeros(Mr,Nsc); for tx=1:Mt data_b=0*round(rand(4,Nsc)); % data data_qam(tx,:)=i*(2*(mod(data_b(1,:)+data_b(2,:),2)+2*dat a_b(1,:))-3)+ 65 2*(mod(data_b(3,:)+data_b(4,:),2)+2*data_b(3,:)) -3; for loop=1:Mt data_qam(tx,pilots+loop -1)=(1+i)*(loop==tx); % pilots end data_time_temp=ifft(data_qam(tx,:)); data_time(tx,:)=[data_ time_temp(end-Ng+1:end) data_time_temp]; end for rx=1:Mr for tx=1:Mt output_temp=conv(data_time(tx,:),h((rx -1)*Mt+tx,:)); output(rx,:)=output_temp(Ng+1:Ng+Nsc)+output(rx,:); end np=(sum(abs(output(rx,:)).^2)/length(output(rx,:)))*sigma2; noise=(randn(size(output(rx,:)))+i*randn(size(output(rx,:))))* sqrt(np); output(rx,:)=output(rx,:)+noise; data_out(rx,:)=fft(output(rx,:)); end %%%%%%%%%%%%%%%%%%%%%% % Channel estimation % %%%%%%%%%%%%%%%%%%%%%% H_act=zeros(Mt*Mr,Nsc); H_est1=zeros(Mt*Mr,Nsc); H_est2=zeros(Mt*Mr,Nsc); i=1; for tx=1:Mt for rx=1:Mr H_est_temp=data_out(rx,pilots+tx -1)./data_qam(tx,pilots+tx -1); %H_est_temp2=aa*abs(H_est_temp1)+bb*abs(H_est2((rx 1)*Mt+tx,:)) ; h_time=ifft(H_est_temp); h_time=[h_time zeros(1,Nsc -length(h_time))]; H_est1((rx-1)*Mt+tx,:)=fft(h_time); H_est2((rx1)*Mt+tx,:)=((aa*abs(H_est1((rx 1)*Mt+tx,:))+bb*abs(H _est2((rx1)*Mt+tx,:))) *H_est1((rx1)*Mt+tx,:))./abs(H_est1 ( rx1)*Mt+tx,:)); if (tx>1) H_est1((rx-1)*Mt+tx,:)=[H_est1((rx-1)*Mt+tx,Nsctx+2:Nsc) H_est1((rx-1)*Mt+tx,1:Nsc-tx+1)]; H_est2((rx-1)*Mt+tx,:)=[H_est2((rx -1)*Mt+tx,Nsctx+2:Nsc) H_est2((rx-1)*Mt+tx,1:Nsc-tx+1)]; 66 end H_act((rx-1)*Mt+tx,:)=fft([h((rx-1)*Mt+tx,:) zeros(1,Nsc CL)]); error1=(abs(H_act((rx -1)*Mt+tx,:)-H_est1((rx-1)*Mt+tx,:)).^2); error2=(abs(H_act((rx -1)*Mt+tx,:)-H_est2((rx-1)*Mt+tx,:)).^2); %error=(abs(H_act((rx1)*Mt+tx,:)H_est((rx1)*Mt+tx,:)).^2)./(ab s(H_act((rx1)*Mt+tx,:)).^2); estimation_error1((rx 1)*Mt+tx,:)=estimation_error1((rx1)*Mt+tx ,:)+error1; estimation_error2((rx 1)*Mt+tx,:)=estimation_error2((rx 1)*Mt+tx ,:)+error2; %subplot(Mt*Mr,3,i),plot([0:Nsc-1],abs(H_act((rx1)*Mt+tx,:))); i=i+1; %subplot(Mt*Mr,3,i),plot([0:Nsc-1],abs(H_est((rx1)*Mt+tx,:))); i=i+1; %subplot(Mt*Mr,3,i),plot([0:Nsc -1],abs(error)); i=i+1; end end end estimation_error1=estimation_error1/iteration_max; estimation_error2=estimation_error2/iteration_max; %estimation_error=min(estimation_error,10*iteration_max*ones(s ize(estimation_error))); %for i=1:Mt*Mr % subplot(Mt*Mr,2,2*i-1),plot([0:Nsc1],estimation_error1(i,:)); % subplot(Mt*Mr,2,2*i),plot([0:Nsc 1],estimation_error2(i,:)); %end MEE1(snrl)=sum(sum(estimation_error1))/(Mt*Mr*Nsc); MEE2(snrl)=sum(sum(estimation_error2))/(Mt*Mr*Nsc); end plot(SNR_dB,10*log10(MEE1)); hold on; plot(SNR_dB,10*log10(MEE2),'r'); H_act=fft([h zeros(1,Nsc -CL)]).'; error1=(abs(H_act-H_est1).^2)./(abs(H_act).^2); error2=(abs(H_act-H_est2).^2)./(abs(H_act).^2); %%%%%%%%% % Plots % %%%%%%%%% %fig=4; %i=1; %subplot(fig,1,i),plot([0:length(H_act) -1],abs(H_act)); i=i+1; 67 %subplot(fig,1,i),plot([0:length(H_est1) -1],abs(H_est1)); i=i+1; %subplot(fig,1,i),plot([0:length(H_est2) -1],abs(H_est2)); i=i+1; %subplot(fig,1,i),plot([0:length(error1) -1],error1); i=i+1; %subplot(fig,1,i),plot([0:length(error2) -1],error2); % -% Xác suất lỗi với điều chế BPSK % Development: Lương Duy Hiếu % -%Dinhnghiathamso EbNodB=0:10; EbNo=10.^(EbNodB./10); sigLen=5*10^5; %TaotinhieuBPSK{+1,-1} s=1-2*(rand(1,sigLen)>=0.5); %TinhtoannangluongbittinhieuEb Eb=norm(s)^2/sigLen; %MatdophoAWGN No=Eb./EbNo; %VonglaptinhtoanBERtheoEb/No fork=1:length(EbNo) %TaoAWGN n=sqrt(No(k)./2).*(randn(1,sigLen)+j*randn(1,sigLen)); %Tinhieuthu y=s+n; %Tachtinhieu shat=sign(real(y)); error=s-shat; noError=length(find(error~=0)); BER(k)=noError/sigLen; End %BERlythuyetcuatruyendanBPSKquakenhAWGN BERtheory=1/2*erfc(sqrt(EbNo)); %Vedothi semilogy(EbNodB,BER,’*’,EbNodB,BERtheory) xlabel(’Eb/No’) ylabel(’BER’) legend(’Bysimulation’,’Bytheory’) title(’BERcuaBPSKquakenhAW GN’) grid % -% Xác suất lỗi với điều chế MQAM 68 % Development: Lương Duy Hiếu % -clear % Xet truc x chay tu - 10 den 40 x = 0:20; t = 10.^(x/10); %u = -1/t; y = 1/2.*(1 - sqrt(t./(1+t))); y2 = 1/2.* erfc(sqrt(2*t)/sqrt(2)); % Lay M =4 y3 = erfc(sqrt(t)*sin(pi/4)); figure semilogy(x,y,'b+-','LineWidth',2); hold on; semilogy(x,y2,'g+-','LineWidth',1); hold on; axis([-0 20 10^-8 1]) grid off xlabel('gamma (dB)'); ylabel('PB'); legend('Rayleigh','AWGN','AWGN.MQAM'); title('Pb cho MQAM'); % -% Xác suất lỗi với điều chế MPAM % Development: Lương Duy Hiếu % -clear % Xet truc x chay tu - 10 den 40 x = 0:20; t = 10.^(x/10); %u = -1/t; y = 1/2.*(1 - sqrt(t./(1+t))); y4 = 1/2 * erfc(sqrt(t/2)); % Dieu che MPAM - voi M =4 figure semilogy(x,y,'b+-','LineWidth',2); hold on; semilogy(x,y4,'g+-','LineWidth',1); hold on; axis([-0 20 10^-8 1]) grid off xlabel('gamma (dB)'); ylabel('PB'); legend('Rayleigh','AWGN','AWGN.MPAM'); 69 title('Pb cho MPAM'); % -% Xác suất lỗi với điều chế QPSK % Development: LD Hi ếu % -clear % Xet truc x chay tu - 10 den 40 x = 0:20; t = 10.^(x/10); %u = -1/t; y = 1/2.*(1 - sqrt(t./(1+t))); y5 = 3/4 * erfc( sqrt(3/15 * t)); % Dieu che BPSK - voi M =4 y6 = erfc( sqrt(1/4 * t)); % Dieu che BFSK Xac suat loi bit y7 = 1/2 * erfc( sqrt(1/2 * t)); % Dieu che BPSK Xac suat loi bit y8 = 1/2 * erfc( sqrt(t)); figure semilogy(x,y,'b+-','LineWidth',2); hold on; semilogy(x,y5,'g+-','LineWidth',1); hold on; axis([-0 20 10^-8 1]) grid off xlabel('gamma (dB)'); ylabel('PB'); legend('Rayleigh','AWGN','QPSK'); title('Pb cho QPSK'); % -% Hàm vẻ đổ thị kênh Rayleigh % Development: LD Hi ếu % -%function plot_Rayleigh(handles) %axes(handles.axes1); delta_ = 1; r = -3:0.01:3; a = -3:0.01:3; [R,A] = meshgrid(r,a); %P = ((R + A)/delta_^2).*exp( -(R.^2 + A.^2)/(2*delta_^2)); P = 1/(sqrt(2*pi)*delta_)*exp(( -R.^2 A.^2)/(2*delta_^2)); mesh(P); 70 title('Hàm phân bố xác suất Gaussian','FontName','.Vn time','FontSize',15,'Color','b'; xlabel('x1','FontSize',14,'Color','b'); ylabel('x2','FontSize',14,'Color','b'); zlabel('P(x1,x2)','FontSize',14,'Color','b'); % -% Hàm vẻ đổ thị trải phổ kênh OFDM thông thường % Development: LD Hi ếu % -%function plot_ofdm_spectrum(handles) global dai_fft so_song_mang do_dai bao_ve %global use_window %axes(handles.axes_hieu_nang); %use_window = 1; dai_fft = 256; so_song_mang = 52; bao_ve = 10; x = randint(1,10000); [x,z] = code_16qam(x); %[sig_out, so_khoi, so_0, giua_0] = phat_cua_so(x); % [sig_out, so_khoi, so_0, giua_0] = phat(x); % N1 = length(sig_out); % m1 =linspace(0,0.5,N1); %w1 = window(@hamming,N); %w1 = window(@blackmanharris,N); %w1 = window(@RECTWIN,N); %sig_out = sig_out.*w1'; % t = 10*log10(abs(fft(sig_out))); % so_song_mang = 1536; % dai_fft = 4096; [sig_out, so_khoi, so_0, giua_0] = phat_cua_so(x); u = 10*log10(abs(fft(sig_out))); %t1 = floor(length(t)/2); N2 = length(sig_out); m2 = linspace(0, 0.5, N2); subplot(1,2,1); %plot(t(1:t1)); %hold on; %subplot(2,1,1); %axes(handles.axes1); %plot(m1,t); % title('Phổ tín hiệu OFDM','FontName','.Vntime','FontSize',14,'Color','b'); 71 % ylabel('Biên độ phổ (dB)','FontName','.Vntime','FontSize',14,' Color','b'); % xlabel('TÇn sè (Hz)','FontName','.Vntime','FontSize',14,'Color','b'); % grid on; % xlim([0 0.5]); % subplot(1,2,2); %xes(handles.axes2); plot(m2,u); title('Phỉ tÝn hiƯu OFDM','FontName','.Vntime','FontSize',14,'Color','b'); xlabel('TÇn số (Hz)','FontName','.Vntime','FontSize',14,'Color','b'); ylabel('Biên độ phổ (dB)','FontName','.Vntime','FontSize',14,'Color','b'); grid on; xlim([0 0.5]); %hold on; %plot(w1,'r'); %ylim([0 5]); %subplot(2,1,1); %plot(sig_out); %[sig_out] = thu(sig_out,so_khoi , so_0, giua_0); % -% Hàm tính BER điều chế 16-QMA % Development: Lương Duy Hiếu % -function so_sanh_ber_16(handles) global so_diem_SNR so_diem_SNR = 80; SNR_min = -50; SNR_max = 50; load ber_sm_16.mat; load ber_cm_16.mat; load ber_lm_16.mat; load ber_no_16.mat; %axes(handles.axes_so_sanh); points = linspace(SNR_min, SNR_max, so_diem_SNR -1); u = 1:so_diem_SNR-1; semilogy(points,ber_no_16(u),' gs'); 72 hold on; semilogy(points,ber_lm_16(u),' bs'); semilogy(points,ber_sm_16(u),' cs'); semilogy(points,ber_cm_16(u),' rs'); grid on; hold off; title('Hiệu BER AOFDM vs 16 QAM','FontName','.Vntime','FontSize',14,'Color','B'); legend('Khong thich ung', 'Thich ung muc dieu che','Thich ung chon loc song mang','Ket hop co che thich ung'); xlabel('SNR (dB)','FontSize',12,'Color','b'); ylabel('BER','FontSize',12,'Color','b'); %ylim([0 0.27]); xlim([5 40]); ... À NỘI TRƯỜNG ĐẠI HỌC CÔNG NGHỆ LƯƠNG DUY HIẾU NÂNG CAO CHẤT LƯỢNG MẠNG TÍNH TỐN DI ĐỘNG THẾ HỆ SAU NGÀNH CƠNG NGHỆ THÔNG TIN CHUYÊN NGÀNH TRUYỀN DỮ LIỆU VÀ MẠNG MÁY TÍNH MÃ SỐ 60 48 15 LUẬN VĂN... thống thông tin di động, hệ gắn liền với cộng nghệ then chốt l àm tàng cho hệ thống Với hệ di động 2G TMDA CDMA Với hệ di động 3G WCDMA CDMA2000 Với hệ 4G cơng nghệ OFDM xem trái tim hệ thống Trong... h ưởng fading Nội dung luận văn đưa đến toán cụ thể sau: Để nâng cao chất lượng hệ thống thông tin di động hệ sau điều chế OFDM từ sử dụng phép biến đổi kỳ dị để biến đổi k ênh từ fading nhiễu