信号与系统第4题.doc

MATLAB实验4 一、实验题目：连续时间信号的时域取样 二、实验目的：掌握连续时间信号的离散化过程，深刻理解时域取样定理。 三、参考书目： [1] 吴湘淇等《信号、系统与信号处理》（修订本），电子工业出版社，2001； [2] 吴湘淇等《信号、系统与信号处理软硬件实现》，电子工业出版社，2002； [3] John R. Buck，《信号与系统Explorations in Signals and Systems》1)信号解析式已知的信号采样； 2)信号解析式未知的信号采样； 3)分析信号采样过程中可能的误差； 4)详细列出连续信号采样的步骤和主要原理。 实验原理及步骤： 实验原理： 时域中的抽样依据是时域抽样定理：若带限信号的最高角频率为wm，则信号f(t)可以用等间隔的抽样值唯一地表示。而抽样间隔必须不大于1/2fm，或者说最低抽样频率为2fm。此时取样信号不会混叠，而且可以复原原信号。 若信号持续时间无限，且无解析表达式，先在图形中找出变化速度最的相邻两个点，并假设两点之间峰与谷的时间间隔td近似等于最高频率的半个周期，即td≈Tm/2, 则有fm=1/Tm=1/2td,选取fs=2.5fm作为取样间隔试验值，即 Tt=1/fs=2/2.5td=0.8td 用反复试探法适当选取参数进行分析。 2、实验步骤： 一般连续非周期信号的处理步骤如下： 1)根据取样定理，确定取样间隔T。 2)确定时域截取长度L(或者窗函数的点数M)。 3）根据时域抽样定理进行取样 a)对信号sa(t)进行取样 sa(t)为连续非周期信号，所以需要对其进行截断，我选取了-15到15的时间范围，另外就是选取抽样间隔，经过实验试探，我发现取Ts=时即能较好地进行取样。 b)持续时间无限且无解析表达式的信号取样 在matlab中选取了一段语音信号，这是没有解析式的，按照上述原理，经过观察，td近似为0.35s，所以Ts=0.7s，此时即能取得较好的取样信号。 实验中出现的现象，分析其原因，得出结论 可以看出，时域信号的取样即是一个离散化的过程。即可在按照取样定理选取的抽样间隔进行取样，即能取得较好的取样波形，采样后取得的包络与原信号大大致相似。不过，由于信号的无限长，需要有截断处理，而且非周期信号不可能完全按照取样定理选取理论的抽样间隔，就会出现混叠误差与截断误差。一般情况下，混叠误差的危害大于截断误差，我们在采样后加一抗混滤波器即可。 回答思考题： 1)在实际工程应用中，如何确定采样速率？ 在实际工程应用中，往往按照fm的3到4倍的速度取样即可。 2)理想采样和实际采样有何区别？ 理想采样是指用理想冲击串δ(t)进行采样 ，而在实际处理中这是不可能的，实际中往往用脉冲进行采样，取得的包络与原信号在误差范围内即可。 3)在连续信号离散化过程中，会出现哪些误差？如何克服或减弱？ 不过，由于信号的无限长，需要有截断处理，而且非周期信号不可能完全按照取样定理选取理论的抽样间隔，就会出现混叠误差与截断误差。一般情况下，混叠误差的危害大于截断误差，我们在采样后加一抗混滤波器即可。 八、附程序及运行结果： a) 编写Matlab程序如下： wm=1; %w信号带宽 wc=wm; %滤波器截止频率 Ts=pi/wm; %采样间隔 ws=2*pi/Ts; %采样角频率 n=-100:100; %时域采样点数 nTs=n.*Ts; %时域采样点 f=sinc(t); t=-15:0.5:15; subplot(2,1,1); plot(t,f); xlabel(t); ylabel(f(t)); title(采样前); %原信号的时域图 t1=-15:0.5:15; f1=sinc(t1/pi); subplot(2,1,2); stem(t1,f1); xlabel(kTs); ylabel(f(kTs)); title(采样后); %采样信号的时域图 所得图如下： b) 编写Matlab程序如下： load mtlb; y=mtlb; Fs=7418; t=(0:length(y)-1)/Fs; subplot(2,1,1); plot(t,y); xlabel(t); ylabel(y(t)); title(采样前); %原信号的时域图 Ts=0.7; n=t./Ts; subplot(2,1,2); stem(n,y); xlabel(n); ylabel(