Phân tích code matlab sim uwb 401

HỌC VIỆN CƠNG NGHỆ BƯU CHÍNH VIỄN THƠNG KHOA VIỄN THÔNG I - - BÀI TẬP VỀ NHÀ CÁC MẠNG THÔNG TIN VÔ TUYẾN Đề bài: “Phân tích code Matlab Sim_UWB_401” Giảng viên : Nguyễn Viết Đảm Hà Nội – 2022 MỤC LỤC Sim_UWB 401 Kết mô 2.1 Tổn hao đường truyền pha đinh phạm vi rộng 2.2 Tổn hao đường truyền mơ hình IEEE 802.16d Sim_UWB 401 Chương trình thực tính tốn biểu diễn tổn hao đường truyền môi trường kênh pha đinh phạm vi rộng tổn hao đường truyền mơ hình kênh IEEE 802.16d cho chương trình Sim_UWB 4.01 Kết chạy chương trình biểu diễn hình 2.1 2.2 phân tích sau Muốn ta phải thiết trị cho tham số đặc trưng cho môi trường kênh pha đinh phạm vi rộng, cụ thể thấy chương trình Dưới ta tóm tắt chương trình dùng Sim_UWB 4.01 Prog 4.01 Chương trình thực tính suy hao đường truyền khơng gian tự Chương trình nhận 04 đầu vào: tần số sóng mang, fc; khoảng cách phát/thu, dist; độ lợi anten phát, Gt; độ lợi anten thu, Gr Kết đầu chương trình tổn thất đường truyền khơng gian tự do, PL Thực chương trình lệnh: [PL] = uwb40101_PL_free(fc, dist, Gt, Gr) Prog 4.02 Chương trình thực tính suy hao đường truyền theo khoảng cách môi trường pha đinh phạm vi rộng che chắn phân bố log-normal Chương trình nhận 05 đầu vào: tần số sáng mang: fc; khoảng cách phát/thu: d; khoảng cách tham chiếu: d0; số mũ tổn hao truyền sóng: n; phương sai: sigma Đầu chương trình suy hao đường truyền môi trường che chắn phân bố log-normal, PL Lệnh thực chương trình là: [PL] = uwb40102_PL_logdist_or_norm(fc, d, d0, n, sigma) Prog 4.03 Chương trình tính tổn thất đường truyền theo mơ hình Hata Chương trình nhận đầu vào: tần số sáng mang: fc; khoảng cách phát/thu: d; độ cao anten phát htx; độ cao anten thu: hrx; loại môi trường ('urban','suburban','open'), Etype Kết đầu chương trình suy hao đường truyền theo mơ hình Hata, PL Lệnh thực chương trình là: [PL] = uwb40104_PL_Hata(fc, d, htx, hrx, Etype) Prog 4.04 Chương trình thực tính tổn hao đường truyền theo mơ hình IEEE 802.16d Chương trình nhận đầu vào: tần số sóng mang: fc; khoảng cách phát/thu: d; độ cao anten phát: htx; độ cao anten thu: hrx; loại môi trường ('A', 'B' 'C' ): Etype; hệ số tương quan: corr_fact; thiết lập loại mơ hình chọn: mod Đầu chương trình suy hao đường truyền theo mơ hình IEEE 802.16d, PL Lệnh thực chương trình là: [PL] = uwb40105_PL_IEEE80216d(fc, d, type, htx, hrx, corr_fact, mod)  Code matlab %Sim_UWB_401 % % % % Prog 4.01 % [PL] = uwb40101_PL_free(fc, dist, Gt, Gr) % % -% % % % Prog 4.02 % [PL] = uwb40102_PL_logdist_or_norm(fc, d, d0, n, sigma) % % -% % % % Prog 4.03 % [PL] = uwb40104_PL_Hata(fc, d, htx, hrx, Etype) % % -% % % % Prog 4.04 % [PL] = uwb40105_PL_IEEE80216d(fc, d, type, htx, hrx, corr_fact, mod) %============================================================== ============ clc; clear all; close all; fc =1.5e9; d0 =100; sigma =3; distance = [1:2:31].^2; Gt = [1 0.5]; Gr = [1 0.5 0.5]; Exp = [2 6]; for k=1:3 y_Free(k,:) = uwb40101_PL_free(fc,distance,Gt(k),Gr(k)); % Pro 4.01 y_logdist(k,:) = uwb40102_PL_logdist_or_norm(fc,distance,d0,Exp(k)); % Pro 4.02 y_lognorm(k,:) = uwb40102_PL_logdist_or_norm(fc,distance,d0,Exp(1),sigma); % end figure(4345) subplot(141) semilogx(distance,y_Free(1,:),'k-o',distance,y_Free(2,:),'b-^',distance,y_Free(3,:),'r-s'); grid on; axis([1 1000 40 110]); title(['Cac mo hinh suy hao khong gian tu do, f_c=',num2str(fc/1e6),'MHz'],'fontname','.vntime','fontsize',10); xlabel('Khoang cach [m]','fontname','.vntime','fontsize',12); ylabel('Suy hao [dB]','fontname','.vntime','fontsize',12); legend('G_t=1, G_r=1','G_t=1, G_r=0.5','G_t=0.5, G_r=0.5'); subplot(142) semilogx(distance,y_logdist(1,:),'k-o',distance,y_logdist(2,:),'b^',distance,y_logdist(3,:),'r-s'); grid on; axis([1 1000 40 110]); title(['Mo hinh PL Log-khoang cach, f_c=',num2str(fc/1e6),'MHz'],'fontname','.vntime','fontsize',10); xlabel('Khoang cach [m]','fontname','.vntime','fontsize',12); ylabel('Suy hao [dB]','fontname','.vntime','fontsize',12); legend('n=2','n=3','n=6') subplot(143) semilogx(distance,y_lognorm(1,:),'k-o',distance,y_lognorm(2,:),'b^',distance,y_lognorm(3,:),'r-s'); grid on; axis([1 1000 40 110]); title(['Mo hinh PL Log-normal, f_c=',num2str(fc/1e6),'MHz, ','\sigma=', num2str(sigma), 'dB', ],'fontname','.vntime','fontsize',10); xlabel('Khoang cach[m]','fontname','.vntime','fontsize',12); ylabel('Suy hao[dB]','fontname','.vntime','fontsize',12); legend('path 1','path 2','path 2'); %%%%%%=================Hata PL model======================================= clear; clc; fc =1.5e9; htx =30; hrx =2; distance = [1:2:31].^2; y_urban = uwb40104_PL_Hata(fc,distance,htx,hrx,'urban'); % Pro 4.03 y_suburban = uwb40104_PL_Hata(fc,distance,htx,hrx,'suburban'); y_open = uwb40104_PL_Hata(fc,distance,htx,hrx,'open'); subplot(144) semilogx(distance,y_urban,'k-s', distance,y_suburban,'k-o', distance,y_open,'k-^'); grid on; axis([1 1000 40 110]); title(['Mo hinh PL Hata, f_c=',num2str(fc/1e6),'MHz'],'fontname','.vntime','fontsize',10); xlabel('Khoang cach[m]','fontname','.vntime','fontsize',12); ylabel('Suy hao [dB]','fontname','.vntime','fontsize',12); legend('urban','suburban','open area'); %%%=============== IEEE 802.16d Path loss Models=========================== clc; clear all; fc=2e9; htx=[30 30]; hrx=[2 10]; distance=[1:1000]; for k=1:2 y_IEEE16d(k,:) = uwb40105_PL_IEEE80216d(fc,distance,'A',htx(k),hrx(k),'atnt'); % Pro 4.04 y_MIEEE16d(k,:) = uwb40105_PL_IEEE80216d(fc,distance,'A',htx(k),hrx(k),'atnt','mod'); end figure(44) subplot(121) semilogx(distance,y_IEEE16d(1,:),'k:','linewidth',1.5), hold on semilogx(distance,y_IEEE16d(2,:),'k-','linewidth',1.5), grid on title(['Mo hinh ton hao IEEE 802.16d , f_c=',num2str(fc/1e6),'MHz'],'fontname','.vntime','fontsize',10); axis([1 1000 10 150]); xlabel('Khoang cach [m]','fontname','.vntime','fontsize',12); ylabel('Suy hao [dB]','fontname','.vntime','fontsize',12); legend('h_{Tx}=30m, h_{Rx}=2m','h_{Tx}=30m, h_{Rx}=10m') subplot(122) semilogx(distance,y_MIEEE16d(1,:),'k:','linewidth',1.5), hold on semilogx(distance,y_MIEEE16d(2,:),'k-','linewidth',1.5), grid on title(['Mo hinh ton bao IEEE 802.16d biÕn ®æi, f_c=', num2str(fc/1e6), 'MHz'], 'fontname','.vntime','fontsize',10); axis([1 1000 10 150]); xlabel('Khoang cach [m]','fontname','.vntime','fontsize',12); ylabel('Ton hao [dB]','fontname','.vntime','fontsize',12); legend('h_{Tx}=30m, h_{Rx}=2m','h_{Tx}=30m, h_{Rx}=10m') Kết mô 2.1 Tổn hao đường truyền pha đinh phạm vi rộng Hình 1: Tổn hao đường truyền pha đinh phạm vi rộng 2.2 Tổn hao đường truyền mơ hình IEEE 802.16d Hình 2: Tổn hao đường truyền mơ hình IEEE 802.16d 10 ... PL Lệnh thực chương trình là: [PL] = uwb4 0105_PL_IEEE80216d(fc, d, type, htx, hrx, corr_fact, mod)  Code matlab %Sim_ UWB_ 401 % % % % Prog 4.01 % [PL] = uwb4 0101_PL_free(fc, dist, Gt, Gr) % %...MỤC LỤC Sim_ UWB 401 Kết mô 2.1 Tổn hao đường truyền pha đinh phạm vi rộng 2.2 Tổn hao đường truyền mơ hình IEEE 802.16d Sim_ UWB 401 Chương trình thực... k=1:3 y_Free(k,:) = uwb4 0101_PL_free(fc,distance,Gt(k),Gr(k)); % Pro 4.01 y_logdist(k,:) = uwb4 0102_PL_logdist_or_norm(fc,distance,d0,Exp(k)); % Pro 4.02 y_lognorm(k,:) = uwb4 0102_PL_logdist_or_norm(fc,distance,d0,Exp(1),sigma);