北邮 通信网实验报告
北邮通信原理软件实验报告
北邮通信原理软件实验报告北邮通信原理软件实验报告一、实验目的本次实验旨在加深对通信原理知识的理解，并通过实际操作掌握常用通信原理技术。
二、实验内容及原理本次实验分为三项实验内容：1、FSK解调实验FSK是一种通信调制技术，在数字信号传输领域中比较常用。
FSK解调实验中，我们使用MATLAB软件编写程序，模拟FSK解调过程，了解解调过程中的基本原理。
FSK是通过改变载波频率来传输信息的调制技术。
在数字信号的传输中，我们一般将数字信号分为两种，0和1，然后分别将它们对应到两个频率上，再将这两个频率进行交错发送，接收方通过检测频率的变化来判断发送方的信息。
在FSK解调实验中，我们使用的解调技术是匹配滤波器法。
解调的过程是将接收到的信号经过低通滤波器，合并成一个信号。
2、QAM解调实验QAM是一种把两路模拟信号叠加的数字调制技术，它是组合了ASK和PSK的数字传输技术。
QAM解调实验中，通过MATLAB软件仿真的方法，模拟QAM解调过程，了解解调过程中的基本原理。
QAM技术是将两路数模合成的模拟信号进行数字化处理，将两路模拟信号进行分别调制成两个独立的数字信号，然后将这两个数字信号通过载波同步合成一个数字信号进行传输。
在QAM解调实验中，我们使用的解调技术是相干解调。
解调的过程是将接收到的信号经过相干解调器解调，得到原始的时域信号，然后通过低通滤波器进行滤波。
3、OFDM调试实验OFDM技术是目前广泛应用于高速数据传输的一种技术，它是通过将信号分成多个子载波进行传输，提高频率利用率，并实现抗多径衰落的效果。
OFDM调试实验中，我们通过软件界面和Matlab代码相结合，模拟OFDM调制和解调过程，了解其中的基本原理。
OFDM技术是通过将原始信号分成多个子信道，每个子信道独立传输，最终将其合并成整个信号。
因此，在OFDM模式下，每个子信道的公共频率就成为可利用的带宽，提高了传输率并减少了所需的带宽。
在OFDM调试实验中，我们使用了MATLAB软件进行调制和解调。
北邮 通信网实验报告
北京邮电大学实验报告通信网理论基础实验报告学院:信息与通信工程学院班级:2013211124学号:姓名:实验一 ErlangB公式计算器一实验内容编写Erlang B公式的图形界面计算器，实现给定任意两个变量求解第三个变量的功能：1）给定到达的呼叫量a和中继线的数目s，求解系统的时间阻塞率B；2）给定系统的时间阻塞率的要求B和到达的呼叫量a，求解中继线的数目s，以实现网络规划；3）给定系统的时间阻塞率要求B以及中继线的数目s，判断该系统能支持的最大的呼叫量a。
二实验描述1 实验思路使用MATLAB GUITOOL设计图形界面，通过单选按钮确定计算的变量，同时通过可编辑文本框输入其他两个已知变量的值，对于不同的变量，通过调用相应的函数进行求解并显示最终的结果。
2 程序界面3 流程图4 主要的函数符号规定如下：b(Blocking)：阻塞率；a(BHT)：到达呼叫量；s(Lines)：中继线数量。
1）已知到达呼叫量a及中继线数量s求阻塞率b 使用迭代算法提高程序效率B(s,a)=a∙B(s−1,a) s+a∙B(s−1,a)代码如下：function b = ErlangB_b(a,s)b =1;for i =1:sb = a * b /(i + a * b);endend2）已知到达呼叫量a及阻塞率b求中继线数量s考虑到s为正整数，因此采用数值逼近的方法。
采用循环的方式，在每次循环中增加s的值，同时调用 B(s,a)函数计算阻塞率并与已知阻塞率比较，当本次误差小于上次误差时，结束循环，得到s值。
代码如下：function s = ErlangB_s(a,b)s =1;Bs = ErlangB_b(a,s);err = abs(b-Bs);err_s = err;while(err_s <= err)err = err_s;s = s +1;Bs = ErlangB_b(a,s);err_s = abs(b - Bs);ends = s -1;end3）已知阻塞率b及中继线数量s求到达呼叫量a考虑到a为有理数，因此采用变步长逼近的方法。
北邮 通信原理硬件实验报告
实验一：双边带抑制载波调幅（DSB-AM）一、实验目的：(1)了解DSB-SC AM信号的产生以及相干解调的原理和实现方法；(2)了解DSB-SC AM信号波形以及振幅频谱特点,并掌握其测量方法；(3)了解在发送DSB-SC AM 信号加导频分量的条件下,收端用锁相环提取载波的原理及其实现方法；(4)掌握锁相环的同步带和捕捉带的测量方法,掌握锁相环提取载波的调试方法。
二、实验系统框图：DSB-SC加导频的产生测量VCO压控灵敏度的框图DSB-SC加导频分量的相干解调及载波提取框图三、实验步骤：SC-DSB 信号的数学表达式为s(t)=Acm(t)cos(Wct)，这个实验产生SC-DSB 的方法很简单，就是用载波跟调制信号直接相乘，其中载波是由主振荡器产生为幅度为1V，频率为100KHZ的正弦波，而调制信号由音频振荡器产生的正弦信号再经缓冲放大器组成，幅度为1V，频率为1KHZ。
1、DSB-SC AM 信号的产生1）按照图连接，将音频振荡器输出的模拟音频信号及主振荡器输出的100KHz模拟载频信号分别用连接线连至乘法器的两个输出端；2）用示波器观看音频输出信号的信号波形的幅度以及振荡频率，调整音频信号的输出频率为10kHz,作为均值为0的调制信号m(t)；3）用示波器观看主振荡器输出信号的幅度以及振幅频谱；4）用示波器观看乘法器的输出波形，并注意已调信号波形的相位翻转与已调信号波形；5）测量已调信号的波形频谱，注意其振幅频谱的特点；6）调整增益G=1:将加法器的B 输出端接地，A 输入端接已调信号，用示波器观看加法器的输出波形以及振幅频谱，使加法器输入与加法器输出幅度一致；7）调整增益g；加法器A 端接已调信号，B 接导频信号。
用频谱仪观看加法器输出信号的振幅频谱，调节增益g 旋钮，使导频信号振幅频谱的幅度为已调信号的边带频谱幅度的0.8倍。
此导频信号功率为已调信号功率的0.32倍。
2、DSB-SC AM 信号的相干解调及载波提取1）锁相环的调试1 单独测量VCO的性能将VCO 模板前面板的频率开关拨到HI 载波频段的位置，VCO 的Vin 输入端暂不接信号。
北邮移动通信实验报告
北邮移动通信实验报告北邮移动通信实验报告一、引言本实验报告旨在总结北邮移动通信实验的实施情况、结果和分析，对实验数据进行归纳和解释，以及提出相应的建议和改进措施。
本实验旨在深入研究移动通信领域的相关技术，并通过实际操作和数据分析，加深对移动通信原理和应用的理解。
二、实验概述1·实验目的本实验的目的是通过模拟移动通信系统的工作原理和性能进行实际操作，熟悉移动通信系统的基本原理、标准和技术，并对系统的性能进行测试和评估。
2·实验设备和软件工具本实验使用的设备和软件工具包括：●移动通信实验设备（包括基站、移动终端、信道仿真器等）●相关的软件平台和工具（如Matlab、C++开发环境等）3·实验步骤本实验的步骤如下：●确定实验需求和目标，设计实验方案●配置实验设备和软件环境●进行实验操作和数据采集●对实验数据进行处理和分析●总结实验结果，提出建议和改进措施三、实验结果与分析1·实验数据收集和处理本实验收集到的数据主要包括移动通信系统的性能参数、信道传输情况、功率消耗等方面的指标。
收集到的数据经过处理和分析后，得出以下结论：●移动通信系统的覆盖范围和容量与基站的功率、天线高度、信道特性等因素相关●数据传输速率与信道带宽、调制方式和信噪比等因素相关2·实验结果分析根据实验数据的分析，可以得出以下结论：●移动通信系统的覆盖范围和容量可以通过调整基站的功率和天线高度来改善●数据传输速率可以通过增加信道带宽、改变调制方式和提高信噪比来提升四、实验总结1·实验成果本实验通过实际操作和数据分析，对移动通信系统的工作原理和性能有了更加深入的认识，对移动通信技术的应用和发展有了更加清晰的了解。
2·实验建议根据本实验的结果和分析，提出以下建议和改进措施：●在设计移动通信系统时，需要充分考虑基站的功率和天线高度对系统覆盖范围和容量的影响●需要注重提升信道传输质量，通过增大信道带宽、改变调制方式和提高信噪比等手段来提高数据传输速率五、附件本文档涉及的附件包括实验数据记录表、实验方案设计图等。
北邮光纤通信实验报告
信息与通信工程学院光纤通信实验报告班级：姓名：学号：实验合作小组：一、OTDR的使用1、实验原理OTDR使用瑞利散射和菲涅尔反射来表征光纤的特性。
瑞利散射是由于光信号沿着光纤产生无规律的散射而形成。
OTDR就测量回到OTDR端口的一部分散射光。
这些背向散射信号就表明了由光纤而导致的衰减（损耗/距离）程度。
形成的轨迹是一条向下的曲线，它说明了背向散射的功率不断减小，这是由于经过一段距离的传输后发射和背向散射的信号都有所损耗。
给定了光纤参数后，瑞利散射的功率就可以标明出来，如果波长已知，它就与信号的脉冲宽度成比例：脉冲宽度越长，背向散射功率就越强。
瑞利散射的功率还与发射信号的波长有关，波长较短则功率较强。
也就是说用1310nm信号产生的轨迹会比1550nm信号所产生的轨迹的瑞利背向散射要高。
在高波长区（超过1500nm），瑞利散射会持续减小，但另外一个叫红外线衰减（或吸收）的现象会出现，增加并导致了全部衰减值的增大。
因此，1550nm是最低的衰减波长；这也说明了为什么它是作为长距离通信的波长。
很自然，这些现象也会影响到OTDR。
作为1550nm 波长的OTDR，它也具有低的衰减性能，因此可以进行长距离的测试。
而作为高衰减的1310nm 或1625nm波长，OTDR的测试距离就必然受到限制，因为测试设备需要在OTDR轨迹中测出一个尖锋，而且这个尖锋的尾端会快速地落入到噪音中。
菲涅尔反射是离散的反射，它是由整条光纤中的个别点而引起的，这些点是由造成反向系数改变的因素组成，例如玻璃与空气的间隙。
在这些点上，会有很强的背向散射光被反射回来。
因此，OTDR就是利用菲涅尔反射的信息来定位连接点，光纤终端或断点。
OTDR的工作原理就类似于一个雷达。
它先对光纤发出一个信号，然后观察从某一点上返回来的是什么信息。
这个过程会重复地进行，然后将这些结果进行平均并以轨迹的形式来显示，这个轨迹就描绘了在整段光纤内信号的强弱。
北邮通信原理软件实验报告
北邮通信原理软件实验报告一、实验目的本次实验的目的是通过使用软件进行通信原理实验，探究数字通信系统的原理和性能。
二、实验内容1.利用软件计算并绘制理想低通滤波器的频率响应曲线。
2.通过软件模拟并比较维纳滤波器与理想低通滤波器的频率响应曲线。
3.仿真带通调制器和解调器在理想信道中的性能。
三、实验步骤1.理想低通滤波器的设计：（1）利用软件，设置滤波器参数，如截止频率和滤波器类型。
（2）计算并绘制理想低通滤波器的频率响应曲线。
2.维纳滤波器与理想低通滤波器的比较：（1）利用软件设置维纳滤波器参数，如截止频率和信噪比。
（2）仿真并比较维纳滤波器与理想低通滤波器的频率响应曲线。
3.带通调制器和解调器的性能仿真：（1）设置带通调制器和解调器的参数，如载波频率和调制系数。
（2）仿真并分析带通调制器和解调器的性能，如频率响应和误码率。
四、实验结果五、实验分析通过本次实验，我们对数字通信系统的原理和性能有了更深入的了解。
首先，理想低通滤波器的频率响应曲线能够更清晰地展现滤波器的特性，帮助我们更好地了解滤波器的设计和应用。
其次，维纳滤波器相对于理想低通滤波器而言，频率响应存在一定的失真，但对于噪声有一定的抑制作用。
在实际应用中，需要根据具体需求选择适合的滤波器类型。
最后，带通调制器和解调器的性能仿真结果能够帮助我们评估系统的性能，如误码率和频率响应。
通过调整调制系数和载波频率，可以使系统在一定范围内具有较好的性能。
六、实验结论通过本次实验的软件仿真，我们探究了数字通信系统的原理和性能。
实验结果显示，理想低通滤波器具有良好的频率响应特性；维纳滤波器能对噪声进行一定的抑制，但频率响应存在一定的失真；带通调制器和解调器在适当的调制系数和载波频率下能够获得较低的误码率和良好的频率响应。
总之，本次实验通过软件仿真，使我们更好地理解了通信原理中的数字通信系统及其性能分析方法，提高了我们的实践能力和理论知识水平。
北邮现代通信技术实验报告
北邮现代通信技术实验报告实验名称：现代通信技术实验实验目的：1. 理解现代通信技术的基本理论和原理。
2. 掌握数字通信系统的基本组成和工作流程。
3. 熟悉通信系统中信号的调制与解调过程。
4. 学会使用通信系统实验设备，进行实验操作和数据分析。
实验原理：现代通信技术主要依赖于数字信号处理技术，通过数字信号的调制与解调实现信息的传输。
在本实验中，我们将学习数字通信系统中的信号调制方法，如幅度键控（ASK）、频率键控（FSK）和相位键控（PSK），以及相应的解调技术。
实验设备与材料：1. 计算机一台，安装有通信仿真软件。
2. 通信原理实验箱一套，包括调制解调模块、信号源模块等。
3. 通信信号发生器。
4. 示波器。
实验步骤：1. 打开通信仿真软件，设置实验参数，如信号频率、调制方式等。
2. 使用通信信号发生器产生模拟信号，输入到通信原理实验箱的信号源模块。
3. 通过实验箱的调制模块对信号进行调制，观察示波器上信号的变化。
4. 将调制后的信号传输至解调模块，观察解调后的信号波形。
5. 记录实验数据，包括调制前后的信号波形、频谱特性等。
实验结果：通过实验，我们得到了以下结果：1. 调制信号与原始信号的波形对比，展示了调制过程中信号的变化。
2. 解调后的信号与原始信号的对比，验证了调制解调技术的准确性。
3. 通过频谱分析，观察到调制信号的频谱特性，理解了调制对信号频谱的影响。
实验分析：在实验过程中，我们发现不同调制方式对信号的影响各有不同。
例如，ASK调制主要改变信号的幅度，而FSK和PSK调制则分别改变信号的频率和相位。
通过解调过程，我们能够从调制信号中恢复出原始信号，验证了通信系统的有效性。
实验结论：通过本次实验，我们深入理解了现代通信技术中的数字信号调制与解调过程。
实验结果表明，通过合理的调制解调技术，可以有效实现信息的传输和恢复。
同时，实验也加深了我们对通信系统基本原理的认识，为进一步学习通信技术打下了坚实的基础。
北邮通信认知实习报告
作为一名电子信息科学与技术专业的学生，了解通信基础知识，掌握通信专业知识的学习方法，明白通信行业最前沿的科学知识对我们未来的发展有重要作用。
通信网概论：通信网是由一定数量的节点(包括终端节点、交换节点)和连接这些节点的传输系统有机地组织在一起的，按约定的信令或协议完成任意用户间信息交换的通信体系。
用户使用它可以克服空间、时间等障碍来进行有效的信息交换。
通信网按功能与用途不同，一般可分为物理网、业务网和支撑网等三种。
物理网是由用户终端、交换系统、传输系统等通信设备所组成的实体结构，是通信网的物质基础，也称装备网。
用户终端是通信网的外围设备，它将用户发送的各种形式的信息转变为电磁信号送入通信网路传送。
交换系统是各种信息的集散中心，是实现信息交换的关键环节。
传输系统是信息传递的通道，它将用户终端与交换系统之间以及交换系统相互之间联接起来，形成网路。
传输系统按传输媒介的不同，可分为有线传输系统和无线传输系统两类。
业务网是疏通电话、电报、传真、数据、图像等各类通信业务的网路，是指通信网的服务功能。
按其业务种类，可分为电话网、移动通信网，数据通信网，智能网等。
电话网是各种业务的基础，数据网可由传输数据信号的电话电路或专用电路构成。
业务网具有等级结构，即在业务中设立不同层次的交换中心，并根据业务流量、流向、技术及经济分析，在交换机之间以一定的方式相互联接。
支撑网是为保证业务网正常运行，增强网路功能，提高全网服务质量而形成的网络。
在支撑管理网中传递的是相应的控制、监测及信令等信号。
按其功能不同，可分为信令网、同步网和管理网。
信令网由信令点、信令转接点、信令链路等组成，旨在为公共信道信令系统的使用者传送信令。
光传送网络：在光通信实验室中，老师向我们介绍了光通信的基本原理及优缺点和光传送网络以及基本原理及优缺点和基于PON的全业务光接入系统。
光传输是在发送方和接收方之间以光信号形态进行传输的技术。
广义的光通信指的一切运用光波作为载体而传送信息的所有通信方式的总称，通常把应用石英光纤的有线光通信称为光通信。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
北京邮电大学实验报告通信网理论基础实验报告学院:信息与通信工程学院班级:2013211124学号:姓名:实验一 ErlangB公式计算器一实验内容编写Erlang B公式的图形界面计算器，实现给定任意两个变量求解第三个变量的功能：1）给定到达的呼叫量a和中继线的数目s，求解系统的时间阻塞率B；2）给定系统的时间阻塞率的要求B和到达的呼叫量a，求解中继线的数目s，以实现网络规划；3）给定系统的时间阻塞率要求B以及中继线的数目s，判断该系统能支持的最大的呼叫量a。
二实验描述1 实验思路使用MATLAB GUITOOL设计图形界面，通过单选按钮确定计算的变量，同时通过可编辑文本框输入其他两个已知变量的值，对于不同的变量，通过调用相应的函数进行求解并显示最终的结果。
2 程序界面3 流程图4 主要的函数符号规定如下：b(Blocking)：阻塞率；a(BHT)：到达呼叫量；s(Lines)：中继线数量。
1）已知到达呼叫量a及中继线数量s求阻塞率b 使用迭代算法提高程序效率B(s,a)=a∙B(s−1,a) s+a∙B(s−1,a)代码如下：function b = ErlangB_b(a,s)b =1;for i =1:sb = a * b /(i + a * b);endend2）已知到达呼叫量a及阻塞率b求中继线数量s考虑到s为正整数，因此采用数值逼近的方法。
采用循环的方式，在每次循环中增加s的值，同时调用 B(s,a)函数计算阻塞率并与已知阻塞率比较，当本次误差小于上次误差时，结束循环，得到s值。
代码如下：function s = ErlangB_s(a,b)s =1;Bs = ErlangB_b(a,s);err = abs(b-Bs);err_s = err;while(err_s <= err)err = err_s;s = s +1;Bs = ErlangB_b(a,s);err_s = abs(b - Bs);ends = s -1;end3）已知阻塞率b及中继线数量s求到达呼叫量a考虑到a为有理数，因此采用变步长逼近的方法。
采用循环的方式，在每次循环中增加a的值（步长为 s/2），同时调用 B(s,a)函数计算阻塞率并与已知阻塞率比较，当本次误差小于预设阈值时，结束循环，得到a值。
代码如下：function a = ErlangB_a(b,s)a =0;step = s /2;b_temp = ErlangB_b(a,s);b_temps = ErlangB_b(a + step,s);while(abs(b - b_temp)>0.000001)if((b - b_temp)*(b - b_temps)<0) step = step /2;elsea = a + step;endb_temp = ErlangB_b(a,s);b_temps = ErlangB_b(a + step,s);endend三结果分析1) 已知到达呼叫量a及中继线数量s求阻塞率b2) 已知到达呼叫量a及阻塞率b求中继线数量s3) 已知阻塞率b及中继线数量s求到达呼叫量a4）分析将程序结果与在线ErlangB公式计算器结果进行比较，结果基本一致，说明程序的正确性。
实验二 M/M/1排队系统一实验内容实现M/M/1单窗口无限排队系统的系统仿真，利用事件调度法实现离散事件系统仿真，并统计平均队列长度以及平均等待时间等值，以与理论分析结果进行对比。
二实验原理1 顾客到达模式设到达过程是一个参数为 λ 的Poisson过程，则长度为 t 的时间内到达 k 个呼叫的概率P k(t)服从Poisson分布，即P k(t)=(λt)kk!∙eλt2 服务模式设每个呼叫的持续时间为 t，则服从参数为 μt 的负指数分布，即P(X>t)=1−e−μt,t≥03 服务规则服务规则为：FIFO4 理论分析设系统到达率为 λ，服务率为 μ，则理论分析如下：理论平均等待时间：E[w]=1μ−λ理论平均排队时间：E[q]=λμ∙(μ−λ)理论系统中平均顾客数：E[N]=λμ−λ理论系统中平均等待队长：E[N q]=λ∙λμ∙(μ−λ)三实验描述1 实验思路仿真时序图示例时序关系图如下：各参数含义如下：b i：第i个事件（到达或离开）发生的时间t i：第i个到达事件发生的时间c i：第i个离开事件发生的时间A i：第（i-1）个顾客与第i个顾客到达时间间隔D i：第i个顾客排队等待时间间隔S i：第i个顾客服务时间长度仿真设计算法1）利用负指数分布与泊松过程的关系，产生符合泊松过程的顾客流；2）分别构建一个顾客到达队列和一个顾客等待队列。
顾客到达后，首先进入到达的队尾排队，并检测是否有顾客等待以及是否有服务台空闲，如果无人等待并且有服务员空闲则进入服务状态，否则顾客将进入等待队列的队尾等待；3）产生符合负指数分布的随机变量作为每个顾客的服务时间；4）当服务员结束一次服务后，就取出等待队列中位于对头的顾客进入服务状态，如果迭代队列为空则服务台空闲等待下一位顾客到来；5）在系统到达稳态时，计算系统的平均等待时间以及平均等待队长等数据。
2 流程图3 主要的函数1) 产生符合泊松过程的顾客流Interval_Arrive=-log(rand(1,SimTotal))/Lambda; %Arrival Time Interval 2) 产生符合负指数分布的服务时间Interval_Serve=-log(rand(1,SimTotal))/Mu;%Service Time3) 计算顾客到达时间%******Arrive Time for each Customer*****t_Arrive(1)=Interval_Arrive(1);ArriveNum(1)=1;for i=2:SimTotalt_Arrive(i)=t_Arrive(i-1)+Interval_Arrive(i);ArriveNum(i)=i;end4) 计算顾客离开时间%%Leave Time for each Customert_Leave(1)=t_Arrive(1)+Interval_Serve(1);LeaveNum(1)=1;for i=2:SimTotalif t_Leave(i-1)<t_Arrive(i)%New customer arrives after the former has been served.t_Leave(i)=t_Arrive(i)+Interval_Serve(i);else%New customer arrives while the former is being served.t_Leave(i)=t_Leave(i-1)+Interval_Serve(i);endLeaveNum(i)=i;end5) 计算平均等待时间t_Wait=t_Leave-t_Arrive;%Waiting Time for each Customert_Wait_avg=mean(t_Wait);%Average Waiting Time6) 计算平均排队时间t_Queue=t_Wait-Interval_Serve;%Queueing Time for each Customert_Queue_avg=mean(t_Queue);%Average Queueing Time7) 计算平均顾客数%%TimelineTimepoint=[t_Arrive,t_Leave];%Record Arrival Time & Leaving Time for each Customer.Timepoint=sort(Timepoint);%Sort all Timepoints.ArriveFlag=zeros(size(Timepoint));%Flag of ArrivalCusNum=zeros(size(Timepoint));temp=2; CusNum(1)=1;for i=2:length(Timepoint)if(temp<=length(t_Arrive))&&(Timepoint(i)==t_Arrive(temp))%if the timepoint is about ArrivalCusNum(i)=CusNum(i-1)+1;%add Number of Customerstemp=temp+1;ArriveFlag(i)=1;else%if the timepoint is about LeavingCusNum(i)=CusNum(i-1)-1;%reduce Number of Customers endend%%Average Number of CustomersTime_interval=zeros(size(Timepoint));Time_interval(1)=t_Arrive(1);for i=2:length(Timepoint)Time_interval(i)=Timepoint(i)-Timepoint(i-1);%Timepoint Interval endCusNum_fromStart=[0 CusNum];%Average Number of CustomersCusNum_avg=sum(CusNum_fromStart.*[Time_interval 0])/Timepoint(end);8) 计算平均队长%%Length of Queue(Excluding Customer being Serviced)QueLength=zeros(size(CusNum));for i=1:length(CusNum)if CusNum(i)>=2QueLength(i)=CusNum(i)-1;%Excluding Customer being Serviced elseQueLength(i)=0;endend%Average Waiting Length of QueueQueLength_avg=sum([0 QueLength].*[Time_interval 0])/Timepoint(end);9) 仿真图：各顾客到达时间与离去时间subplot(2,2,1);stairs([0 ArriveNum],[0 t_Arrive],'b');hold on;stairs([0 LeaveNum],[0 t_Leave],'r');legend('Arrive Time','Leave Time');title('Arrive Time and Leave Time for each customer'); xlabel('Customer ID');ylabel('Arrive/Leave Time(Accumulated)');10) 仿真图：系统等待队长分布subplot(2,2,2);stairs(Timepoint,CusNum,'b')axis([0 Timepoint(end)0 max(CusNum)])title('Distribution of System Customer Number');xlabel('Time(Accumulated)');ylabel('Customer Number');11) 仿真图：各顾客在系统中的排队时间和等待时间subplot(2,2,3);stairs([0 ArriveNum],[0 t_Queue],'b');hold on;stairs([0 LeaveNum],[0 t_Wait],'r');hold off;axis([0 max(max(ArriveNum),max(LeaveNum))0max(max(t_Queue),max(t_Wait))]);title('Queueing Time and Waiting Time for each customer'); xlabel('Customer ID');ylabel('Queueing/Waiting Time');legend('Queueing Time','Waiting Time');四结果分析为了检验算法的正确性，我们进行了多次实验。