通信原理课程设计----基带信号眼图仿真

1课程设1需求分析1．设基带传输系统响应是α＝1的升余弦滚降系统，画出在接收端的基带数字信号波形及其眼图。2．设定二进制数字基带信号an∈｛＋1,－1｝，g(t)=10≤t≤Ts；t为其他值时g(t)=0。系统加性高斯白噪声的双边功率谱密度为0。画出：(1)经过理想低通H(f)=1│f│≤5/(2Ts)后的眼图。(2)经过理想低通H(f)=1│f│≤1/Ts后的眼图。(3)比较分析上面图形。在该部分中叙述：对题目中要求的功能进行的简单的叙述分析，把题目内容给介绍一下，还需要介绍分工情况。2概要设计1、基带传输特性基带系统的分析模型如图3-1所示，要获得良好的基带传输系统，就应该2nsnatnT()HnsnahtnT基带传输抽样判决图2-1基带系统的分析模型抑制码间干扰。设输入的基带信号为nsnatnT，sT为基带信号的码元周期，则经过基带传输系统后的输出码元为nsnahtnT。其中1()()2jthtHed（3-1）理论上要达到无码间干扰，依照奈奎斯特第一准则，基带传输系统在时域应满足：10()0,skhkTk，为其他整数(3-2)频域应满足：()0,ssTTH，其他(3-3)0sT()HsTsT图2-2理想基带传输特性此时频带利用率为2/BaudHz,这是在抽样值无失真条件下，所能达到的最高频率利用率。由于理想的低通滤波器不容易实现，而且时域波形的拖尾衰减太慢，因此在得不到严格定时时，码间干扰就可能较大。在一般情况下，只要满足：222(),sissssiHHHHTTTTT(3-4)3基带信号就可实现无码间干扰传输。这种滤波器克服了拖尾太慢的问题。从实际的滤波器的实现来考虑，采用具有升余弦频谱特性()H时是适宜的。(1)(1)1sin(),2(1)()1,0(1)0,ssssssTTTTHTT(3-5)这里称为滚降系数，01。所对应的其冲激响应为：222sincos()()14sssstTtThtttTT(3-6)此时频带利用率降为2/(1)Baud/Hz，这同样是在抽样值无失真条件下，所能达到的最高频率利用率。换言之，若输入码元速率'1/ssRT，则该基带传输系统输出码元会产生码间干扰。2、眼图所谓眼图就是将接收滤波器输出的，未经再生的信号，用位定时以及倍数作为同步信号在示波器上重复扫描所显示的波形（因传输二进制信号时，类似人的眼睛）。干扰和失真所产生的畸变可以很清楚的从眼图中看出。眼图反映了系统的最佳抽样时间，定时的灵敏度，噪音容限，信号幅度的畸变范围以及判决门限电平，因此通常用眼图来观察基带传输系统的好坏。4图2-3眼图示意图部分的算法设计说明（可以是描述算法的流程图）；每个程序中使用的存储结构设计说明（如果题目已经指定了数据存储的，按照指定的设计，并且写出该存储结构的定义）。也就是简单介绍一下你程序中自定义的函数（函数名，参数，以及功能）等。3运行环境Windows74开发工具和编程语言MATLABV6.05详细设计clearall;closeall;N_eye=7;%眼图中码元个数N0=0;%高斯白噪声功率谱设为0Ts=2;%码元宽度设为1N_sample=10;%每个码元抽样10个点dt=Ts/N_sample;gt=ones(1,N_sample);%矩形不归零脉冲N_code=1000;%码元个数5t=0:dt:N_code*N_sample*dt-dt;data=sign(randn(1,N_code));data1=sigexpand(data,N_sample);st=conv(data1,gt);%产生NRZ基带信号h1t=5/Ts*sinc(5*(t-5)/Ts);%滤波器1的冲激响应h2t=2/Ts*sinc(2*(t-5)/Ts);%滤波器2的冲激响应h3t=1/Ts*sinc(1*(t-5)/Ts);%滤波器3的冲激响应h4t=0.5/Ts*sinc(0.5*(t-5)/Ts);r1t=conv(h1t,st);r2t=conv(h2t,st);r3t=conv(h3t,st);r4t=conv(h4t,st);r1t=[0r1t*dt];r2t=[0r2t*dt];r3t=[0r3t*dt];r4t=[0r4t*dt];sigma1=sqrt(2.5*N0/Ts);sigma2=sqrt(N0/Ts);sigma3=sqrt(0.5*N0/Ts);sigma4=sqrt(0.25*N0/Ts);n1=sigma1*randn(1,length(r1t));n2=sigma2*randn(1,length(r2t));n3=sigma2*randn(1,length(r3t));%加入白噪声，该白噪声是通过滤波器后的带限白噪声6n4=sigma2*randn(1,length(r4t));r1t=r1t+n1;r2t=r2t+n2;r3t=r3t+n3;r3t=r3t+n3;%至此得到基带信号分别通过滤波器1和滤波器2滤波器3滤波器4的输出波形ttt=0:dt:N_eye*N_sample*dt-dt;figure(1)subplot(411)%画出滤波器1输出的眼图forn=9:1:100ss=r1t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss);holdon;endylabel('B=2.5Hz');subplot(412)%画出滤波器2输出的眼图forn=9:1:100ss=r2t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss)holdon;endylabel('B=1Hz')subplot(413)%画出滤波器3输出的眼图forn=9:1:100ss=r3t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;7plot(ttt,ss)holdon;endylabel('B=0.5Hz')subplot(414)%画出滤波器4输出的眼图forn=9:1:100ss=r4t(n*N_sample*N_eye+1:(n+1)*N_sample*N_eye);drawnow;plot(ttt,ss)holdon;endylabel('B=0.25Hz')6调试分析刚开始的时候，程序出现了问题???Undefinedcommand/function'sigexpand'.之后添加了一个函数8function[out]=sigexpand(d,M)%将输入的序列扩展成间隔为N-1个0的序列N=length(d);out=zeros(M,N);%例如，M=8，N=1000，产生一个8行1000列的矩阵out(1,:)=d;%d的每一个元素赋值给out矩阵每列的第一个元素，因为d中共1000个元素（0，-1或者1），且out矩阵共有1000列，所以当然成立out=reshape(out,1,M*N);%把out这个矩阵重新变成一个1行M*N列的矩阵，即1行8000列的有8000个元素的矩阵，其中每隔7个0就是一个原来的d中的元素之后再调用sigexpand()函数的到时候就就没有问题了。7测试结果9程序运行结果如图7-1所示图7-1程序运行结果图比较上图中的1和3可以发现，3的性能比1的性能。1中的失真比较严重。参考文献[1]孙屹，MATLAB通信仿真开发手册，北京：国防工业出版社[2]李建新，现代通信系统分析与仿真－MATLAB通信工具箱，西安：西安电子科技大学出版社[3]曹志刚，现代通信原理，北京：清华大学出版社[4]教学用“通信原理”教材10心得体会通过这次的课程设计，我更加理解了眼图的用途，通过眼图可以迅速的定性评价系统的性能。可以观察到麻将干扰和信道噪声等影响因素的情况虫儿快速的判断出系统性能的优劣程度。最佳抽烟时刻是‘眼睛’张开最大的时刻，定时误差灵敏度是眼图斜边的斜率，斜率越大，对应的定时误差越敏感。图的阴影区的垂直高度便是抽样时刻上信号受噪音干扰的畸变程度。图中央的横轴位置对应于判决门限电平，抽样时刻时上下两阴影区的间隔距离之半为噪声容限。此外，又学习了用matlab仿真通信原理实验的有关方法，及常用的相关函数。通过对程序错误的调试，更加熟练的应用软件。通过仿真，更容易理解了通信原理的相关内容。通信原理是一门很有用的课程，但是内容繁杂，需要记忆的很多，所以我们应该加深的理解而不是单纯的记忆，另外如果增加实践的内容，我们将学得更11多，更快！路漫漫其修远兮，吾将上下而求索！