实验一低通采样定理和内插与抽取实现a.docxVIP

实验一：低通采样定理和内插与抽取实现 实验目的 连续信号和系统的表示方法，以及坊真方法。 用MATLAB实现连续信号采用与重构的方法， 采样信号的插值和抽取等重采样实现方法。 用时域采样信号重构连续时域信号的原理和方法。 用MATLAB绘图函数表示信号的基本方法，实验数据的可视化表示。 原理 1、时域抽样定理 令连续信号xa(t)的傅里叶变换为Xa(ji」),抽样脉冲序列p(t)傅里叶变换为 Pj」),抽样后的信号xA(t)的傅里叶变换为XV1)若采用均匀抽样,抽样周期Ts, 抽样频率为门s=2二fs,由前面分析可知：抽样的过程可以通过抽样脉冲序列 p(t) 与连续信号xa(t)相乘来完成，即满足:xA(t)=xa(t) p(t)，又周期信号f(t)傅里叶变换 为： m Ilf ( i) / = 二心可\( Wf jn= — -oa 故可以推得p(t)的傅里叶变换为： P( 5^) = 2n c- d H J=- — 其中： p n =月 z p D l at 根据卷积定理可知： X(g = p(jo) 得到抽样信号x(t)的傅里叶变换为： W X(j Q) = 。一 ?Q j / n= — e 其表明:信号在时域被抽样后，他的频谱X(j门)是连续信号频谱X(j门)的形状 以抽样频率门为间隔周期重复而得到，在重复过程中幅度被p(t)的傅里叶级数Pn 加权。因为Pn只是n的函数，所以Xj」)在重复的过程中不会使其形状发生变化。 假定信号x(t)的频谱限制在 亠m~+「m的范围内,若以间隔Ts对xa(t)进行 抽样,可知抽样信号XA(t)的频谱XW)是以门s为周期重复。显然，若在抽样的 过程中门s2i m,则XT门)将发生频谱混叠现象，只有在抽样的过程中满足 11 s=2「m条件,XW门)才不会产生频谱的混叠，接收端完全可以由xA(t)恢复原 连续信号xa(t),这就是低通信号抽样定理的核心内容。 2、信号的重建 从频域看，设信号最高频率不超过折叠频率 Xa(j 门)=Xa(j 门) 『s/2 Xa『」)=0 p | 门 s/2 则理想取样后的频谱就不会产生混叠，故有： Ar ( j = J, 2Z A： . j o- j?t o. j tij I I ￡^/2 让取样信号xA(t)通过一带宽等于折叠频率的理想低通滤波器 H(j )=T |门 s/2 H(j「)=0 |门 |「s/2 滤波器只允许通过基带频谱，即原信号频谱，故： Y(j 门)=XA(j 门)H(j 门)=Xa(j 门) 因此在滤波器的输出得到了恢复的原模拟信号： y(t)=xa(t) 从时域上看，上述理想的低通滤波器的脉冲响应为 , JT mu 于* □n: 根据卷积公式可求得理想低通滤波器的输出为 y(t) — nT) ha(t — nT ) 由上式显然可得 ha(t- nT) = sin(nT)/(^/T)(t- nT) JT JT 上式表明只要满足取样频率高于两倍信号最高频率，连续时间函数xa(t)就可 用他的取样值xa(nT)来表达而不损失任何信息，这时只要把每一个取样瞬时值与 内插函数式相乘求和即可得出xa(t),在每一取样点上，由于只有该取样值所对应的 内插函数式不为零，所以各个取样点上的信号值不变。 内容 内容 1?设计连续时间信号线性滤波器分离信号组份 99 已知信号x t [=為(m 1)cos 2二(100m 50)t ，试设计滤波器，分离出如下信号: m ^0 m =1,2,3……50 m =51,52,53……100 m =40,41,42……60 m =1,2,3……40,61,62……100 据以下采样频率：(a) fs = 20000Hz (b) f^ 10000Hz (c) f^30000Hz 求信号频谱及相应的滤波器。 参考程序如下: 设计一个Butterworth 模拟带通滤波器，设计指标为：通带频率： 1000- 2000Hz,两侧过渡带宽500Hz,通带波纹1dB,阻带衰减100dB假设一个信号, 其中f1=100Hz,f2=1500Hz,f3=2900Hz 。信号的采样频率为10000Hz=试将原信号 与通过该滤波器的模拟信号进行比较。参考程序如下： wp=[1000 2000]*2*pi;ws=[500 2500]*2*pi;Rp=1;Rs=100; % 滤波器设计参数，对于给 定Hz应乘以2 [N,Wn]=buttord(wp,ws,Rp,Rs,s); % 求得滤波器的最小阶数和截止频率 w=li nspace(1,3000,1000)*2*pi; %设置绘制频率响应的频率点 [b,a]=butter(N,Wn,s); %设计模拟 Butterworth 滤波器 H=freqs(b,a,w); %计算给定频率点的