话务工程程序清单.doc

%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