Phân tích code matlab sim uwb 402

1

HỌC VIỆN CÔNG NGHỆ BƯU CHÍNH VIỄN THÔNG

KHOA VIỄN THÔNG I

- -

Đề bài:

BÀI TẬP VỀ NHÀ CÁC MẠNG THÔNG TIN VÔ TUYẾN

“Phân tích code Matlab Sim_UWB_402”

Hà Nội – 2022

Mục lục

1 Sim_UWB_402a 2

 Code Matlab: 7

 Kết quả mô phỏng 8

2 Sim_UWB_402b 8

 Code Matlab 11

2

1 Sim_UWB_402a

4

6

 Code Matlab:

% function Sim_UWB_402a

clc;

0, 1 hàng N cột

0, 2 hàng N cột

% Rayleigh model

Pro 4.05

dùng để xây dựng biểu đồ mật độ xác suất dựa trên mẫu có sẵn

% Rician model

Pro 4.06

[temp x]=hist(abs(Rician_ch(i,:)),level);

end

','fontname','.vntime','fontsize',14);

Rician',

'fontname','.vntime','fontsize',14);

8

2 Sim_UWB_402b

10

 Code Matlab

%===========================================================

===============

clc;

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;

% 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));

12

[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));

mu_l = sum((C1.*cos(2*pi*F1*f_max*(n*T_s + t_0) +

(sum((C2.*cos(2*pi*F2*f_max*(n*T_s + t_0) +

RHO.*exp(j*2*pi*F_RHO*f_max*(n*T_s + t_0));

h_t(n+1) = sum(mu_l);

end

H_t = abs(h_t);

figure(49)

subplot(2,2,1);

H_t = H_t/max(H_t);

plot(t,20*log10(H_t));

[s]','FontName','.VnTime','FontSize',12);

[dB]','FontName','.VnTime','FontSize',14);

subplot(2,2,3);

th1 = rand(size(tau_l))*2*pi;

e_l=exp(-j*(2*pi*k*f_s*tau_l + th1 ));

s_f = abs(s_f);

S_f = abs(s_f);

S_f = S_f/max(S_f);

[dB]','FontName','.VnTime','FontSize',14);

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));

[s]','FontName','.VnTime','FontSize',12);

[dB]','FontName','.VnTime','FontSize',14);