MATLAB通信仿真

南京工程学院课程设计说明书（论文）题目带限数字基带传输系统的仿真课程名称Matlab通信仿真设计院(系、部、中心)通信工程学院专业无线通信班级无线通信131学生姓名陈明星学号208130115设计地点信息楼C322指导教师刘伟伟设计起止时间：2015年12月14日至2014年12月18日成绩一、课程设计目的和要求1.综合运用信号与线性系统，通信原理，以及matlab语言等综合应用知识。2.基本熟悉软件开发流程、要点以及注意事项。3.掌握基带数字传输系统的仿真方法。4.熟悉基带传输系统的基本结构。5.掌握带限信道的仿真以及性能分析。6.通过观察眼图判断信号的传输质量。7.设计和编程，以加深对通信系统基础知识的理解，提高学生利用matlab在通信领域进行实践的应用能力、分析解决问题的能力；增强学生对实践应用流程模式的理解。二：课设设备1.硬件：PC机2.软件：MATLAB7.1，Windows7操作系统三：课程设计原理及背景知识在实际通信中传输信道的带宽是有限的，这样的信道称为带限信道。带限信道的冲激响应在时间上是无限的，因此一个时隙内的代表数据的波形经过带限信道后将在邻近的其他时隙上形成非零值，称为波形的拖尾。拖尾和邻近其他时隙上的传输波形相互叠加后，形成传输数据之间的混叠，造成符号间干扰，也称为码间干扰。接收机中，在每个传输时隙中的某时间点上，通过对时域混叠后的波形进行采样，然后对样值进行判决来恢复接收数据。在采样时间位置上符号间的干扰应最小化（该采样时刻称为最佳采样时刻），并以适当的判决门限来恢复接收数据，使误码率最小（该门限称为最佳判决门限）。在工程上，为了便于观察接收波形中的码间干扰情况，可在采样判决设备的输入端以恢复的采样时钟作为同步，用示波器观察该端口的接收波形。利用示波管显示的暂时记忆特性，在示波管上将显示多个时隙内接收信号的重叠波形图案，称为眼图。对于传输符号为等概的双极性码，最佳判决门限为0，最佳采样时刻为眼图开口最大处，因为该时刻上的码间干扰最小。当无码间干扰时，在最佳采样时刻上眼图波形将会聚为一点。显然，只要带限信号冲击响应的拖尾波在时隙周期整数倍上的值为0，那么在采样时刻就没有码间干扰，例如抽样函数xxxSasin)(。然而，抽样函数的频谱时门函数，物理不可实现，即使近似实现也十分困难。还存在一类无码间干扰的时域函数，具有升余弦频率特性，幅频响应是缓变的，在工程上易于实现，其冲激响应为:2/41/cos/t/Tsin22scosssstrTtTtTth其中，sT为码元传输时隙宽度，10为滚降系数。α=0时，退化为矩形门函数；α=1时，为全升余弦。其傅里叶变换cosrH即为相应的系统函数。设发送滤波器为TG，物理信道的传递函数为C，接收滤波器为RG，则带限信号总的传递函数为RTGCGH。对于物理信道是AWGN信道的情况，可以证明，当发送滤波器与接收滤波器相互匹配时，即RTGG*，系统误码率最小。对于理想的物理信道（C=1），收发滤波器相互匹配时有2*TTGGGHR由此求得收发滤波器传递函数的实数解为HGGRT无码间干扰情况下，信道总的传递函数是滚升余弦的，匹配的收发滤波器称为平方根升余弦滤波器，有rcosHRTGG工程上，滚升余弦滤波器和平方根滚升余弦滤波器通常用FIR滤波器来近似实现。FIR滤波器的分母系数为1，分子系数向量等于冲激响应的采样序列。Matlab通信工具箱中提供了设计升余弦滤波器的函数rcosine。函数rcosine用于计算FIR滤波器时的用法如下。Num=rcosine(Fd,Fs,’fir/normal’,r,delay)‘fir/normal’用于FIR滚升余弦滤波器设计num=rcosine(Fd,Fs,’fir/sqrt’,r,delay)‘fir/sqrt’用于FIR平方根滚升余弦滤波器设计r是滚降系数；Fd是输入数字序列的采样率即码元速率；Fs为滤波器的采样率，Fs必须是Fd的整数倍；delay是输入到响应峰值之间的延迟（单位是码元时隙数）。四：课程设计任务任务一：升余弦波形及其频谱研究用MATLAB通信工具箱提供的rcosine函数作出一组滚升余弦滤波器的冲击响应，滚降系数为0，0.5，0.75，1，并通过FFT求出其幅频特性。码元时隙为1ms，在一个码元时隙内采样10次，滤波器延时为5个码元时隙。任务二：滚升余弦滤波器设计设计一个滚升余弦滤波器，滚降系数为0.75。输入为4元双极性数字序列，符号速率为1000波特，设滤波器采样率为10000次/s，即在一个符号间隔中有10个采样点。请建立simulin仿真模型观察升余弦滤波器的输出波形，眼图及功率谱。任务三：带限基带传输模型建立一个基带传输模型，发送数据为二进制双极性不归零码，发送滤波器为平方根升余弦滤波器，滚降系数为0.5，信道为AWGN信道，接收滤波器与发送滤波器相匹配。发送数据率为1000bps，要求观察接收信号眼图，并设计接收机采样判决部分，对比发送数据与恢复数据波形，统计误码率。五：课程设计成果内容及相关截图任务一：MATLAB编程如下：clc;Fd=1e3;Fs=Fd*10;delay=5;forr=[0,0.5,0.75,1]num=rcosine(Fd,Fs,'fir/normal',r,delay)t=0:1/Fs:1/Fs*(length(num)-1);figure(1);plot(t,num);axis([00.01-0.21.1]);holdon;Hw=abs(fft(num,1000));f=(1:Fs/1000:Fs)-1;figure(2);plot(f,Hw);axis([01500012]);holdon;end升余弦波形图：升余弦波形图对应的幅频特性曲线：任务二：滚升余弦滤波器和眼图测试模型模块分析：系统仿真步进设为1e-4s,采用RandomIntegerGenerator产生采样间隔为1e-3的4元整数（0，1，2，3），并用UnipolartoBipolarConverter模块将其转换为双极性（-3，-1，1，3）。通过升速率模块Upsample将基带数据的采样速率升高为10000次每秒，其输出为冲激脉冲形式的数据序列。滚升余弦FIR滤波器以DiscreteFilter模块实现，其分母系数设置为1，分子系数通过rcosine函数计算。这样，就得到了滚降系数为0.75的滚升余弦滤波器，滤波器延时时间为3个数据时隙，即30个滤波器采样间隔。滤波器输出通过Downsample模块降低4倍采样速率，使送入频谱仪的采样率为2500次/s，这样频谱仪显示的频谱范围是0~1250HZ。同时，滤波输出送入通信模块库中的眼图显示模块Discrete-TimeEyeDiagramScope显示眼图。在眼图显示模块中需要设置：（1）每个数据的采样点数，设为10；（2）每次扫描显示的符号个数设为2，这样眼图将显示2个符号时间宽度；（3）显示所保留的扫描波形轨迹图，可使用默认值；（4）每次显示的新轨迹数，也可使用默认值；（5）Discrete-TimeEyeDiagramScope模块可同时显示同相支路和正交支路上的波形眼图，本题只有一个支路，可选择In-phaseOnly选项。仿真结果截图任务三：高斯信道下基带传输系统测试模型模块分析：设计系统仿真采样率为1e4Hz,滤波器采样速率等于系统仿真采样率。数字信号速率为1000bps，故在进入发送滤波器之前需要10倍升速率，接收解码后再以10倍降速率来恢复信号传输比特率。仿真模型如图所示，系统分为二进制信源、发送滤波器、高斯信道、接收匹配滤波器、接收采样、判决恢复以及信号测量等7部分。二进制信源输出双极性不归零码，并想接收端提供原始数据以便对比和统计误码率。发送滤波器和接收滤波器均为平方根升余弦滤波器，高斯信道采用随机数发生器和加法器实现。接收定时假设为理想的，采用脉冲发生器实现1000HZ的矩形脉冲作为恢复定时脉冲，以乘法器实现在最佳采样时刻对接收滤波器输出的采样。然后对抽样结果进行门限判决，最佳判决门限为零，判决输出结果在一个传输码元时隙内不变，最后以10倍降速率采样得出采样率为1000HZ的恢复数据。TriggeredSubsystem保持子系统模型高斯信道下基带传输系统测试眼图仿真结果测试误码率部分截图（错误码元为9个，共发送10000个码元）发送和接收的二进制双极性不归零码六：课程设计心得与体会在这次课程设计中，我有如下几点总结与思考：首先：我认识了MATLAB中的SIMULINK仿真，通过对照课设指导书，能独立地完成课设任务，通过参考相关MATLAB方面的书籍，自己能设计出一个SIMULINK的例子。其次：在此次课程设计中，需要用到通信原理方面的相关知识，然而通信原理是上学期学的，其中相关的知识都已经忘得差不多了，我不得不重新拿起通信原理的教材，简单的看一看，所以说以前学过的知识要复习，要温故。第三点：不管做啥事都要有计划性，要按照计划，逐个攻克，不能“盲，忙”。课设前两天，我无计划的做着任务，一遇到问题，就跳过，接着做下一个，结果两天下来，一个任务也没有完成。于是我不得不停下来，调整思路，制定自己的计划。最后一点：通过此次课程设计，我不仅复习了通信原理的知识，而且学习了matlab的SIMULINK仿真（不仅复习了旧知识，而且学习了新知识，虽然这新知识对我以后所从事的工作无关，但毕竟开阔了眼界），如果让我自己去设计，我相信有了这次课程设计的基础，我在设计的时候应该会更容易，更有信心去完成好设计。