ch12时域及采样.ppt

图1.5.7 理想恢复 + + + 理想恢复 Matlab Implementation t=t1:Dt:t2; n=n1:n2; Fs=1/Ts; nTs = n*Ts t=t1:Dt:t2; n=n1:n2; Fs=1/Ts; nTs = n*Ts Xa = x*sinc(Fs*(ones(length(n),1)*t- nTs’*ones(1,length(t)))) Sinc function : sinc(x) Description : sinc computes the sinc function of an input vector or array, where the sinc function is Xa = x*sinc(Fs*(ones(length(n),1)*t- nTs’*ones(1,length(t)))) 图1.5.8 实际D/AC方框图 实际数模转换器 Practical D/A converters Zero-order-hold (ZOH) interpolation: In this interpolation a given sample value is held for the sample interval until the next sample is received. 由时域离散信号xa(nT)恢复模拟信号的过程是在采样点内插的过程。理想低通滤波的方法是用g(t)函数作内插函数，还可以用一阶线性函数作内插。零阶保持器是将前一个采样值进行保持，一直到下一个采样值来到，再跳到新的采样值并保持，因此相当于进行常数内插。零阶保持器的单位冲激函数h(t)以及输出波形如图1.5.9所示。对h(t)进行傅里叶变换，得到其传输函数： 图1.5.9 零阶保持器的输出波形 (1.5.10) 图1.5.10 零阶保持器的频率特 性? 由该图看到零阶保持器是一个低通滤波器，能够起到将时域离散信号恢复成模拟信号的作用。图中虚线表示理想低通滤波器的幅度特性。零阶保持器的幅度特性与其有明显的差别，主要是在|Ω|π/T区域有较多的高频分量，表现在时域上，就是恢复出的模拟信号是台阶形的。因此需要在D/AC之后加平滑低通滤波器，滤除多余的高频分量，对时间波形起平滑作用，这也就是在图1.5.1模拟信号数字处理框中，最后加平滑滤波的原因。虽然这种零阶保持器恢复的模拟信号有些失真，但简单、易实现，是经常使用的方法。 First-order-hold (FOH) interpolation In this case the adjacent samples are joined by straight lines. Matlab Implementation: ZOH: stairs FOH: plot Examples----Reconstruction of signal Let 1、Determine and plot its Fourier transform Xa (j?); 2、Sample xa(t) at Fs=5000 sam/sec to obtain x1(n). Determine and plot X1(e jω). 3、Sample xa(t) at Fs=1000 sam/sec to obtain x2(n). Determine and plot X2(e jω). % 用sinc函数重构信号 % 模拟信号 Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t)); subplot(1,1,1) subplot(2,1,1);plot(t*1000,xa); a=axis;text(a(2)+0.2,a(3),t 毫秒); ylabel(xa(t)) title(模拟信号) % 离散时间信号 x1(n) Ts = 0.0002; Fs = 1/Ts; n = -25:1:25; nTs = n*Ts; x = exp(-1000*abs(nTs)); % 模拟信号重构 t = -0.005:Dt:0.005; xa = x * sinc(Fs*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); % 校验 error = max(abs(xa - exp(-1000*abs(t)))) subplot(2,1,2);plot(t*1000,xa); xlabel(t 毫秒); ylabel(xa(