реф. Исследование методов повышения пропускной способности в сетях UMTS

title('UL Load Factor vs. number of users');

% Расчет выигрыша от обработки

Gop = 10*log10(W./R);

% Расчет предельной пропускной способности при различных типах соты

figure (1)

Npr1 = (1+(Gop./(ebn01.*v))).*(1+i1);

Npr12 = (1+(Gop./(ebn01.*v))).*(1+i2);

Npr13 = (1+(Gop./(ebn01.*v))).*(1+i3); %plot(ebn01,Npr1,'b*-','LineWidth',2)

plot(ebn01,Npr1,'r*-',ebn01,Npr12,'b*-',ebn01,Npr13,'g*- ','LineWidth',2)

legend('Picocell','Microcell', 'Macrocell'); grid

%legend('received power','throughput'); title('Number of connections vs. Eb/N0'); xlabel('Eb/N0 , dB');

ylabel('Number of connections');

94

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

Оценка пропускной способности соты при изменении состава оборудования в нисходящем канале

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

Исходные данные

ebn0 = 5; % отношение Eb/N0

W = 3840; % чиповая скорость в сети UMTS

R = 12.2; % скорость передачи речи AMR

N = 0:5:60; % число пользователей в соте

v = 0.67; % коэффициент активности абонентов i = .65; % отношение помех для в соте

f1 = .825; % коэффициент ортогональности при открыт. Мест-ти f2 = .65; % коэффициент ортогональности в пригороде

f3 = .525; % коэффициент ортогональности в городе

f4 = .4; % коэффициент ортогональности при плот. Застройке

% Расчет коэффициента нагрузки и увеличения помех при открытой Местности figure(1);

qul1 = (ebn0./(W./R)).*N.*(1-f1+i); NR1 = 1./(1-qul1);

NRdb1 = -10*log10(1-qul1);

plot (N,NRdb1,'*g-','LineWidth',1) hold on;

% Расчет коэффициента нагрузки и увеличения помех в пригороде qul2 = (ebn0./(W./R)).*N.*(1-f2+i);

NR2 = 1./(1-qul2);

NRdb2 = -10*log10(1-qul2);

plot (N,NRdb2,'*r-','LineWidth',1) hold on;

% Расчет коэффициента нагрузки и увеличения помех в городе qul3 = (ebn0./(W./R)).*N.*(1-f3+i);

NR3 = 1./(1-qul3);

NRdb3 = -10*log10(1-qul3);

plot (N,NRdb3,'*b-','LineWidth',1) hold on;

% Расчет коэффициента нагрузки и увеличение помех при плот. Застройке qul4 = (ebn0./(W./R)).*N.*(1-f4+i);

NR4 = 1./(1-qul4);

NRdb4 = -10*log10(1-qul4);

plot (N,NRdb4,'*c-','LineWidth',1)

95

grid on;

legend('Open country','Suburb', 'City','Dense buildings');

%Расчет выигрыша от обработки для различных типов местности

Gop = (W./R);

%Расчет предельной пропускной способности для различных типов местности

figure (1)

Npr1 = (1+(Gop./(ebn01.*v))).*(1-f1+i);

Npr12 = (1+(Gop./(ebn01.*v))).*(1-f2+i);

Npr13 = (1+(Gop./(ebn01.*v))).*(1-f3+i);

Npr14 = (1+(Gop./(ebn01.*v))).*(1-f4+i);

%plot(ebn01,Npr1,'b*-','LineWidth',2) plot(ebn01,Npr1,'r*',ebn01,Npr12,'b*-',ebn01,Npr13,ebn01,Npr14,'c*-

…','LineWidth',2)

legend('Open country','Suburb', 'City','Dense buildings'); grid

%legend('received power','throughput'); title('Number of connections vs. Eb/N0'); xlabel('Eb/N0 , dB');

ylabel('Number of connections');

96

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

Оценка относительной загрузки соты при использовании AMR-кодека в восходящем каналах

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

clc

format short

% Исходные данные

ebn01=0:1:15;% отношение Eb/N0

W = 3840; % чиповая скорость в сети UMTS

R1 = 12.2; % скорость передачи речи AMR 12.2

R2 = 7.4; % скорость передачи речи AMR 7.4

R3 = 4.75; % скорость передачи речи AMR 4.75

N = 0:5:50; % число пользователей в соте

v = 0.67; % коэффициент активности абонентов i = .55; % отношение помех от других сот

f = .67; % коэффициент ортогональности

T = .02; % нагрузка создаваемая одним абонентом

N = 140; % число пользователей в соте

% Расчет выигрыша от обработки при различных скоростях AMR

Gop1 = 10*log10(W./R1); % выигрыша от обработки AMR 12.2

Gop2 = 10*log10(W./R2); % выигрыша от обработки AMR 7.4

Gop3 = 10*log10(W./R3); % выигрыша от обработки AMR 4.75

% Расчет относительной загрузки в восходящем канале для AMR 12.2 figure (1)

Npr1 = (1+(Gop1./(ebn01.*v))).*(1+i); qdl1 = (T.*N)./Npr1;

% Расчет относительной загрузки в восходящем канале для AMR 7.4

Npr12 = (1+(Gop2./(ebn01.*v))).*(1+i); qdl2 = (T.*N)./Npr12;

% Расчет относительной загрузки в восходящем канале для AMR 4.75 Npr13 = (1+(Gop3./(ebn01.*v))).*(1+i);

qdl3 = (T.*N)./Npr13; %plot(ebn01,Npr1,'b*-','LineWidth',2)

plot(ebn01,qdl1,'r*-',ebn01,qdl2,'b*-',ebn01,qdl3,'g*-','LineWidth',2) legend('AMR 12.2','AMR 7.4', 'AMR 4.75');

grid

%legend('received power','throughput'); title('Fractional load UL-AMR vs. Eb/N0'); xlabel('Eb/N0 , dB');

ylabel('Fractional load');

97

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

Оценка относительной загрузки сот при использовании AMR-кодека в нисходящем каналах

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

clc

format short

% Исходные данные

ebn01=0:1:15;% отношение Eb/N0

W = 3840; % чиповая скорость в сети UMTS

R1 = 12.2; % скорость передачи речи AMR 12.2

R2 = 7.4; % скорость передачи речи AMR 7.4

R3 = 4.75; % скорость передачи речи AMR 4.75

N = 0:5:50; % число пользователей в соте

v = 0.67; % коэффициент активности абонентов i = .55; % отношение помех от других сот

f = .67; % коэффициент ортогональности

T = .02; % нагрузка создаваемая одним абонентом

N = 140; % число пользователей в соте

% Расчет выигрыша от обработки при различных скоростях AMR

Gop1 = 10*log10(W./R1); % выигрыша от обработки AMR 12.2

Gop2 = 10*log10(W./R2); % выигрыша от обработки AMR 7.4

Gop3 = 10*log10(W./R3); % выигрыша от обработки AMR 4.75

% Расчет относительной загрузки в нисходящем канале для AMR 12.2 figure (1)

Npr1 = (1+(Gop1./(ebn01.*v))).*(1-f+i); qdl1 = (T.*N)./Npr1;

% Расчет относительной загрузки в нисходящем канале для AMR 7.4

Npr12 = (1+(Gop2./(ebn01.*v))).*(1-f+i); qdl2 = (T.*N)./Npr12;

% Расчет относительной загрузки в нисходящем канале для AMR 4.75

Npr13 = (1+(Gop3./(ebn01.*v))).*(1-f+i); qdl3 = (T.*N)./Npr13; %plot(ebn01,Npr1,'b*-','LineWidth',2)

plot(ebn01,qdl1,'r*-',ebn01,qdl2,'b*-',ebn01,qdl3,'g*-','LineWidth',2) legend('AMR 12.2','AMR 7.4', 'AMR 4.75');

grid

title('Fractional load DL-AMR vs. Eb/N0'); xlabel('Eb/N0 , dB');

ylabel('Fractional load');

98

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

Оценка относительной загрузки соты при изменении состава оборудования в нисходящем каналах

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

clc

format short

% Исходные данные

ebn01 =0:1:15;% отношение Eb/N0

W = 3840; % чиповая скорость в сети UMTS

R = 12.2; % скорость передачи речи

N = 0:5:50; % число пользователей в соте

v = 0.67; % коэффициент активности абонентов i1 = .15; % отношение помех в пикосоте

i2 = .3; % отношение помех в микросоте

i3 = .45; % отношение помех в макросоте f = .67; % коэффициент ортогональности

T = .02; % нагрузка создаваемая одним абонентом

N = 140; % число пользователей в соте

% Расчет выигрыша от обработки

Gop = 10*log10(W./R);

% Расчет относительной загрузки в восходящем канале для пикосоты figure (1)

Npr1 = (1+(Gop./(ebn01.*v))).*(1+i1); qdl1 = (T.*N)./Npr1;

% Расчет относительной загрузки в восходящем канале для микросоты

Npr12 = (1+(Gop./(ebn01.*v))).*(1+i2); qdl2 = (T.*N)./Npr12;

% Расчет относительной загрузки в восходящем канале для макросоты

Npr13 = (1+(Gop./(ebn01.*v))).*(1+i3); qdl3 = (T.*N)./Npr13; %plot(ebn01,Npr1,'b*-','LineWidth',2)

plot(ebn01,qdl1,'r*-',ebn01,qdl2,'b*-',ebn01,qdl3,'g*-','LineWidth',2) legend('Picocell','Microcell', 'Macrocell');

grid

%legend('received power','throughput'); title('Fractional load vs. Eb/N0'); xlabel('Eb/N0 , dB'); ylabel('Fractional load');

99

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

Оценка относительной загрузки соты при изменении состава оборудования в нисходящем каналах

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

clc

format short

% Исходные данные

ebn01 =0:1:15;% отношение Eb/N0

W = 3840; % чиповая скорость в сети UMTS

R = 12.2; % скорость передачи речи

N = 0:5:50; % число пользователей в соте

v = 0.67; % коэффициент активности абонентов i = .15; % отношение помех в соте

f = .825; % коэффициент ортогональностив откыт-ой местн-сти f2 = .65; % коэффициент ортогональности в пригороде

f3 = .525; % коэффициент ортогональности в городе

f4 = .4;% коэффициент ортогональности в плотной застройке

T = .02; % нагрузка создаваемая одним абонентом

N = 140; % число пользователей в соте

% Расчет выигрыша от обработки

Gop = 10*log10(W./R);

% Расчет относительной загрузки в нисходящем канале при откыт-ой местн-сти figure (1)

Npr1 = (1+(Gop./(ebn01.*v))).*(1-f1+i); qdl1 = (T.*N)./Npr1;

% Расчет относительной загрузки в нисходящем канале для пригорода

Npr12 = (1+(Gop./(ebn01.*v))).*(1-f2+i); qdl2 = (T.*N)./Npr12;

% Расчет относительной загрузки в нисходящем канале для города

Npr13 = (1+(Gop./(ebn01.*v))).*(1-f3+i); qdl3 = (T.*N)./Npr13;

% Расчет относительной загрузки в нисходящем канале при плотной застройке

Npr14 = (1+(Gop./(ebn01.*v))).*(1-f4+i); qdl4 = (T.*N)./Npr14;

plot(ebn01,qdl1,'r*-',ebn01,qdl2,'b*-',ebn01,qdl3,'c*-',ebn01,qdl4,'g*-

…','LineWidth',2)

legend('Open country','Suburb', 'City','Dense buildings'); grid

%legend('received power','throughput'); title('Fractional load vs. Eb/N0'); xlabel('Eb/N0 , dB'); ylabel('Fractional load');

100