信号与系统实验三 四信号与系统实验三 四.doc

实验三、 连续信号的采样与恢复 一、实验目的 加深理解采样对信号的时域和频域特性的影响； 加深对采样定理的理解和掌握，以及对信号恢复的必要性； 掌握对连续信号在时域的采样与重构的方法。 二、实验原理 (1) 信号的采样 信号的采样原理图如下图所示，其数学模型表示为： = 其中的f(t)为原始信号，δTs(t)为理想的开关信号（冲激采样信号）δTs(t) =，fs(t)为采样后得到的信号称为采样信号。由此可见，采样信号在时域的表示为无穷多冲激函数的线性组合，其权值为原始信号在对应采样时刻的定义值。 令原始信号f(t)的复里叶变换为F(ejw)=FT(f(t))，则采样信号fs(t) 的复里叶变换Fs(ejw)=FT(fs(t))=。由此可见，采样信号fs(t)的频谱就是将原始信号f(t)的频谱在频率轴上以ws为周期进行周期延拓后的结果（幅度为原频谱的1/Ts）。如果原始信号为有限带宽的信号（即当|w||wm|时，有F(ejw)=0），则有：如果取样频率ws≥2wm时，频谱不发生混叠；否则会出现频谱混叠。 (2) 信号的重构 设信号f(t)被采样后形成的采样信号为fs(t)，信号的重构是指由fs(t)经过内插处理后，恢复出原来的信号f(t)的过程。因此又称为信号恢复。 由前面的介绍可知，在采样频率ws≥2wm的条件下，采样信号的频谱Fs(ejw)是以ws为周期的谱线。选择一个理想低通滤波器，使其频率特性H(jw)满足： H(jw)= 式中的wc称为滤波器的截止频率，满足wm≤wc≤ws/2。将采样信号通过选择的理想低通滤波器，输出信号的频谱将与原信号的频谱相同。根据信号的时域表示与频率表示的一一对应关系可得，经过理想滤波器还原得到的信号即为原信号本身。信号重构的原理图见下图。 通过以上分析，得到如下的时域采样定理：一个带宽为wm的带限信号f(t)，可唯一地由它的均匀取样信号fs(nTs)确定，其中，取样间隔Tsπ/wm, 该取样间隔又称为奈奎斯特（Nyquist）间隔。 根据时域卷积定理，求出信号重构的数学表达式为： = = = 式中的抽样函数Sa(wct)起着内插函数的作用，信号的恢复可以视为将抽样函数进行不同时刻移位后加权求和的结果，其加权的权值为采样信号在相应时刻的定义值。利用MATLAB中的抽样函数Sinc(t)=sin(πt)/(πt)来表示Sa(t)，有Sa(t)=Sinc(t/π)，于是，信号重构的内插公式也可表示为： = (3) 模拟低通滤波器的设计 在任何滤波器的设计中，第一步是确定滤波器阶数N及适当的截止频率Ωc。对于巴特沃斯滤波器，可使用MATLAB命令buttord来确定这些参数，设计滤波器的函数为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)为被采样信号。利用MATLAB实现对信号f(t)的采样，显示原信号与采样信号的时域和频域波形。 因为门信号并非严格意义上的有限带宽信号，但是，由于其频率f1/τ的分量所具有的能量占有很少的比重，所以一般定义fm=1/τ为门信号的截止频率。其中的τ为门信号在时域的宽度。在本例中选取fm=0.5，临界采样频率为fs=1，过采样频率为fs1（为了保证精度，可以将其值提高到该值的50倍），欠采样频率为fs1。 % 显示信号及其Fourier变换示例 R=0.01;t=-4:R:4; f=rectpuls(t,2) w1=2*pi*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变换 F=abs(F);wk=[-fliplr(wk),wk(2:1001)]; F=[fliplr(F),F(2:1001)]; % 计算对应负频率的频谱 figure; subplot(2,1,1); plot(t,f);xlabel(t); ylabel(f(t));title(f(t)=u(t+1)-u(t-1)); subplot(2,1,2); plot(wk,F)