采样的时频域分析

电子科技大学实验报告一、实验室名称：数字信号处理实验室二、实验项目名称：采样的时域及频域分析三、实验原理：1、采样的概念：采样是将连续信号变化为离散信号的过程。1.A、理想采样：即将被采样信号与周期脉冲信号相乘B、实际采样：将被采样信号与周期门信号相乘，当周期门信号的宽度很小，可近似为周期脉冲串。根据傅里叶变换性质0000()()()()ˆˆ()()()()()()(())FTFTaaTnnFTaaTaTaannxtXjTjxtxtTxnTtnTXjXjn式中T代表采样间隔，01T由上式可知：采样后信号的频谱是原信号频谱以0为周期的搬移叠加结论：时域离散化，频域周期化；频谱周期化可能造成频谱混迭。)(txa)(tT^)(txa)(txa()Tpt^)(txa()axt)(jXaC、低通采样和Nyquist采样定理设()()aaxtXj且()0,2aMMXjf当，即为带限信号。则当采样频率满足2/22sMMff时,可以从采样后的^()()()aassnxtxnTtnT信号无失真地恢复()axt。称2Mf为奈奎斯特频率，12NMTf为奈奎斯特间隔。注意：实际应用中，被采信号的频谱是未知的，可以在ADC前加一个滤波器（防混迭滤波器）。2、低通采样中的临界采样、欠采样、过采样的时域及频域变化情况。低通采样中的临界采样是指在低通采样时采样频率2sMff低通采样中的欠采样是指在低通采样时采样频率2sMff低通采样中的欠采样是指在低通采样时采样频率2sMff设一带限信号的频谱如下：ˆ()axt)(ˆjXa()aGj0mm()Pj0TT2Tm（1）临界采样（2）过采样（3）欠采样由上图可知，当为临界采样和过采样时，理论上可以无失真的恢复采样信号，但是实际在临界采样时，由于实际滤波器的性能限制，无法无失真的恢复，在欠采样时只能部分恢复原信号的频谱特性。因此过采样时使用最为广泛的采样方式，当需要注意的是对临界采样和欠采样由于采样频率可以降低，在不需要恢复出信号的全部频谱特征时，则往往使用这两种采样方式。随着信号处理技术的发展，信号的频率越来越高，这两种方式也有着广泛的应用前景。()Pj0TT2TmTT()Pj0TT2TmTT()PGj()PGj()PGjTT1/T1/T1/T在理论分析中使用的带限信号在实际应用是并不存在的，因为要求该信号在时域上是无限长的，因此无论采样频率有多大，实际采样的信号都是会发生混叠的，如下图所示：在实际应用中，我们只需使采样频率满足能够恢复出我们需要的信号即可。3、带通采样过程及带通采样定理。带通采样是对于带通信号进行采样的过程。LHHL0||,称为带通信号的带宽。此时采样频率为2()1221HLsfm其中m是当采样频率满足122sf时最大的正整数。此时信号可以被无失真的恢复，这就是带通采样定理。原理：采样后的带通信号同样是原信号的周期搬移叠加，但由于带通信号在某个频带不存在信号分量，采样后得到信号频谱存在间隔，当采样频率满足一定条件（不满足底通采样定理）时，同样可以无失真的恢复。示意图如下：（1）当最高频率H是带宽的整数倍，即()HM，而选择的抽样频率22()HTM，此时有()Pj0TT2TmTT()PGj1/T()aGj0mm从图中可以看出，当把该采样信号通过一个理想带通滤波器时，可以恢复出原信号。（2）当最高频率H不是带宽的整数倍，我们可以认为的扩展带宽，使得该带通信号的()HM，而选择的抽样频率22()HTM，此时有()aGj0HLHL()Pj0TT2Tm()pGj0HLHLHLL00从上图可以看出同样能无失真的恢复出原带通信号(拓展知识):4、变采样率的数字信号处理A、降采样率（整数倍抽取）的实现原理，时域和频域的变化情况。降采样率是指每次抽样保留输入序列中的第M个样本，而除去中间的M-1个样本：[][]ynxnM用框图表示为可以得到11/01()()MMkMkYzXzWM，以2倍下抽样器为例，即L=2，可得/2/2/2(2)/211(){()()}{()()}22jjjjjYeXeXeXeXe，如下图所示0H()Pj0TT2Tm()pGj0HLHL00()jXe220/2()jXe220/2()jXeMx[n]y[n]可以知道，在降采样率时，()jXe的原形状会丢失，即发生混叠现象。M倍下抽样器的输出和输入之间傅氏变换的关系为：1(2)/01()()MjjkMkYeXeM在下抽样以前，为了避免引起混叠，信号需要通过一个低通滤波器来带限到||/M即：B、升采样率（整数倍内插）的实现原理，时域和频域的变化情况。升采样率是指通过在对原离散信号的两个连续样本间插入L-1个等距的样本值（不一定为零），亦即抽样因子为L的上抽样。上抽样后的序列长度为原来的L倍：[/],0,,2[]0,uxnLnLLxnotherwise，框图表示为可以得到：()()LuXzXz，()()jjLuXeXe，对于L=2时，可得下图：x[n]xu[n]L混叠()jXe220()jXe220Mx[n]y[n]()Hz如图，2倍的抽样率扩展导致频谱的2倍重复，表明傅里叶变换以2倍压缩。因此可得输入频谱的一个额外镜像，这个过程也叫做映射。上采样后不必要的镜像必须用一个称为内插滤波器的低通滤波器H(z)来消除，即：C、分数倍变采样率的实现原理，时域和频域的变化情况。采样率的分数转换可以用M倍抽取器和L倍内插器级联而成，其中M和L都是正整数。这样级联有两种可能的形式四、实验目的：深刻理解低通采样中的临界采样的时域及频域变化情况。深刻理解低通采样中的欠采样的时域及频域变化情况。深刻理解低通采样中的过采样的时域及频域变化情况深刻理解带通采样过程及带通采样定理。五、实验内容：本实验要求学生运用MATLAB编程完成可变采样率采样（抽取）程序，并对提供的离散时间信号分别进行临界采样、过采样、欠采样时信号时域和频域的信号变化情况，以加深对相关教学内容的深刻理解。进而拓展到可变采样率信号处理的基本方法的MATLAB实现，得到信号的时频域变化情况，使学有余力的同学进一步加深对变采样率信号处理相关知识的理解。六、实验器材（设备、元器件）：Pc机,DSP试验箱七、实验步骤：1、在MATLAB中设计完成可变采样率采样（抽取）程序。2、对比观察、分析各种采样（临界采样、过采样、欠采样）时域频域的情况。八、实验数据及结果分析：x[n]xu[n]()HzLy[n]M()dHz()HzL()uHzM()dHzL()uHzML1.可变采样率采样（抽取）程序:clearclcf=input('thesampingfrequencyf=');N=input('采样长度N=');n=0:0.1:N-1;t=0:0.01:10;x=sin(20*pi*t);subplot(2,2,1)plot(t,x,'r');axis([00.2-1.51.5]);gridxlabel('t'),ylabel('margin')title('origialsingal')y=sin(20*pi*n/f);holdonstem(n/f,y);axis([00.5-1.51.5]);gridtitle('thesampedsingal')holdofftitle('timesingalofbeforesampeandsamped')h1=fft(x,1024);h=fft(y,1024);subplot(2,2,2)plot(N*[-511:512]/1024,abs(h1));gridaxis([0501max(abs(h1))])title('Frequencydomainwaveformoriginalsingal')subplot(2,2,3)plot(f*[-511:512]/1024,abs(h));gridtitle('Frequencydomainwaveformsmapedsingal')由以上程序可以知道函数信号x的频率为10Hz，以下是对该信号的三种采样方式：（1）临界采样，有奈切斯特定理有临界频率为20Hz，输入f=20，采样长度20，得到以下波形由上图可以看出每隔0.1采样一个点，原函数的频谱图在右上，主要是集中在3、27的附近，当采样后是对频谱的搬移，搬移20个单位得到左下的图，从-3搬到了17，而3的搬到了-17处，27的搬到了47，为了保证分析的有效性该图未显示后两种。（2）欠采样，输入f=5.N=20；得到以下图形由上图知道当搬移5个单位后，可以看到一些频谱混叠再来一起。（3）过采样，输入f=30,N=20，得到以下波形从以上图形知道在-2的点搬到了28处，搬移后无混叠现象，能够很好地恢复出波形。2.降采样的程序：九、实验结论：经过matlab仿真的程序波形与理论的比较的吻合，临界采样可以刚好的恢复出原始的波形，欠采样产生了混叠不能有效地恢复波形，过采样也可以很好地恢复原始波形。十、总结及心得体会：通过本实验，对于采样有了一个比较全面的认识。报告评分：指导教师签字：