%erlangb.m%Erlang-B公式,計算阻塞概率%調用語句為:erlangb(業務量,中繼線數(整數)function erlangb(a,n)global bk=0;b=1;b1=(a*b)/(a*b+k);while k=n b=b1; b1=(a*b)/(a*b+k); k=k+1;endb=b1fprintf(a=%d,n=%d,b=%dn,a,n,b)function b=p2(a,n)%Erlang-B公式,計算阻塞概率%調用語句為:p2(業務量,中繼線數(整數)t=1;b=t;k=1;t=t*(n-k+1)/a;b=b+t;while kn k=k+1; t=t*(n-k+1)/a; b=b+t;endb=1/bfprintf(a=%d,n=%d,b=%dn,a,n,b)%Erlang-B公式,計算阻塞概率%調用語句為:p3(業務量,中繼線數(實數)function b=p3(a,x)n=fix(x);x0=x-n;c0=1;c1=-(2+a)/(1+3*a+a*a);c2=1/(1+4*a+4*a*a+a*a*a);b=c0+c1*x0+c2*x0*x0;for k=1:n b=a*b/(a*b+x0+k);endbfprintf(a=%d,x=%d,b=%dn,a,x,b)%Erlang-B公式,計算阻塞概率%調用語句為:p4(業務量,中繼線數(實數)function b=p4(a,x)n=fix(x);x0=x-n;b=(2-x0)*a+a*a)/(x0+2*a+a*a);for k=1:n b=a*b/(a*b+x0+k);endbfprintf(a=%d,x=%d,b=%dn,a,x,b)%Erlang-B公式,計算阻塞概率%調用語句為:p5(業務量,中繼線數(實數),損失率)function p5(a,x,b);n=fix(x);h=x-n;b=1;for j=1:n b=j/a*b+1;endb1=(n+1)/a*b+1;b2=(n+2)/a*b1+1;y=b(1-h)*b1h;y=y*(b1*b1/b/b2)(h*(1-h)/2);b=1/ym0=a*b;z0=1-m0+a/(x+1+m0-a);m1=fix(100000*m0+0.5)/100000z1=fix(100000*z0+0.5)/100000c=fix(100000*b+0.5)/100000fprintf(a=%d,x=%d,b=%dn,a,x,b)fprintf(溢出業務量=%d,迂回路由占線數的方差與溢出業務比=%d,n,m1,z1)%Erlang-B公式,由阻塞概率,中繼線數計算業務量 %調用語句為:p6(阻塞概率,中繼線數(整數)或在定義全局變量a的情況下p6(b,n,a)function p6(b,n,a)a=n/(1-b)-1;for k=1:5 p=1; for j=1:n p=j/a*p+1; end f=1-b*p; f1=n/a-1+1/p; a=a-f/f1enda=(100000*a+0.5)/100000fprintf(b=%d,n=%d,a=%dn,b,n,a) %Erlang-B公式,由業務量,阻塞概率計算中繼線數 %調用語句為:p7(業務量,阻塞概率) 或在定義全局變量c的情況下p7(a,p,c)function p7(a,p,c)x=(1-p)*(1+a);for k=1:10 n=fix(x); h=x-n; b=1; for j=1:n b=j/a*b+1; end b1=(n+1)/a*b+1; b2=(n+2)/a*b1+1; y=b(1-h)*b1h; y=y*(b1*b1/b/b2)(h*(1-h)/2); g=(x+1)/a+1/y; x=x+(1-p*y)/log(g);endc=fix(100000*x+0.5)/100000fprintf(a=%d,p=%d,c=%dn,a,p,c) %Engset公式,計算阻塞概率(方法1)%調用語句為:p8(每用戶的業務量,中繼線數(整數),用戶數)%或在設定全局變量b的前提下p8(a,n,m,b)function b=p8(a,n,m)b=1;k=1;b=(m-k)*a*b/(k+(m-k)*a*b);while k=n b=(m-k)*a*b/(k+(m-k)*a*b); k=k+1;endbfprintf(a=%d,n=%d,m=%d,b=%dn,a,n,m,b) %Engset公式,計算阻塞概率(方法2)%調用語句為:p9(每用戶的業務量,中繼線數(整數),用戶數)function b=p9(a,n,m)t=1;k=1;t=1+t*k/(m-k)*a);while k=n t=1+t*k/(m-k)*a); k=k+1;endb=1/tfprintf(a=%d,n=%d,m=%d,b=%dn,a,n,m,b) %Engset公式,計算阻塞概率(方法3)%調用語句為:p10(每用戶的業務量,中繼線數(整數),用戶數)%或在設定全局變量b的前提下p10(a,n,m,b)function b=p10(a,n,m)t=1;b=t;k=1;t=t*(n-k+1)/(m-n+k-1)*a);b=b+t;while kb1)&(bb(k+1)b1)&(bb(k+1)b1)&(bb(k+1)0,均等待時間function erlangc(lam,mu,n)r=lam/mu;a=r;k=0;b=1;b1