话务工程程序清单.docVIP

%erlangb.m %Erlang-B公式,计算阻塞概率 %调用语句为:erlangb(业务量,中继线数(整数)) function erlangb(a,n) global b k=0; b=1; b1=(a*b)/(a*b+k); while k=n b=b1; b1=(a*b)/(a*b+k); k=k+1; end b=b1 fprintf(a=%d,n=%d,b=%d\n,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; end b=1/b fprintf(a=%d,n=%d,b=%d\n,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); end b fprintf(a=%d,x=%d,b=%d\n,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); end b fprintf(a=%d,x=%d,b=%d\n,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; end b1=(n+1)/a*b+1; b2=(n+2)/a*b1+1; y=b^(1-h)*b1^h; y=y*(b1*b1/b/b2)^(h*(1-h)/2); b=1/y m0=a*b; z0=1-m0+a/(x+1+m0-a); m1=fix(100000*m0+0.5)/100000 z1=fix(100000*z0+0.5)/100000 c=fix(100000*b+0.5)/100000 fprintf(a=%d,x=%d,b=%d\n,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/f1 end a=(100000*a+0.5)/100000 fprintf(b=%d,n=%d,a=%d\n,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)*b1^h; y=y*(b1*b1/b/b2)^(h*(1-h)/2); g=(x+1)/a+1/y; x=x+(1-p*y)/log(g); end c=fix(100000*x+0.5)/100000 fprintf(a=%d,p=%d,c=%d\n,a,p,c) %Engset公式,计算阻塞概率(方法1) %调用语句为:p8(每用户的业务量,中继线数(整数),用户数) %或在设定全局变量b的前提下p8(a,n,m,b) function b=p8(a,n,m) b=1;