带通抽样定理.doc

该【带通抽样定理 】是由【wawa】上传分享，文档一共【10】页，该文档可以免费在线阅读，需要了解更多关于【带通抽样定理 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。《信号与系统A(2)》课程自学报告实施报告题目:带通采样定理与软件无线电带通抽样定理实际中遇到的很多信号是带通型信号,这种信号的带宽往往远小于信号中心频率。若带通信号的上截止频率为,下截止频率为,这时并不须要抽样频率高于两倍上截止频率,可依据带通抽样定理确定抽样频率。[定理]带通抽样定理:一个频带限制在内的时间连续信号,信号带宽,令,这里为不大于的最大正整数。假如抽样频率满意条件,(-9)则可以由抽样序列无失真的重建原始信号。对信号以频率抽样后,得到的采样信号的频谱是的频谱经过周期延拓而成,延拓周期为,如图3-3所示。为了能够由抽样序列无失真的重建原始信号,必需选择合适的延拓周期(也就是选择采样频率),使得位于和的频带重量不会和延拓重量出现混叠,这样运用带通滤波器就可以由采样序列重建原始信号。由于正负频率重量的对称性,我们仅考虑的频带重量不会出现混叠的条件。在抽样信号的频谱中,在频带的两边,有着两个延拓频谱重量:和。为了避开混叠,延拓后的频带重量应满意(-10)(-11)综合式(-10)和式(-11)并整理得到(-12)这里是大于等于零的一个正数。假如取零,则上述条件化为(-13)这时事实上是把带通信号看作低通信号进行采样。取得越大,则符合式(-12)的采样频率会越低。但是有一个上限,因为,而为了避开混叠,延拓周期要大于两倍的信号带宽,即。因此(-14)由于为不大于的最大正整数,因此不大于的最大正整数为,故有综上所述,要无失真的复原原始信号,采样频率应满意,(-15)图3-3带通采样信号的频谱带通抽样定理在频分多路信号的编码、数字接收机的中频采样数字化中有重要的应用。作为一个特例,我们考虑()的状况,即上截止频率为带宽的整数倍。若按低通抽样定理,则要求抽样频率,抽样后信号各段频谱间不重叠,采纳低通滤波器或带通滤波器均能无失真的复原原始信号。依据带通抽样,若将抽样频率取为(值取为),抽样后信号各段频谱之间仍不会发生混叠。采纳带通滤波器仍可无失真地复原原始信号,但此时抽样频率远低于低通抽样定理的要求。图3-4所示为,时抽样信号的频谱。图3-4,时的抽样频谱在带通抽样定理中,由于,带通抽样信号的抽样频率在到之间改变,如图3-5所示。图3-5带通抽样定理由以上探讨可知,低通信号的抽样和复原比起带通信号来要简洁。通常,当带通信号的带宽大于信号的最低频率时,在抽样时把信号当作低通信号处理,运用低通抽样定理,而在不满意上述条件时则运用带通抽样定理。模拟电话信号经限带后的频率范围为300Hz~3400Hz,在抽样时按低通抽样定理,抽样频率至少为6800Hz。由于在实际实现时滤波器均有肯定宽度的过渡带,抽样前的限带滤波器不能对3400Hz以上频率重量完全予以抑制,在复原信号时也不行能运用志向的低通滤波器,所以对语音信号的抽样频率取为8kHz。这样,在抽样信号的频谱之间便可形成肯定间隔的爱护带,既防止频谱的混叠,又放松了对低通滤波器的要求。这种以适当高于奈奎斯特频率进行抽样的方法在实际应用中是很常见的。软件无线电所覆盖的频率范围一般都比较宽,只有宽频段才能具有广泛的适应性,这样宽的频段号的频谱 Xs(ω) 可表示为:+∞Xs(ω)=1Tsn=-∑∞X(ω- nωs)用Nyquist低通采样[1]是不现实的,而事实上软件无线电中的各种无线电信号的瞬时信号带宽都比较窄,这样采纳带通采样将使采样速率大大降低。本文重点是通过对带通采样定理进行严格的数学推导和论证,利用一系列不等式组提出了带通采样理论上能达到的最小采样速率。用matlab对定理进行仿真fs=100000;%抽样频率N=2^16;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);%频率为250HZfori=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));gridon;title('滤波前信号y');xlabel('时间/s');%原始信号subplot(2,1,2);plot(W(1:N/64),HW(1:N/64));%查看信号频谱gridon;title('滤波前信号频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');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));gridon;title('抽样信号chyang');xlabel('时间/s');subplot(2,1,2);plot(W(1:N/64),HWCH(1:N/64));agridon;title('抽样信号频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');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));gridon;title('采样之后');xlabel('时间/s');subplot(2,1,2);plot(W(1:N/64),HW2(1:N/64));%设置带通滤波器ws=[1130/fs*2,1370/fs*2];%ws1ws2阻带截止频率wp=[1180/fs*2,1320/fs*2];%wp1wp2通带截止频率[n,wn]=buttord(wp,ws,3,15);%求滤波器阶数n与截止频率wn[b,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,,-,]);gridon;title('滤波后信号');xlabel('时间/s');subplot(2,1,2);plot(W(1:N/64),HW3(1:N/64));gridon;title('滤波后信号频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');信号还原完成