带通抽样定理.doc

_信号与系统A（2）课程自学报告实施报告题 目：带通采样定理与软件无线电带通抽样定理实际中遇到的许多信号是带通型信号，这种信号的带宽往往远小于信号中心频率。若带通信号的上截止频率为，下截止频率为，这时并不需要抽样频率高于两倍上截止频率，可按照带通抽样定理确定抽样频率。定理 带通抽样定理：一个频带限制在内的时间连续信号，信号带宽，令，这里为不大于的最大正整数。如果抽样频率满足条件，（3.1-9）则可以由抽样序列无失真的重建原始信号。对信号以频率抽样后，得到的采样信号的频谱是的频谱经过周期延拓而成，延拓周期为，如图3-3所示。为了能够由抽样序列无失真的重建原始信号，必须选择合适的延拓周期（也就是选择采样频率），使得位于和的频带分量不会和延拓分量出现混叠，这样使用带通滤波器就可以由采样序列重建原始信号。由于正负频率分量的对称性，我们仅考虑的频带分量不会出现混叠的条件。在抽样信号的频谱中，在频带的两边，有着两个延拓频谱分量：和。为了避免混叠，延拓后的频带分量应满足（3.1-10）（3.1-11）综合式（3.1-10）和式（3.1-11）并整理得到（3.1-12）这里是大于等于零的一个正数。如果取零，则上述条件化为（3.1-13）这时实际上是把带通信号看作低通信号进行采样。取得越大，则符合式（3.1-12）的采样频率会越低。但是有一个上限，因为，而为了避免混叠，延拓周期要大于两倍的信号带宽，即。因此（3.1-14）由于为不大于的最大正整数，因此不大于的最大正整数为，故有综上所述，要无失真的恢复原始信号，采样频率应满足， （3.1-15）图3-3 带通采样信号的频谱带通抽样定理在频分多路信号的编码、数字接收机的中频采样数字化中有重要的应用。作为一个特例，我们考虑（）的情况，即上截止频率为带宽的整数倍。若按低通抽样定理，则要求抽样频率，抽样后信号各段频谱间不重叠，采用低通滤波器或带通滤波器均能无失真的恢复原始信号。根据带通抽样，若将抽样频率取为（值取为），抽样后信号各段频谱之间仍不会发生混叠。采用带通滤波器仍可无失真地恢复原始信号，但此时抽样频率远低于低通抽样定理的要求。图3-4所示为，时抽样信号的频谱。图3-4 ，时的抽样频谱在带通抽样定理中，由于，带通抽样信号的抽样频率在到之间变化，如图3-5所示。图3-5 带通抽样定理由以上讨论可知，低通信号的抽样和恢复比起带通信号来要简单。通常，当带通信号的带宽大于信号的最低频率时，在抽样时把信号当作低通信号处理，使用低通抽样定理，而在不满足上述条件时则使用带通抽样定理。模拟电话信号经限带后的频率范围为300Hz3400Hz，在抽样时按低通抽样定理，抽样频率至少为6800Hz。由于在实际实现时滤波器均有一定宽度的过渡带，抽样前的限带滤波器不能对3400Hz以上频率分量完全予以抑制，在恢复信号时也不可能使用理想的低通滤波器，所以对语音信号的抽样频率取为8kHz。这样，在抽样信号的频谱之间便可形成一定间隔的保护带，既防止频谱的混叠，又放松了对低通滤波器的要求。这种以适当高于奈奎斯特频率进行抽样的方法在实际应用中是很常见的。软件无线电所覆盖的频率范围一般都比较宽 ,只有宽频段才能具有广泛的适应性 ,这样宽的频段号的频谱Xs()可表示为 :+ Xs() =1Ts n = -X( -ns)用 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(fs/250); %频率为250HZ for i=1:N if(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); % 查看信号频谱grid on;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(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,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(时