资源描述：
《连续信号的采样与恢复实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验六、连续信号的采样与恢复一、实验目的1.加深理解采样对信号的时域和频域特性的影响；2.加深对采样定理的理解和掌握，以及对信号恢复的必要性；3.掌握对连续信号在时域的采样与重构的方法。二、实验原理(1)信号的采样信号的采样原理图如下图所示，其数学模型表示为：=其中的f(t)为原始信号，为理想的开关信号（冲激采样信号）δTs(t)=，fs(t)为采样后得到的信号称为采样信号。由此可见，采样信号在时域的表示为无穷多冲激函数的线性组合，其权值为原始信号在对应采样时刻的定义值。令原始信号f(t)的傅立叶变换为F(j
2、w)=FT(f(t))，则采样信号fs(t)的傅立叶变换Fs(jw)=FT(fs(t))=。由此可见，采样信号fs(t)的频谱就是将原始信号f(t)的频谱在频率轴上以采样角频率ws为周期进行周期延拓后的结果（幅度为原频谱的1/Ts）。如果原始信号为有限带宽的信号，即当
3、w
4、>
5、wm
6、时，有F(jw)=0，则有：如果取样频率ws≥2wm时，频谱不发生混叠；否则会出现频谱混叠。(2)信号的重构设信号f(t)被采样后形成的采样信号为fs(t)，信号的重构是指由fs(t)经过内插处理后，恢复出原来的信号f(t)的过程
7、。因此又称为信号恢复。由前面的介绍可知，在采样频率ws≥2wm的条件下，采样信号的频谱Fs(jw)是以ws为周期的谱线。选择一个理想低通滤波器，使其频率特性H(jw)满足：H(jw)=式中的wc称为滤波器的截止频率，满足wm≤wc≤ws/2。将采样信号通过该理想低通滤波器，输出信号的频谱将与原信号的频谱相同。因此，经过理想滤波器还原得到的信号即为原信号本身。信号重构的原理图见下图。通过以上分析，得到如下的时域采样定理：一个带宽为wm的带限信号f(t)，可唯一地由它的均匀取样信号fs(nTs)确定，其中，取样间
8、隔Ts<π/wm,该取样间隔又称为奈奎斯特（Nyquist）间隔。根据时域卷积定理，求出信号重构的数学表达式为：式中的抽样函数Sa(wct)起着内插函数的作用，信号的恢复可以视为将抽样函数进行不同时刻移位后加权求和的结果，其加权的权值为采样信号在相应时刻的定义值。利用MATLAB中的抽样函数来表示Sa(t)，有，，于是，信号重构的内插公式也可表示为：(3)模拟低通滤波器的设计在任何滤波器的设计中，第一步是确定滤波器阶数N及适当的截止频率Ωc。对于巴特沃斯滤波器，可使用MATLAB命令buttord来确定这些参
9、数，设计滤波器的函数为butter。其调用形式为[N,wn]=buttord(wp,ws,rp,rs,'s')[b,a]=butter[N,wn,'s']其中，wp、ws、rp、rs为待设计滤波器的技术指标，分别代表通带截止频率、阻带截止频率、通带最大衰减和阻带最小衰减；'s'表示设计滤波器的类型为模拟滤波器；N、wn为设计得到的滤波器的阶数和3dB截止频率，b、a为滤波器系统函数的分子和分母多项式的系数矢量，假定系统函数的有理分式表示为：三、程序示例示例1：选取门信号f(t)=g2(t)为被采样信号。利用M
10、ATLAB实现对信号f(t)的采样，显示原信号与采样信号的时域和频域波形。因为门信号并非严格意义上的有限带宽信号，但是，由于其频率f>1/τ的分量所具有的能量占有很少的比重，所以一般定义fm=1/τ为门信号的截止频率。其中的τ为门信号在时域的宽度。在本例中选取fm=0.5，临界采样频率为fs=1，过采样频率为fs>1（为了保证精度，可以将其值提高到该值的50倍），欠采样频率为fs<1。%显示原信号及其Fourier变换示例R=0.01;%采样周期t=-4:R:4;f=rectpuls(t,2);w1=2*pi
11、*10;%显示从-20*pi到20*pi频率范围内的频谱N=1000;%计算出2*1000+1个频率点的值k=0:N;wk=k*w1/N;F=f*exp(-j*t'*wk)*R;%利用数值计算求连续信号的Fourier变换，详细原理见附录Fudu=abs(F);%计算频谱的幅度wk=[-fliplr(wk),wk(2:1001)];Fudu=[fliplr(Fudu),Fudu(2:1001)];%计算对应负频率的频谱figure;subplot(2,1,1);plot(t,f);xlabel('t');yl
12、abel('f(t)');title('f(t)=u(t+1)-u(t-1)');subplot(2,1,2);plot(wk,Fudu);xlabel('w');ylabel('F(jw)');title('f(t)的Fourier变换');程序运行后的结果见下图。%显示采样信号及其Fourier变换示例R=0.25;%可视为过采样t=-4:R:4;f=rectpuls(t,2);w1=2*pi*
