《采样定理(软件)分析》由会员分享，可在线阅读，更多相关《采样定理(软件)分析（10页珍藏版）》请在金锄头文库上搜索。
1、、实验项目名称：采样的时域及频域分析二、实验原理:1、采样的概念：采样是将连续信号变化为离散信号的过程1. A、理想采样：即将被采样信号与周期脉冲信号相乘B、实际采样：将被采样信号与周期门信号相乘, 可近似为周期脉冲用。.当周期门信号的宽度很小,PT（t）根据傅里叶变换性质Xa(t)FT Xa(j )nxa(t)Xa(t) T(T)nT(T) FT 0 0(j )Xa(nT) T(t nT) FT 兄(j )n 0)式中T代表采样间隔,由上式可知：采样后信号的频谱是原信号频谱以为周期的搬移叠加结论：时域离散化，频域周期化；频谱周期化可能造成频谱混迭xa(t)设 Xa(t)Xa(j )且*,。)
2、 0,当门 M 2 fM ,即为带限信号。则当采样频率满足fs 2 M/22 fM时，可以从采样后的Xa(t)Xa(nTs) (t nTs)信号无失真地恢复Xa(t)。称2fM为奈奎斯特频率,n1Tn -fM为奈奎斯特间隔。2一、/汪忠：实际应用中，被采信号的频谱是未知的，可以在ADC前加一个滤波器(防混迭滤波器)。2、低通采样中的临界采样、欠采样、过采样的时域及频域变化情况低通采样中的临界采样是指在低通采样时采样频率fs 2 fM低通采样中的欠采样是指在低通采样时采样频率fs 2 fM低通采样中的欠采样是指在低通采样时采样频率fs 2 fMA P(j )(1)临界采样fs2 fMP(j )(
3、2)过采样fs 2 fMP(j )(3)欠采样fs 2 fM由上图可知，当为临界采样和过采样时，理论上可以无失真的恢复采样信号， 但是实际在临界采样时，由于实际滤波器的性能限制，无法无失真的恢复，在欠 采样时只能部分恢复原信号的频谱特性。因此过采样时使用最为广泛的采样方式， 当需要注意的是对临界采样和欠采样由于采样频率可以降低，在不需要恢复出信号的全部频谱特征时，则往往使用这两种采样方式。随着信号处理技术的发展， 信号的频率越来越高，这两种方式也有着广泛的应用前景。在理论分析中使用的带限信号在实际应用是并不存在的，因为要求该信号在 时域上是无限长的，因此无论采样频率有多大，实际采样的信号都是会
4、发生混叠 的，如下图所示：Ga(j )m 0mT 2 mP(j ) 鼻鼻小在实际应用中，我们只需使采样频率满足能够恢复出我们需要的信号即可3、带通采样过程及带通采样定理。带通采样是对于带通信号进行采样的过程。0 L | | H, HL称为带通信号的带宽。此时采样频率为22_HD其中m是当采样频率满足fs 2时最大的正整数。此2 2m 12时信号可以被无失真的恢复，这就是 带通采样定理。原理：采样后的带通信号同样是原信号的周期搬移叠加， 但由于带通信号在某个 频带不存在信号分量，采样后得到信号频谱存在间隔，当采样频率满足一定条件（不满足底通采样定理）时，同样可以无失真的恢复。示意图如下：H是带宽
5、的整数倍，即HM (),而选择的抽样频率(1)当最图频率卜 P(j )IJT2 mT 0 T从图中可以看出，当把该采样信号通过一个理想带通滤波器时，可以恢复出原信号。(2)当最高频率 h不是带宽的整数倍，我们可以认为的扩展带宽，使得该带一，、一一 2 .通信号的H M(),而选择的抽样频率 T 2( ) H ,此时有MJkiP(j )T2 mGa(j )Gp(j )T 0 T从上图可以看出同样能无失真的恢复出原带通信号4、变采样率的数字信号处理(拓展知识)A、降采样率(整数倍抽取)的实现原理，时域和频域的变化情况。降采样率是指每次抽样保留输入序列中的第M个样本，而除去中间的M-1个样本：yn
6、xnM用框图表示为 IMxnynM 1可以得到Y(z) X(zMWMk),以2倍下抽样器为例，即L=2,可得M k 0Y(ej )2X(ej /2) X(ej/2)$X(ej/2) X(ej( 2)/2),如下图所示202可以知道，在降采样率时，X(ej )的原形状会丢失，即发生混叠现象。 M倍下M 1抽样器的输出和输入之间傅氏变换的关系为：Y(ej ) X(ej(2k)/M)M k 0在下抽样以前，为了避免引起混叠，信号需要通过一个低通滤波器来带限到| | /M 即：B、开采样率(整数倍内插)的实现原理，时域和频域的变化情况。开采样率是指通过在对原离散信号的两个连续样本间插入L-1个等距的样
7、本值(不一定为零)，亦即抽样因子为L的上抽样。上抽样后的序列长度为原来的L,、xn/L,n 0, L, 2L倍： xun,框图表小为0, otherwisexnxun202如图，2倍的抽样率扩展导致频谱的2倍重复，表明傅里叶变换以2倍压缩。因 此可得输入频谱的一个额外镜像， 这个过程也叫做映射。上采样后不必要的镜像 必须用一个称为内插滤波器的低通滤波器 H（z）来消除，即：C、分数倍变采样率的实现原理，时域和频域的变化情况。采样率的分数转换可以用 M倍抽取器和L倍内插器级联而成，其中M和L都是 正整数。这样级联有两种可能的形式三、实验目的：深刻理解低通采样中的临界采样的时域及频域变化情况。深刻
8、理解低通采样中的欠采样的时域及频域变化情况。深刻理解低通采样中的过采样的时域及频域变化情况深刻理解带通采样过程及带通采样定理。（拓展内容）理解降采样率（整数倍抽取）的实现原理，时域和频域的变化情况。理解开采样率（整数倍内插）的实现原理，时域和频域的变化情况。理解分数倍变样率的实现原理，时域和频域的变化情况。通过具体实践，理解在高倍数变采样率的情况中，应当采用多级实现方案。学习设计用于抽取和内插的滤波器。四、实验内容：本实验要求学生运用MATLAB!程完成可变采样率采样（抽取）程序，并对 提供的离散时间信号分别进行临界采样、 过采样、欠采样时信号时域和频域的信 号变化情况，以加深对相关教学内容的
9、深刻理解。进而拓展到可变采样率信号处 理的基本方法的MATLA殴现，得到信号的时频域变化情况，使学有余力的同学 进一步加深对变采样率信号处理相关知识的理解。五、实验步骤：1、在MATLAB 中设计完成可变采样率采样(抽取)程序。2、对比观察、分析各种采样(临界采样、过采样、欠采样)时域频域的情况。3、 (拓展要求)设计完成整数倍内插的MATLAB 程序，观察时域频域的变化情况，提出相应滤波器设计要求。4、 (拓展要求)设计分数倍变采样率的MATLAB 程序，观察时域频域的变化情况，提出相应滤波器设计要求。5、(拓展要求)通过硬件(DSP实验箱演示上述信号的采样时域(示波器)波 形及频域波形(计
10、算结果)。并与MATLA勰序作比较对照。六、实验数据及结果分析：新建文件为caiyang.m,在该文件中输入程序如下：function fz=caiyang(fy,fs)%第一个输入变量是原信号函数，信号函数fy 以字符串的格式输入%第二个输入变量是采样频率fs0=10000; tp=0.1;t=-tp:1/fs0:tp;k1=0:999; k2=-999:-1;m1=length(k1); m2=length(k2);f=fs0*k2/m2,fs0*k1/m1;% 设置原信号的频率数组w=-2*pi*k2/m2,2*pi*k1/m1;fx1=eval(fy);FX1=fx1*exp(-j*1
11、:length(fx1)*w);%求原信号的离散时间傅里叶变换figure% 画原信号波形subplot(2,1,1),plot(t,fx1,r)title(原信号),xlabel(时间 t (s)axis(min(t),max(t),min(fx1),max(fx1) % 画原信号幅度频谱subplot(2,1,2),plot(f,abs(FX1),r)title(原信号幅度频谱),xlabel(频率f (Hz)axis(-100,100,0,max(abs(FX1)+5)% 对信号进行采样Ts=1/fs;%采样周期t1=-tp:Ts:tp;%采样时间序列f1=fs*k2/m2,fs*k1/
12、m1;%设置采样信号的频率数组t=t1;%变量替换fz=eval(fy);%获取采样序列FZ=fz*exp(-j*1:length(fz)*w);%采样信号的离散时间傅里叶变换figure% 画采样序列波形subplot(2,1,1),stem(t,fz,.),title(取样信号)，xlabel(时间 t (s)line(min(t),max(t),0,0)% 画采样信号幅度频谱subplot(2,1,2),plot(f1,abs(FZ),m)title(取样信号幅度频谱),xlabel(频率f (Hz)%主函数f1=cos(2*pi*50*t) ;%输入一个信号fs0=caiyang(f1,80); %频率fs2fmax ，即 欠采样fs1=caiyang(f1,100);%S率 fs 2 fmax,临 界采样fs2=caiyang(f1,150);%率 fs2 fmax,即 过采样七、实验结论：
