带通采样定理证明

推荐精选推荐精选带通信号的采样与重建一、带通采样定理的理论基础基带采样左理只讨论了其频谱分布在（0,/〃）的基带信号的采样问题。作为接收机的模数转换来说：接收信号大多为已调制的射频信号。射频信号相应的频率上限远高于基带信号的频率上限。这时如果想采用基带采样就需要非常髙的采样速率！这是现实中的A/D难以实现的。这时，低通采样能理已经不能满足实际中的使用要求。带通采样立理是适用于这样的带通信号的采样理论基础，下面给出左理。带通采样立理：设一个频率带限信号兀⑴其频带限制在（厶,几）内，如果其采样速率人满足式：「迭屮 （2-1）式中，〃取能满足fx>2（fH-fL）的最大整数（0,1,2…），则用£进行等间隔釆样所得到的信号采样值双川7；）能准确的确左原信号x（/）o带通采样左理使用的前提条件是：只允许在苴中一个频带上存在信号，而不允许在不同的频带同时存在信号，否则将会引起信号混叠⑴。如图2.3所示，为满足这一条件的一种方案，采用跟踪滤波器的办法来解决，即在采样前先进行滤波⑴，也就是当需要对位于某一个中心频率的带通信号进行采样时，就先把跟踪滤波器调到与之对应的中心频率几上，滤出所感兴趣的带通信号俎（/）,然后再进行采样，以防止信号混叠。这样的跟踪滤波器称之为抗混叠滤波器。•X(t)X(n)F0=(2n-l)B/2Fs=2B•X(t)X(n)F0=(2n-l)B/2Fs=2B图2.3带通信号采样式（2J）用带通信号的中心频率九和频带宽度B也可用式（2-2）表示:f.=^fx （2-2）式中，九=（£+九）/2，"取能满足fs>2B（B为频带宽度）的最大正整数。当频带宽帯B一立时，为了能用最低采样速率即两倍频带宽度的采样速率（fs=2B）,带通信号的中心频率必须满足f严。也即信号的最髙或最低频率是信号的整数倍。2带通采样理论的应用大大降低了所需的射频采样频率，为后而的实时处理奠左了基础。但是从软件无线电的要求来看，带通采样的带宽应是越宽越好，这样对不同基带带宽的信号会有更好的适应性，在相同的工作频率范囤内所需要的“盲区”采样频率数量减少，有利于简化系统设讣。另外，当对于一个频率很高的射频信号采样时，如果采样频率设的太低，对提高采样量化的信噪比是不利的。所以在可能的情况下，带通采样频率应该尽可能选的髙一些，使瞬时采样带宽尽可能宽。但是随着采样速率的提髙带来的一个问题是采样后的数据流速率很髙。因此一个实际的无线电通信带宽一般为几千赫兹到几百赫兹。实际对单信号采样时采样率是不高的。所以对这种窄带信号的采样数据流降速是完全可能的。多速率信号处理技术为这种降速处理实现提供了理论依拯。二、带通采样过程待采样信号为中频是100MHz,带宽为2MHz的带通信号:fc0=100e6;//中频频率fcl=99e6;//信号-的频率fc2=101e6;//信号二的频率fsl=3e6;//欠采样的采样频率fs2=4e6;//临界采样采样频率fs3=6e6;//大于2倍带宽的采样速率f0=250e6//用以模拟连续信号的离散信号采样速率（远人于nyquist率)t=0:l/f0:le-6;xt0=2*cos(2*pi*fcl*t)+4*cos(2*pi*fc2*t);待采样信号的波形和频谱如图所示：待采样信号波形待采样信号频谱图1原信号波形及频谱按照如上的三种采样频率对待采样信号进行采样，得到的三个信号及其幅度谱如图所示:1）欠采样条件下得到的釆样信号：欠采样信号波形欠采样信号频谱图2欠采样信号的波形及幅度谱2）临界釆样:临界采样信号波形x10*4临界采样信号频谱150010005000.513.5大于临界采样信号波形图34x10-%1.5 2 2.6数字频率临界采样信号的波形及幅度谱3）满足釆样条件的采样信号大于临界采祥信号频谱0数字频率響1000遞豐500x1016000 2 4 6 8 10 12 14Frequency(Hz) x100 2 4 6 8 10 12 14Frequency(Hz) x1070>pn三dll<图4正常采样信号的波形及幅度谱三、信号重建1)设讣通带为99〜101MHz的4阶巴特沃兹一型IIR滤波器来作为信号恢复的模拟滤波器，滤波器的设计过程及幅频响应特性如图所示：fs=250e6;N=4;figure(6);Wn=[0.7920.808];[b,a]=butter(N,Wnz*bandpass*);[h,w]=freqz(b,a)；plot(w/pi*fs/2,abs(h));grid;title(*AmplitudeResponse1);xlabel('Frequency(Hz)1);ylabel('Amplitude*);AmplitudeResponse141210.80.60.40.20图1滤波器的幅度响应2)将三种采样信号的到的数据流通过信号恢复滤波器，即上述的巴特沃兹一型滤波器，得到的结果如下：0.108

nW06040200204006080■06

o08o■020040060-040.108

nW06040200204006080■06

o08o■020040060-04002o欠采样恢复60 80 100 120 140 160图2欠采样恢复临界采样恢复100 150t200250图3临界采样恢复正常采样恢复图3正常采样恢复四、结果分析从结果可见：欠采样信号进行恢复时，从信号完整角度进行观察便可以发现信号的不同,缺少了原信号的部分信息。而临界采样信号的恢复.在信号周期上有一上的变化，比原信号的信息有所丢失。正常采样下的信号可以恢复岀原信号一个完整周期的频谱。但由于所选取信号长度的问题，没能恢复全部信号，但这已经能够证明带通采样定理的正确性，即采样后的信号没有丢失原信号的信息。附：%z0l"-ENuSI'AipANedn%%z0l'O99A5；|xAEiafc=100MHz,9A%EpAEIafl=59MHz°ff2=101MHz%%-ENuEUAE・O±6Ia3MHz;<?4MHz;$6MHz%%%%%B=2MHz,AiAUEIp±2ENu5；pAE'工卩込雯Q士贸z0^iE±£o9A°ApA.Ep^x：»»a»ip^j£%%%%clearall;clc;closeall;fc0=100e6;fcl=99e6;%9A°Al)iA5；pAEfc2=101e6;%BA°A2pA3：pAEfsl=3e6;%pf6uSt>±5'o^ipA-ENu^pAEfs2=4e6;%AU^c-ENuZpAEfs3=6e6;%z66u5t>±5zu^2.pAcENuZpAEf0=250e6;%±iE^O-DA°ApA-ENuSpAE±仓逛玄：总舖1臣泊菱色>B0»,oEyA^M^I%%%pII"-ENu5'A£pA2ENux^iao-A乞鸟飭人。航沁t=0:l/f0:le-6;N=le-6*f0;xt0=2*cos(2*pi*fcl*t)+4*cos(2*pi*fc2*t);figure(1);subplot(2z1z1);plot(t,xtO);title('zy-ENu9A°A-DI1);xlabml(1x1);ylabml('xt1);yjwO=fft(xtOzN);absyO=abs(yjwO);f=(0:N-l)*fO/N;subplot(2,1,2);plot(f,absyO);title(1'y-ENuDA%E^x•)xlabel('EyxO.EyxAE');ylabel「・u^E^x1)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%z0l■■-ENuSI'AiIApA^・-ENu-t=0:l/fsl:5e-5;N=5e-5*fsi;xtsl=2*cos(2*pi*fcl*t)+4*cos(2*pi*fc2*t);figure(2);subplot(2,1z1);plot(tzxtsl);title(‘g・-ENuQA°A-DI1);xlabel('x1);ylabel(1xt1);yjwsl=fft(xtslzN);absysl=abs(yjwsl);f=(O:N-l)*fsl/N;subplot(2,1,2);plot(f,absysl);xlabel('EyxO.EyxAE');ylabel「・u^E^x1)title('9•-ENu9A%2p5；x•)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%z0l'-ENu^I'AilAAU^9-ENu)iA-t=0:l/fs2:5e-5;N=5e-5*fs2;xts2=2*cos(2*pi*fcl*t)+4*cos(2*pi*fc2*t);figure(3);subplot(2,1,1);plot(t,xts2);title(1AU^9£ENuDA°A-DI');xlabel('x1);ylabel(1xt1);yjws2=fft(xts2,N);absys2=abs(yjws2);f=(0:N-l)*fs2/N;subplot(2,1,2);plot(f,absys2);xlabel('EyxO^pAE');ylabel「・u^lE.Ex')title(1AU^92ENu9A%E^Ex•);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%zoOUAU^92ENu^pAEpAz0l'2ENuyiA-'91^°疋p臣x%%%t=0:1/fs3:5e-5;N=5e-5*fs3;xts3=2*cos(2*pi*fcl*t)+4*cos(2*pi*fc2*t);figure(4);subplot(2z1z1);plot(t,xts3);title('zo6uAU^9cENuDAoA-DI1);xlabel(•x1);ylabel('xt1);yjws3=fft(xts3zN);absys3=abs(yjws3);f=(0:N-l)*fs3/N;subplot(2,1,2);plot(fzabsys3);xlabel(1EyxO^yxAE');ylabel(1・uSlEaEx•)title(1zoOUAU^9-ENu9A%Ep5；x1);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fs=250e6;N=4;figure(6);Wn=[0.7920.808];[b,a]=butter(N,Wn,1bandpass*);[h,w]=freqz(b,a)；plot(w/pi*fs/2,abs(h));grid;title(*AmplitudeResponse*);xlabel(*Frequency(Hz)*);ylabel(*Amplitude*);%%%9・-ENuDA°ApA»d/%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%figure(7);yretsl=filter(b,azxtsl);plot(yretsl);title('g・£ENu»0>'')%%%%%AU^9-ENuBA°ApA»0z%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%figure(8);yrets2=filter(