Cellular Codes in Matlab
Cellular Codes in Matlab
clc;
clear all;
close all;
bw=33000;
sim_ch_bw =25;
disp('channel Bandwidth..');
dup_ch_bw=2*sim_ch_bw;
t_ch=(bw/dup_ch_bw);
disp(dup_ch_bw);
disp('Total available channel');
disp(t_ch);
cc_bw=1000;
t_cc=cc_bw/dup_ch_bw;
disp('Total control channel');
disp(t_cc);
for N=[4 7 12]
%n=size(N);
%for i=1:1:n
ch=(t_ch/N);
ch_per_cell=round(ch);
disp('ch per cell');
disp(N);
disp(ch_per_cell);
c=(t_cc/N);
cc=round(c);
vc=(ch_per_cell-cc);
disp('cc & vc are....');
%disp(N);
disp(cc);
disp(vc);
end
OUTPUT:
channel Bandwidth.. 50
Total available channel 660
Total control channel 20
ch per cell 4 165
cc & vc are.... 5 160
ch per cell 7 94
cc & vc are.... 3 91
ch per cell 12 55
cc & vc are....2 53
Input:
r_si=15;
io=6;
for n=[4 3]
N=7
q=sqrt(3*N);
disp('n');
disp(n);
disp('q');
disp(q);
si=10*(log10((q^n)/io));
disp('si');
disp(si);
if(si<r_si)
i=2;j=2;
N1=(i*i)+(i*j)+(j*j);
q=sqrt(3*N1);
disp('n');
disp('q');
disp(q);
si1=10*(log10((q^n)/io));
disp('si1');
disp(si1);
end
end
Output:
Q 4.5826
si 18.6629
N = 7
N 3
Q 4.5826
Si 12.0518
N
Q 6
si1
15.5630
Input:
clc;
close all;
clear all;
gos=0.5/100;
au=0.1;
a=[0.005 1.13 3.96 11.10 80.9];
c=[1 5 10 20 100];
disp('blocking probability');
disp(gos);
disp('traffic intensity per user');
disp(au);
disp('traffic intensity');
disp(a);
disp('channel');
disp(c);
U=(a/au);
u=round(U);
disp('number of users');
disp(u);
Output:
blocking probability
0.0050
traffic intensity
0.0050 1.1300 3.9600 11.1000 80.9000
channel
1 5 10 20 100
number of users
0 11 40 111 809
Input:
clc;
clear all;
blocking_probability=2/100;
au=(2/60)*3;
disp('For system A');
number_A_channel_per_cell=19;
a=12
u=a/au;
A=u*394;
percentage_market_penetration_for_A=(A/2000000)*100;
disp('Total number of subscriber of A= ');
disp(A);
disp('percentage_market_penetration_for_A');
disp(percentage_market_penetration_for_A);
disp('For system B');
number_A_channel_per_cell=57;
a=45
u=a/au;
B=u*98;
percentage_market_penetration_for_B=(B/2000000)*100;
disp('Total number of subscriber of B= ');
disp(B);
disp('percentage_market_penetration_for_B');
disp(percentage_market_penetration_for_B);
disp('For system C');
number_A_channel_per_cell=100;
a=88
u=a/au;
C=u*49;
percentage_market_penetration_for_C=(C/2000000)*100;
disp('Total number of subscriber of C= ');
disp(C);
disp('percentage_market_penetration_for_C');
disp(percentage_market_penetration_for_C);
t=A+B+C;
disp('Total number of subscriber supported by A,B,C= ');
disp(t);
percentage_market_penetration_for_three_system=(t/2000000)*100;
disp('percentage_market_penetration_for_three_system');
disp(percentage_market_penetration_for_three_system);
Output:
For system A
a =
12
Total number of subscriber of A=
47280
percentage_market_penetration_for_A
2.3640
For system B
a =
45
percentage_market_penetration_for_B
2.2050
For system C
a =
88
percentage_market_penetration_for_C
2.1560
percentage_market_penetration_for_three_system 6.7250
Input:
clc;
clear all;
close all;
total_city_coverage_area=1300;
radius=4;
each_cell_cover=floor(2.5981*radius^2);
disp('(a)')
number_of_cells=floor(total_city_coverage_area/each_cell_cover)
disp('(b)')
number_of_channel_per_cell=floor(40000/(60*7))
disp('(c)')
traffic_intensity_per_cell=84
disp('(d)')
maximum_carried_traffic= number_of_cells*traffic_intensity_per_cell
disp('(e)')
total_number_of_user=maximum_carried_traffic/0.03
disp('(f)')
number_of_mobile_per_channel=floor(total_number_of_user/
(number_of_channel_per_cell*7))
disp('(g)')
theorital_maximum_number_of_User_that_could_be_served=number_of_channe
l_per_cell*number_of_cells
Output:
a) number_of_cells =
31
(b)number_of_channel_per_cell =
95
(c)traffic_intensity_per_cell =
84
(d)maximum_carried_traffic =
2604
(e)total_number_of_user =
86800
(f)number_of_mobile_per_channel =
130
(g)theorital_maximum_number_of_User_that_could_be_served = 2945
Input:
clc;
clear all;
close all;
r=1.387;
n=4;
N=60;
A=2.598*r^2;
c=N/4;
disp('(a)')
traffic_intensity =9;
number_of_user=floor(traffic_intensity/(0.029*A))
disp('(b)')
lemda=1;
h=(0.029/lemda)*60*60;
the_probability_to_wait= (exp(-(c-traffic_intensity)*10)/h)*100
disp('(c)')
the_probability_of_delay=0.05*the_probability_to_wait
Output:
(a)
number_of_user =
62
(b)
the_probability_to_wait =
8.3875e-27
(c)
the_probability_of_delay =
4.1937e-28
Input:
clc;
clear all;
close all;
pt=50;
f=900;
gt=1;
gr=1;
d=100;
disp('(a)')
transmitter_power_in_dbm =ceil(10*log10(50*1000))
disp('(b)')
transmitter_power_in_dBw =ceil(10*log10(50*1))
pr_mw=((pt*gt*gr*(300/900)^2)/((4*3.1416)^2*d^2*1)*1000)
received_power_in_dBm =10*log10(pr_mw)
pr_10km = received_power_in_dBm + (20*log10(100/10000))
Output:
(a)transmitter_power_in_dbm =
47
(b)
transmitter_power_in_dBw =
17
pr_mw =
0.0035
received_power_in_dBm =
-24.5369
pr_10km = -64.5369
Input:
clc;
clear all;
close all;
t_r_field=5;
E_field=10^-3;
d =5*10^3;
f=900;
lemda=300/900;
disp('(a)')
length_of_antenna= lemda/4
gain= (10^(2.55/10))
disp('(b)')
er_d=(2*E_field*1000*2*3.1416*50*1.5)/(lemda*d^2)
ae=(gain*lemda^2)/(4*3.1416);
pr_d=(er_d^2/377)*ae
received_power_at_5_km_distance= 10*log10(pr_d)
OutPut:
(a)
length_of_antenna =
0.0833
gain =
1.7989
(b)
er_d =
1.1310e+08
pr_d = 5.3965e+11
received_power_at_5_km_distance =
-122.6789
Input:
clc;
clear all;
close all;
hm=2;
hb=100;
f=900;
d=4;
a=(3.2*(log10(11.75*hm))^2)-4.97
l=69.55+26.16*log10(f)-13.82*log10(hb)-a+(44.9-
6.55*log10(hb))*log10(d);
disp('loss path')
disp(l)
Output:
a =
1.0454
loss path
137.2930
Input:
clc;
clear all;
close all;
f=1.8;
hb=20;
d=((20^2+30^2)^0.5)/1000;
path_loss=135.41+(12.49*log10(f))-(4.99*log10(hb))+((46.84-
2.34*log10(hb))*log10(d))
Output:
path_loss =
68.9079