带通信号的采样与重建(优选)甄选.docx
带通信号的采样与重建#(优选.)带通信号的采样与重建一、带通采样定理的理论基础基带采样定理只讨论了其频谱分布在(O,f)的基带信号的采样问题。作为接收机的模数转换H来说：接收信号大多为已调制的射频信号。射频信号相应的频率上限远高于基带信号的频率上限。这时如果想采用基带采样就需要非常高的采样速率！这是现实中的A/D难以实现的。这时，低通采样定理已经不能满足实际中的使用要求。带通采样定理是适用于这样的带通信号的采样理论基础，下面给出定理。带通采样定理：设一个频率带限信号Mf)其频带限制在(f7)内，如果其采样速率/满足式：1Hsf=2(f+f)(2-1)J/1H式中，取能满足f2(/一7)的最大整数(0,1,2)，则用/进行等间隔采样所得到5H1s的信号采样值x(nT)能准确的确定原信号x(t)。S带通采样定理使用的前提条件是：只允许在其中一个频带上存在信号，而不允许在不同的频带同时存在信号，否则将会引起信号混叠。如图2.3所示，为满足这一条件的一种方案，采用跟踪滤波器的办法来解决，即在采样前先进行滤波，也就是当需要对位于某一个中心频率的带通信号进行采样时，就先把跟踪滤波器调到与之对应的中心频率,上，滤出所感兴趣的带通信号XW，Onn然后再进行采样，以防止信号混叠。这样的跟踪滤波器称之为抗混叠滤波器。X(t)跟踪滤波器采样X(n)FO=(2n+DB/2Fs=2B图2.3带通信号采样式(2-1)用带通信号的中心频率/和频带宽度B也可用式(2-2)表示：0f=22+1/(2-2)04s式中，f=（/+f）/2,取能满足/2B（B为频带宽度）的最大正O1H'S整数。当频带宽带B一定时，为了能用最低采样速率即两倍频带宽度的采样速率（/=2B）,带通2n+1信号的中心频率必须满足九一。也即信号的最高或最低频率是信号的整数倍。带通采样理论的应用大大降低了所需的射频采样频率，为后面的实时处理奠定了基础。但是从软件无线电的要求来看，带通采样的带宽应是越宽越好，这样对不同基带带宽的信号会有更好的适应性,在相同的工作频率范围内所需要的“盲区”采样频率数量减少，有利于简化系统设计。另外，当对于一个的。所以在随着采样速一般为几千据流降速是频率很高的射频信号采样时，如果采样频率设的太低，对提高采样量化的信噪比是不利可能的情况下，带通采样频率应该尽可能选的高一些，使瞬时采样带宽尽可能宽。但是率的提高带来的一个问题是采样后的数据流速率很高。因此一个实际的无线电通信带宽赫兹到几百赫兹。实际对单信号采样时采样率是不高的。所以对这种窄带信号的采样数完全可能的。多速率信号处理技术为这种降速处理实现提供了理论依据。二、带通采样过程待采样信号为中频是IoOMHz,带宽为2MHz的带通信号:fc0=100e6;中频频率fc1=99e;fc2=101e6;信号一的频率信号二的频率fs1=3e;为采样的采样频率fs2=4e6;I1今界采样采样频率fs3=6e6;大于2倍带宽的采样速率f0=250e6加以模拟连续信号的离散信号采样速率（远大于nyquist率）t=0:1/f0:1e-6;xt0=2*cos(2*pi*fc1*t)+4*cos(2*pi*fc2*t);待采样信号的波形和频谱如图所示：艾600o0.511.5数字颉平2.5XIO8图1原信号波形及频谱按照如上的三种采样频率对待采样信号进行采样，得到的三个信号及其幅度谱如图所示:1）欠采样条件下得到的采样信号：Iii20O!111.40.60_811.2×XKM欠采样信号频诣欠采样信号浣形80C600烹娶400I.11,1200I1031523-3数字赁幸yin图2欠采样信号的波形及幅度谱2）临界采样:啦界采样信号敕形50书XX10.4I倡界采样信号颌谱1500Otr1.522.5数字颉字3.54白X1O图3临界采样信号的波形及幅度谱3）满足采样条件的采样信号大于临界采样信号波形大于临界采样信号频谱III0123456数字频率X1O6入IV图4正常采样信号的波形及幅度谱三、信号重建设计通带为99IOIMHZ的4阶巴特沃兹一型HR滤波器来作为信号恢复的模拟滤波器，滤波器的设计过程及幅频响应特性如图所示：fs=250e;N=4;figure(6);Wn=0.7920.808;bza=butter(N,Wn,'bandpass,);h,w=freqz(bza);p1ot(wpi*fs2,abs(h);grid;tit1e('Amp1itudeResponse');x1abe1('Frequency(Hz),);y1abe1('Amp1itude');AmpitudeRespons仑1-''T.J-r-1.-111：1-»1<I>(1-*".1，-tz1Inso*U-,.ftI1_.一.Z.I；_-.一，IiiiPIII-0200?"15A1214FreqruencyHiz)图1滤波器的幅度响应将三种采样信号的到的数据流通过信号恢复滤波器，即上述的巴特沃兹一型滤波器，得到的结果如下：欠采样恢复OJ0.080.06i004-0.02.0.04-006-0.0,OI,八02040-6080100120140160图2欠采样恢复临界采样恢复0.060.040.02-004-0.06-oo50100150200250图3临界采样恢复正常采样恢复图3正常采样恢复四、结果分析从结果可见：欠采样信号进行恢复时，从信号完整角度进行观察便可以发现信号的不同，缺少了原信号的部分信息。而临界采样信号的恢复，在信号周期上有一定的变化，比原信号的信息有所丢失。正常采样下的信号可以恢复出原信号一个完整周期的频谱。但由于所选取信号长度的问题，没能恢复全部信号，但这已经能够证明带通采样定理的正确性，即采样后的信号没有丢失原信号的信息。附：%z0t26rJu,AiAt1eO0%z0i-ODDAAetafC=IOOMHz,DAoAAetaf1=99MHzotf2=101MHz%26UufiOAO±ta3MHzi¢4MHZj6MHz%B=2MHz,AiA06I±26f3uffiAziS1f±SIz0iifi±£oDAoAAJEffi×2>>>>a>>ip;£%c1eara11;c1c;c1osea11;fc0=100e6;fc1=99e6;fc2=101e6;%AoA1A1Afi%DAoA2AJEAfifs1=3e6;如f60q±q'0"演2网山即Acfs2=4e6;fs3=6e;%A026b)uffiAfi%'66(Hp±9'uiiA26b1uffiAfif0=250e;%±i¾-DAoAA2fcb1uffAfi1-.±JEa26f)uffiAt1d>>,fiyi3I%tt,26fiu1AiX26f1u×÷ta0-AiAaDAoA%t=0:1/f0:1e-6;N=1e-6*f0;xt0=2*cos(2*pi*fc1*t)+4*cos(2*pi*fc2*t);figure(1);subp1ot(2,1,1);p1ot(t,xt);tit1e(''y26fiuDAoA2Dt,);x1abe1(,x,);y1abe1(,xt,);yjw=fft(xt,N);absy=abs(yjw);f=(0:N-1)*fON;subp1ot(2,1,2);p1ot(f,absy);tit1e('zy26fiuDAoAffi×,)x1abe1(,Ey×OAE,);y1abe1(,uR×,)%IAA22正口龙X%t=0:1/fs1:5e-5;N=5e-5*fs1;xts1=2*cos(2*pi*fc1*t)+4*cos(2*pi*fc2*t);figure(2);subp1ot(2,1,1);p1ot(t,xts1);tit1e(,QdJ)；x1abe1(,x,);y1abe1(,xt');yjws1=fft(xts1,N);absys1=abs(yjws1);f=(O:N-1)*fs1N;subp1ot(2,1,2);p1ot(f,absys1);x1abe1(,fiy×OJEAfi,);y1abe1(,Mffi×,)tit1e(,Q26f1uDAoA1×,)%,0I"2E1AiiAAUbQ2ERuyA2-DI¼o×%t=0:1fs2:5e-5;N=5e-5*fs2;xts2=2*cos(2*pi*fc1*t)+4*cos(2*pi*fc2*t);figure(3);subp1ot(2z1,1);p1ot(tzxts2);tit1e(,A0926RuDAoA2-Di,);x1abe1('x');y1abe1(,xt,);yjws2=fft(xts2,N);absys2=abs(yjws2);f=(0:N-1)*fs2N;subp1ot(2,1,2);p1ot(fzabsys2);x1abe1(,fiy×OffiAfi');y1abe1(',)tit1e('AU2ESuDAoA×,);%,OUAU2ESuAEAz01"2EuA2"DI¼0比无%t=0:1/fs3:5e-5;N=5e-5*fs3;xts3=2*cos(2*pi*fc1*t)+4*cos(2*pi*fc2*t);figure(4);subp1ot(2,1,1);p1ot(tzxts3);tit1e(,60UAU2ERuDA0A2"DI,);x1abe1('x,);y1abe1(,xt');yjws3=fft(xts3,N);absys3=abs(yjws3);f=(0:N-1)*fs3N;subp1ot(2,1,2);p1ot(fzabsys3);x1abe1(,Ey×OSAE');y1abe1(,1Effi×')tit1e(
