《通信原理实验评测报告》由会员分享，可在线阅读，更多相关《通信原理实验评测报告（26页珍藏版）》请在金锄头文库上搜索。
1、通信原理实验报告实验名称： 实验一数字基带传输系统的MATLAB方真实验二模拟信号幅度调制仿真实验班级：10通信项目三班学 号： 2018550920姓 名： 彭龙龙指导老师： 王仕果实验一数字基带传输系统的MATLA仿真一、实验目的1、 熟悉和掌握常用的用于通信原理时域仿真分析的MATLAB函数；2、掌握连续时间和离散时间信号的 MATLAB产生；3、牢固掌握冲激函数和阶跃函数等函数的概念，掌握卷积表达式及其物理意义，掌握 卷积的计算方法、卷积的基本性质；4、 掌握利用MATLAB计算卷积的编程方法，并利用所编写的MATLAB程序验证卷积 的常用基本性质；5、掌握MATLAB描述通信系统中不
2、同波形的常用方法及有关函数，并学会利用 MATLAB求解系统功率谱，绘制相应曲线。基本要求：掌握用 MATLAB描述连续时间信号和离散时间信号的方法，能够编写 MATLAB程序，实现各种常用信号的MATLA实现，并且以图形的方式再现各种信号的波形。二、实验内容1、编写MATLAB程序产生离散随机信号2、编写MATLAB程序生成连续时间信号3、编写MATLAB程序实现常见特殊信号三、实验原理从通信的角度来看，通信的过程就是消息的交换和传递的过程。而从数学的角度来 看，信息从一地传送到另一地的整个过程或者各个环节不外乎是一些码或信号的交换过 程。例如信源压缩编码、纠错编码、AMI编码、扰码等属于码
3、层次上的变换，而基带成形、滤波、调制等则是信号层坎上的处理。码的变换是易于用软件来仿真的。要仿真信号 的变换，必须解决信号与信号系统在软件中表示的问题。3.1信号及系统在计算机中的表示时域取样及频域取样一般来说，任意信号 s（t是定义在时间区间-+R）上的连续函数，但所有计算机段。为此将把S(t按区间 T, T 截短为ST(t，再对ST(t按时间间隔 t均匀取样，得IL 22到取样点数为：TNt(3-1At仿真时用这个样值集合来表示信号s(t。显然 t反映了仿真系统对信号波形的分辨率， t越小则仿真的精确度越高。据通信原理所学，信号被取样以后，对应的频谱时频率1的周期函数，其重复周期是 。如果
4、信号的最高频率为At证不发生频域混叠失真。设fH,那么必须有Bs12：t则称Bs为仿真系统的系统带宽。如果在仿真程序中设定的采样间隔是1fH三一才能保2.：t(3-2 t，那么不能用此仿真程序来研究带宽大于 Bs的信号或系统。此外，任意信号s(t的频谱S(f通常来说也是定义在时间区间-8, + 8)上的连续函数，所以仿真频域特性时，也必须把S(f截短并取样。考虑到系统带宽为截短区间设计为-Bs, Bs然后再按间隔 f均匀取样，得到取样点数为：Nf2Bs :fBs，便把频谱的(3-3将式3-2)代入式3-3 )得Nf1 ：t：f(3-41同样，信号在频域被离散后，对应到时域也是一个周期信号，其周
5、期为。如果时Af域截短时间为 T，那么必须 T 1/df才能保证不发生频域混叠失真。也就是说，如果仿真程序中设定得频域采样间隔是 f，那么就不能仿真截短时间超过的信号。所以，可以(3-5把频域的取样间隔设计为：将式3-5)代入式3-1 )得Nt(3-61N =。要提高仿也也f这样一来，时域的总取样点数及频域的总取样点数都相等，为真的精度，就必须降低时域取样间隔 t及频域取样间隔 f，也就是要加大总取样点数N。这说明仿真的精度与仿真系统的运算量直接有关。为了处理上的方便，我们今后规定采样点数N为2的整幕。举例来说，例如设计要求2Bs的系统带宽为1MHz，频域最好分辨率为 10kHz，那么据此可求
6、得 Nf s= 200，则取AfN = 256。对应的其 他参数为Bs = 1MHz , f二绝 =7.8kHz , T二丄=128S ,NAfit0.51 s oN频域分析限于篇幅，在此不介绍MATLAB中关于傅立叶变换的有关函数。为了方便仿真，我们利用MATLAB提供的函数编写了两个函数t2f和f2t t2f的功能是做傅立叶变换，f2t的功能是做傅立叶反变换，它们的引用格式分别为X=t2f(x及x=f2t(X，其中x是时域信号x(tj截短并采样所得的取样值矢量，X是对x(t的傅立叶变换 X(f截短并采样所得的取样值矢量。这两个函数分别如附录所示。我们关心的另一个指标是信号的功率谱密度，任意
7、信号s(t的功率谱的定义是PS(f) TimT_SC，其中ST(f)是s(t 截短后所得信号sr(t)的傅立叶变化,Sr(f)2是st (t)的能量谱,是sr(t)在截短时间T内的功率谱。对于MATLAB 仿真系统，若 x是时域取样值矢量,X是对应的傅立叶变换，那么x的功率谱便为矢量P =(X.conj(X)/T o3.2与随机信号产生相关的指令高斯噪声的产生Bs,并且各样点的值互不因为函数randn(1,N产生N个互不相关的、均值为零、方差为 1的高斯随机数，所以 可用它来产生高斯白噪声。设仿真系统的取样点数是N ,系统带宽为Bs,矢量x二n0* Bs* randn(1,N)的总功率为n0*
8、 Bs，最高频率分量为相关，故它代表双边功率谱密度为W/Hz )的白噪声。23.2.2 随机码序列的产生语句round(rand(1,M产生M个取值1、0等概的随机码。函数round表示四舍五入。函数rand产生均匀分布于区间0, 1的随机数。语句sign(rand(1,M产生M个取值土 1等概的随机码。函数 sign(x对矢量x的元素取 正负号，而高斯数 randn取正负数的概率是相等的。产生数字随机信号的一般方法一般来说，随机数字信号可以直接或间接表示成PAM信号请参阅通信原理教材)。PAM信号是指所有形如s(t)二aig(t -iTs)的信号。它可以用如下的等效模型来表示：i =0wim
9、p(f)二 丫叫灭_码)Gif故对所有不同的数字信号，都可以用相同的方法来产生。1) 产生随机序列矢量 aa=round(rand(1,M。为了方便起见，一般规定a的长度M是2的整幕。2) 产生冲击序列信号imp(timp=zeros(1,N。imp(1:L:M=a/dt 。矢量imp代表信号imp(t。L是每码元内的采样点其中，N是imp的矢量长度，M是码元矢量a中的码元数, 数。现规定M、N都是2的整幕，于是L自然也是2的整幕。3) 产生PAM数字信号s(ts(t)二 imp(t) ： g(t)g(t)代表脉冲波形，例如，升余弦滚降要求的波形有:g(t)=sin 二t/Ts t /Tsco
10、s(a 二t/Ts)1 -(2a 二 t/Ts)2设g(t)的傅立叶变换为 G(f)，s(t)的傅立叶变换为S(f)，则有：S( f)二t2 f (imp). * G%t2f为傅立叶变换函数s(t) = f 2t(S( f) %f2t为傅立叶反变换函数设矢量s代表数字信号s(t)，矢量g代表脉冲波形 g(t)，矢量 G代表其频谱G( f)。那么s的产生方法是:s=conv(imp,g。其中，函数 conv表示卷积。卷积后 s的长度是length(imp+length(g-1。扣除延迟 时间及拖尾时间后，数字信号为：Ii=find(g=max(g 。s=s(1:N+ii(1。也可用频域的方法产生
11、数字信号s(t):s=f2t(t2f(imp.*G。注意，此时imp的点数应与g或G相同。若g的宽度小于imp，则应用零补齐。3.3信号的仿真连续时间信号的仿真如前所述，MATLAB有很多内部数学函数可以用来产生这样的数字序列，例如sin(、cos(、exp(等函数可以直接产生一个按照正弦、余弦或指数规律变化的数字序列。例 如，运行如下程序%program1_1% This program is used to gen eratea sinu soidalsig nal and draw itsplotclear, % Clear all variablesclose all, % Close
12、 all figure win dowsdt = 0.01。 % Specify the step of time variablet = -2:dt:2。 % Specify the interval of timex = sin( 2*pi*t。 % Gen erate the sig nalplot(t,x % Open a figure window and draw the plot of x(t title(Sinusoidal signal x(txlabel(Time t (sec在通信原理课程中，单位阶跃信号u(t和单位冲激信号 S (t是二个非常有用的信号。它们的定义如下：
13、J 方(M = 1r=-iE= 0.心 0t0/dt = 0.01。y = (u(t-u(t-dt/dt 。产生单位阶跃信号的扩展函数为：% Un it step functionfunction y = u(ty = (t=0。% y = 1 for t 0, else y = 0请将这二个MATLAB函数分别以delta和u为文件名保存在 work文件夹中，以后，就可 以像教材中的方法使用单位冲激信号S (t和单位阶跃信号u(t。离散时间信号的仿真程序Program1_2用来产生离散时间信号xn=sin(0.2n n% Program1_2% This program is used to
14、 gen erate a discrete-time sinu soidal sig nal and draw its plotclear, % Clear all variablesclose all, % Close all figure win dowsn = -10:10。 % Specify the in terval of timex = sin( 0.2*pi* n。 % Gen erate the sig nalstem (n,x % Open a figure window and draw the plot of xntitle (Sinusoidal signal xnxlabel (Time index n请仔细阅读该程序，比较程序Program1_1和Program1_2中的不同之处，以便自己编程时能够正确使用这种方法方针连续时间信号和离散时间信号。程序Program1_3用来仿真下
