%Functions
%cis.m-main function
%place_ms.m-for placing MS called from main function
%place_ms_in.m-function called from place_ms.m for internal calculation
%sinr1.m- function for calculation of SINR of FRF(7,3,1)
%sinr2.m- function for calculation of SINR of FRF(7,4,2)
%sinr11.m- function for calculation of SINR of FRF(11,5,2)
%sinr13.m- function for calculation of SINR of FRF(13,4,1)
%frf1.m-function for calculation of SINR of FRF 1 and FRF 3
%cell_throghput.m-function for calculation of different FRF's
clear all
clc
N=6;%No. of cell in neighbourhood of a cell
n=100;% No. of mobile station per cell
r=1500;% radius of cell
a=sqrt(3);
h(1)=5000;%X co-ordinate of central cell
m=1000;
k(1)=5000;%Y co-ordinate of central cell
% Calculation of centers of neighbour cells
for i=1:N
h(i+1)=h(1)+a*r*cos((2*i-1)*pi/6);
k(i+1)=k(1)+a*r*sin((2*i-1)*pi/6);
end
%Calculation of centers of 2nd tier cells
h(8)=h(1)+3*r;
k(8)=k(1);
i=8;
for j=1:N
i=i+1;
h(i)=h(1)+2*a*r*cos((2*j-1)*pi/6);
k(i)=k(1)+2*a*r*sin((2*j-1)*pi/6);
i=i+1;
h(i)=h(1)+3*r*cos(j*pi/3);
k(i)=k(1)+3*r*sin(j*pi/3);
end
%Calculation of centers of 3rd tier cells
for j=1:N
i=i+1;
h(i)=h(1)+sqrt(21)*r*cos(0.19+(j-1)*pi/3);
k(i)=k(1)+sqrt(21)*r*sin(0.19+(j-1)*pi/3);
i=i+1;
h(i)=h(1)+3*a*r*cos(pi/6+(j-1)*pi/3);
k(i)=k(1)+3*a*r*sin(pi/6+(j-1)*pi/3);
i=i+1;
h(i)=h(1)+sqrt(21)*r*cos(0.857+(j-1)*pi/3);
k(i)=k(1)+sqrt(21)*r*sin(0.857+(j-1)*pi/3);
end
%calculation of Averaged sinr for m times
[p p1]=place_ms(h,k,n,N,r);
%for 13,4,1
frf=13/4;
c1=[h(1),k(1);h(3),k(3);h(5),k(5);h(7),k(7);h(8),k(8);h(11),k(11);h(14),k(14)];
c2=[h(1),k(1);h(6),k(6);h(17),k(17);h(4),k(4);h(15),k(15)];
c3=[h(1),k(1);h(16),k(16);h(10),k(10);h(2),k(2);h(19),k(19)];
c4=[h(1),k(1);h(9),k(9);h(12),k(12);h(13),k(13);h(18),k(18)];
[y13 mean13 sd13]=sinr13(c1,c2,c3,c4,p1,n);
[ce13 p0]=cell_throghput(y13,n,frf);
cdfs13=cdf('normal',y13,mean13,sd13);
outageprobablityoffrf13=p0
%for FRF (11,5,2)
%frf=11/5;
%c1=[h(1),k(1);h(18),k(18);h(15),k(15);h(5),k(5);h(19),k(19);h(3),k(3);h(4),k(4);h(8),k(8);h(20),k(20);h(21),k(21);h(24),k(24);h(28),k(28);h(27),k(27);h(28),k(28);h(32),k(32);h(35),k(35);h(37),k(37)];
%c2=[h(1),k(1);h(11),k(11);h(16),k(16);h(12),k(12);h(9),k(9);h(15),k(15);h(18),k(18);h(7),k(7);h(14),k(14);h(37),k(37);h(28),k(28);h(26),k(26);h(30),k(30);h(31),k(31);h(34),k(34);h(21),k(21);h(35),k(35);h(20),k(20);h(27),k(27);h(32),k(32)];
%c3=[h(1),k(1);h(4),k(4);h(6),k(6);h(8),k(8);h(12),k(12);h(9),k(9);h(17),k(17);h(10),k(10);h(37),k(37);h(28),k(28);h(23),k(23);h(29),k(29);h(36),k(36);h(33),k(33);h(20),k(20);h(27),k(27);h(32),k(32)];
%c4=[h(1),k(1);h(2),k(2);h(3),k(3);h(6),k(6);h(11),k(11);h(13),k(13);h(16),k(16);h(37),k(37);h(28),k(28);h(23),k(23);h(29),k(29);h(36),k(36);h(33),k(33);h(20),k(20);h(27),k(27);h(32),k(32)];
%c5=[h(1),k(1);h(2),k(2);h(5),k(5);h(7),k(7);h(10),k(10);h(13),k(13);h(14),k(14);h(17),k(17);h(19),k(19);h(37),k(37);h(28),k(28);h(23),k(23);h(29),k(29);h(36),k(36);h(22),k(22);h(24),k(24)];
%[y111 mean111 sd111]=sinr111(c1,c2,c3,c4,c5,p1,n);
%ce111=cell_throghput(y111,n,frf);
%cdfs111=cdf('normal',y111,mean111,sd111);
%for FRF (11,5,2)
frf=11/5;
c1=[h(1),k(1);h(18),k(18);h(15),k(15);h(5),k(5);h(19),k(19);h(3),k(3);h(4),k(4);h(8),k(8)];
c2=[h(1),k(1);h(11),k(11);h(16),k(16);h(12),k(12);h(9),k(9);h(15),k(15);h(18),k(18);h(7),k(7);h(14),k(14)];
c3=[h(1),k(1);h(4),k(4);h(6),k(6);h(8),k(8);h(12),k(12);h(9),k(9);h(17),k(17);h(10),k(10)];
c4=[h(1),k(1);h(2),k(2);h(3),k(3);h(6),k(6);h(11),k(11);h(13),k(13);h(16),k(16)];
c5=[h(1),k(1);h(2),k(2);h(5),k(5);h(7),k(7);h(10),k(10);h(13),k(13);h(14),k(14);h(17),k(17);h(19),k(19)];
[y11 mean11 sd11]=sinr11(c1,c2,c3,c4,c5,p1,n);
[ce11 p0]=cell_throghput(y11,n,frf);
cdfs11=cdf('normal',y11,mean11,sd11);
outageprobablityoffrf11=p0
% for FRF 7/3
n1=7;
%Cells having subband 1
c1=[h(1),k(1);h(4),k(4);h(6),k(6);h(8),k(8);h(9),k(9);h(12),k(12);h(17),k(17)];
%Cells having subband 2
c2=[h(1),k(1);h(3),k(3);h(5),k(5);h(10),k(10);h(15),k(15);h(18),k(18);h(19),k(19)];
%Cells having subband 3
c3=[h(1),k(1);h(2),k(2);h(7),k(7);h(13),k(13);h(11),k(11);h(14),k(14);h(16),k(16)];
[y73 mean1 sd1]=sinr1(c1,c2,c3,p1,n,n1);%Average SINR for each M.S.
frf=7/3;
[ce73 p0]=cell_throghput(y73,n,frf);
cdfs73=cdf('normal',y73,mean1,sd1);%cdf of sinr for FRF 7/3
outageprobablityoffrf73=p0
%for FRF 7/4
c1=[h(1),k(1);h(4),k(4);h(5),k(5);h(7),k(7);h(8),k(8);h(10),k(10);h(11),k(11);h(14),k(14);h(17),k(17);h(19),k(19)];
c2=[h(1),k(1);h(3),k(3);h(4),k(4);h(6),k(6);h(8),k(8);h(9),k(9);h(12),k(12);h(15),k(15);h(17),k(17);h(18),k(18)];
c3=[h(1),k(1);h(2),k(2);h(3),k(3);h(5),k(5);h(10),k(10);h(13),k(13);h(15),k(15);h(16),k(16);h(18),k(18);h(19),k(19)];
c4=[h(1),k(1);h(2),k(2);h(6),k(6);h(7),k(7);h(9),k(9);h(11),k(11);h(12),k(12);h(13),k(13);h(14),k(14);h(16),k(16)];
n1=10;
[y74 mean4 sd4]=sinr2(c1,c2,c3,c4,p1,n,n1);
frf=7/4;
[ce74 p0]=cell_throghput(y74,n,frf);
cdfs74=cdf('normal',y74,mean4,sd4);%cdf of sinr for FRF 7/4
outageprobablityoffrf74=p0
%for FRF 1
n1=19;
[y1 mean1 sd1]=frf1(h,k,p1,n,n1);
frf=1;
[ce1 p0]=cell_throghput(y1,n,frf);
cdfs1=cdf('normal',y1,mean1,sd1);%cdf of sinr for FRF 1
outageprobablityoffrf1=p0
%for FRF 3
c1=[h(1),h(9),h(11),h(13),h(15),h(17),h(19)];
c2=[k(1),k(9),k(11),k(13),k(15),k(17),k(19)];
n1=7;
[y3 mean3 sd3]=frf1(c1,c2,p1,n,n1);
frf=3;
[ce3 p0]=cell_throghput(y3,n,frf);
cdfs3=cdf('normal',y3,mean3,sd3);%cdf of sinr for FRF 3
outageprobablityoffrf3=p0
%plot different figures
ce=[ce1,ce3,ce73,ce74,ce11,ce13];
x={'frf1','frf3','frf7/3','frf7/4','frf11/5','frf13/4'};
%MS in different cellular location
figure(1)
plot(h,k,'*',p(:,1),p(:,2),'.r')
hold on
for j=1:38
t=(0:N)*2*pi/N;
plot( r*cos(t)+h(j), r*sin(t)+k(j));
end
xlabel('Cell site')
hold off
% CDF of Different FRF schemes
figure(2)
plot(y1,cdfs1,'*-g',y3,cdfs3,'*-r',y74,cdfs74,'s-b',y11,cdfs11,'*-y',y13,cdfs13,'s-g');
grid
title('Different FRF Schemes')
xlabel('SINR (dB)')
ylabel('Cumulative density function')
hold off
%CDF of Integer FRF's
figure(3)
plot(y11,cdfs11,'*-r',y74,cdfs74,'.-g',y73,cdfs73,'s-b')
title('Integer FRF')
grid
xlabel('SINR (dB)')
ylabel('Cumulative density function')
hold off
%CDF of Non-Integer FRF's
%figure(4)
%plot(y73,cdfs73,'.-y',y74,cdfs74,'s-m')
%title('Non-Integer FRF')
%grid
%xlabel('SINR (dB)')
%ylabel('Cumulative density function')
%hold off
% Plotting of Cell -Throughput of Different Schemes
figure(5)
pareto(ce,x)
ylabel('Cell Throughput (kbps)')
axis([0 8 0 5000])
评论0