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_402” Giảng viên : Nguyễn Viết Đảm Hà Nội – 2022 Mục lục Sim_UWB_402a Code Matlab: Kết mô Sim_UWB_402b Code Matlab .11 1 Sim_UWB_402a Code Matlab: % function Sim_UWB_402a clc; clear all; close all; N = 200000; gán giá trị N= 200000 level = 30; gán giá trị level =30 K_dB = [-40 15]; gán ma trận Rayleigh_ch = zeros(1,N); tạo ma trận có gtri 0, hàng N cột Rician_ch = zeros(2,N); tạo ma trận có gtri 0, hàng N cột color line marker = ['k']; = ['-']; = ['s','o','^']; % Rayleigh model Rayleigh_ch = uwb40201_Ray_model(N); % Pro 4.05 [temp,x] = hist(abs(Rayleigh_ch(1,:)),level); lệnh dùng để xây dựng biểu đồ mật độ xác suất dựa mẫu có sẵn plot(x,temp,['k-' marker(2)],'linewidth',2,'color','b'); hold on; % Rician model for i=1:length(K_dB); Rician_ch(i,:)=uwb40202_Ric_model(K_dB(i),N); % Pro 4.06 [temp x]=hist(abs(Rician_ch(i,:)),level); plot(x,temp,['k-' marker(i+1)],'linewidth',1.5,'color','r'); end xlabel('x','fontname','.vntime','fontsize',12); ylabel('Mật độ xác suất ','fontname','.vntime','fontsize',14); PT = legend('Rayleigh','Rician, K=-40dB','Rician, K=15dB'); set(PT,'FontName','.VnTime','FontSize',14); grid on; title('Ph©n bè cđa kênh pha đinh Rayleigh kênh pha đinh Rician', 'fontname','.vntime','fontsize',14); Kết mô Sim_UWB_402b 10 Code Matlab %=========================================================== =============== clc; clear all; close all; T_s = 0.1E-5; T_sim = 0.5E-1; t = 0:T_s:T_sim; m_s = 10; B = 1E+6; f_s = 2E+3; f = f_s:f_s:B; AREA = 'bu'; % input('Seclect area for simulation : [ra;tu;bu;ht]= '); %['ra';'tu';'bu';'ht']; f_max = 90; t_0 = 0.01; N_1 = 10; %=========================================================== =============== %===== Generate & plot time-domain input signal of the channel simulator %=========================================================== =============== x_t = ones(1,length(t)); 11 [C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,q_l,T,tau_l] = uwb40203_F_S_K_pcal2(N_1,AREA,f_max,T_s); %%%%%%%%%%%%%%%%%%%%%%%%%%% %%Initialization; %%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%% mu_1=zeros(size(q_l)); e_l = zeros(size(q_l)); h_t = zeros(1,length(t)); for n =0:length(t)-1, if (n/m_s) == fix(n/m_s), mu_l = sum((C1.*cos(2*pi*F1*f_max*(n*T_s + t_0) + TH1)).').* exp(-j*2*pi*F01*f_max*(n*T_s + t_0)) + j* (sum((C2.*cos(2*pi*F2*f_max*(n*T_s + t_0) + TH2)).').* exp(-j*2*pi*F02*f_max*(n*T_s + t_0))) + RHO.*exp(j*2*pi*F_RHO*f_max*(n*T_s + t_0)); end h_t(n+1) = sum(mu_l); end H_t = abs(h_t); save Sim_UWB_402b.mat; figure(49) subplot(2,2,1); H_t = H_t/max(H_t); plot(t,20*log10(H_t)); grid on; xlabel('Thời gian [s]','FontName','.VnTime','FontSize',12); ylabel(' Biên độ [dB]','FontName','.VnTime','FontSize',14); title('Đáp ứng kênh CIR miÒn thêi gian ', 'FontName','.VnTime','color','b','FontSize',16); subplot(2,2,3); s_f = zeros(1,length(f));% th1 = rand(size(tau_l))*2*pi; for k=1:length(f), e_l=exp(-j*(2*pi*k*f_s*tau_l + th1 )); s_f(k) = sum(e_l);% end s_f = abs(s_f); S_f = abs(s_f); 12 S_f = S_f/max(S_f); plot(f,20*log10(S_f),'r'); grid on; xlabel('TÇn sè [Hz]','FontName','.VnTime','FontSize',12) ylabel('Biên độ [dB]','FontName','.VnTime','FontSize',14); title('Đáp ứng kênh miền tần sè', 'FontName','.VnTime','color','b','FontSize',16); subplot(2,2,[2,4]); size(H_t) size(S_f) [H,S] = meshgrid(H_t,S_f); H_ft = H.*S; [x,y] = size(H_ft); t=linspace(0,T_sim,y); f=linspace(0,B,x); mesh(t,f,20*log10(H_ft)); xlabel('Thêi gian [s]','FontName','.VnTime','FontSize',12); ylabel('Tần số [Hz]','FontName','.VnTime','FontSize',12) zlabel('Biên độ [dB]','FontName','.VnTime','FontSize',14); title('Đáp ứng kênh miền thời gian tần số ', 'FontName','.nTimeh','color','b','FontSize',16); 13 ... lục Sim_ UWB_ 402a Code Matlab: Kết mô Sim_ UWB_ 402b Code Matlab .11 1 Sim_ UWB_ 402a Code Matlab: % function Sim_ UWB_ 402a... set(PT,''FontName'',''.VnTime'',''FontSize'',14); grid on; title( ''Phân bố kênh pha đinh Rayleigh kênh pha đinh Rician'', ''fontname'',''.vntime'',''fontsize'',14); Kt qu mô Sim_ UWB_ 402b 10 Code Matlab %===========================================================... clear all; close all; T_s = 0.1E-5; T _sim = 0.5E-1; t = 0:T_s:T _sim; m_s = 10; B = 1E+6; f_s = 2E+3; f = f_s:f_s:B; AREA = ''bu''; % input(''Seclect area for simulation : [ra;tu;bu;ht]= ''); %[''ra'';''tu'';''bu'';''ht''];