ErlandB公式实验报告.docVIP

《通信网性能分析基础》 实 验 报 告 实验二 Erland B公式计算器 实验目的 Erland B公式表达了电话交换系统的时间阻塞率B(s,a)与到达的呼叫量a以及中继线数目s之间的关系，为电话网络的规划和中继线容量配置奠定了基础。 本次实验要求用计算机语言编写Erland B公式计算器，实现给定任意两个变量求解第三个变量的功能，使学生掌握Erland B公式的计算方法，培养数值分析计算与算法设计能力。 实验原理 Erlang即时拒绝系统用于建模电话交换系统。设电话呼叫流的到来服从Poisson过程，λ为呼叫的到达率，呼叫的持续时间服从参数μ的负指数分布。系统有s条中继线，呼叫到来时，如系统中有空闲的中继线，则呼叫可到达任意一个空闲的中继线，如果没有空闲的中继线，就拒绝该呼叫。在这样的情况下，该交换系统的排队系统模型为M/M/S(s)。 利用生灭过程稳态分布的结论，可求解M/M/s(s)系统的时间阻塞率为。该公式给出了B，s，a三者之间的关系。 实验要求 本次实验要求编写Erlang B公式的计算器，实现给定任意两个变量求解第三个变量的功能： 　　(1) 给定到达的呼叫量a和中继线的数目s，求解系统的时间阻塞率B； 　　(2) 给定系统的时间阻塞率的要求B和到达的呼叫量a，求解中继线的数目s，以实现网络规划； (3) 给定系统的时间阻塞率要求B以及中继线的数目s，判断该系统能支持的最大的呼叫量a。 实验内容 采用的语言 本次实验采用的软件是VC6.0，所采用的语言是C++。 数据结构 由于本次实验算法与函数间的关系较为简单，只运用了普通的double型、int型等简单变量，并没有设计指针、树、链表等较为复杂的数据关系。 主要的函数及算法流程图 1）已知a和s，求B的函数 void CErland_BDlg::OnBUTTONcalB() { int i; UpdateData(TRUE); double B=1; for(i=0;i=m_s1;i++) { B=(m_a1*B)/((double)i+m_a1*B); } m_B1=B;//求出B(s,a) UpdateData(FALSE); } 在该函数中，我们并没有运用Erland B公式，而是推出另一个与Erland B公式相等的递推公式，利用这种算法，可使程序的实现较为简便。下面是程序的流程图： 已知B和s，求a的函数 void CErland_BDlg::OnBUTTONcala() { UpdateData(TRUE); int i; double a,equ,B; a=0;equ=1;B=1; while ((equ0.0001)||(equ-0.0001)) { a=a+0.001; for(i=0;i=m_s2;i++) { B=B*a/(i+a*B); } equ=m_B2-B; } m_a2=a;//求出a UpdateData(FALSE); } 在该函数中，我们运用了穷举的算法。先设定一个从0开始的a值，利用该a值带入公式中，算出B(s,a)值并与输入的B值比较，若二者之差不在精度范围（本程序默认精度为0.0001）内，则增加a的值，直至算出的B(s,a)与实际B值之差相差在精度范围内。此时则默认该a值为所求。具体实现流程图如下： 已知B和a，求s的函数 void CErland_BDlg::OnBUTTONcals() { UpdateData(TRUE); int i,s; double B,equ; s=0;B=1;equ=1; while ((equ0.01)||(equ-0.01)) { s=s+1; for(i=0;i=s;i++) { B=m_a3*B/(i+m_a3*B); } equ=B-m_B3; } m_s3=s;//求出s UpdateData(FALSE); } 在本函数中，所采用的算法的基本思想与上述函数类似，同样也是采取穷举的方法。具体实现流程图如下： 清零 void CErland_BDlg::OnBUTTONzero() { UpdateData(TRUE); m_a1=0;m_s1=0;m_B1=0; m_a2=0;m_s2