移动通信

Path Loss Using HATA for CDMA
clc;close all; clear all;
%for_big_cities
fc=1900;
hb=50;
hm=1.5;
dis_big_city=1.675;
dis_medium_city=3.84;
dis_rural_area=15.13;
total_area_big_city=100;
total_area_medium_city=200;
total_area_rural_area=700;
if(300<=fc)
ch_b=(8.29*(log10(1.54hm))^2) -1.1;
disp(‘r’)
elseif(300>=fc)
disp(‘b’)
ch_b=(3.2(log10(11.75hm))^2) -4.97;
end
Path_loss_big_cities=69.55+26.16log10(fc)-13.82log10(hb)-ch_b+(44.9-6.55log10(hb))log10(dis_big_city)
%small_cities
ch=.8+(1.1log10(fc)-.7)hm-1.56log10(fc);
LU=69.55+26.16log10(fc)-13.82log10(hb)-ch+(44.9-6.55*log10(hb))log10(dis_medium_city)
Path_loss_small_cities=LU-(2(log10(fc./28)).^2) -5.4
%open or rural areas

LU=69.55+26.16log10(fc)-13.82log10(hb)-ch+(44.9-6.55log10(hb))log10(dis_rural_area)
Path_loss_rural_areas=LU-(4.78(log10(fc))^2) +18.33log10(fc) -40.94
no_of_cell_media_city=total_area_big_city/dis_big_city
no_of_cell_medium_city=total_area_medium_city/dis_medium_city
no_of_cell_rural=total_area_rural_area/dis_rural_area
Path Loss Using HATA for LTE
clc;close all; clear all;
%for_big_cities
fc=1900;
hb=30;
hm=1.5;
dis_big_city=6.45;
dis_medium_city=14.3;
dis_rural_area=53;
total_area_big_city=100;
total_area_medium_city=200;
total_area_rural_area=700;
if(300<=fc)
ch_b=(8.29*(log10(1.54hm))^2) -1.1;
disp(‘r’)
elseif(300>=fc)
disp(‘b’)
ch_b=(3.2(log10(11.75*hm))^2) -4.97;
end

Path_loss_big_cities=69.55+26.16log10(fc)-13.82log10(hb)-ch_b+(44.9-6.55log10(hb))log10(dis_big_city)
%small_cities
ch=.8+(1.1log10(fc)-.7)hm-1.56log10(fc);
LU=69.55+26.16log10(fc)-13.82log10(hb)-ch+(44.9-6.55log10(hb))log10(dis_medium_city)
Path_loss_small_cities=LU-(2(log10(fc./28)).^2) -5.4
%open or rural areas
LU=69.55+26.16log10(fc)-13.82log10(hb)-ch+(44.9-6.55log10(hb))log10(dis_rural_area)
Path_loss_rural_areas=LU-(4.78(log10(fc))^2) +18.33log10(fc) -40.94
no_of_cell_media_city=total_area_big_city/dis_big_city
no_of_cell_medium_city=total_area_medium_city/dis_medium_city
no_of_cell_rural=total_area_rural_area/dis_rural_area
Bandwidth
clc;close all;clear all
R=14.410^3;
s_r=5.1;
b=0.85;
v=0.4
%big_city
N_b=98;
B_bc=(N_b-1)(s_r*(1+b)v)R
%sub_urban
N_s=72;
B_su=(N_s-1)(s_r(1+b)v)R
%sub_rural
N_r=58;
B_rural=(N_r-1)(s_r(1+b)*v)*R

Erlang
clc;close all;clear all;
area_media_city=100;
area_sub_urban=200;
area_rural=1000-200-100;
user_density_per_km_mc=500;
user_density_per_km_su=100;
user_density_per_km_rural=5;
T_subscriber_media_city=(area_media_cityuser_density_per_km_mc)
T_subscriber_sub_urban=(area_sub_urbanuser_density_per_km_su)
T_subscriber_rural=(area_ruraluser_density_per_km_rural)
total_subscriber=(area_media_cityuser_density_per_km_mc)+(area_sub_urbanuser_density_per_km_su)+(area_ruraluser_de
nsity_per_km_rural)
calls_per_hour=5;
holding_time_in_hours=2/60
total_traffic_during_the_busy_hour=total_subscribercalls_per_hourholding_time_in_hours
sector=3;
total_cell_media_city=60
erlang_media_city=total_traffic_during_the_busy_hour./(total_cell_media_citysector)
%sub_urban
calls_per_hour=3;
holding_time_in_hours=2/60
total_traffic_during_the_busy_hour=total_subscribercalls_per_hourholding_time_in_hours
total_cell_sub_city=53
erlang_sub_city=total_traffic_during_the_busy_hour./(total_cell_sub_citysector)
%rural_area
calls_per_hour=2;
holding_time_in_hours=2/60

total_traffic_during_the_busy_hour=total_subscribercalls_per_hourholding_time_in_hours
total_cell_rural_area=47
erlang_rural_area=total_traffic_during_the_busy_hour./(total_cell_rural_areasector)
BER UL
clc;close all;clear all;
N=98;
B= 5.2715e+06;
R=14.410^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
figure;
loglog((Eb_no),(pb))
hold on
N=98;
B= 5.2715e+06;
R=14.410^3;
R=R/2;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2*Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),‘r’)
legend(‘BPSK’,‘QPSK’)
xlabel(‘E_b/N_o’)
xlabel(‘P_b’)
title(‘N=98’)

N=72;
B= 3.8585e+06;
R=14.410^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
figure;
loglog((Eb_no),(pb))
hold on
N=72;
B= 3.8585e+06;
R=14.410^3;
R=R/2;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),‘r’)
legend(‘BPSK’,‘QPSK’)
xlabel(‘E_b/N_o’)
xlabel(‘P_b’)
title(‘N=72’)
N=58;
B= 3.0977e+06;
R=14.4*10^3;
G=B/R;

Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
figure;
loglog((Eb_no),(pb))
hold on
N=58;
B= 3.0977e+06;
R=14.410^3;
R=R/2
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),‘r’)
legend(‘BPSK’,‘QPSK’)
xlabel(‘E_b/N_o’)
xlabel(‘P_b’)
title(‘N=58’)
BER UL and DL
clc;close all;clear all;
%uplink
N=98;
B= 5.2715e+06;
R=14.410^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))

figure;
loglog((Eb_no),(pb),’+blue’)
hold on
N=98;
B= 5.2715e+06;
R=14.410^3;
R=R/2;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),’-r’)
%downlink
N=98;
B= 5.2715e+06;
R=25610^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),‘oblack’)
hold on
N=98;
B= 5.2715e+06;
R=25610^3;
R=R/2;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2*Eb_no))
pb=qfunc(1./sqrt(hh))

loglog((Eb_no),(pb),‘green’)
legend(‘uplink_BPSK’,‘uplink_QPSK’,‘Downlink_BPSK’,‘Downlink_QPSK’)
xlabel(‘E_b/N_o’)
xlabel(‘P_b’)
title(‘N=98’)
N=72;
B= 3.8585e+06;
R=14.410^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
figure;
loglog((Eb_no),(pb),’+blue’)
hold on
N=72;
B= 3.8585e+06;
R=14.410^3;
R=R/2;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),’-r’)
%downlink
N=72;
B= 3.8585e+06;
R=25610^3;
G=B/R;
Eb_no=0:.001:7;

hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),‘oblack’)
hold on
N=72;
B= 3.8585e+06;
R=25610^3;
R=R/2;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),'green’)
legend(‘uplink_BPSK’,‘uplink_QPSK’,‘Downlink_BPSK’,‘Downlink_QPSK’)
xlabel(‘E_b/N_o’)
xlabel(‘P_b’)
title(‘N=72’)
N=58;
B= 3.0977e+06;
R=14.410^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
figure;
loglog((Eb_no),(pb),’+blue’)
hold on
N=58;
B= 3.0977e+06;
R=14.410^3;

R=R/2
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),’-r’)
N=58;
B= 3.0977e+06;
R=25610^3;
G=B/R;
Eb_no=0:.001:7;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),‘oblack’)
N=58;
B= 3.0977e+06;
R=25610^3;
R=R/2
G=B/R;
Eb_no=0:.001:14;
hh=((N-1)./(3G))+(1./(2Eb_no))
pb=qfunc(1./sqrt(hh))
loglog((Eb_no),(pb),’*green’)
legend(‘uplink_BPSK’,‘uplink_QPSK’,‘Downlink_BPSK’,‘Downlink_QPSK’)
xlabel(‘E_b/N_o’)
xlabel(‘P_b’)
title(‘N=58’)

BER
clc;close all;clear all;
%uplink
N=98;
B= 5.2715e+06;
R=14.410^3;
G=B/R;
Eb_no=10^.7;
hh=((N-1)./(3G))+(1./(2Eb_no));
pb=qfunc(1./sqrt(hh))
N=98;
B= 5.2715e+06;
R=14.410^3;
R=R/2;
G=B/R;
Eb_no=7;
hh=((N-1)./(3G))+(1./(2Eb_no));
pb=qfunc(1./sqrt(hh))
%downlink
N=98;
B= 5.2715e+06;
R=25610^3;
G=B/R;
Eb_no=6;
hh=((N-1)./(3G))+(1./(2*Eb_no));
pb=qfunc(1./sqrt(hh))
hold on
N=98;
B= 5.2715e+06;

R=25610^3;
R=R/2;
G=B/R;
Eb_no=6;
hh=((N-1)./(3G))+(1./(2Eb_no));
pb=qfunc(1./sqrt(hh))
LTE BER
clc;close all;clear all;
snr=0:3.1623e+03;
Pb_qpsk=.5erfc(sqrt(snr/2));
Pb_16_qam=(3/4)erfc(sqrt(snr/10));
Pb_m_qam=2(1-(1./sqrt(64)))erfc(sqrt((3./(263)snr)))-(1-(2./sqrt(64))+(1/64))(erfc(sqrt((3./(2*63)*snr)))).^2;
loglog(snr,Pb_qpsk)
hold on
loglog(snr,Pb_16_qam,‘r’)
loglog(snr,Pb_m_qam,‘black’)