实验报告连续时间信号的数字处理.docVIP

实 验 报 告 专业班级 姓名 学号 实验日期 实验名称 连续时间信号的数字处理 一、实验目的： 1、连续信号抽样为离散时间信号在时域和频域的效果。 2、模拟滤波器的设计。 二、实验内容： 1、时域中的抽样过程 连续正弦信号xa(t)的抽样。 % Program P5_1 % Illustration of the Sampling Process % in the Time-Domain clf; t = 0:0.0005:1; f = 13; xa = cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);grid xlabel(Time, msec);ylabel(Amplitude); title(Continuous-time signal x_{a}(t)); axis([0 1 -1.2 1.2]) subplot(2,1,2); T = 0.1; n = 0:T:1; xs = cos(2*pi*f*n); k = 0:length(n)-1; stem(k,xs);grid; xlabel(Time index n);ylabel(Amplitude); title(Discrete-time signal x[n]); axis([0 (length(n)-1) -1.2 1.2])； 运行结果： 将正弦信号频率改为3HZ。 运行结果： 将正弦信号频率改为7HZ。 运行结果： 由图可知，在正弦信号3种不同频率下相应的等效离散时间信号间没有差别。因为xs=cos(2πfn), xs(f=13)=cos(0.6πn’); xs(f=3)=cos(0.6πn’); xs(f=7)=cos(0.6πn’);xs的表达式相同，故等效离散时间信号相同。 2、时域中的混叠效果 （1）、从P5.1所产生的x[n]中产生重构的信号ya(t)。 % Program P5_2 % Illustration of Aliasing Effect in the Time-Domain % Program adapted from [Kra94] with permission from % The Mathworks, Inc., Natick, MA. clf; T = 0.1;f = 13; n = (0:T:1); xs = cos(2*pi*f*n); t = linspace(-0.5,1.5,500); ya = sinc((1/T)*t(:,ones(size(n))) - (1/T)*n(:,ones(size(t))))*xs; plot(n,xs,o,t,ya);grid; xlabel(Time, msec);ylabel(Amplitude); title(Reconstructed continuous-time signal y_{a}(t)); axis([0 1 -1.2 1.2]); 运行结果： （2）、在程序P5.2中，t的范围-0.5至1.5，时间增量的值为2/（500-1）=0.004008。在图中，t的范围是0至1。 将图中t的范围改成-0.5至1.5， % Program P5_2 % Illustration of Aliasing Effect in the Time-Domain % Program adapted from [Kra94] with permission from % The Mathworks, Inc., Natick, MA. clf; T = 0.1;f = 13; n = (0:T:1); xs = cos(2*pi*f*n); t = linspace(-0.5,1.5,500); ya = sinc((1/T)*t(:,ones(size(n))) - (1/T)*n(:,ones(size(t))))*xs; plot(n,xs,o,t,ya);grid; xlabel(Time, msec);ylabel(Amplitude); title(Reconstructed continuous-time signal y_{a}(t)); axis([-0.5 1.5 -1.2 1.2]); 运行结果： 这种改变后产生的曲线在时间轴0至1上逼近重构的连续时间信号ya(t),在这范围外差别大。？？？？？ （3）、恢复原始信号显示范围并通过分别改变正弦信号的频率到3