通信原理课程设计——模拟信号的抽样

1、模拟信号的抽样过程仿真、基本原理1.1抽样定理抽样时时间上连续的模拟信号变成一系列时间上离散的抽样序列的过程。抽样定理 要解决的事，能否由此抽样序列无失真的恢复出模拟信号。对一个频带受限的、时间连续的模拟信号抽样，当抽样速度达到一定的数值时，那 么根据它的抽样值就能无失真恢复原模拟信号。也就是说，若要传输模拟信号，不一定 要传输模拟信号本身，只需要传输由抽样得到的抽样即可。因此，抽样定理是模拟信号 数字化的理论依据。抽样的过程是将输入的模拟信号与抽样信号相乘，通常抽样信号时一个周期为t的周期脉冲信号，抽样后得到的信号称为抽样信号。理想的抽样为6mt)= a t-nts in j其中 6(t )

2、 = 0：:，fs =工称为抽样速率。因此抽样后信号为tsqqxs t = x tt 八 x nts 二 t -ntsk 二二二1.2低通抽样定理一个频带限制在(0，fh)内的连续信号x(n,如果以ts=1(2fh)秒的时间间隔对 它进行均匀抽样，则x(t %各被所得到的抽样值完全确定，可以由抽样值序列无失真地重 建原始彳t号。ts=1 (2fh”抽样的最大间隔，称为奈奎斯特间隔。低通信号的抽样可以从频域来理解，抽样的时域、频域对照如图4-3所示，根据频域卷积定理，x(t)的频域表达式为iiixs = - x 产. = 一 x- -nsx -ns2 二ts il n 品，,tsn=二由上式可见

3、，在 露的整数倍(n=1,2,)处存在x )的复制谱。如图4-3 (c) 所示，抽样后信号的频谱是原信号频谱平移nfs后叠加而成，因此如果不发生频谱重叠，可以通过低通滤出原信号。如果抽样频率8s 2eh ,即抽样间隔ts a142 fh ),则抽样信号的频谱会发生混叠 现象，此时不可能无失真地重建原始信号。x 。滤波器的作用等效将抽样后的信号 xs (切)通过介质频率为6 h的低通滤波器，只允许低于缶h的频率分量通过，滤除过高的频率分量，从而恢复出原来被抽样的信号于用一门函数与 xs1、斗目乘。低通滤波器的特性如图43 (c)上虚线表示。在时域上就是与冲激响应ht)作运算。即1 .，xt =h

4、 t xs tf nts sa-h t - nts 1tsn口式中，抽样信号sat)= nt vt就是h(t 也就是h(8)的傅里叶逆变换。图43 (d)从几何意义上来说,以每个抽样值为峰值画一个sa函数的波形，则合成的波形就是x(t )。图43低通抽样的时域、频域对照1.3、 带通抽样定理一个带通信号 m (t),其频率限制在fl与fh之间，带宽为 b=fh-fl ,如果最小抽样 速率f=2fh/m , m是一个不超过fh /b的最大整数，那么m (t)可以完全由其抽样值确定。下面分两种情况加以说明：(1)若最高频率fh为带宽的整数倍，即 fh =nb。此时fh/b=n是整数，m=n所以 抽

5、样速率&=2fh/m=2b。图6-7画出了 fh=5b时的频谱图。图中抽样信号的频谱 ms3 )既没有混叠也没有留空隙，而且包含有m(t)的频谱m(co),如图6-7 (c)中两对虚线所框的部分所示。这样，采用带通滤波器就能无失真恢复原信号，且此时抽样速率(2b)远低于按低通抽样定理时 =10b的要求。显然，若fs再减小，即fs2b时必然会出现混叠失真。由此可知：当 fh =nb时，能重建原信号 m(t) 的最小抽样频率为 fs=2b 。(2)若最高频率fh不为带宽的整数倍，即 fh =nb+kb,0k1此时, 定理知m是一个不超过n+k的最大整数，显然，m=n,所以能恢复出原信号fh /b=

6、n+k ,由m(t)的最小抽样速率为叼。m式中，n是一个不超过nfh/b的最大整数，0kb。由于带通信号一般为窄带信号，容易满足flb,因此带通信号通常可按2b速率抽样。、设计过程2.1、 matlab 编程简介：matlab各高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而 被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作，而且 利用matlaet品的开放式结构matlab 特点 ：(1)高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来；(2)具有完备的图形处理功能，实现计算结果和编程的可视化 ；(3)友好的用户界面及接近数学表达式的

7、自然化语言，使学者易于学习和掌握；(4)功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等)，为用户提供了大量 方便实用的处理工具。2.2、 计思路模拟信号首先被抽样，通常抽样是按照等时间间隔进行的，模拟信号被抽样后，成为抽样信号，它在时间上是离散的但是其取值仍然是连续的，所以是离散模拟信号。在理论上抽样过程可以看成是用周期性单位冲激脉冲和此模拟信号相乘。抽样结果得到的是一系列周期性的冲激脉冲，其面积和模拟信号的取值成正比。下面举例说明：给出某一低通信号如：x(t)=0.5+0.1cos0.15 兀+1.5sin2.5 兀 t+0.5cos4 % t(1)画出该低通信号波形；(2)画出抽样速

8、率分别为2hz、4hz、6hz时的抽样序列；(3)从抽样序列恢复出原始信号。分析：此次实验是根据要求先用matlab产生一个的x(t)=0.5+0.1cos0.15 兀+1.5sin2.5兀t+0.5cos4 nt的连续信号，然后对该信号进行抽样，产生频谱。抽样频率不同，即单位时 间内的抽样点数不同，会产生不同的频谱， 再设计一个低通滤波器，使抽样后的信号经过此低通滤波器，恢复出原始信号，由于抽样频率不同会恢复出不同的原始信号。2.3、 流程图2.4、 仿真程序第一步：产生低通信号波形t=0:pi/10:30;x2=0.5+0.1*cos(0.15*pi)+1.5*sin(2.5*pi*t)+

9、0.5*cos(4*pi*t);figureplot(t,x2);xlabel(t);ylabel(x(t);title(原时域连续信号);grid第二步：分别以2hz,4hz,6hz的频率对信号进行抽样n1=input(请输入采样点数n:);n=0:n1;zb=size(n);figuret=30*n/zb(2);sinf=0.5+0.1*cos(0.15*pi*t)+1.5*sin(2.5*pi*t)+0.5*cos(4*pi*t);stem(n,sinf,.);xlabel(n);ylabel(x(n);title(采样后的时域信号);gridn的值为采样点数。n=60时为20hz,n=

10、120为4hz,n=180为6hz第三步：设计低通滤波器b,a=butter(8,350/500);h,w=freqz(b,a,512,2000);figure;plot(w*2000/(2*pi),abs(h);xlabel(hz);ylabel(频率响应幅度)；title(低通滤波器);grid；第四步：根据抽样后的信号恢复原信号n1=input(请输入采样点数 n:)n=0:n1;zb=size(n);figuret=30*n/zb(2);sinf=0.5+0.1*cos(0.15*pi*t)+1.5*sin(2.5*pi*t)+0.5*cos(4*pi*t);subplot(2,1,1

11、);stem(n,sinf,.);xlabel(n);ylabel(x(n);title(采样后的时域信号);gridy=filter(b,a,sinf);subplot(2,1,2);plot(y);xlabel(t);ylabel(x(t);title(恢复后的连续信号);grid;输入n的值，输出抽样n的抽样波形和对应的经过低通滤波器的恢复波形。2hz时信号失真严重，4hz时能基本恢复波形，6hz时能很好的恢复波形，分别见图别见图时，6hz时能三、仿真结果第一步：n=60n=120bq fiqurci 3，=，回 .钾n=180第三步:第四步:izzin=60目n=120re zlh f

12、igure zt-r回 3n=180四、仿真结果分析(1)第一步：产生要求的信号波形，横坐标表示时间t,长度为30,描点间隔为pi/10 ,纵坐标表示幅度x2,幅度与时间的关系为：x2=0.5+0.1cos(0.15*pi)+1.5*sin(0.5*pi*t)+0.5*cos(4*pi*t)(2)第二步；分别以2hz,4hz,6hz的频率对信号进行抽样，程序中 n为程序运行后手动输入，n的值为采样点数。n=6时为2hz,n=120为4hz,n=180为6hz.(3)第三步：设计出低通滤波器，利用 b,a=butter(n,wn,ftype)一个阶数为n频率为 wn的低通滤波器其中参数ftype

13、的形式确定了滤波器的形式当它为high时得到高通滤波器若wn一个含有两个元素向量w1 w2则返回的a,b所构成的滤波器是阶数为2n的带通滤波器滤波器的通带范围是w1wnw2.h,w=freqz(b,a,n)b 和a分别为离散系统的系统函数 分子，分母多项式的系统向量，返回量 h则包含了离散系统频响在 0pi范围内n个频率等 分点的值(其中n为正整数)，w则包含了范围内n个频率等分点。调用默认的 n时，其值 是 512.(4)第四步：根据抽样信号后的信号恢复原信号，输入n的值，输出抽样n的抽样波形和对应的经过低通滤波器输出的恢复波形 2hz时信号失真严重，4hz时基本能恢复波形，6hz 时能很好地恢复波形。五、总结这次设计中，程序的调试遇到很多问题，从基本的函数图形的产生，到后来的抽样产生频谱，再到后来的设计低通滤波器，到最后的信号恢复，从无知到熟知，这个过程不仅巩固了 matlab的应用，更让我对信号的产生，抽样处理，低通滤波器的设计实现，及信号的 恢复有了