Đánh giá chất lượng hệ thống thông tin vô tuyến sử dụng phần mềm matlab

Tr-ờng đại học vinh khoa điện tử viễn thông === === đồ án tốt nghiệp đại học Đề tài: đánh giá chất l-ợng hệ thống thông tin vô tuyến sư dơng phÇn mỊm matlab Ngườ ng d n : ThS Phạm mạnh toàn n v nt c n : lê nhân toàn p : 49K - ĐTVT Mó số sinh viên : 0851080311 nghÖ an - 01/2013 BỘ GIÁO DỤC VÀ ĐÀO TẠO CỘNG HỊA XÃ HƠI CHỦ NGHĨA VIỆT NAM TRƢỜNG ĐẠI HỌC VINH Độc lập - Tự - Hạnh phúc NHIỆM VỤ ĐỒ ÁN TỐT NGHIỆP Họ tên sinh viên: Lê Nhân Toàn Mã số sinh viên: 0851080311 Ngành: Điện tử - Viễn thông Khoá: 49 Đầu đề đồ án: Các số liệu liệu ban đầu: Nội dung phần thuyết minh tính tốn: Các vẽ, đồ thị (ghi rõ loại kích thƣớc vẽ): Họ tên giảng viên hướng dẫn: Ngày giao nhiệm vụ đồ án: Ngày hoàn thành đồ án: TRƯỞNG BỘ MƠN ThS Phạm Mạnh Tồn / /20 / /20 Ngày tháng năm 2013 NGƯỜI HƯỚNG DẪN Sinh viên hoàn thành nộp đồ án tốt nghiệp ngày tháng năm 2013 CÁN BỘ PHẢN BIỆN BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC VINH - BẢN NHẬN XÉT ĐỒ ÁN TỐT NGHIỆP Họ tên sinh viên: Ngành: Giảng viên hướng dẫn: Cán phản biện: Lê Nhân Toàn Điện tử - Viễn thơng Mã số sinh viên: 0851080311 Khố: 49 ThS Phạm Mạnh Toàn Nội dung thiết kế tốt nghiệp: Nhận xét cán phản biện: Ngày tháng năm Cán phản biện (Ký, ghi rõ họ tên) MỤC LỤC Trang LỜI NÓI ĐẦU i TÓM TẮT ĐỒ ÁN ii DANH MỤC HÌNH VẼ iii CÁC THUẬT NGỮ VIẾT TẮT v CHƢƠNG I TỔNG QUAN VỀ HỆ THỐNG THÔNG TIN VÔ TUYẾN 1.1 Giới thiệu chương 1.2 Khái niệm 1.3 Lịch sử phát triển hệ thống thông tin vô tuyến 1.4 Khái niệm kênh truyền vô tuyến 1.5 Các đặc tính kênh vơ tuyến 1.6 Truyền dẫn băng tần sở truyền dẫn băng thông 1.7 Mộ số kỹ thuật xử lý tín hiệu .6 1.8 Điều chế 10 1.8.1 Khái niệm điều chế số 10 1.8.2 Tại phải dùng điều chế số 11 1.8.3 Điều chế tín hiệu nhiều mức nhằm nâng cao hiệu phổ 11 1.8.4 Lựa chọn tối ưu tập tín hiệu 13 1.9 Sự phân chia tài nguyên vô tuyến 14 1.10 Phân loại hệ thống thông tin vô tuyến 16 1.11 Kết luận chương 18 CHƢƠNG II TRUYỀN DẪN TRÊN KÊNH VÔ TUYẾN SỐ 19 2.1 Giới thiệu chương 19 2.2 Lý thuyết kênh vô tuyến 19 2.2.1 Truyền dẫn phân tập đa đường 19 2.2.2 Kênh không phụ thuộc thời gian 19 2.2.3 Hiệu ứng Doppler kênh phụ thuộc thời gian 20 2.2.4 Bề rộng độ định thời gian kênh 22 2.2.5 Các mơ hình kênh 22 2.3 Nhiễu thông tin vô tuyến số 22 2.3.1 Tạp âm cộng trắng chuẩn 22 2.3.2 Nhiễu xuyên kênh 23 2.3.3 Nhiễu đồng kênh 23 2.3.4 Nhiễu đa truy nhập 24 2.4 Méo tuyến tính 25 2.4.1 Khái niệm 25 2.4.2 Các biện pháp khắc phục 26 2.5 Méo phi tuyến 28 2.5.1 Khái niệm 28 2.5.2 Các biện pháp khắc phục 29 2.6 Fading 31 2.6.1 Khái niệm 31 2.6.2 Phân loại fading 31 2.6.3 Mơ hình tốn học fading 32 2.6.4 Sự ảnh hưởng chuyển động MS 33 2.6.5 Hậu truyền sóng fading đa đường 35 2.6.6 Các loại kênh fading 37 2.6.7 Các biện pháp khắc phục ảnh hưởng fading 40 2.7 Hiện tượng xuyên nhiễu kí hiệu (ISI: Inter Synbol Interference) 44 2.7.1 Khái niệm 44 2.7.2 Ảnh hưởng 44 2.7.3 Điều kiện truyền ISI 44 2.7.4 Giảm nhiễu ISI sử dụng phương pháp lọc 45 2.8 Kết luận chương .52 CHƢƠNG III MÔ PHỎNG VÀ ĐÁNH GIÁ CHẤT LƢỢNG HỆ THỐNG THÔNG TIN VÔ TUYẾN BẰNG PHẦN MỀM MATLAB 53 3.1 Giới thiệu chương 53 3.2 Sơ lược phần mềm Matlab 53 3.3 Vai trị mơ 54 3.4 Đánh giá chất lượng hệ thống thông tin vô tuyến mô 55 3.5 Kênh truyền AWGN 56 3.5.1 BPSK qua kênh truyền AWGN 57 3.5.2 QPSK qua kênh AWGN 60 3.6 Kênh Fading 65 3.6.1 M-PSK qua kênh fading Rayleigh 66 3.6.2 QAM qua kênh fading Rayleigh 71 3.6.3 Mô so sánh AWGN với Fading Rayleigh kỹ thuật điều chế 74 3.7 Kết luận chương .78 KẾT LUẬN 79 TÀI LIỆU THAM KHẢO 80 PHẦN PHỤ LỤC 81 LỜI NÓI ĐẦU So với lĩnh vực truyền thơng khác, thơng tin vơ tuyến có tăng trưởng nhanh chóng năm gần Kênh truyền vơ tuyến chịu nhiều tác động phức tạp như: AWGN, fading, tượng Doppler, ISI, méo Nên việc phân tích, tính tốn, đánh giá chất lượng hệ thống trở nên khó khăn so với kênh truyền hữu tuyến Đề tài “Đán g c ất lượng h thống thông tin vô tuyến sử dụng phần mềm Matlab” phương tiện biểu đạt mới, sử dụng mô thay hệ thống thực, cho phép người học có nhìn trực quan sâu sắc vấn đề kỹ thuật phức tạp Nội dung đồ án chia thành chương: Chương 1: Tổng quan hệ thống thông tin vô tuyến Chương 2: Truyền dẫn kênh vô tuyến số Chương 3: Mô đánh giá chất lượng hệ thống thông tin vô tuyến sử dụng phần mềm matlab Xin trân trọng cảm ơn ThS Phạm Mạnh Toàn giới thiệu, cung cấp tài liệu, tận tình hướng dẫn nội dung phương pháp, giúp tơi hồn thành đồ án Xin chân trọng cảm ơn thầy cô giáo khoa Điện tử Viễn thông trường Đại học Vinh giúp đỡ suốt thời gian học tập hồn thành chương trình đào tạo Sinh viên Lê Nhân Tồn i TĨM TẮT ĐỒ ÁN Đồ án tìm hiểu, lich sử phát triển, cấu trúc, đặc tính kênh truyền hệ thống thơng tin vơ tuyến Phân tích tác động kênh truyền vơ tuyến chất lượng hệ thống Để đơn giản q trình tính tốn thiết kế hệ thống Các giả định lý tưởng đặt ra: Băng tần truyền dẫn vơ hạn, kênh truyền tuyến tính chiệu tác động tạp âm cộng trắng chuẩn tác động phức tạp khác fading, ISI, hiệu ứng Doppler, méo tín hiệu phân tích xem xét, từ bổ xung thiết bị, phương pháp khắc phục ảnh hưởng tượng Đồ án sử dụng phần mềm Matlab mô đánh giá chất lượng hệ thống trường hợp kênh chiệu tác động AWGN kênh chiệu tác động fading Rayleigh ABSTRACT The project explores the history of development, structure, channel characteristics of wireless communication system Analysis of the impact of radio channels for system quality To simplify the calculation process of system design Ideal assumptions set out: Transmission band is infinite, channel is linear and is affected by parasitic noise plus standard white and other complex effects such as fading, ISI, Doppler effect, distortion is analyzed to consider, from there additional devices, methods to overcome the effects of the phenomenon that The project uses simulation Matlab software and assess the quality of the system in case only channel affected by AWGN and channels affected by Rayleigh fading ii DANH MỤC HÌNH VẼ Hình 1.1 Trang Sơ đồ khối đơn giản hệ thống thông tin vô tuyến Hình 1.2 Sơ đồ khối hệ thống thơng tin vô tuyến số Hình 2.1 Hàm truyền đạt kênh 21 Hình 2.2 Sự lan truyền đường l tới trạm MS 33 Hình 2.3 Đáp ứng xung lọc xung có chiều dài hữu hạn 36 Hình 2.4 Hàm phân bố Rayleigh với   38 Hình 2.5 Hàm phân bố Rice cho giá trị khác K với Ap = 39 Hình 2.6 Phân tập không gian sử dụng anten 42 Hình 2.7 Phân tập khơng gian tần số sử dụng anten 43 Hình 2.8 Chuyển mạch bảo vệ kênh dự phịng 43 Hình 2.9 Mơ hình hệ thống băng gốc với tín hiệu xung đơn vị 46 Hình 2.10 Dạng xung qua lọc lý tưởng 47 Hình 2.11 Hàm truyền lọc tổng cộng 48 Hình 3.1 Sơ đồ mơ truyền dẫn kênh AWGN 57 Hình 3.2 Giao diện bắt đầu vào mô kỹ thuật điều chế 59 Hình 3.3 Giao diện chương trình „chuongtrinh‟ 59 Hình 3.4 Mơ BER điều chế BPSK qua kênh AWGN 60 Hình 3.5 Mô BER điều chế QPSK qua kênh AWGN 62 Hình 3.6 Mơ BER điều chế 4-QAM qua kênh AWGN 64 Hình 3.7 Mơ BER điều chế 16-QAM qua kênh AWGN 65 Hình 3.8 Mơ hình truyền sóng đa đường 66 Hình 3.9 Sơ đồ mơ truyền dẫn MPSK qua kênh fading sử dụng tách tín hiệu đồng (coherent detection) 66 Hình 3.10 Sơ đồ phân bố tín hiệu (signai constellation) tín hiệu 67 Hình 3.11 Mô BER điều chế BPSK qua kênh Fading 69 Hình 3.12 Mơ BER điều chế QPSK qua kênh Fading 70 Hình 3.13 Sơ đồ mơ truyền dẫn MQAM qua kênh fading sử dụng tách tín hiệu đồng (coherent detection) 71 Hình 3.14 Mô BER điều chế 4-QAM qua kênh Fading 73 iii Hình 3.15 Mơ BER điều chế 16-QAM qua kênh Fading 74 Hình 3.17 Mơ BER điều chế BPSK 75 Hình 3.18 Mô BER điều chế QPSK 76 Hình 3.19 Mơ BER điều chế 4-QAM 76 Hình 3.20 Mơ BER điều chế 16-QAM 77 Hình 3.21 Mơ BER điều chế kỹ thuật điều chế qua kênh AWGN 77 Hình 3.22 Mơ BER điều chế kỹ thuật điều chế qua kênh Fading 78 iv sigma=sqrt (Es./ (2*m.*EbNo)); % Iteration for estimating BER for k=1:length (EbNo) % Generate AWGN with variance sigma^2 nk=sigma (k)* (randn (1,N)+j*randn (1,N)); % Rx signal yk=sk.*gk+nk; ykHat=conj (gk).*yk./abs (gk).^2; % Coherent detection % Decision bkHat=qamdemod (ykHat,M); % Error calculation noSymErr=sum ((bkHat-bk)~=0); % Number of symbol errors SER (k)=noSymErr/N; % Symbol Error rate [noBitErr,BER (k)]=biterr (bkHat,bk,m); end % BER Ly thuyet cua truyen dan 4QAM qua kenh Fadinh Rayleigh Lythuyet_BER=1/2* (1-1./sqrt (1+1./EbNo)); % Ve thi close all figure semilogy (EbNodB,BER,'bp-','LineWidth',2) hold on semilogy (EbNodB,Lythuyet_BER,'rp','LineWidth',2); axis ( [0 30 10^-6 1]) xlabel ('Eb/No') ylabel ('BER') legend ('Dieu che-4QAM','Ly thuyet-4QAM') title ('BER dieu che 4QAM qua kenh Fadinh Rayleigh') grid ========================================================== muoisauQAMquaFading.m 89 clear all M=16; m=log2 (M); % Number of bits/symbol N=5*10^5; % So symbol gui di EbNodB=0:2:30; EbNo=10.^ (EbNodB./10); % Generate random integers ranging from 0:M-1 b=randint (1,N, [0 M-1]); % Random intergers from 0:M-1 sk=qammod (b,M); Es=norm (sk).^2/length (sk); % Generate faded signal according to channel type g=1/sqrt (2)* (randn (1,N)+j*randn (1,N)); s=sk.*g; % Faded signal % Noise variance of AWGN sigma=sqrt (Es./ (2*m.*EbNo)); % Iteration for estimating BER for k=1:length (EbNo) % Generate AWGN with variance sigma^2 n=sigma (k)* (randn (1,N)+j*randn (1,N)); % Rx signal yk=s+n; ykHat=conj (g).*yk./abs (g).^2; % Coherent detection % Decision bkHat=qamdemod (ykHat,M); % Error calculation noSymErr=sum ((bkHat-b)~=0); % Number of symbol errors SER (k)=noSymErr/N; % Symbol Error rate [noBitErr,BER (k)]=biterr (bkHat,b,m);% BER w/o Gray encoding end % BER Ly thuyet cua truyen dan 16QAM qua kenh Fadinh Rayleigh Lythuyet_BER=3/8* (1-1./sqrt (1+5./ (2*EbNo))); 90 % Ve thi close all figure semilogy (EbNodB,BER,'bp-','LineWidth',2) hold on semilogy (EbNodB,Lythuyet_BER,'rp','LineWidth',2); axis ( [0 30 10^-6 1]) xlabel ('Eb/No') ylabel ('BER') legend ('Dieu che-16QAM','Ly thuyet-16QAM') title ('BER dieu che 16QAM qua kenh Fadinh Rayleigh') grid ========================================================== BPSK.m clear; EbNodB=0:2:30; EbNo=10.^ (EbNodB./10); N=5*10^5; % Tao tin hieu BPSK {+1,-1} bk=rand (1,N)>0.5; %awgn sk=1-2*bk; %fading s = 2*bk-1; % Tinh toan nang luong bit tin hieu Eb Eb=norm (sk)^2/N; % Mat AWGN No=Eb./EbNo; % Vong lap tinh toan BER theo Eb/No for k=1:length (EbNo) nk=sqrt (No (k)./2).* (randn (1,N)+j*randn (1,N)); 91 % Tao AWGN % Tin hieu thu yk=sk+nk; % Tach tin hieu shat=sign (real (yk)); error=sk-shat; noError=length (find (error~=0)); BER (k)=noError/N; %tao fading g = 1/sqrt (2)* [randn (1,N) + j*randn (1,N)]; % Rayleigh channel % Channel and noise Noise addition y = s.*g + nk; % equalization yHat = y./g; % receiver - hard decision decoding ipHat = real (yHat)>0; % counting the errors nErr (k) = size (find ( [bk- ipHat]),2); ber = nErr/N; end % BER ly thuyet cua truyen dan BPSK qua kenh AWGN Lythuyet_BER=1/2*erfc (sqrt (EbNo)); %BER ly thuyet cua truyen dan BPSK qua kenh Fading Lt_BER = 0.5.* (1-sqrt (EbNo./ (EbNo+1))); % Ve thi close all figure close all semilogy(EbNodB,ber,'bp-',EbNodB,BER,'rp',EbNodB,Lt_BER,'mp',EbNodB,Lythuyet_BER,'kp'); axis ( [0 30 10^-6 1]) 92 legend ('Sim qua fading', 'Sim qua awgn','Theory qua fading','Theory qua awgn'); xlabel ('Es/No, dB') ylabel ('BER') title ('BER dieu che BPSK') grid ========================================================== QPSK.m clear all; close all; N=5*10^5; EbNodB=0:2:20; EbNo=10.^ (EbNodB/10); %tao tin hieu cho fading M=4; m=log2 (M); % Number of bits/symbol Es=1; % Symbol energy b=randint (1,N, [0 M-1]); % Random intergers from 0:M-1 skk=sqrt (Es)*exp (j* (2*pi/M*b+pi/M)); % Tao Rayleigh g=1/sqrt (2)* (randn (1,N)+j*randn (1,N)); for n=1:length (EbNodB) % tinh ber qua awgn %tao tin hieu cho fading si=2* (round (rand (1,N))-0.5); sq=2* (round (rand (1,N))-0.5); sk=si+j*sq; nk= (1/sqrt (2*EbNo (n)))* (randn (1,N)+j*randn (1,N)); yk=sk+nk; si_=sign (real (yk)); sq_=sign (imag (yk)); ber1= (N-sum (si==si_))/N; 93 ber2= (N-sum (sq==sq_))/N; ber (n)=mean ( [ber1 ber2]); %tinh ber qua fading g=1/sqrt (2)* (randn (1,N)+j*randn (1,N)); ykk=skk.*g+nk; % Coherent detection ykHat=conj (g).*ykk; % Decision thetak=angle (ykHat); thetak=mod (thetak+2*pi,2*pi); % Convert from [0,+/-pi] > [0,2*pi) bkHat=floor (M/ (2*pi)*thetak); % Error calculation noSymErr=sum ((bkHat-b)~=0); % Number of symbol errors SER (n)=noSymErr/N; % Symbol Error rate [noBitErr,BER (n)]=biterr (bkHat,b,m); end % BER ly thuyet cua truyen dan QPSK qua kenh AWGN Lythuyet_BER = 1/2*erfc (sqrt (EbNo*2)*sin (pi/4)) % BER ly thuyet cua truyen dan QPSK qua kenh Fading Lt_BER=1/2* (1-1./sqrt (1+1./EbNo)); % ve thi close all figure close all semilogy(EbNodB,BER,'bp-',EbNodB,ber,'rp',EbNodB,Lt_BER,'mp',EbNodB,Lythuyet_BER,'kp'); axis ( [0 20 10^-6 1]) legend ('Sim qua fading', 'Sim qua awgn','Theory qua fading','Theory qua awgn'); xlabel ('Es/No, dB') ylabel ('BER') title ('BER dieu che QPSK') 94 grid ========================================================== bonQAM.m clear all M=16; % Define parameters m=log2 (M); % Number of bits/symbol N=5*10^5; % So symbol gui di EbNodB=0:2:30; EbNo=10.^ (EbNodB./10); % Generate random integers ranging from 0:M-1 bk=randint (1,N, [0 M-1]); % Random intergers from 0:M-1 sk=qammod (bk,M); Es=norm (sk).^2/length (sk); % Noise variance of AWGN sigma=sqrt (Es./ (2*m.*EbNo)); % Iteration for estimating BER for k=1:length (EbNo) % Generate AWGN with variance sigma^2 nk=sigma (k)* (randn (1,N)+j*randn (1,N)); %tao awgn % Rx signal yk=sk+nk; bkHat=qamdemod (yk,M); % Error calculation noSymErr=sum ((bkHat-bk)~=0); % Number of symbol errors SER (k)=noSymErr/N; % Symbol Error rate [noBitErr,BER (k)]=biterr (bkHat,bk,m);% %tao fading gk=1/sqrt (2)* (randn (1,N)+j*randn (1,N)); ykk=sk.*gk+nk; 95 ykkHat=conj (gk).*ykk./abs (gk).^2; % Coherent detection % Decision bkkHat=qamdemod (ykkHat,M); % Error calculation nErr=sum ((bkkHat-bk)~=0); % Number of symbol errors ser (k)=nErr/N; % Symbol Error rate [noBitErr,ber (k)]=biterr (bkkHat,bk,m);% end % BER ly thuyet cua truyen dan 4QAM qua kenh AWGN Lythuyet_BER=1/2*erfc (sqrt (EbNo)); % BER ly thuyet cua truyen dan 4QAM qua kenh Fading Lythuyet_BER=3/8*erfc (sqrt (2/5*EbNo))-9/64* (erfc (2/5*EbNo)).^2; % ve thi close all figure close all semilogy (EbNodB,ber,'bp-',EbNodB,BER,'rp- ',EbNodB,Lt_BER,'mp',EbNodB,Lythuyet_BER,'kp'); axis ( [0 30 10^-6 1]) legend ('Sim qua fading', 'Sim qua awgn','Theory qua fading','Theory qua awgn'); xlabel ('Es/No, dB') ylabel ('BER') title ('BER dieu che 4QAM') grid ========================================================== muoisauQAM.m clear all M=16; % Define parameters m=log2 (M); % Number of bits/symbol N=10^5; % So symbol gui di 96 EbNodB=0:2:30; EbNo=10.^ (EbNodB./10); % Generate random integers ranging from 0:M-1 b=randint (1,N, [0 M-1]); % Random intergers from 0:M-1 sk=qammod (b,M); Es=norm (sk).^2/length (sk); % Noise variance of AWGN sigma=sqrt (Es./ (2*m.*EbNo)); % Iteration for estimating BER for k=1:length (EbNo) % Generate AWGN with variance sigma^2 n=sigma (k)* (randn (1,N)+j*randn (1,N)); %tao awgn % Rx signal yk=sk+n; bkHat=qamdemod (yk,M); % Error calculation noSymErr=sum ((bkHat-b)~=0); % Number of symbol errors SER (k)=noSymErr/N; % Symbol Error rate [noBitErr,BER (k)]=biterr (bkHat,b,m);% %tao fading g=1/sqrt (2)* (randn (1,N)+j*randn (1,N)); s=sk.*g; ykk=s+n; ykkHat=conj (g).*ykk./abs (g).^2; % Coherent detection % Decision bkkHat=qamdemod (ykkHat,M); % Error calculation nErr=sum ((bkkHat-b)~=0); % Number of symbol errors ser (k)=nErr/N; % Symbol Error rate [noBitErr,ber (k)]=biterr (bkkHat,b,m);% 97 end % BER ly thuyet cua truyen dan 16QAM qua kenh AWGN Lythuyet_BER=3/8*erfc (sqrt (2/5*EbNo))-9/64* (erfc (2/5*EbNo)).^2; % BER ly thuyet cua truyen dan 16QAM qua kenh Fading Lt_BER=3/8* (1-1./sqrt (1+5./ (2*EbNo))); % ve thi close all figure close all semilogy (EbNodB,ber,'bp ',EbNodB,BER,'rp ',EbNodB,Lt_BER,'mp ',EbNodB,Lythuyet_BER,'kp '); axis ( [0 30 10^-6 1]) legend ('Sim qua fading', 'Sim qua awgn','Theory qua fading','Theory qua awgn'); xlabel ('Es/No, dB') ylabel ('BER') title ('BER dieu che 16QAM') grid ========================================================== AWGNdieuche.m EbNodB=0:2:20; EbNo=10.^ (EbNodB./10); N=5*10^5; for k=1:length (EbNodB) % Tinh BER BPSK % Tao tin hieu BPSK {+1,-1} b=rand (1,N)>0.5; s=1-2*b; % Tinh toan nang luong bit tin hieu Eb Eb=norm (s)^2/N; % Mat AWGN No=Eb./EbNo; 98 nk=sqrt (No (k)./2).* (randn (1,N)+j*randn (1,N)); % Tin hieu thu y=s+nk; % Tach tin hieu shat=sign (real (y)); error=s-shat; nError=length (find (error~=0)); ber (k)=nError/N; % Tinh BER QPSK si=2* (round (rand (1,N))-0.5); sq=2* (round (rand (1,N))-0.5); Eb=norm (s)^2/N; % Mat AWGN No=Eb./EbNo; nk=sqrt (No (k)./2).* (randn (1,N)+j*randn (1,N)); sk=si+j*sq; yk=sk+nk; si_=sign (real (yk)); sq_=sign (imag (yk)); ber1= (N-sum (si==si_))/N; ber2= (N-sum (sq==sq_))/N; berk (k)=mean ( [ber1 ber2]); % Tinh BER 4QAM bkk=randint (1,N, [0 4-1]); % Random intergers from 0:M-1 skk=qammod (bkk,4); Es=norm (skk).^2/length (skk); sigma=sqrt (Es./ (4.*EbNo)); nkk=sigma (k)* (randn (1,N)+j*randn (1,N)); % Rx signal ykkk=skk+nkk; bkkkHat=qamdemod (ykkk,4); 99 % Error calculation noSymErr=sum ((bkkkHat-bkk)~=0); % Number of symbol errors serk (k)=noSymErr/N; % Symbol Error rate [nErr,berkk (k)]=biterr (bkkkHat,bkk,2); % Tinh BER 16QAM bkkk=randint (1,N, [0 16-1]); % Random intergers from 0:M-1 skkk=qammod (bkkk,16); Es=norm (skkk).^2/length (skkk); % Noise variance of AWGN sigma=sqrt (Es./ (8.*EbNo)); nkkk=sigma (k)* (randn (1,N)+j*randn (1,N)); % Rx signal ykkkk=skkk+nkkk; bkkkkHat=qamdemod (ykkkk,16); % Error calculation nSymErr=sum ((bkkkkHat-bkkk)~=0); % Number of symbol errors serkk (k)=nSymErr/N; % Symbol Error rate [nkErr,berkkk (k)]=biterr (bkkkkHat,bkkk,4);% end % Plot BER close all figure close all semilogy(EbNodB,ber,'bp',EbNodB,berk,'rp',EbNodB,berkk,'k.',EbNodB,berkkk,'mp-'); axis ( [0 20 10^-6 1]) legend ('BPSK', 'QPSK','4-QAM','16-QAM'); xlabel ('Es/No, dB') ylabel ('BER') title ('BER dieu che qua kenh AWGN') grid 100 ==========================================================F adingDieuche.m EbNodB=0:2:30; EbNo=10.^ (EbNodB./10); N=10^5; for k=1:length (EbNodB) % Tinh BER BPSK b = rand (1,N)>0.5; % generating 0,1 with equal probability s = 2*b-1; % BPSK modulation -> -1; -> n = 1/sqrt (2)* [randn (1,N) + j*randn (1,N)]; % white gaussian noise, 0dB variance h = 1/sqrt (2)* [randn (1,N) + j*randn (1,N)]; % Rayleigh channel % Channel and noise Noise addition y = h.*s + 10^ (-EbNodB (k)/20)*n; % equalization yHat = y./h; % receiver - hard decision decoding ipHat = real (yHat)>0; % counting the errors Err (k) = size (find ( [b- ipHat]),2); ber = Err/N; % Tinh BER QPSK Es=1; % Tao tin hieu QPSK bk=randint (1,N, [0 4-1]); % Random intergers from 0:M-1 sk=sqrt (Es)*exp (j* (2*pi/4*bk+pi/4)); % Tao Rayleigh fading g=1/sqrt (2)* (randn (1,N)+j*randn (1,N)); % Mat AWGN sigma=sqrt (Es./ (4.*EbNo)); % Tao AWGN voi variance sigma^2 101 nk=sigma (k)* (randn (1,N)+j*randn (1,N)); % Tin hieu thu yk=sk.*g+nk; % Coherent detection ykHat=conj (g).*yk; % Decision thetak=angle (ykHat); thetak=mod (thetak+2*pi,2*pi); % Convert from [0,+/-pi] > [0,2*pi) bkHat=floor (4/ (2*pi)*thetak); % Error calculation nErr=sum ((bkHat-bk)~=0); % Number of symbol errors ser (k)=nErr/N; % Symbol Error rate [nErr,berk (k)]=biterr (bkHat,bk,2); % Tinh BER 4QAM bkk=randint (1,N, [0 4-1]); % Random intergers from 0:M-1 skk=qammod (bkk,4); Esk=norm (skk).^2/length (skk); skkk=skk.*g; % Faded signal % Noise variance of AWGN sigmak=sqrt (Esk./ (4.*EbNo)); % Generate AWGN with variance sigma^2 nkk=sigmak (k)* (randn (1,N)+j*randn (1,N)); % Rx signal ykk=skkk+nkk; ykkHat=conj (g).*ykk./abs (g).^2; % Coherent detection % Decision bkkHat=qamdemod (ykkHat,4); % Error calculation nkErr=sum ((bkkHat-bkk)~=0); % Number of symbol errors serk (k)=nkErr/N; % Symbol Error rate [nkErr,berkk (k)]=biterr (bkkHat,bkk,2); 102 % Tinh BER 16QAM bkkk=randint (1,N, [0 16-1]); % Random intergers from 0:M-1 skkk=qammod (bkkk,16); Eskk=norm (skkk).^2/length (skkk); skkkk=skkk.*g; % Faded signal % Noise variance of AWGN sigmakk=sqrt (Eskk./ (8.*EbNo)); % Generate AWGN with variance sigma^2 nkkk=sigmakk (k)* (randn (1,N)+j*randn (1,N)); % Rx signal ykkk=skkkk+nkkk; ykkkHat=conj (g).*ykkk./abs (g).^2; % Coherent detection % Decision bkkkHat=qamdemod (ykkkHat,16); % Error calculation nkkErr=sum ((bkkkHat-bkkk)~=0); % Number of symbol errors serkk (k)=nkkErr/N; % Symbol Error rate [nkkErr,berkkk (k)]=biterr (bkkkHat,bkkk,4); end % Plot BER close all figure close all semilogy(EbNodB,ber,'bp-',EbNodB,berk,'rp-',EbNodB,berkk,'k.',EbNodB,berkkk,'mp-'); axis ( [0 30 10^-6 1]) legend ('BPSK', 'QPSK','4-QAM','16-QAM'); xlabel ('Es/No, dB') ylabel ('BER') title ('BER dieu che qua kenh Fading Rayleigh') grid 103 ... tốn, đánh giá chất lượng hệ thống trở nên khó khăn so với kênh truyền hữu tuyến Đề tài “Đán g c ất lượng h thống thông tin vô tuyến sử dụng phần mềm Matlab? ?? phương tiện biểu đạt mới, sử dụng. .. Tóm lại, thơng tin viba bao gồm: - Thông tin vệ tinh hệ thống thông tin vô tuyến chuyển tiếp, trạm chuyển tiếp đặt vệ tinh ngồi Trái đất - Thơng tin di động hệ thống thông tin vô tuyến chuyển tiếp... thống thông tin vô tuyến Chương 2: Truyền dẫn kênh vô tuyến số Chương 3: Mô đánh giá chất lượng hệ thống thông tin vô tuyến sử dụng phần mềm matlab Xin trân trọng cảm ơn ThS Phạm Mạnh Tồn giới