实验三 信号的采样与恢复
()t s
实验三 信号的采样与恢复
一、实验目的
1、了解电信号的采样方法与过程以及信号恢复的方法。
2、验证抽样定理。
二、实验内容
1、观察抽样脉冲、抽样信号、抽样恢复信号。
2、观察抽样过程中，发生混叠和非混叠时的波形。
三、实验原理
1、离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。抽样信号
()t f s 可以看成连续信号()t f 和一组开关函数()t s 的乘积。()t s 是一组周期性窄脉冲，见图
2-1-1，T S 称为抽样周期，其倒数S s T f 1
=称抽样频率。
τ S T
图 2-1-1矩形抽样脉冲
对抽样信号进行傅里叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的原信号频率。平移的频率等于抽样频率s f 及其谐波频率s f 2、s f 3……。当抽样信号是周期性窄脉冲时，平移后的频率幅度按
()
x x sin 规律衰减。
抽样信号的频谱是原信号频谱周
期的延拓，它占有的频带要比原信号频谱宽得多。
2、正如测得了足够的实验数据以后，我们可以在坐标纸上把一系列数据点连起来，得到一条光滑的曲线一样，抽样信号在一定条件下也可以恢复到原信号。只要用一截止频率等于原信号频谱中最高频率f n 的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器输出可以得到恢复后的原信号。
3、但原信号得以恢复的条件是B f s 2≥，其中s f 为抽样频率，B 为原信号占有的频带宽度。而B f 2min =为最低抽样频率又称“奈奎斯特抽样率”。当B f s 2<时，抽样信号的频谱会发生混迭，从发生混迭后的频谱中我们无法用低通滤波器获得原信号频谱的全部内容。在实际使用中，仅包含有限频率的信号是极少的。因此即使B f s 2=，恢复后的信号失真还是难免的。图2-2-2画出了当抽样频率B f s 2≥（不混叠时）及当抽样频率B f s 2<（混
t
叠时）两种情况下冲激抽样信号的频谱。
1
(a) 连续信号的频谱
1 0
（b ） 高抽样频率时的抽样信号及频谱（不混叠）
1
（c ） 低抽样频率时的抽样信号及频谱（混叠）
图2-2-2 抽样过程中出现的三种情况
4、为了实现对连续信号的抽样和抽样信号的复原，可用实验原理框图2-2-3的方案。除选用足够高的抽样频率外，常采用前置低通滤波器来防止原信号频谱宽而造成抽样后信号频谱的混叠。但这也会造成失真。如实验选用的信号频带较窄，则可不设前置低通滤波器。本实验就是如此。
图 2-2-3 抽样定理实验方框图 四、实验说明
1、把系统时域与频域分析模块插在主板上，用导线接通此模块“电源接入”和主板上的电源（看清标识，防止接错，带保护电路），并打开此模块的电源开关。
2、将函数信号发生器产生一正弦波幅度（峰－峰值）为2V 左右，为便于观察，被抽
S
T 1
m ω- m ω
ω
()ωF
t
()t f
m ω- m ω
s ω
s ω-
ω
()ωs F 0 s T
t
()t f s
S
T 1
m ω- m ω
s ω
s ω- ω
()ωs F 0 s T
t
()t f s
信 号
输 入 抽样门 低 通
滤 波
抽 样
脉 冲
样信号频率一般选择50HZ ～400HZ 的范围，此时可以选择150Hz ，将其送入抽样器，即用导线将函数信号发生器的输出端与本实验模块的输入端相连；抽样脉冲的频段由开关SK1000进行选择，有“高”“中”“低”档，此时可以选择低频，频率则是通过电位器“频率调节”来调节的，调节抽样脉冲的频率为800Hz 左右，抽样脉冲的脉宽则是由电位器“脉宽调节”进行调节的，调节抽样信号脉宽为30％以上；用示波器测试“抽样信号”的波形，观察经抽样后的正弦波。
3、改变抽样脉冲的频率为B f s 2≥和B f s 2<，用导线将“抽样信号”和“低通输入”相连，用示波器测试测试钩“抽样恢复”，观察复原后的信号，比较其失真程度。
4、（对于要求高的学生可以进行以下实验）设计一定截至频率的低通滤波器，用于信号的抽样恢复。（可以参考第三章的实验六“模拟滤波器的设计”）
五、实验报告
1、整理并绘出原信号、抽样信号以及复原信号的波形，你能得出什么结论？
2、整理在三种不同抽样频率情况下，()t f s 波形，比较后得出结论。
3、实验调试中的体会。
六、思考题
1、若连续时间信号为50HZ 的正弦波，开关函数为T S =0.5ms 的窄脉冲，试求抽样后信号()t f s 。
2、设计一个二阶RC 低通滤波器，截止频率为1.5KHZ 左右。
3、若连续时间信号取频率为200HZ~300HZ 的方波和三角波，计算其有效的频带宽度。该信号经频率为s f 的周期脉冲抽样后，若希望通过低通滤波后的信号失真较小，则抽样频率和低通滤波器的截止频率应取多大，试设计一满足上述要求的低通滤波器。
语音信号处理实验指导书
语音信号处理实验指导书 实验一 语音信号采集与简单处理 一、 实验目的、要求 （1）掌握语音信号采集的方法 （2）掌握一种语音信号基音周期提取方法 （3）掌握短时过零率计算方法 （4）了解Matlab 的编程方法 二、 实验原理 基本概念： （a ）短时过零率： 短时内，信号跨越横轴的情况，对于连续信号，观察语音时域波形通过横轴的情况；对于离散信号，相邻的采样值具有不同的代数符号，也就是样点改变符号的次数。 对于语音信号，是宽带非平稳信号，应考察其短时平均过零率。 其中sgn[.]为符号函数 ?? ?? ?<=>=0 x(n)-1sgn(x(n))0 x(n)1sgn(x(n)) 短时平均过零的作用 1.区分清/浊音： 浊音平均过零率低，集中在低频端； 清音平均过零率高，集中在高频端。 2.从背景噪声中找出是否有语音，以及语音的起点。 （b ）基音周期 基音是发浊音时声带震动所引起的周期性，而基音周期是指声带震动频率的倒数。基音周期是语音信号的重要的参数之一，它描述语音激励源的一个重要特征，基音周期信息在多个领域有着广泛的应用，如语音识别、说话人识别、语音分析与综合以及低码率语音编码，发音系统疾病诊断、听觉残障者的语音指导等。因为汉语是一种有调语言，基音的变化模式称为声调，它携带着非常重要的具有辨意作用的信息，有区别意义的功能，所以，基音的提取和估计对汉语更是一个十分重要的问题。 ∑--= -=1 )]1(sgn[)](sgn[21N m n n n m x m x Z
由于人的声道的易变性及其声道持征的因人而异，而基音周期的范围又很宽，而同—个人在不同情态下发音的基音周期也不同，加之基音周期还受到单词发音音调的影响，因而基音周期的精确检测实际上是一件比较困难的事情。基音提取的主要困难反映在：①声门激励信号并不是一个完全周期的序列，在语音的头、尾部并不具有声带振动那样的周期性，有些清音和浊音的过渡帧是很难准确地判断是周期性还是非周期性的。②声道共振峰有时会严重影响激励信号的谐波结构，所以，从语音信号中直接取出仅和声带振动有关的激励信号的信息并不容 易。③语音信号本身是准周期性的(即音调是有变化的)，而且其波形的峰值点或过零点受共振峰的结构、噪声等的影响。④基音周期变化范围大，从老年男性的50Hz 到儿童和女性的450Hz ，接近三个倍频程，给基音检测带来了一定的困难。由于这些困难，所以迄今为止尚未找到一个完善的方法可以对于各类人群(包括男、女、儿童及不向语种)、各类应用领域和各种环境条件情况下都能获得满意的检测结果。 尽管基音检测有许多困难，但因为它的重要性，基音的检测提取一直是一个研究的课题，为此提出了各种各样的基音检测算法，如自相关函数(ACF)法、峰值提取算法(PPA)、平均幅度差函数(AMDF)法、并行处理技术、倒谱法、SIFT 、谱图法、小波法等等。 三、使用仪器、材料 微机（带声卡）、耳机，话筒。 四、 实验步骤 （1）语音信号的采集 利用Windows 语音采集工具采集语音信号，将数据保存wav 格式。 采集一组浊音信号和一组清音信号，信号的长度大于3s 。 （2）采用短时相关函数计算语音信号浊音基音周期，考虑窗长度对基音周期计算的影响。采用倒谱法求语音信号基音周期。 （3）计算短时过零率，清音和浊音的短时过零率有何区别。 五、实验过程原始记录（数据，图表，计算） 短时过零率 短时相关函数 P j j n s n s j R N j n n n n ,,1) ()()(1 =-=∑-= ∑--=-=10 )]1(sgn[)](sgn[21N m n n n m x m x Z
数字信号处理实验一
实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
（2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结
实验五 信号的采样与恢复
信号与系统实验报告 【实验原理】 1、离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。抽样信号f s (t )可以看成连续信号f (t )和一组开关函数s (t )的乘积。s (t )是一组周期性窄脉冲，见图1，T s 称为抽样周期，其倒数T s =1T S ?称抽样频率。 图1矩形抽样脉冲 对抽样信号进行傅里叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的原信号频率。平移的频率等于抽样频率f s 及其谐波频率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 ，恢复后的信号失真还是难免的。图2画出了当抽样频率f s ≥2B （不混叠时）及当抽样频率f s <2B （混叠时）两种情况下冲激抽样信号的频谱。 (a)连续信号的频谱
数字信号处理实验一
一、实验目的 1. 通过本次实验回忆并熟悉MATLAB这个软件。 2. 通过本次实验学会如何利用MATLAB进行序列的简单运算。 3. 通过本次实验深刻理解理论课上的数字信号处理的一个常见方法——对时刻n的样本附近的一些样本求平均，产生所需的输出信号。 3. 通过振幅调制信号的产生来理解载波信号与调制信号之间的关系。 二、实验内容 1. 编写程序在MATLAB中实现从被加性噪声污染的信号中移除噪声的算法，本次试验采用三点滑动平均算法，可直接输入程序P1.5。 2. 通过运行程序得出的结果回答习题Q1.31-Q1.33的问题，加深对算法思想的理解。 3. 编写程序在MATLAB中实现振幅调制信号产生的算法，可直接输入程序P1.6。 4. 通过运行程序得出的结果回答习题Q1.34-Q1.35的问题，加深对算法思想的理解。 三、主要算法与程序 1. 三点滑动平均算法的核心程序： %程序P1.5 %通过平均的信号平滑 clf; R=51; d=0.8*(rand(R,1)-0.5);%产生随噪声 m=0:R-1; s=2*m.*(0.9.^m);%产生为污染的信号 x=s+d';%产生被噪音污染的信号 subplot(2,1,1); plot(m,d','r-',m,s,'g--',m,x,'b-.');
xlabel('时间序号n');ylabel('振幅'); legend('d[n]','s[n]','x[n]'); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2); plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]'); xlabel('时间序号n');ylabel('振幅'); 2. 振幅调制信号的产生核心程序：（由于要几个结果，因此利用subplot函数画图） %程序P1.6 %振幅调制信号的产生 n=0:100; m=0.1;fH=0.1;fL=0.01; m1=0.3;fH1=0.3;fL1=0.03; xH=sin(2*pi*fH*n); xL=sin(2*pi*fL*n); y=(1+m*xL).*xH; xH1=sin(2*pi*fH1*n); xL1=sin(2*pi*fL1*n); y1=(1+m1*xL).*xH; y2=(1+m*xL).*xH1; y3=(1+m*xL1).*xH; subplot(2,2,1); stem(n,y); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.1;fH=0.1;fL=0.01;'); subplot(2,2,2); stem(n,y1); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.1;fL=0.01;'); subplot(2,2,3); stem(n,y2); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.3;fL=0.01;'); subplot(2,2,4); stem(n,y3); grid;
数字信号处理实验
实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π=
1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2）
语音信号处理实验报告
语音信号处理实验 班级： 学号： 姓名： 实验一基于MATLAB的语音信号时域特征分析（2学时）
1）短时能量 （1）加矩形窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=linspace(1,1,2.^(i-2)*N);%形成一个矩形窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if(i==2) ,legend('N=32'); elseif(i==3), legend('N=64'); elseif(i==4) ,legend('N=128'); elseif(i==5) ,legend('N=256'); elseif(i==6) ,legend('N=512'); end end
00.51 1.52 2.5 3 x 10 4 -1 1 x 10 4 024 x 10 4 05 x 10 4 0510 x 10 4 01020 x 10 4 02040 （2）加汉明窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=hanning(2.^(i-2)*N);%形成一个汉明窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if(i==2), legend('N=32'); elseif(i==3), legend('N=64'); elseif(i==4) ,legend('N=128');
数字信号处理实验答案完整版
数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
实验一熟悉Matlab环境 一、实验目的 1.熟悉MATLAB的主要操作命令。 2.学会简单的矩阵输入和数据读写。 3.掌握简单的绘图命令。 4.用MATLAB编程并学会创建函数。 5.观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉了MATLAB基本命令的基础上，完成以下实验。 上机实验内容： （1）数组的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。 clear all; a=[1 2 3 4]; b=[3 4 5 6]; c=a+b; d=a-b; e=a.*b; f=a./b; g=a.^b; n=1:4; subplot(4,2,1);stem(n,a); xlabel('n');xlim([0 5]);ylabel('A'); subplot(4,2,2);stem(n,b); xlabel('n');xlim([0 5]);ylabel('B'); subplot(4,2,3);stem(n,c); xlabel('n');xlim([0 5]);ylabel('C'); subplot(4,2,4);stem(n,d); xlabel('n');xlim([0 5]);ylabel('D'); subplot(4,2,5);stem(n,e); xlabel('n');xlim([0 5]);ylabel('E'); subplot(4,2,6);stem(n,f); xlabel('n');xlim([0 5]);ylabel('F'); subplot(4,2,7);stem(n,g); xlabel('n');xlim([0 5]);ylabel('G'); （2）用MATLAB实现下列序列： a) x(n)= 0≤n≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。
信号采样与重建的编程实现
课程设计任务书 学生：凯鑫专业班级：电信1203班 指导教师：阙大顺，王虹工作单位：信息工程学院 题目: 信号采集与重建的编程实现 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“数字信号处理原理与实现”、“Matlab及 在电子信息课程中的应用”等； 3.先修课程：信号与系统、数字信号处理、Matlab应用实践及信号处理类课程等。 要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计时间：1周（课实践）； 2.课程设计容：信号采样与重建的编程实现，具体包括：连续信号的时域采样、频谱混叠分析、 由离散序列恢复模拟信号等； 3.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具 体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结； 4.课程设计说明书按学校“课程设计工作规”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献； ⑦其它必要容等。 时间安排： 1）第1-2天，查阅相关资料，学习设计原理。 2）第3-4天，方案选择和电路设计仿真。 3）第4-5天，电路调试和设计说明书撰写。 4）第6天，上交课程设计成果及报告，同时进行答辩。
数字信号处理实验三
实验三 离散付里叶变换（DFT ） 一、实验目的： 1． 掌握离散付里叶级数 2． 掌握DFT 变换。 3． 掌握DFT 特性。 4． 掌握利用DFT 计算线性卷积。 5． 掌握快速付里叶变换(FFT)。 二、实验原理： 1．离散付里叶级数（DFS ） )(~n x 为周期序列，其频率为基本频率（N /2π ）的倍数（或谐波）。其离散付叶级 数（DFS ）为：；∑ -=-±== 1 2,,1,0,)(~)(~N n kn N j k e n x k X π IDFS 为：∑ -== 1 2)(~ 1)(~N k kn N j e k X N n x π 2．离散付里叶变换（DFT ） )(n x 为长度N 的有限长序列，其DFT 为：21 1()()N j kn N k x n X k e N π--== ∑ IDFT 为：21 ()(),0,1,N j kn N n X k x n e k π-== =±∑ 3.DFT 的特性： （1） 线性性：)]([)]([)]()([2121n x bDFT n x aDFT n bx n ax DFT +=+ （2） 循环折叠（圆周对称）性：?? ? -≤≤-==-1 1) (0) 0())((N n n N x n x n x N （3） 共轭性: N k X n x DFT ))(()]([* *-= （4） 实序列的对称性（圆周共轭对称性）：N k X k X ))(()(* -= （5） 序列的圆周移位：N m n x m n x ))(()(~ -=-
（6） 频域中的圆周移位：)())(()]([|ln k R l k X n x W DFT N N N -=- （7） 时域循环卷积：)()()]()([2121k X k X n x n x DFT =? （8） 频域循环卷积（乘法性）：)()(1)]()([2121k X k X N n x n x DFT ?= （9） 帕塞瓦尔（Parseval ）定理：∑ ∑ -=-== = 1 2 1 2 ) (1)(N k N n x k X N n x E 4.用DFT 计算线性卷积: 设)(1n x 为1N 点序列，)(2n x 为2N 点序列，)(3n x 为)(1n x 和)(2n x 的线性卷积，其为121-+N N 点序列，)(4n x 为)(1n x 和)(2n x 的圆卷积，其长度为 N ，当121-+=N N N 时，)()(43n x n x =。实际中，采用分段卷积法，即重 叠保留法和重叠相加法。需要对数据流进行分块处理，这时直接采用DFT 计算线性卷积会产生一些问题，而应该将)(n x 通过重复前M-1个取样进行分块，这样可得到正确结果。 5.快速付里叶变换（FFT ）： 掌握基2-时域抽取FFT(DIT-FFT)和基2-频域抽取FFT(DIF-FFT)。MATLAB 提供fft 函数来计算x 的DFT 。fft 函数是用机器语言写的，采用混合基法，其调用形式为：),(N x fft X =。如N 为2的幂，则得到高速的基2-FFT 算法；若N 不是2的乘方，则将N 分解成质数，得到较慢的混合基FFT 算法；最后，若N 为质数，则fft 函数采用的是原始的DFT 算法。 三、实验步骤： 1.离散付里叶级数（DFS ） （1）自已动手:编写实现离散付里叶级数和逆离散付里叶级数的函数。 （2）已知周期性序列如下所示：}3,2,1,0,3,2,1,0,3,2,1,0{)(~ ↑ =n x 求其离散付里叶级数。 2.离散付里叶变换（DFT ） (1) 编写实现DFT 和IDFT 的函数。 (2) 已知)(n x 是一个六点序列，如下所示： ?? ?≤≤=e ls e n n x 0 501 )( 要求计算该序列的离散时间的付里叶变换和离散付里叶变换，并绘出它们的幅度和相
数字信号处理实验三
实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期
x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response
信号实验：连续信号的采样和恢复
电子科技大学 实 验 报 告 学生姓名： 学号: 指导老师： 日期：2016年 12月 10日
一、实验室名称： 连续信号的采样和恢复 二、实验项目名称： 实验项目四：连续信号的采样和恢复 三、实验原理： 实际采样和恢复系统如图3.4-1所示。可以证明，奈奎斯特采样定理仍然成立。 ? ) x t ) (t P T ) 图3.4-1 实际采样和恢复系统 采样脉冲： 其中，T s πω2=， 2/)2/sin(τωτωτs s k k k T a =，T <<τ。 采样后的信号： ∑∞ -∞ =-=?→←k s S F S k j X T j X t x ) ((1)()(ωωω 当采样频率大于信号最高频率两倍，可以用低通滤波器)(ωj H r 由采样后的 ()()2() F T T k s k p t P j a k ωπδωω+∞ =-∞ ←?→= -∑
信号)(t x S 恢复原始信号)(t x 。 目的：1、使学生通过采样保持电路理解采样原理。 2、使学生理解采样信号的恢复。 任务：记录观察到的波形与频谱；从理论上分析实验中信号的采样保持与恢 复的波形与频谱，并与观察结果比较。 四、实验内容 实验内容（一）、采样定理验证 实验内容（二）、采样产生频谱交迭的验证 五、项目需用仪器设备名称：数字信号处理实验箱、信号与系统实验板的低通滤 波器模块U11和U22、采样保持器模块U43、PC 机端信号与系统实验软件、＋5V 电源 六、实验步骤： 打开PC 机端软件SSP.EXE ，在下拉菜单“实验选择”中选择“实验六”；使用串口电缆连接计算机串口和实验箱串口，打开实验箱电源。 实验内容（一）、采样定理验证 实验步骤： 1、连接接口区的“输入信号1”和“输出信号”，如图3.4-2所示。 图3.4-2 观察原始信号的连线示意图 2、信号选择：按“3”选择“正弦波”，再按“＋”或“－”设置正弦波频率为“2.6kHz ”。 按“F4”键把采样脉冲设为10kHz 。 七、实验数据及结果分析：
语音信号处理实验报告
通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程1502班 指导教师： 设计时间：2018/10/22-2018/11/23 评语： 通信与信息工程学院 二〇一八年 实验题目：语音信号分析与处理 一、实验内容 1. 设计内容 利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分
（1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。 （4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分 （5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。 （7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。 二、实验原理 1.设计原理分析 本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。 首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。 对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析 在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。
数字信号处理实验4
数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线
（2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：
数字信号处理实验六-时域采样与信号的重建
实验目的： 1.了解用MATLAB语言进行时域抽样与信号重建的方法 2.进一步加深对时域信号抽样与恢复的基本原理的理解 3.掌握采样频率的确定方法和内插公式的编程方法。 二．实验内容 1认真阅读并输入实验原理与方法中介绍的例子，观察输出波形曲线，理解每一条语句的含义。. 2．已知一个连续时间信号f(t)=sinc(t)。取最高有限带宽频率fm=1Hz。（1）分别显示原连续时间信号波形和Fm=fm、Fm=2fm、Fm=3fm三种情况下抽样信号的波形。 实验程序： dt=0.1; f0=1; T0=1/f0; fm=f0; Tm=1/fm; t=-2:dt:2; f=sinc(t); subplot(4,1,1),plot(t,f,'k'); 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=-2:Ts:2; f=sinc(n); subplot(4,1,i+1),stem(n,f,'filled','k'); axis([min(n) max(n) 1.1*min(f) 1.1*max(f)]); end 实验截图： （2）求解原连续信号波形和抽样信号所对应的幅度谱。实验程序： dt=0.1;t=-4:dt:4;
N=length(t);f=sinc(t);Tm=1;fm=1/Tm; wm=2*pi*fm;k=1:N; w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt; subplot(4,1,1),plot(w1/(2*pi),abs(F1));grid axis([0 max(4*fm) 1.1*min(F1) 1.1*max(F1)]); for i=1:3; if i<= 2 c=0 ,else c=0.2,end fs=(4-i+c)*fm; Ts=1/fs; n=-4:Ts:4; f=sinc(n); N=length(n); wm=2*pi*fs; k=1:N; w=k*wm/N; F=f*exp(-j*n'*w)*Ts; subplot(4,1,5-i),plot(w/(2*pi),abs(F),'k');grid axis([0 max(4*fm) 1.1*min(F) 1.1*max(F)]); end 实验截图：
大学本科语音信号处理实验讲义8学时
语音信号处理实验讲义 时间：2011-12
目录 实验一语音信号生成模型分析 (3) 实验二语音信号时域特征分析 (7) 实验三语音信号频域特征分析 (12) 实验四语音信号的同态处理和倒谱分析 (16)
实验一 语音信号生成模型分析 一、实验目的 1、了解语音信号的生成机理，了解由声门产生的激励函数、由声道产生的调制函数和由嘴唇产生的辐射函数。 2、编程实现声门激励波函数波形及频谱，与理论值进行比较。 3、编程实现已知语音信号的语谱图，区分浊音信号和清音信号在语谱图上的差别。 二、实验原理 语音生成系统包含三部分：由声门产生的激励函数()G z 、由声道产生的调制函数()V z 和由嘴唇产生的辐射函数()R z 。语音生成系统的传递函数由这三个函数级联而成，即 ()()()()H z G z V z R z = 1、激励模型 发浊音时，由于声门不断开启和关闭，产生间隙的脉冲。经仪器测试它类似于斜三角波的脉冲。也就是说，这时的激励波是一个以基音周期为周期的斜三角脉冲串。单个斜三角波的频谱表现出一个低通滤波器的特性。可以把它表示成z 变换的全极点形式 12 1()(1) cT G z e z --= -? 这里c 是一个常数，T 是脉冲持续时间。周期的三角波脉冲还得跟单位脉冲串的z 变换相乘： 112 1 ()()()1(1)v cT A U z E z G z z e z ---=?= ?--? 这就是整个激励模型，v A 是一个幅值因子。 2、声道模型 当声波通过声道时，受到声腔共振的影响，在某些频率附近形成谐振。反映在信号频谱图上，在谐振频率处其谱线包络产生峰值，把它称为共振峰。 一个二阶谐振器的传输函数可以写成 12 ()1i i i i A V z B z C z --= -- 实践表明，用前3个共振峰代表一个元音足够了。对于较复杂的辅音或鼻音共振峰要到5个以上。多个()i V z 叠加可以得到声道的共振峰模型 12 1 11 ()()11R r r M M i r i N k i i i i k k b z A V z V z B z C z a z -=---======---∑∑∑ ∑ 3、辐射模型 从声道模型输出的是速度波，而语音信号是声压波。二者倒比称为辐射阻抗，它表征了
信号与系统——信号的采样与恢复实验
实验六 信号与系统实验 1.信号的采样与恢复实验 1.1实验目的 (1)熟悉信号的采样与恢复的过程 (2)学习和掌握采样定理 (3)了解采样频率对信号恢复的影响 1.2实验原理及内容 (1)采样定理 采样定理论述了在一定条件下，一个连续时间信号完全可以用该信号等时间间隔上瞬时值表示，这些值包含该信号全部信息，利用这些值可以恢复原信号。采样定理是连续时间信号与离散时间信号的桥梁。 采样定理：对于一个具有有限频谱且最高频率为max w 的连续信号进行采样，当采样频率s w >=2max w 时，采样函数能够无失真地恢复出原信号。 (2)采样信号的频谱 连续周期信号经过周期矩形脉冲抽样后，抽样信号的频谱为 )]([)2 ( )(s n s s nw w j F nw Sa T A jw F -= ∑ +∞ -∞ =τ τ 它包含了原信号频谱以及重复周期为s w 的原信号频谱的搬移，且幅度按 )2 (ττ s nw Sa T A 规律变化。所以抽样信号的频谱便是原信号频谱的周期性拓延。 （3）采样信号的恢复 将采样信号恢复成原信号，可以是用低通滤波器。低通滤波器的截止频率c f 应当满足 max max f f f f x c -≤≤。实验中采用的低通滤波器的截止频率固定为 Hz RC f 8021≈=π （4）单元构成 本实验电路由脉冲采样电路和滤波器两部分构成，滤波器部分不再赘述，其中采样保持部分电路由一片CD4052完成。此电路有两个输入端，其中IN1端输入被采样信号，Pu 端输入采样脉冲。 1.3实验步骤 本实验在脉冲与恢复单元完成。 (1)信号的采样 1)使波形发生器第一路输出幅值3V 、频率10Hz 的三角波信号；第二路输出幅值5V 、频率100Hz 、占空比50%的脉冲信号，将第一路信号接入IN1端；作为输入信号，第二路信号接入Pu 端，作为采样脉冲。 2)用示波器分别测量IN1端和OUT1端，观察采样前后波形的差异。 3)增加采样脉冲的频率为200、500、800等值。观察OUT1端波形的变化。解释现象产生的原因。
语音信号处理实验报告实验二
通信工程学院12级1班 罗恒 2012101032 实验二 基于MATLAB 的语音信号频域特征分析 一、 实验要求 要求根据已有语音信号，自己设计程序，给出其倒谱、语谱图的分析结果，并根据频域分析方法检测所分析语音信号的基音周期或共振峰。 二、 实验目的 信号的傅立叶表示在信号的分析与处理中起着重要的作用。因为对于线性系统来说，可以很方便地确定其对正弦或复指数和的响应，所以傅立叶分析方法能完善地解决许多信号分析和处理问题。另外，傅立叶表示使信号的某些特性变得更明显，因此，它能更深入地说明信号的各项红物理现象。 由于语音信号是随着时间变化的，通常认为，语音是一个受准周期脉冲或随机噪声源激励的线性系统的输出。输出频谱是声道系统频率响应与激励源频谱的乘积。声道系统的频率响应及激励源都是随时间变化的，因此一般标准的傅立叶表示虽然适用于周期及平稳随机信号的表示，但不能直接用于语音信号。由于语音信号可以认为在短时间内，近似不变，因而可以采用短时分析法。 三、 实验设备 1.PC 机； 2.MATLAB 软件环境； 四、 实验内容 1.上机前用Matlab 语言完成程序编写工作。 2.程序应具有加窗（分帧）、绘制曲线等功能。 3.上机实验时先调试程序，通过后进行信号处理。 4.对录入的语音数据进行处理，并显示运行结果。 5.依次给出其倒谱、语谱图的分析结果。 6. 根据频域分析方法检测所分析语音信号的基音周期或共振峰。 五、 实验原理及方法 1、短时傅立叶变换 由于语音信号是短时平稳的随机信号，某一语音信号帧的短时傅立叶变换的定义为： 其中w(n -m)是实窗口函数序列，n 表示某一语音信号帧。令n -m=k'，则得到 ()()()jw jwm n m X e x m w n m e ∞-=-∞= -∑
数字信号处理实验报告实验三(DOC)
物理与电子信息工程学院 实验报告 实验课程名称：数字信号处理 实验名称：用FFT对信号作频谱分析班级：1012341 姓名：严娅 学号：101234153 成绩：_______ 实验时间：2012年12月6日
一、实验目的 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便正确应用FFT。 二、实验原理 用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D和分析误差。频谱分辨率直接和FFT的变换区间N有关，因为FFT能够实现的频率分辨率是N/ 2π，因此要求D 2π。可以根据此式选择FFT的变换区间N。误差主要来自于/ N≤ 用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时离散谱的包络才能逼近于连续谱，因此N要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。 对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。 三、实验步骤及内容 （1）对以下序列进行谱分析。
?????≤≤-≤≤-=?? ???≤≤-≤≤+==其它n n n n n n x 其它n n n n n n x n R n x ,074, 330,4)(,074, 830, 1)() ()(3241 这些都是时域离散非周期信号，选择FFT 的变换区间N 为8和16 两种情况进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。 （2）对以下周期序列进行谱分析。 4()cos 4x n n π= 5()cos(/4)cos(/8)x n n n ππ=+ 这些是时域离散周期信号，选择FFT 的变换区间N 为8和16 两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。 （3）对模拟周期信号进行谱分析 6()cos8cos16cos20x t t t t πππ=++ 这是时域连续周期信号，选择采样频率Hz F s 64=，变换区间N=16,32,64 三种情况进行谱分析。分别打印其幅频特性，并进行分析和讨论。 四、实验程序清单 %第10章实验3程序exp3.m % 用FFT 对信号作频谱分析