带通抽样定理

1、信号与系统 A（2）课程自学报告实施报告题目：带通采样定理与软件无线电带通抽样定理实际中遇到的许多信号是带通型信号，这种信号的带宽往往远小于信号中心频率。若带通信号的上截止频率为f H ，下截止频率为 f L ，这时并不需要抽样频率高于两倍上截止频率f H ，可按照带通抽样定理确定抽样频率。定理 带通抽样定理 ：一个频带限制在( fL , f H ) 内的时间连续信号x(t ) ，信号带宽 B f Hf L ，令 M f H / BN ，这里 N 为不大于 f H / B 的最大正整数。如果抽样频率f s 满足条件2 f Hf s2 f L ， 0 m N1 （3.1-9）m 1m则可以由抽样

2、序列无失真的重建原始信号x(t ) 。对信号 x(t ) 以频率 f s 抽样后，得到的采样信号 x(nTs ) 的频谱是 x( t ) 的频谱经过周期延拓而成，延拓周期为 f s ，如图 3-3 所示。为了能够由抽样序列无失真的重建原始信号 x( t) ，必须选择合适的延拓周期（也就是选择采样频率），使得位于 ( f L , f H ) 和 ( f H , f L ) 的频带分量不会和延拓分量出现混叠，这样使用 带通滤波器 就可以由采样序列重建原始信号。由于正负频率分量的对称性， 我们仅考虑 ( f L , f H ) 的频带分量不会出现混叠的条件。在抽样信号的频谱中， 在 ( f L ,

3、f H ) 频带的两边，有着两个延拓频谱分量：( f H mf s, f L mf s ) 和 ( f H(m1) f s , f L ( m 1) f s ) 。为了避免混叠， 延拓后的频带分量应满足f Lmf sf L （3.1-10）f H (m1) f sf H （3.1-11）综合式（ 3.1-10）和式（ 3.1-11）并整理得到2 f Hf s2 f Lm1m （3.1-12）这里 m 是大于等于零的一个正数。如果 m 取零，则上述条件化为f s2 f H （3.1-13）这时实际上是把带通信号看作低通信号进行采样。m 取得越大，则符合式（ 3.1-12）的采样频率会越低。但是

4、m 有一个上限，因为 fs 2 f L ，而为了避免混叠，延拓周期要大于两倍的m信号带宽，即fs2B 。因此2 f L2 f Lf Lm2BB （3.1-14）f s由于 N 为不大于 f H / B 的最大正整数，因此不大于f L / B 的最大正整数为 N1 ，故有 0 m N 1综上所述，要无失真的恢复原始信号x(t ) ，采样频率f s 应满足2 f Hfs2 f L ， 0m N 1 （3.1-15）m 1mX ( f )ff Hf L0fLf HX s ( f )f H(m1) f sf Lmfsf Lmf sf H(m1) f sf sf sffH fL0fLf H图3-3 带通

5、采样信号的频谱带通抽样定理在频分多路信号的编码、 数字接收机的中频采样数字化中有重要的应用。作为一个特例，我们考虑f HNB （ N1 ）的情况，即上截止频率为带宽的整数倍。若按低通抽样定理，则要求抽样频率 f s 2 NB ，抽样后信号各段频谱间不重叠， 采用低通滤波器或带通滤波器均能无失真的恢复原始信号。根据带通抽样，若将抽样频率取为 f s 2B （ m 值取为N 1 ），抽样后信号各段频谱之间仍不会发生混叠。采用带通滤波器仍可无失真地恢复原始信号，但此时抽样频率远低于低通抽样定理fs2NB 的要求。图3-4 所示为 f H3B ， f s2B 时抽样信号的频谱。X ( f )f-3 B

6、 -2B02B3BX( ffs)f-B04B 5BX ( ffs)f-5B-4B0BX s( f )f-B0B图 3-4 f H3B ， f s2 B 时的抽样频谱在带通抽样定理中，由于 0M1 ，带通抽样信号的抽样频率在 2B到 4B之间变化，如图3-5 所示。f s5BM/N 131 / 2B1/ 31/ 41 / 52BB0234B56f HBBBBB图 3-5带通抽样定理由以上讨论可知，低通信号的抽样和恢复比起带通信号来要简单。通常，当带通信号的带宽B 大于信号的最低频率f L 时，在抽样时把信号当作低通信号处理，使用低通抽样定理，而在不满足上述条件时则使用带通抽样定理。模拟电话信号经

7、限带后的频率范围为300Hz 3400Hz，在抽样时按低通抽样定理，抽样频率至少为6800Hz。由于在实际实现时滤波器均有一定宽度的过渡带，抽样前的限带滤波器不能对3400Hz以上频率分量完全予以抑制，在恢复信号时也不可能使用理想的低通滤波器，所以对语音信号的抽样频率取为8kHz。这样，在抽样信号的频谱之间便可形成一定间隔的保护带，既防止频谱的混叠，又放松了对低通滤波器的要求。这种以适当高于奈奎斯特频率进行抽样的方法在实际应用中是很常见的。软件无线电所覆盖的频率范围一般都比较宽,只有宽频段才能具有广泛的适应性,这样宽的频段号的频谱Xs()可表示为:+ (1Xs) =X(-Tsn = -n s)

8、用 Nyquist 低通采样 1 是不现实的 ,而实际上软件无线电中的各种无线电信号的瞬时信号带宽都比较窄 ,这样采用带通采样将使采样速率大大降低 。本文重点是通过对带通采样定理进行严格的数学推导和论证 ,利用一系列不等式组提出了带通采样理论上能达到的最小采样速率 。用 matlab 对定理进行仿真fs=100000;%抽样频率N=216;t=(0:N-1)/fs;%t 时间W=(0:N-1)*fs/N;%W 频率y=cos(1200*2*pi*t)+cos(1250*2*pi*t)+cos(1300*2*pi*t);% y 原始信号chyang=zeros(1,N);%抽样信号T=round

9、(fs/250);%频率为 250HZfor i=1:Nif(mod(i,T)=0)chyang(i)=1;endendy2=y.*chyang;%y2 为抽样后的信号hw=fft(y,N);%快速傅里叶变换HW=abs(hw);%取幅度值W=(0:N-1)*fs/N;%频率值figure(1);subplot(2,1,1);plot(t(1:8000),y(1:8000);grid on;title(' 滤波前信号 y');xlabel('时间 /s');% 原始信号subplot(2,1,2);plot(W(1:N/64),HW(1:N/64);% 查看信号

10、频谱gridon;title(' 滤 波 前 信 号 频 谱 图 ');xlabel(' 频 率 /Hz');ylabel(' 振 幅|H(ejw)|');hwch=fft(chyang,N);HWCH=abs(hwch);%取幅度值W=(0:N-1)*fs/N;%频率值figure(2);subplot(2,1,1);stem(t(1:10000),chyang(1:10000);grid on;title(' 抽样信号 chyang');xlabel('时间 /s');subplot(2,1,2);plot(W

11、(1:N/64),HWCH(1:N/64);agrid on;title(' 抽 样 信 号 频 谱 图 ');xlabel(' 频 率 /Hz');ylabel(' 振 幅 |H(ejw)|');hw2=fft(y2,N);%快速傅里叶变换HW2=abs(hw2);%取幅度值W=(0:N-1)*fs/N;%频率值figure(3);subplot(2,1,1);plot(t(1:8000),y2(1:8000);grid on;title(' 采样之后 ');xlabel(' 时间 /s');subplot(2,

12、1,2);plot(W(1:N/64),HW2(1:N/64);%设置带通滤波器ws=1130/fs*2,1370/fs*2;%ws1 ws2 阻带截止频率wp=1180/fs*2,1320/fs*2;%wp1 wp2 通带截止频率n,wn=buttord(wp,ws,3,15);%求滤波器阶数 n 与截止频率 wnb,a=butter(n,wn);%求滤波器 H(s)的表达式的分子分母系数 b，ay3=filter(b,a,y2);%y2 为 y 经过滤波器之后的函数hw3=fft(y3,N);HW3=abs(hw3);figure(4);subplot(2,1,1);plot(t(1:8000),y3(1:8000);%axis(0,0.1,-0.02,0.02);grid on;title(' 滤波后信号 ');xlabel(' 时间 /