dnuysig=zeros(9,5); for j=1:9
for k=1:5 s1=0; s2=0;
for q=1:5 s1=s1+f4(j,q)*nuysig(j,q);
s2=s2+f4(j,q)*nuyc(j,q)*nuysig(j,q); end;
dnuysig(j,k)=-(d-y5)*muy(j)*(s1*(x3(k)*nuyc(j,k)*eps(j,k)+x3(k)*...
nuyc(j,k)*nuysig(j,k)*depnuysig(j,k))-s2*(x3(k)*eps(j,k)+x3(k)*... nuysig(j,k)*depnuysig(j,k)))/(s1^2);
nuysig42(j,k)=nuysig(j,k)+eta*dnuysig(j,k); end
end;
function nuyc41=BP41(d,y5,muy,f4,nuyc,nuysig,x3,depnuyc,eta) nuyc41=zeros(9,5);
dnuyc=zeros(9,5);
for j=1:9
for k=1:5 s1=0; s2=0;
for q=1:5 s1=s1+f4(j,q)*nuysig(j,q);
s2=s2+f4(j,q)*nuyc(j,q)*nuysig(j,q); end;
dnuyc(j,k)=-(d- y5)*muy(j)*(s1*(nuyc(j,k)*nuysig(j,k)*x3(k)*depnuyc(j,k)+...
f4(j,k)*nuysig(j,k))-s2*nuysig(j,k)*x3(k)*depnuyc(j,k))/(s1^2); nuyc41(j,k)=nuyc(j,k)+eta*dnuyc(j,k);
end;
end;
function ck3=BP31(d,y5,muy,f4,nuyc,nuysig,f3,x3,ck,sigk,depck,eps,eta) dck=zeros(1,5);
ck3=zeros(1,5); for k=1:5
s3=0;
for i=1:9
s1=0; s2=0;
for q=1:5 s1=s1+f4(i,q)*nuysig(i,q);
s2=s2+f4(i,q)*nuyc(i,q)*nuysig(i,q); end;
s3=s3+muy(i)*(nuyc(i,k)*nuysig(i,k)*s1-nuysig(i,k)*s2)*... (x3(k)*depck(i,k)+(2*(f3(k)-ck(k))/(sigk(k)^2))*...
gaussmf(f3(k),[ck(k) sigk(k)])*eps(i,k))/(s1^2);
end;
dck(k)=-(d-y5)*s3; ck3(k)=ck(k)+eta*dck(k); end;
function sig3=BP32(d,y5,muy,f4,nuyc,nuysig,f3,x3,ck,sigk,depck,eps,eta) dsigk=zeros(1,5);
sig3=zeros(1,5); for k=1:5
s3=0;
for i=1:9
s1=0; s2=0;
for q=1:5 s1=s1+f4(i,q)*nuysig(i,q);
s2=s2+f4(i,q)*nuyc(i,q)*nuysig(i,q); end;
s3=s3+muy(i)*(nuyc(i,k)*nuysig(i,k)*s1-nuysig(i,k)*s2)*... (x3(k)*depck(i,k)+(2*(f3(k)-ck(k))/(sigk(k)^2))*...
gaussmf(f3(k),[ck(k) sigk(k)])*eps(i,k))/(s1^2);
end;
dsigk(k)=-(d-y5)*s3; sig3(k)=sigk(k)+eta*dsigk(k); end;
function c21new=BP21c(d,y5,muy,f4,nuyc,nuysig,eps,f3,ck,sigk,f21,ome1,c21,sig21,eta, x)
c21new=zeros(1,5); dc21=zeros(1,5); for j=1:5
s5=0;
for i=1:9
s4=0;
for l=1:5 s1=0; s2=0; s3=0;
for q=1:5 s1=s1+f4(i,q)*nuysig(i,q);
s2=s2+f4(i,q)*nuyc(i,q)*nuysig(i,q); end;
for k=1:5
s3=s3+ome1(k)*exp(f21(j))*(2*(x-c21(j))/(sig21(j)^2)); end;
s4=s4+((nuyc(i,l)*nuysig(i,l)*s1-nuysig(i,l)*s2)*eps(i,l)*...
exp(-((f3(l)-ck(l))/sigk(l))^2)*(-2*(f3(l)- ck(l))/(sigk(l)^2))*s3)/(s1^2);
end; s5=s5+muy(i)*s4;
end;
end;
dc21(j)=-(d-y5)*s5; c21new(j)=c21(j)+eta*dc21(j);
function sig21new=BP21sig(d,y5,muy,f4,nuyc,nuysig,eps,f3,ck,sigk,f21,ome1,c21,sig21, eta,x)
sig21new=zeros(1,5); dsig21=zeros(1,5); for j=1:5
s5=0;
for i=1:9
s4=0;
for l=1:5 s1=0; s2=0; s3=0;
for q=1:5 s1=s1+f4(i,q)*nuysig(i,q);
s2=s2+f4(i,q)*nuyc(i,q)*nuysig(i,q); end;
for k=1:5
s3=s3+ome1(k)*exp(f21(j))*(((x-c21(j))^2)/(sig21(j)^3)); end;
s4=s4+((nuyc(i,l)*nuysig(i,l)*s1-nuysig(i,l)*s2)*eps(i,l)*...
exp(-(f3(l)-ck(l))/sigk(l))^2)*(-2*(f3(l)-ck(l))/(sigk(l)^2)*s3)/(s1^2); end;
s5=s5+muy(i)*s4; end;
end;
dsig21(j)=-(d-y5)*s5; sig21new(j)=sig21(j)+eta*dsig21(j);
function ome1new=BPOme1(d,y5,muy,f4,nuyc,nuysig,f3,ck,sigk,ome1,eps,eta,x1) ome1new=zeros(1,5);
dome1new=zeros(1,5); for k=1:5
s3=0;
for i=1:9
s1=0; s2=0;
for q=1:5 s1=s1+f4(i,q)*nuyc(i,q);
s2=s2+ f4(i,q)*nuyc(i,q)*nuysig(i,q); end;
s3=s3+muy(i)*(nuyc(i,k)*nuysig(i,k)*s1-nuysig(i,k)*s2)*... gaussmf(f3(k),[sigk(k) ck(k)])*eps(i,k)*(-2*(f3(k)-
ck(k)/(sigk(k)^2))*x1(k)); end;
end;
dome1new(k)=-(d-y5)*s3; ome1new(k)=ome1(k)+eta*dome1new(k);
2. Một số modul chương trình mô phỏng mạng di động tế bào và DDBNFS
function [meshnum, meshposition] = cellmesh()
Fineness = 50; % Tham so luoi cac diem vi tri k = 1;
j = 0;
ds = sqrt(3.0) / Fineness; xmesh = -0.5 * sqrt(3.0); while xmesh < 0.5 * sqrt(3.0)
xmesh = (k - 1) * ds - 0.5 * sqrt(3.0); if xmesh < 0.0
ymin = -xmesh/sqrt(3.0) - ds - 1.0; ymax = xmesh/sqrt(3.0) + 1.0;
elseif xmesh == 0.0 ymin = -1.0;
ymax = 1.0; else
ymin = xmesh/sqrt(3.0) - ds -1.0; ymax = -xmesh/sqrt(3.0) + 1.0;
end
k= k + 1; ymesh = ymin;
while ymesh < ymax ymesh = ymesh + ds; if ymesh > ymax
break
end
end
end
j = j + 1; posi(j,1) = xmesh; posi(j,2) = ymesh;
meshnum = j; meshposition = posi;
% Các cell nhom compact N=19 function [out] = basest()
1) | = | 0.0; | |
baseinfo(1, | 2) | = | 0.0; |
baseinfo(2, | 1) | = | -0.5*sqrt(3.0); |
baseinfo(2, | 2) | = | 1.5; |
baseinfo(3, | 1) | = | -sqrt(3.0); |
baseinfo(3, | 2) | = | 0.0; |
baseinfo(4, | 1) | = | -0.5*sqrt(3.0); |
baseinfo(4, | 2) | = | -1.5; |
baseinfo(5, | 1) | = | 0.5*sqrt(3.0); |
baseinfo(5, | 2) | = | -1.5; |
baseinfo(6, | 1) | = | sqrt(3.0); |
baseinfo(6, | 2) | = | 0.0; |
baseinfo(7, | 1) | = | 0.5*sqrt(3.0); |
baseinfo(7, | 2) | = | 1.5; |
baseinfo(8, | 1) | = | 0.0; |
baseinfo(8, | 2) | = | 3.0; |
baseinfo(9, | 1) | = | -sqrt(3.0); |
baseinfo(9, | 2) | = | 3.0; |
Có thể bạn quan tâm!
- S. W. Halpern (1983), “Reuse Partitioning In Cellular Systems”, Proceeding Of The 1983 33Rd Ieee Vehicular Technology Conference, Pp. 322-327, New York, Ny, Usa.
- Alexe E. Leu And Brian L. Mark (2002), “Modeling And Analysis Of Fast Handoff Algorithms For Microcellular Networks”, Proceeding Of The 10Th Ieee Mascots’2002, Pp. 321-328.
- Về một phương pháp mới điều khiển mượn, khóa kênh tần số mạng di động tế bào - 17
Xem toàn bộ 149 trang tài liệu này.
1) | = | -1.5*sqrt(3.0); | |
baseinfo(10, | 2) | = | 1.5; |
baseinfo(11, | 1) | = | -2.0*sqrt(3.0); |
baseinfo(11, | 2) | = | 0.0; |
baseinfo(12, | 1) | = | -1.5*sqrt(3.0); |
baseinfo(12, | 2) | = | -1.5; |
baseinfo(13, | 1) | = | -sqrt(3.0); |
baseinfo(13, | 2) | = | -3.0; |
baseinfo(14, | 1) | = | 0.0; |
baseinfo(14, | 2) | = | -3.0; |
baseinfo(15, | 1) | = | sqrt(3.0); |
baseinfo(15, | 2) | = | -3.0; |
baseinfo(16, | 1) | = | 1.5*sqrt(3.0); |
baseinfo(16, | 2) | = | -1.5; |
baseinfo(17, | 1) | = | 2.0*sqrt(3.0); |
baseinfo(17, | 2) | = | 0.0; |
baseinfo(18, | 1) | = | 1.5*sqrt(3.0); |
baseinfo(18, | 2) | = | 1.5; |
baseinfo(19, | 1) | = | sqrt(3.0); |
baseinfo(19, | 2) | = | 3.0; |
out = baseinfo;
%Nhom compact N=19
function [out] = wrap(inputmat)
= = | 1:19; [ 2, 9,10, 3, 1, 7, | 8,14,13,17,16,11,12, 4, 5, | |||||
6,18,19,15]; inputmat( 3,1:19) | = | [ 3,10,11,12, 4, 1, | 2, 9,17,16,15,19,18,13,14, 5, 6, 7, | ||||
8]; inputmat( 4,1:19) | = | [ 4, 3,12,13,14, 5, | 1, 2,10,11,19,18,17, 9, 8,15,16, 6, | ||||
7]; inputmat( 5,1:19) | = | [ 5, 1, 4,14,15,16, | 6, 7, 2, 3,12,13, 9, | ||||
8,19,11,10,17,18]; inputmat( 6,1:19) | = | [ | 6, | 7, | 1, | 5,16,17,18,19, 8, | 2, 3, 4,14,15,11,10, |
9,13,12]; inputmat( 7,1:19) | = | [ | 7, | 8, | 2, | 1, 6,18,19,15,14, | 9,10, 3, 4, |
5,16,17,13,12,11]; |
inputmat( 8,1:19) = [ 8,14, 9, 2, 7,19,15, 5, 4,13,17,10, 3, 1,
6,18,12,11,16];
inputmat( 9,1:19) = [ 9,13,17,10, 2, 8,14, 4,12,18, 6,16,11, 3, 1, 7,19,15,
5];
= = | [10,17,16,11, 3, 2, 9,13,18, 6, 5,15,19,12, 4, 1, 7, [11,16,15,19,12, 3,10,17, 6, 5,14, 8, 7,18,13, 4, 1, | 2, | |
9]; | |||
inputmat(12,1:19) | = | [12,11,19,18,13, 4, 3,10,16,15, 8, 7, 6,17, 9,14, 5, | 1, |
2]; | |||
inputmat(13,1:19) | = | [13,12,18,17, 9,14, 4, 3,11,19, 7, 6,16,10, 2, 8,15, | 5, |
1]; |
= | [14, | 4,13, 9, 8,15, 5, 1, | 3,12,18,17,10, 2, 7,19,11,16, | ||||
inputmat(15,1:19) 3,10,17]; inputmat(16,1:19) | = = | [15, [16, | 5,14, 8,19,11,16, 6, 6, 5,15,11,10,17,18, | 1, 4,13, 9, 2, 7,18,12, 7, 1, 4,14, 8,19,12, 3, 2, | |||
9,13]; | |||||||
inputmat(17,1:19) 4]; inputmat(18,1:19) | = = | [17,18, [18,19, | 6,16,10, 9,13,12,19, 7, 6,17,13,12,11,15, | 7, 8, | 1, 2, | 5,15,11, 3, 1, 5,16,10, | 2, 8,14, 9,14, 4, |
3]; |
inputmat(19,1:19) = [19,15, 8, 7,18,12,11,16, 5,14, 9, 2, 1, 6,17,13, 4,
3,10];
out = inputmat;