信号与系统实验信号抽样与内插
信号与系统实验四实验报告
实验四 时域抽样与频域抽样一、实验目的加深理解连续时间信号的离散化过程中的数学概念和物理概念，掌握时域抽样定理的基本内容。
掌握由抽样序列重建原连续信号的基本原理与实现方法，理解其工程概念。
加深理解频谱离散化过程中的数学概念和物理概念，掌握频域抽样定理的基本内容。
二、 实验原理时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件：对于基带信号，信号抽样频率sam f 大于等于2倍的信号最高频率m f ，即m sam f f 2≥。
时域抽样是把连续信号x (t )变成适于数字系统处理的离散信号x [k ] ；信号重建是将离散信号x [k ]转换为连续时间信号x (t )。
非周期离散信号的频谱是连续的周期谱。
计算机在分析离散信号的频谱时，必须将其连续频谱离散化。
频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件。
三．实验内容1. 为了观察连续信号时域抽样时抽样频率对抽样过程的影响，在[0，0.1]区间上以50Hz 的抽样频率对下列3个信号分别进行抽样，试画出抽样后序列的波形，并分析产生不同波形的原因，提出改进措施。
)102cos()(1t t x ⨯=π答: 函数代码为: t0 = 0:0.001:0.1;x0 =cos(2*pi*10*t0);plot(t0,x0,'r')hold onFs =50;t=0:1/Fs:0.1;x=cos(2*pi*10*t); stem(t,x); hold offtitle('连续信号及其抽样信号')函数图像为:)502cos()(2t t x ⨯=π同理,函数图像为:)0102cos()(3t t x ⨯=π同理,函数图像为:由以上的三图可知,第一个图的离散序列,基本可以显示出原来信号,可以通过低通滤波恢复,因为信号的频率为20HZ,而采样频率为50>2*20,故可以恢复,但是第二个和第三个信号的评论分别为50和100HZ,因此理论上是不能够恢复的,需要增大采样频率,解决的方案为,第二个信号的采样频率改为400HZ,而第三个的采样频率改为1000HZ,这样可以很好的采样,如下图所示:2. 产生幅度调制信号)200cos()2cos()(t t t x ππ=，推导其频率特性，确定抽样频率，并绘制波形。
信号与系统采样实验报告
实验5采样采样定理给定了一些条件，在这些条件之下，一个带限的连续时间信号能够完全用它的离散样本表示。
所得到的离散时间信号)(][nT x n x c =包含了在连续时间信号中的全部信息。
只要这个连续时间信号是充分在频率上带限的，即T j X c π≥Ω=Ω,0)(。
当满足这一条件时，原连续时间信号能够完全用样本][n x 之间的内插予以重建。
如果][n x 满足采样定理，就有可能完全在离散时间域中处理][n x 而得到另一个序列，这个序列本该以不同的采样率对)(t x c 采样而得到。
这个处理称为采样率转换。
离散时间系统的灵活性对于连续时间LTI 系统的实现提供了一种强有力的手段，这就是连续时间信号的离散时间系统处理。
在这一技术中，一个带限的连续时间输入被采样，用一个离散时间系统所得到的样本，然后将这个离散时间系统的输出样本进行内插，给出连续时间输出信号。
本章练习将研究涉及信号采样和重建中的许多问题。
注意，该章用Ω代表连续时间频率变量，而用ω代表离散时间频率变量。
§5.1由欠采样引起的混叠目的这个练习讨论信号经采样后其频谱的变化以及由于欠采样而在而在带限内插重建信号上引起的混叠效果。
相关知识如果一个连续时间信号)(t x 每隔T 秒采样一次，那么信号的样本就形成了离散时间序列)(][nT x n x =。
奈奎斯特采样定理说的是，如果)(t x 的带宽小于s π=Ω2，即2,0)(s c j X Ω≥Ω=Ω，那么)(t x 就完全可以由它的样本)(nT x 予以重建。
带限内插或信号重建是最容易将)(t x 首先乘以冲激串后而看出来的 ∑∞-∞=-=n p nT t nT x t x )()()(δ 用一个截止频率2s Ω的理想低通滤波器对)(t x p 滤波，就能从)(t x p 中将)(t x 恢复出来。
定义)(t x r 为低通过滤)(t x p 而得到的重建信号。
若)(t x 的带宽大于2s Ω，那么样本)(nT x 就不能完全确定)(t x ，)(t x r 一般说来不等于)(t x 。
信号与系统实验四-信号的采样及恢复
实验四 信号的采样及恢复一、实验目的1、加深理解连续时间信号离散化过程中的数学概念和物理概念；2、掌握对连续时间信号进行抽样和恢复的基本方法；3、通过实验验证抽样定理。
二、实验内容1、为了观察连续信号时域抽样时，抽样频率对抽样过程的影响，在[0，0.1]区间上以50Hz 的抽样频率对下列3个信号分别进行抽样，试画出抽样后序列的波形，并分析产生不同波形的原因，提出改进措施。
（1）)102cos()(1t t x ⨯=π（2）)502cos()(2t t x ⨯=π （3）)1002cos()(3t t x ⨯=π2、产生幅度调制信号)200cos()2cos()(t t t x ππ=，推导其频率特性，确定抽样频率，并绘出波形。
3、对连续信号)4cos()(t t x π=进行抽样以得到离散序列，并进行重建。
（1）生成信号)(t x ，时间t=0:0.001:4，画出)(t x 的波形。
（2）以10=sam f Hz 对信号进行抽样，画出在10≤≤t 范围内的抽样序列)(k x ；利用抽样内插函数)/1()(sam r f T T t Sa t h =⎪⎭⎫⎝⎛=π恢复连续信号，画出重建信号)(t x r 的波形。
)(t x 与)(t x r 是否相同，为什么？ （3）将抽样频率改为3=sam f Hz ，重做（2）。
4、利用MATLAB 编程实现采样函数Sa 的采样与重构。
三、实验仪器及环境计算机1台，MATLAB7.0软件。
四、实验原理对连续时间信号进行抽样可获得离散时间信号，其原理如图8-1。
采样信号)()()(t s t f t f s ∙=，)(t s 是周期为s T 的冲激函数序列，即)()()(∑∞-∞=-==n sT nT t t t s sδδ则该过程为理想冲激抽样。
其中s T 称为采样周期，ss T f 1=称为抽样频率， ss s T f ππω22==称为抽样角频率。
信号和系统实验_信号抽样和内插
武汉大学教学实验报告电子信息学院电子信息学院专业 2014 年 11 月 27 日实验名称信号的抽样与内插指导教师姓名 ** 年级大三学号 201230******* 成绩6）画出各信号的频谱图，程序代码如下：N=length(time); %离散点的个数Ts=(time(N) - time(1))/N; %抽样周期m=floor(N/2); %因为DFT是对称的，只需要取一半 Ws=2*pi/Ts;W=Ws*(0:m)/N;F=fft(z1,N); FF=F(1:m+1); F11=abs(FF);plot(W,F14,'b',-W,F14,'b');%axis([-500,500,0,750])title('恢复后信号的幅频特性');xlabel('频率（Rad/s）');1.输入信号为1Hz的正弦波，波形如左下图，经过1.1Hz低通滤波后的波形如右下图在临界抽样频率2Hz下的到的抽样波形如左下图，经过1.1Hz滤波后得到的波形如右下图临界抽样频率下的频谱如下图2.当抽样频率变为1.5Hz时，为欠抽样，抽样波形如左下图，经过1.1Hz低通滤波后得到的波形如右下图1.5Hz欠抽样情况下得到的频谱如下图3.当抽样频率变为5Hz时，为过抽样，抽样波形如下图经过1.1Hz低通滤波后得到的波形如右下图过采样的频谱如下图4.当信号波形为方波时，左下图为1Hz方波信号波形，右下图为经过10Hz滤波后方波的波形，经过预滤波后，方波中高频成分被滤除，波形变得有些失真，但是和最后抽样滤波的结果更加近似左下图为没有经过滤波的方波被20Hz抽样信号抽样，经过10Hz低通滤波器后的波形如右下图，抽样滤波时在从抽样信号中提取方波信号时方波中的高平成分也被滤除，因而最终结果更加接近于预滤波后的方波波形。
波形的频谱如下图5.当信号波形为锯齿波时，左下图为1Hz锯齿波信号波形，右下图为经过10Hz滤波后锯齿波的波形，经过预滤波后，方波中高频成分被滤除，波形变得有些失真，但是和最后抽样滤波的结果更加近似下图1 为没有经过滤波的锯齿波被20Hz抽样信号抽样，经过10Hz低通滤波器后的波形如下图2 ，抽样滤波时在从抽样信号中提取方波信号时锯齿波中的高平成分也被滤除，因而最终结果更加接近于预滤波后的锯齿波波形。
《信号与系统实验》信号的采样与恢复(抽样定理)实验一
《信号与系统实验》信号的采样与恢复（抽样定理）实验一、实验目的1、了解电信号的采样方法与过程以及信号恢复的方法。
2、验证抽样定理。
二、实验设备1、信号与系统实验箱2、双踪示波器三、原理说明1、离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。
抽样信号f s(t)可以看成连续f(t)和一组开关函数s (t)的乘积。
s (t)是一组周期性窄脉冲，见实验图5-1，T s(t)称为抽样周期，其倒数f s(t)= 1/T s称为抽样频率。
图5-1 矩形抽样脉冲对抽样信号进行傅立叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的信号频率。
平移的频率等于抽样频率f s(t)及其谐波频率2f s、3f s》》》》》》。
当抽样信号是周期性窄脉冲时，平移后的频率幅度(sinx)/x规律衰减。
抽样信号的频谱是原信号频谱周期的延拓，它占有的频带要比原信号频谱宽得多。
2、正如测得了足够的实验数据以后，我们可以在坐标纸上把一系列数据点连起来，得到一条光滑的曲线一样，抽样信号在一定条件下也可以恢复到原信号。
只要用一截止频率等于原信号频谱中最高频率f n的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器输出可以得到恢复后的原信号。
3、但原信号得以恢复的条件是f s 2B，其中f s为抽样频率，B为原信号占有的频带宽度。
而f min=2B为最低抽样频率又称“柰奎斯特抽样率”。
当f s<2B时，抽样信号的频谱会发生混迭，从发生混迭后的频谱中我们无法用低通滤波器获得原信号频谱的全部内容。
在实际使用中，仅包含有限频率的信号是及少的，因此即使f s=2B，恢复后的信号失真还是难免的。
图5-2画出了当抽样频率f s>2B（不混叠时）f s<2B（混叠时）两种情况下冲激抽样信号的频谱。
t f(t)0F()t 0m ωm ω-(a)连续信号的频谱Ts t 0f s (t)F()t0m ωm ω-s ω-s ω（）(b)高抽样频率时的抽样信号及频谱 不混叠图5-2 冲激抽样信号的频谱实验中f s >2B 、f s =2B 、f s <2B 三种抽样频率对连续信号进行抽样，以验证抽样定理——要使信号采样后能不失真地还原，抽样频率f s 必须大于信号频率中最高频率的两倍。
信号与系统实验报告总结
信号与系统实验报告总结信号与系统实验报告实验信号抽样及信号重建⼀、实验⽬的1、进⼀步理解信号的抽样及抽样定理;2、进⼀步掌握抽样信号的频谱分析;3、掌握和理解信号抽样以及信号重建的原理;⼆、实验内容及步骤练习1、什么是抽样定理,信号采样后重建的步骤，抽样频率如何设置, 答:(1).抽样,Sampling,:就是从连续时间信号中抽取⼀系列的信号样本，从⽽，得到⼀个离散时间序列(Discrete-time sequence).抽样定理:设时间连续信号f(t)，其最⾼截⽌频率为f m ，如果⽤时间间隔为T<=1/2fm的开关信号对f(t)进⾏抽样时，则f(t)就可被样值信号唯⼀地表⽰。
(2).步骤:从频域看,设信号最⾼频率不超过折叠频率:,,,,Xa(j)=Xa(j) ||,,,Xa(j)=0 ||>s/2则理想取样后的频谱就不会产⽣混叠,故有:让取样信号x^(t)通过⼀带宽等于折叠频率的理想低通滤波器:,,,H(j)=T ||,,,H(j)=0 ||>s/2滤波器只允许通过基带频谱,即原信号频谱,故:,,,,Y(j)=X^(j)H(j)=Xa(j)因此在滤波器的输出得到了恢复的原模拟信号:y(t)=xa(t)从时域上看,上述理想的低通滤波器的脉冲响应为:根据卷积公式可求得理想低通滤波器的输出为:由上式显然可得:则:上式表明只要满⾜取样频率⾼于两倍信号最⾼频率,连续时间函数xa(t)就可⽤他的取样值xa(nT)来表达⽽不损失任何信息，这时只要把每⼀个取样瞬时值与内插函数式相乘求和即可得出xa(t),在每⼀取样点上,由于只有该取样值所对应的内插函数式不为零,所以各个取样点上的信号值不变。
(3).频率设置:根据抽样定理 ws/wm的值必须⼤于或等于2练习2、给范例程序Program4_1加注释。
% Programclear, close all,tmax = 4; dt = 0.01;t = 0:dt:tmax;Ts = 1/10; % Sampling periodws = 2*pi/Ts; % Sampling frequencyw0 = 20*pi; dw = 0.1; % The frequency of x(t)w = -w0:dw:w0;n = 0:1:tmax/Ts; % Make the time variable to be discretex = exp(-4*t).*u(t);xn = exp(-4*n*Ts); % The sampled version of x(t)subplot(221) % Plot the original signal x(t)plot(t,x), title('A continuous-time signal x(t)'),xlabel('Time t'), axis([0,tmax,0,1]), grid onsubplot(223) % Plot xnstem(n,xn,'.'), title('The sampled version x[n] of x(t)'), xlabel('Time index n'), axis([0,tmax/Ts,0,1]), grid on Xa = x*exp(-j*t'*w)*dt;X = 0;for k = -8:8; % Periodically extend X to form a periodic signal X = X + x*exp(-j*t'*(w-k*ws))*dt;endsubplot(222) % Plot xaplot(w,abs(Xa))title('Magnitude spectrum of x(t)'), grid onaxis([-60,60,0,1.8*max(abs(Xa))])subplot(224)plot(w,abs(X))title('Magnitude spectrum of x[n]'), xlabel('Frequency inradians/s'),grid on axis([-60,60,0,1.8*max(abs(Xa))])练习3、分别进⾏设置ws/wm= 2，ws/wm= 1，ws/wm= 3，并运⾏抽样信号重建程序，并根据抽样定理及重建条件分析三种设置情况下的结果。
信号与系统实验报告——采样(含程序)
电 子 科 技 大 学实 验 报 告学生姓名：苏晓菁 学 号： 2804301026 指导教师：张鹰 一、实验室名称：信号与系统实验室 二、实验项目名称：连续信号的采样和恢复 三、实验原理：实际采样和恢复系统如图3.6-1所示。
可以证明，奈奎斯特采样定理仍然成立。
x )(t P T)图3.6-1 实际采样和恢复系统采样脉冲：Tsπω2=，其中，2/)2/sin(τωτωτs s k k k Ta =，T <<τ。
采样后的信号：∑∞-∞=-=−→←k s SFS k j X Tj Xt x )((1)()(ωωω当采样频率大于信号最高频率两倍，可以用低通滤波器)(ωj Hr由采样后的()()2()FT T k s k p t P j a k ωπδωω+∞=-∞←−→=-∑信号)(txS 恢复原始信号)(tx。
四、实验目的：目的：使学生通过采样保持电路理解采样原理。
使学生理解采样信号的恢复。
五、实验内容：实验内容（一）、采样定理验证实验内容（二）、采样产生频谱交迭的验证六、实验器材（设备、元器件）：数字信号处理实验箱、信号与系统实验板的低通滤波器模块U11、高通滤波器模块U21、PC机端信号与系统实验软件、＋5V电源七、实验步骤：实验内容（一）、采样定理验证1、连接接口区的“输入信号1”和“输出信号”，如图1所示。
图12、信号选择：按“3”选择“正弦波”，再按“＋”或“－”设置正弦波频率为“2.6kHz”。
按“F4”键把采样脉冲设为10kHz。
3、观察原始正弦波。
4、按图2的模块连线示意图连接各模块。
5、观察采样后的波形。
6、用截止频率为3kHz的低通滤波器U11恢复采样后的信号。
按图3的模块连线示意图连接各模块。
图37、观察恢复后的波形。
实验内容（二）、采样产生频谱交迭的验证 重复实验内容（一）的实验步骤；注意在第2步中正弦波的频率仍设为“2.6kHz ”后，按“F4”键把采样脉冲频率设为“5kHz ”；在第6步中用3kHz 的恢复滤波器（U11）。
1_1信号的抽取与内插
内插矩阵[2]的列
第0列 第1列 第n列
[k] [k-2] [k-2n]
内插矩阵[2] k行n列
[ 2]k ,n [k 2n]
抽取矩阵与内插矩阵
1 0 0 [ 2] 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0
利用MATLAB实现序列内插
1 0.5 0 -0.5 -1 0
2
4
6
8
10
12
14
16
18
20
1
0.5
0 -0.5 -1 0
10
20
30
40
50
60
抽取的矩阵描述
x0 x 1 x 1 0 0 0 0 0 0 x 0 0 1 0 0 0 x2 [ 2] x 2 x 3 x 0 0 0 0 1 0 4 x 4 x5
X D ( z) x[kM ]z k
k
1 X D (z) M
M 1
l0
X (z W )
1 M
l M
X D (e
j
1 ) M
M 1
l0
X (e
j
2 πl M
)
M倍抽取后频谱的变换规律
X D(e )
X (e
j
j
1 M
M 1
l 0
X (e
j
2πl
k
k k是L的整数倍
x[k / L]z k
x[n]z nL
n
X I(z) X (z )
信号与系统 抽样定理实验
信号与系统实验报告实验六抽样定理实验六抽样定理一、实验内容:(60分)1、阅读并输入实验原理中介绍的例题程序,观察输出的数据与图形,结合基本原理理解每一条语句的含义。
2、已知一个连续时间信号f(t)=sinc(t),取最高有限带宽频率f m=1Hz。
(1)分别显示原连续信号波形与F s=f m、F s=2f m、F s=3f m三种情况下抽样信号的波形;程序如下:dt=0、1;f0=0、2;T0=1/f0;fm=5*f0;Tm=1/fm;t=-10:dt:10;f=sinc(t);subplot(4,1,1);plot(t,f);axis([min(t),max(t),1、1*min(f),1、1*max(f)]);title('ÔÁ¬ÐøÐźźͳéÑùÐźÅ');for i=1:3;fs=i*fm;Ts=1/fs;n=-10:Ts:10;f=sinc(n);subplot(4,1,i+1);stem(n,f,'filled');axis([min(n),max(n),1、1*min(f),1、1*max(f)]); end运行结果如下:（2）求解原连续信号与抽样信号的幅度谱;程序: dt=0、1;fm=1;t=-8:dt:8;N=length(t);f=sinc(t);wm=2*pi*fm;k=0:N-1;w1=k*wm/N;F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1、1*min(abs(F1)),1、1*max(abs(F1))]);for i=1:3;if i<=2 c=0;else c=1;endfs=(i+c)*fm;Ts=1/fs;n=-6:Ts:6;N=length(n);f=sinc(n);wm=2*pi*fs;k=0:N-1;w=k*wm/N;F=f*exp(-1i*n'*w)*Ts;subplot(4,1,i+1);plot(w/(2*pi),abs(F));axis([0,max(4*fm),0、5*min(abs(F)),1、1*max(abs(F))]); end波形如下:（3）用时域卷积的方法(内插公式)重建信号。
数字信号处理讲义-信号的抽取与内插
j2πl
X(e M
)
12
M倍抽取后频谱的变换规律
XD(ej)M 1M l01
2πl
j
X(e M )
X (e j
)
扩 M 倍
X
j
(e M
)
周 期 化 2π为
1 M1
2πl
j
X(e M )
M l0
13
证明
~M[k]
M1
1 kl WM
M l0
XD(z)x[kM ]zk
n
x[n]z M
k
n是M的整数倍
1X (ej( )
13 X D (ej )
序列抽取不混叠的条件 X(ej)=0，||>/M15
1 X(ej)
X(ej) 1
1
X(ej()
2XD(ej) 1
2倍抽取产生的频谱混叠
16
抽取和内插的变换域描述
(b) L倍内插
XI(z) xI[k]zk
Ml0
H(z)M1
M l0
1
X(zMWM l )
20
内插等式
x[k ] L
H (z L ) y3[k]
x[k ] H (z)
y4[k] L
Y3(z)X(zL)H(zL) Y4(z)X(z)H(z)LX(zL)H(zL)
21
基本单元的连接
x[k ]
L v1[k] M y1[k] ?
x[k] M v2[k] L
0
3
6
9
k
xD[k]
k
0
1
2
3
5
例： M倍抽取是时变系统。
x[k ]
xD [k], M 2
信号抽样与内插
3）分别在正常采样与欠采样条件下，配置各模块的参数（如信号源的频率，抽样脉冲的间隔，低通滤波器的截止频率等）。
4）在模型文件的菜单中选择Simulation－>Start，运行在正常采样、与欠采样条件下的仿真模型；
5）仿真结束后，打开示波器，观察在正常采样与欠采样条件下的仿真结果。
F=fft(z3,N); FF=F(1:m+1); F13=abs(FF);
F=fft(z4,N); FF=F(1:m+1); F14=abs(FF);
subplot(221);
plot(W,F11,'b',-W,F11,'b');
title('输入信号的幅频特性');
xlabel('频率（Rad/s）');
subplot(222);
plot(W,F12,'b',-W,F12,'b');
title('滤波后信号的幅频特性');
xlabel('频率（Rad/s）');
subplot(223);
plot(W,F13,'b',-W,F13,'b');
title('抽样后信号的幅频特性');
xlabel('频率（Rad/s）');
2、因输入信号是周期的，所以频谱都是离散的，以方波过抽样为例，基波为2*pi，频谱间隔为4*pi，抽样后将频谱搬移200*pi的整数倍。
3、抽样后信号包络是Sa函数，是因为实验采用方波抽样；另外，占空比越大，抽样后信号幅度越大。
抽取和内插
多速率信号处理及抽取和内插一：多速率信号处理1、在信号处理系统中有时需要不同的抽样率，这样做的目的有时是为了适应不同系统之间的级联，以利于信号的处理、编码、传输和存储，有时则是为了节省计算工作量。
数据速率的转换两种途径：1）数字信号→数模转换→模拟信号→模数转换→另一抽样率抽样2）数字信号处理→数字信号处理基本方法→抽样率转换目的：改变原有数字信号的频率方法：抽取和内插，低通滤波。
低通滤波：抽取和内插的前提条件是信号频带内没有频谱混叠，实现这一点需要用到低通滤波。
2、多速率滤波器-->具有线性相位的FIR滤波器。
常用的多速率滤波器：多速率FIR滤波器，积分梳状滤波器（CIC）和半带滤波器(HB);3、常用多速率信号处理结构第一级：CIC滤波器。
用于实现抽取和低通滤波第二级：fir实现的半带滤波器优点：工作在较低频率下，且滤波器参数得到优化，更容易以较低阶数实现，达到节省资源，降低功耗的目的。
二：抽取概念：使抽样率降低的转换。
1、整数倍抽取当信号的抽取数据量太大时，为了减少数据量以便于处理和计算，我们把抽样数据每隔（D-1）个取一个，这里D是一个整数。
这样的抽取称为整数抽取，D称为抽取因子。
2、抽取后结果：信号的频谱：信号的频谱周期降低1/D;信号的时域：信号的时域每D个少了（D-1）信号。
3、抗混叠滤波：在抽取前，对信号进行低通滤波，把信号的频带限制在抽样后频率的一半以下，这样，整数倍抽取的的问题就变成了一个低通滤波的问题。
信号时域图信号频域图程序运行后所得到的滤波前后信号的时域图，滤波器的频率响应图如上图。
从图中可以看出，经半带滤波器滤波后的信号，与原信号相比，波形没有改变，但抽样速率降低了一半；半带滤波器通阻带容限相同，具有严格线性相位。
三：内插概念：使抽样率升高的转换。
1、整数倍内插：在已知的相邻抽样点之间等间隔插入（I-1）个零值点。
然后进行低通滤波，即可求得I倍内插的结果。
2、内插后结果：信号的时域：已知抽样序列的两相邻抽样点之间等间隔多了I-1个值信号的频谱：信号的频谱周期增加了I倍。
ch7_1信号的抽取与内插
Y2 ( z) X ( z)H ( z M ) M
1
M 1
1
1
X ( z M WMl ) H (( z M WMl ) M )
M l0
H (z) M 1 M l0
1
X (z M WMl )
基本单元
内插等式
x[k ] L
H (z L ) y3[k]
x[k ] H (z)
y4[k] L
n
X I (z) X (z L ) XI(ej)= X(ejL)
基本单元
XI(ej)= X(ejL)
L=5时内插序列的频谱
1 X(ej)
镜像
1 XI(ej)
镜像
基本单元的连接
M
N
y[k]
x1[k]
1
x2[k]
2
M
y[k]
基本单元
基本单元的连接
x[k]
L v1[k] M y1[k] ?
x[k]
M
v2 [k ] L
y2[k]
如M和L互素，即M和L无公因子，则上述两种级联等价。
V1(ej ) X(ejL )
V2(ej )
1 M
M 1 k0
2πk j
X(e M )
Y1(ej )
Y3 (z) X (z L )H (z L )
Y4(z) X (z)H(z) L X (z L )H (z L )
基本单元
基本单元的连接
x[k]
L v1[k] M y1[k] ?
x[k]
M
v2 [k ] L
y2[k]
例: L=M=2
信号与系统实验报告2015解析
信号与系统实验报告 13级通信四班赵豆学号：20131870实验一 连续时间信号的采样一、 实验目的1、进一步加深对采样定理和连续信号傅立叶变换的理解。
2、了解MATLAB 提供的在相邻点间内插的几种方法。
3、使用sinc 函数并且研究时域中的混叠问题二、 实验步骤1．复习采样定理和采样信号的频谱采样定理如果采样频率s F 大于有限带宽信号)(t x a 带宽0F 的两倍，即02F F s > （1）则该信号可以由它的采样值)()(s a nT x n x =重构。
否则就会在)(n x 中产生混叠。
该有限带宽模拟信号的02F 被称为乃魁斯特频率。
必须注意，在)(t x a 被采样以后，)(n x 表示的最高模拟频率为2/s F Hz （或πω=）。
2．熟悉如何用MATLAB 语言实现模拟信号表示严格地说，除了用符号处理工具箱(Symbolics)外，不可能用MATLAB 来分析模拟信号。
然而如果用时间增量足够小的很密的网格对()a x t 采样，就可得到一根平滑的曲线和足够长的最大时间来显示所有的模态。
这样就可以进行近似分析。
令t ∆是栅网的间隔且s t T ∆<<，则()()G a x m x m t ∆=∆ （2）可以用一个数组来仿真一个模拟信号。
不要混淆采样周期s T 和栅网间隔t ∆，因为后者是MATLAB 中严格地用来表示模拟信号的。
类似地，付利叶变换关系也可根据（2）近似为：∑∑∆Ω-∆Ω-∆=∆≈Ωmt m j G m t m j G a e m x t t em x j X )()()( （3） 现在，如果)(t x a （也就是)(m x G ）是有限长度的。
则公式（3）与离散付利叶变换关系相似，因而可以用同样的方式以MATLAB 来实现，以便分析采样现象。
3．根据提供的例子程序，按照要求编写实验用程序；三、实验内容（1）通过例一熟悉用MATLAB 语言实现描绘连续信号的频谱的过程，并在MATLAB 语言环境中验证例1的结果；例1 令t a e t x 1000)(-=，求出并绘制其付利叶变换。
信号与系统实验报告模板6
信号与系统实验报告模板6武汉大学教学实验报告电子信息学院电子信息工程专业 2015 年 9 月 25 日实验名称信号的抽样与内插指导教师姓名年级 2013级学号成绩一、预习部分1.实验目的2.实验基本原理3.主要仪器设备（含必要的元器件、工具）一、实验目的1.熟悉信号的抽样与恢复过程2.观察欠采样与过采样时信号频谱的变化3. 掌握采样频率的确定方法二、实验基本原理由时域抽样定理可知，若有限带宽的连续时间信号f (t)的最高角频率为Wm，信号f (t)可以用等间隔的抽样值唯一表示，且抽样间隔T s必须不大于1/2 f m , 或者说抽样频率ωs ≥2ωm。
下图所示为信号抽样与恢复示意图，其中(a)中为抽样前带限信号f (t)，其频谱F(ω)为图(b)所示，最高频率为ωm。
当该信号被抽样间隔为T s的冲激序列抽样时，若s T 大于1/2 f m（过采样），则抽样后信号f (t) s 的频谱为图(f)所示，频谱没有产生混迭现象。
将抽样后信号f s(t)通过一个低通滤波器，能恢复原信号f (t)。
若T s小于1/2 f m（欠采样），则抽样后信号f (t) s 的频谱将产生混迭现象，不能从抽样后信号f (t) s 中恢复原信号f (t).三、主要仪器设备（含必要的元器件、工具）要用到的matlab函数及工具箱1. Simulink 仿真利用Simulink 完成信号的抽样与内插实验仿真设计。
2. fft 函数功能：离散傅里叶变换。
调用格式：y = fft(x, n)3. abs 函数功能：求绝对值和复数的模。
调用格式：y = abs(x)二、实验操作部分1.实验数据、表格及数据处理2.实验操作过程（可用图表示）3.实验结论一、实验数据、表格及数据处理1、方波的抽样与恢复（过抽样）方波频率1Hz，振幅为1，第一级低通滤波器截止频率50Hz 抽样信号为120Hz,占空比为50%的矩形波方波信号源的波形图（f=1Hz amplitude=1）方波抽样后波形恢复的方波波形方波频谱图2、三角波的抽样与恢复（欠抽样）三角波频率1Hz，振幅为1，滤波器截止频率为25Hz，抽样信号为频率40Hz，占空比为50%的矩形波。
实验三 信号的抽样和内插
∫
p(t)e− jnωst dt
抽样后信号：fs (t) = f (t) p(t)
∞ 1 ∴Fs (ω) = F(ω) * P(ω) = ∑P F(ω − nωs ) n 2π n=−∞
冲激抽样
f (t )
FT
(1)
0 P (t )
1
∞
F (ω )
t
FT
0
P (ω ) P(ω) = 2π ∑Pδ (ω − nωs ) n
信号在时域被冲激抽样后，其频谱Fs(ω)是 ω 原信号的频谱F(ω)以抽样频率ωs为间隔周期 重复而得到的。 在重复过程中，F(ω)幅度被抽样脉冲p(t)的 傅立叶系数所加权，加权系数取决于抽样脉冲 序列的形状。由于冲激序列的傅立叶系数为常 数，所以冲激抽样后的频谱是将原信号的频谱 F(ω)以ωs为周期等幅地重复。
实验三 信号的抽样与内插
实验目的： 实验目的
熟悉信号采样过程，并通过本实验观察 欠采样时信号频谱的混迭现象，了解采 样前后信号频谱的变化，加深对采样定 理的理解，掌握采样频率的确定方法。
实验原理
设：f (t) ↔ F(ω), p(t) ↔ P(ω), fs (t) ↔ Fs (ω)
p(t)为 期 号 周 为 s， ： 周 信 ， 期 T 有 1 P(ω) = 2π ∑Pδ (ω − nωs ), 其 P = 中n n Ts n=−∞
∞
f s (t )
抽样间隔T 抽样间隔 s
0
Fs (ω ) 1 Ts
ω s = 2ωm ωc = ω m
t
−ωs
0
ωs
H(ω)
ω
h(t)
Ts
Ts
t
信号与系统实验(MATLAB 西电版)实验14 信号的抽样与恢复
实验14 信号的抽样与恢复 图 14.5 信号采样的性质
实验14 信号的抽样与恢复
4) 模拟低通滤波器设计 MATLAB clf; Fp=3500; Fs=4500; Wp=2*pi*Fp; Ws=2*pi*Fs; ［N，Wn］=buttord(Wp，Ws，0.5，30，′s′); ［b，a］=butter(N，Wn，′s′); wa=0：(3*Ws)/511：3*Ws;
实验14 信号的抽样与恢复
11) MATLAB % clf; L=input(′过采样因子=′); M=input(′欠采样因子=′); n=0：29; x=sin(2*pi*0.43*n) + sin(2*pi*0.31*n); y=resample(x，L，M); subplot(2，1，1);
实验14 信号的抽样与恢复 图 14.3 正弦信号的采样
实验14 信号的抽样与恢复
2) MATLAB clf; T=0.1; f=13; n=(0：T：1)′; xs=cos(2*pi*f*n);
实验14 信号的抽样与恢复
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(′时间，msec′); ylabel(′幅值′); title(′重构连续信号y_{a}(t)′); axis(［0 1 -1.2 1.2］); 正弦信号的采样与重构结果如图14.4
实验14 信号的抽样与恢复
plot(w/pi，abs(Xz)); axis(［0 1 0 1］); grid title(′输入谱′); subplot(2，1，2); L=input(′过采样因子 =′); y=zeros(1，L*length(x)); y(［1： L： length(y)］)=x; ［Yz，w］=freqz(y，1，512); plot(w/pi，abs(Yz)); axis(［0 1 0 1］); grid title(′输出谱′); 信号的频域过采样结果如图14.9
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Ｎ=ｌenｇtｈ(time); %离散点的个数
Ts=(tｉｍe(N） － ｔｉｍe(１））／Ｎ; %抽样周期
ｍ=flｏｏr(N/2); %因为DFT是对称的,只需要取一半
Ws=2*pｉ/Ｔs;
W=Ｗs*（０:ｍ)/N;
Ｆ=fft(z1,Ｎ); FF=F（1:m+1）； F1１=ａｂs(FF);
四、教师评语
指导教师年月日
三、实验效果分析（包括仪器设备等使用效果）
思考题
1.说明采样频率变化对信号时域和频域特性的影响。
答：对于连续频谱信号，采样频率Ｆｓ大于信号最大频率的2倍时,信号抽样把不会产生混迭，才能不失真恢复原波形
2.分析采样与内插仿真模型中两个低通滤波器的作用。
答：前者预滤波，滤除信号中的高频噪声，后者用于恢复抽样后的波形
１/（２fm)，或者说抽样频率ωs≥2ωm。
三、 涉及的 MＡTLAB 相关内容
1． Sｉmulｉnk 仿真
利用 Siｍulink 完成信号的抽样与内插实验仿真设计。
2. fft 函数
功能:离散傅里叶变换。
调用格式:y = fｆt（x， n)
3． abs 函数
功能:求绝对值和复数的模。
调用格式:ｙ ＝ abｓ（x)
二、实验操作部分
1.实验数据、表格及数据处理
2.实验操作过程(可用图表示）
3.实验结论
ﻬ
四、实验内容与方法
设计信号x（t)＝sin(２πft），f＝1Hz的抽样与恢复的实验，实验步骤如下：
1）在MATLAB命令窗口中输入“simulink”，启动Simuｌink LｉbｒａryBｒowｓｅｒ;
2)在Simｕlinｋ Libｒary Brｏｗser中，新建一个模型文件,编辑模型文件，
%FFt的结果为复数，求aｂ1:m＋1); F1２=abs(ＦF)；
F=fft(ｚ３，N); FF=F(１：m+1); F１3=abs(ＦF）;
F=fｆt(ｚ4,N）; ＦF＝Ｆ(1：m+1); Ｆ14＝ａbs(ＦF)；
sｕbplot(２,2,１)；
plot(W,F1１,＇b',-W,F11,＇b');
一、预习部分
1.实验目的
2.实验基本原理
3.主要仪器设备(含必要的元器件、工具)
一、实验目的
１.熟悉信号的抽样与恢复过程;
2.观察欠采样与过采样时信号频谱的变化;
3.掌握采样频率的确定方法。
二、实验基本原理
由时域抽样定理可知，若有限带宽的连续时间信号ｆ(t)的最高角频率为
ωm,则信号f(t)可以用等间隔的抽样值唯一表示，且抽样间隔Ts必须不大于
subplot(２，2,3);
plｏｔ(Ｗ，F13,＇ｂ',－Ｗ,F1３,'ｂ');
%axis([－500,5００,0,750])
titlｅ('抽样后信号的幅频特性');
xlabeｌ('频率(Rad/s）');
suｂｐlｏt(2,2,4);
plｏｔ(W,F14,＇b',－W，F14,'b');
%axｉs([－500，５00,0，750]）
1．5Ｈz欠抽样情况下得到的频谱如下图
3.当抽样频率变为５Hz时，为过抽样，抽样波形如下图
经过1.１Hz低通滤波后得到的波形如右下图
过采样的频谱如下图
4.当信号波形为方波时，左下图为1Hz方波信号波形,右下图为经过１0Ｈz滤波后方波的波形,经过预滤波后,方波中高频成分被滤除,波形变得有些失真，但是和最后抽样滤波的结果更加近似
信号与系统实验－信号抽样与内插
———————————————————————————————— 作者:
———————————————————————————————— 日期：
武汉大学教学实验报告
电子信息学院电子信息学院专业２０14年11月27日
实验名称信号的抽样与内插指导教师
姓名*＊年级大三学号2０12３０＊**＊＊**成绩
建立如图2所示的抽样与内插的仿真模型，并保存为saｍple.mdl；
3)分别在欠采样与过采样条件下，配置各模块的参数（如信号源的频率,抽样脉冲的间隔,低通滤波器的截止频率等)。
4）在模型文件的菜单中选择Sｉmulatｉoｎ－>Ｓtarｔ,运行在欠采样、与过采样条件下的仿真模型；
5）仿真结束后，打开示波器，观察在欠采样与过采样条件下的仿真结果。
%aｘｉｓ([－500，5０0,0,１500］)
titｌｅ（'输入信号的幅频特性＇）；
xlａbel('频率(Rad/s)'）;
subplot（2，2，2);
ｐloｔ(W,F12,'b＇,-W,F12,'b');
%aｘis([-50０,500,0，1500］）
tｉtle('滤波后信号的幅频特性＇）;
xｌaｂel('频率(Rad/s）')；
titlｅ(＇恢复后信号的幅频特性');
xｌａｂel('频率(Rad/s）');
1.输入信号为1Hz的正弦波，波形如左下图，经过１.1Ｈz低通滤波后的波形如右下图
在临界抽样频率2Hz下的到的抽样波形如左下图，经过1.1Hｚ滤波后得到的波形如右下图
临界抽样频率下的频谱如下图
2.当抽样频率变为1．5Hｚ时,为欠抽样,抽样波形如左下图，经过1．１Hz低通滤波后得到的波形如右下图
左下图为没有经过滤波的方波被2０Hz抽样信号抽样,经过１0Ｈｚ低通滤波器后的波形如右下图,抽样滤波时在从抽样信号中提取方波信号时方波中的高平成分也被滤除，因而最终结果更加接近于预滤波后的方波波形。
波形的频谱如下图
5.当信号波形为锯齿波时,左下图为1Hｚ锯齿波信号波形，右下图为经过１０Ｈｚ滤波后锯齿波的波形,经过预滤波后，方波中高频成分被滤除,波形变得有些失真,但是和最后抽样滤波的结果更加近似
下图１为没有经过滤波的锯齿波被20Hz抽样信号抽样,经过１０Hｚ低通滤波器后的波形如下图2 ，抽样滤波时在从抽样信号中提取方波信号时锯齿波中的高平成分也被滤除，因而最终结果更加接近于预滤波后的锯齿波波形。
锯齿波经过20Hz抽样的频谱如下图
锯齿波和方波相似,频谱分布广，抽样时不可避免会发生频谱混叠,最后滤波结果会出现失真。