目录 摘要 (2) 第一章绪论 (3) 1.1 语音课设的意义 (3) 1.2 语音课设的目的与要求 (3) 1.3 语音课设的基本步骤 (3) 第二章设计方案论证 (5) 2.1 设计理论依据 (5) 2.1.1 采样定理 (5) 2.1.2 采样频率 (5) 2.1.3 采样位数与采样频率 (5) 2.2 语音信号的分析及处理方法 (6) 2.2.1 语音的录入与打开 (6) 2.2.2 时域信号的FFT分析 (6) 2.2.3 数字滤波器设计原理 (7) 2.2.4 数字滤波器的设计步骤 (7) 2.2.5 IIR滤波器与FIR滤波器的性能比较 (7) 第三章图形用户界面设计 (8) 3.1 图形用户界面概念 (8) 3.2 图形用户界面设计 (8) 3.3 图形用户界面模块调试 (9) 3.3.1 语音信号的读入与打开 (9) 3.3.2 语音信号的定点分析 (9) 3.3.3 N阶高通滤波器 (11) 3.3.4 N阶低通滤波器 (12) 3.3.5 2N阶带通滤波器 (13) 3.3.6 2N阶带阻滤波器 (14) 3.4 图形用户界面制作 (15) 第四章总结 (18) 附录 (19) 参考文献 (24)
摘要 数字信号处理是将信号以数字方式表示并处理的理论和技术。数字信号处理与模拟信号处理是信号处理的子集。 数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。因此在进行数字信号处理之前需要将信号从模拟域转换到数字域，这通常通过模数转换器实现。而数字信号处理的输出经常也要变换到模拟域，这是通过数模转换器实现的。 数字信号处理的算法需要利用计算机或专用处理设备如数字信号处理器（DSP）和专用集成电路（ASIC）等。数字信号处理技术及设备具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等突出优点，这些都是模拟信号处理技术与设备所无法比拟的。 数字信号处理的核心算法是离散傅立叶变换(DFT)，是DFT使信号在数字域和频域都实现了离散化，从而可以用通用计算机处理离散信号。而使数字信号处理从理论走向实用的是快速傅立叶变换(FFT)，FFT的出现大大减少了DFT的运算量，使实时的数字信号处理成为可能、极大促进了该学科的发展。 MATLAB是矩阵实验室（Matrix Laboratory）的简称，和Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。
应用Matlab对含噪声的语音信号进行频谱分析及滤波 一、实验内容 录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；在语音信号中增加正弦噪声信号(自己设置几个频率的正弦信号)，对加入噪声信号后的语音信号进行频谱分析；给定滤波器的性能指标，采用窗函数法和双线性变换设计数字滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比试听，分析信号的变化。 二、实现步骤 1．语音信号的采集 利用Windows下的录音机，录制一段自己的话音，时间在1 s内。然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，（可用默认的采样频率或者自己设定采样频率）。 2．语音信号的频谱分析 要求首先画出语音信号的时域波形；然后对语音号进行快速傅里叶变换，得到信号的频谱特性。 在采集得到的语音信号中加入正弦噪声信号，然后对加入噪声信号后的语音号进行快速傅里叶变换，得到信号的频谱特性。并利用sound试听前后语音信号的不同。
分别设计IIR和FIR滤波器，对加入噪声信号的语音信号进行去噪，画出并分析去噪后的语音信号的频谱，并进行前后试听对比。 3.数字滤波器设计 给出数字低通滤波器性能指标:如，通带截止频率fp＝10000 Hz，阻带截止频率fs＝12000 Hz（可根据自己所加入噪声信号的频率进行阻带截止频率设置），阻带最小衰减Rs＝50 dB，通带最大衰减Rp＝3 dB（也可自己设置），采样频率根据自己语音信号采样频率设定。
报告内容 一、实验原理 含噪声语音信号通过低通滤波器，高频的噪声信号会被过滤掉，得到清晰的无噪声语音信号。 二、实验内容 录制一段个人自己的语音信号，并对录制的信号进行采样；画出采样后语音信号的时域波形和频谱图；在语音信号中增加正弦噪声信号(自己设置几个频率的正弦信号)，对加入噪声信号后的语音信号进行频谱分析；给定滤波器的性能指标，采用窗函数法和双线性变换设计数字滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采集的信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比试听，分析信号的变化。给出数字低通滤波器性能指标:如，通带截止频率fp＝10000 Hz，阻带截止频率fs＝12000 Hz （可根据自己所加入噪声信号的频率进行阻带截止频率设置），阻带最小衰减Rs＝50 dB，通带最大衰减Rp＝3 dB（也可自己设置），采样频率根据自己语音信号采样频率设定。 三、实验程序 1、原始信号采集和分析 clc;clear;close all; fs=10000; %语音信号采样频率为10000 x1=wavread('C:\Users\acer\Desktop\voice.wav'); %读取语音信号的数据，赋给x1 sound(x1,40000); %播放语音信号 y1=fft(x1,10240); %对信号做1024点FFT变换 f=fs*(0:1999)/1024; figure(1); plot(x1) %做原始语音信号的时域图形 title('原始语音信号'); xlabel('time n'); ylabel('fuzhi n'); figure(2); plot(f,abs(y1(1:2000))); %做原始语音信号的频谱图形 title('原始语音信号频谱') xlabel('Hz'); ylabel('fuzhi');
1-1画出下列各信号的波形【式中)()(t t t r ε=】为斜升函数。 （2）∞<<-∞=-t e t f t ,)( （3）)()sin()(t t t f επ= （4）)(sin )(t t f ε= （5）)(sin )(t r t f = （7）)(2)(k t f k ε= （10）)(])1(1[)(k k f k ε-+= 解：各信号波形为 （2）∞<<-∞=-t e t f t ,)( （3）)()sin()(t t t f επ= （4）)(sin )(t t f ε= （5）)(sin )(t r t f = （7）)(2)(k t f k ε= （10）)(])1(1[)(k k f k ε-+= 1-2 画出下列各信号的波形[式中)()(t t t r ε=为斜升函数]。 （1）)2()1(3)1(2)(-+--+=t t t t f εεε （2）)2()1(2)()(-+--=t r t r t r t f （5）)2()2()(t t r t f -=ε （8）)]5()([)(--=k k k k f εε （11） )]7()()[6sin()(--=k k k k f εεπ （12）)]()3([2)(k k k f k ---=εε 解：各信号波形为 （1）)2()1(3)1(2)(-+--+=t t t t f εεε （2）)2()1(2)()(-+--=t r t r t r t f
（5） )2()2()(t t r t f -=ε （8）)]5()([)(--=k k k k f εε （11） )]7()()[6sin()(--=k k k k f εεπ （12）)]()3([2)(k k k f k ---=εε 1-3 写出图1-3所示各波形的表达式。 1-4 写出图1-4所示各序列的闭合形式表达式。 1-5 判别下列各序列是否为周期性的。如果是，确定其周期。 （2）)6 3cos()443cos()(2ππππ+++=k k k f （5）)sin(2cos 3)(5t t t f π+= 解： 1-6 已知信号)(t f 的波形如图1-5所示，画出下列各函数的波形。 （1）)()1(t t f ε- （2）)1()1(--t t f ε （5） )21(t f - （6）)25.0(-t f （7）dt t df ) ( （8）dx x f t ?∞-)( 解：各信号波形为 （1）)()1(t t f ε- （2） )1()1(--t t f ε （5）)21(t f -
信号与系统期末考试试题 一、选择题（共10题，每题3分 ，共30分，每题给出四个答案，其中只有一个正确的） 1、 卷积f 1(k+5)*f 2(k-3) 等于 。 （A ）f 1(k)*f 2(k) （B ）f 1(k)*f 2(k-8)（C ）f 1(k)*f 2(k+8)（D ）f 1(k+3)*f 2(k-3) 2、 积分 dt t t ? ∞ ∞ --+)21()2(δ等于 。 （A ）1.25（B ）2.5（C ）3（D ）5 3、 序列f(k)=-u(-k)的z 变换等于 。 （A ） 1-z z （B ）-1-z z （C ）11-z （D ）1 1--z 4、 若y(t)=f(t)*h(t),则f(2t)*h(2t)等于 。 （A ） )2(41t y （B ）)2(21t y （C ）)4(41t y （D ）)4(2 1 t y 5、 已知一个线性时不变系统的阶跃相应g(t)=2e -2t u(t)+)(t δ,当输入f(t)=3e —t u(t)时，系 统的零状态响应y f (t)等于 （A ）(-9e -t +12e -2t )u(t) （B ）(3-9e -t +12e -2t )u(t) （C ）)(t δ+(-6e -t +8e -2t )u(t) （D ）3)(t δ +(-9e -t +12e -2t )u(t) 6、 连续周期信号的频谱具有 （A ） 连续性、周期性 （B ）连续性、收敛性 （C ）离散性、周期性 （D ）离散性、收敛性 7、 周期序列2)455.1(0 +k COS π的 周期N 等于 （A ） 1（B ）2（C ）3（D ）4 8、序列和 ()∑∞ -∞ =-k k 1δ等于 （A ）1 (B) ∞ (C) ()1-k u (D) ()1-k ku 9、单边拉普拉斯变换()s e s s s F 22 12-+= 的愿函数等于 ()()t tu A ()()2-t tu B ()()()t u t C 2- ()()()22--t u t D 10、信号()()23-=-t u te t f t 的单边拉氏变换()s F 等于 ()A ()()()232372+++-s e s s ()() 2 23+-s e B s
2020/3/27 2009-2010 学年第二学期 通信工程专业《数字信号处理》（课程）参考答案及评分标准 一、 选择题 （每空 1 分，共 20 分） 1．序列 x( n) cos n sin n 的周期为（ A ）。 4 6 A ． 24 B ． 2 C ． 8 D ．不是周期的 2．有一连续信号 x a (t) cos(40 t) ，用采样间隔 T 0.02s 对 x a (t) 进行采样，则采样所得的时域离散信 号 x(n) 的周期为（ C ） A ． 20 B ． 2 C ． 5 D ．不是周期的 3．某线性移不变离散系统的单位抽样响应为h(n) 3n u( n) ，该系统是（ B ）系统。 A ．因果稳定 B ．因果不稳定 C ．非因果稳定 D ．非因果不稳定 4．已知采样信号的采样频率为 f s ，采样周期为 T s ，采样信号的频谱是原模拟信号频谱的周期函数，周 期为（ A ），折叠频率为（ C ）。 A ． f s B ． T s C ． f s / 2 D ． f s / 4 5．以下关于序列的傅里叶变换 X ( e j ) 说法中，正确的是（ B ）。 A ． X ( e B ． X ( e C ． X (e D ． X (e j j j j ) 关于 是周期的，周期为 ) 关于 是周期的，周期为 2 ) 关于 是非周期的 ) 关于 可能是周期的也可能是非周期的 6．已知序列 x(n) 2 (n 1) (n)(n 1) ，则 j X (e ) 的值为（ ）。 C
2020/3/27 A ． 0 B ． 1 C ． 2 D ． 3 N 1 7．某序列的 DFT 表达式为 X (k ) x(n)W M nk ，由此可看出，该序列的时域长度是（ A ），变换后数字域 n 0 上相邻两个频率样点之间的间隔（ C ）。 A ． N B ． M C ．2 /M D ． 2 / N 8．设实连续信号 x(t) 中含有频率 40 Hz 的余弦信号，现用 f s 120 Hz 的采样频率对其进行采样，并利 用 N 1024 点 DFT 分析信号的频谱，得到频谱的谱峰出现在第（ B ）条谱线附近。 A ． 40 B ． 341 C ． 682 D ．1024 9．已知 x( n) 1，2,3,4 ，则 x ( ) R 6 ( ) （ ）， x ( n 1) R 6 (n) （ ） n 6 n 6 A C A ． 1,0,0,4,3,2 B ． 2,1,0,0,4,3 C ． 2,3,4,0,0,1 D ． 0,1,2,3,4，0 10．下列表示错误的是（ B ）。 A ． W N nk W N ( N k) n B ． (W N nk ) * W N nk C ． W N nk W N (N n) k D ． W N N /2 1 11．对于 N 2L 点的按频率抽取基 2FFT 算法，共需要（ A ）级蝶形运算，每级需要（ C ）个蝶形运算。 A ． L B ． L N 2 C ． N D ． N L 2 12．在 IIR 滤波器中，（ C ）型结构可以灵活控制零极点特性。 A ．直接Ⅰ B ．直接Ⅱ C ．级联 D ．并联 13．考虑到频率混叠现象，用冲激响应不变法设计 IIR 数字滤波器不适合于（ B ）。 A ．低通滤波器 B ．高通、带阻滤波器 C ．带通滤波器 D ．任何滤波器
系（院）物理与电子工程学院专业电子信息工程题目语音信号的处理与分析 学生姓名 指导教师 班级 学号 完成日期：2013 年5 月 目录 1 绪论.............................................................................................................. 错误!未定义书签。 1.1课题背景及意义................................................................................. 错误!未定义书签。 1.2国内外研究现状................................................................................. 错误!未定义书签。 1.3本课题的研究内容和方法................................................................. 错误!未定义书签。 1.3.1 研究内容................................................................................ 错误!未定义书签。 1.3.2 开发环境................................................................................ 错误!未定义书签。 2 语音信号处理的总体方案............................................................................ 错误!未定义书签。 2.1 系统基本概述.................................................................................... 错误!未定义书签。 2.2 系统基本要求与目的........................................................................ 错误!未定义书签。 2.3 系统框架及实现................................................................................ 错误!未定义书签。 2.3.1 语音信号的采样.................................................................... 错误!未定义书签。 2.3.2 语音信号的频谱分析............................................................ 错误!未定义书签。 2.3.3 音乐信号的抽取.................................................................... 错误!未定义书签。 2.3.4 音乐信号的AM调制.............................................................. 错误!未定义书签。 2.3.5 AM调制音乐信号的同步解调............................................... 错误!未定义书签。 2.4系统设计流程图................................................................................. 错误!未定义书签。 3 语音信号处理基本知识................................................................................ 错误!未定义书签。 3.1语音的录入与打开............................................................................. 错误!未定义书签。 3.2采样位数和采样频率......................................................................... 错误!未定义书签。 3.3时域信号的FFT分析......................................................................... 错误!未定义书签。 3.4切比雪夫滤波器................................................................................. 错误!未定义书签。 3.5数字滤波器设计原理......................................................................... 错误!未定义书签。 4 语音信号实例处理设计................................................................................ 错误!未定义书签。 4.1语音信号的采集................................................................................. 错误!未定义书签。
数字信号处理实验二：语音信号分析与处理 学号 姓名 注：1）此次实验作为《数字信号处理》课程实验成绩的重要依据，请同学们认真、独立完成，不得抄袭。 2）请在授课教师规定的时间内完成； 3）完成作业后，请以word 格式保存，文件名为：学号+姓名 4）请通读全文，依据第2及第3 两部分内容，认真填写第4部分所需的实验数据，并给出程序内容。 1. 实验目的 (1) 学会MATLAB 的使用，掌握MATLAB 的程序设计方法 (2) 掌握在windows 环境下语音信号采集的方法 (3) 掌握MATLAB 设计FIR 和IIR 滤波器的方法及应用 (4) 学会用MATLAB 对语音信号的分析与处理方法 2. 实验内容 录制一段自己的语音信号，对录制的语音信号进行采样，画出采样后语音信号的时域波形和频谱图，确定语音信号的频带范围；使用MATLAB 产生白噪声信号模拟语音信号在处理过程中的加性噪声并与语音信号进行叠加，画出受污染语音信号的时域波形和频谱图；采用双线性法设计出IIR 滤波器和窗函数法设计出FIR 滤波器，画出滤波器的频响特性图；用自己设计的这两种滤波器分别对受污染的语音信号进行滤波，画出滤波后语音信号的时域波形和频谱图；对滤波前后的语音信号进行时域波形和频谱图的对比，分析信号的变化；回放语音信号，感觉与原始语音的不同。 3. 实验步骤 1）语音信号的采集与回放 利用windows 下的录音机或其他软件录制一段自己的语音（规定：语音内容为自己的名字，以wav 格式保存，如wql.wav ），时间控制在2秒之内，利用MATLAB 提供的函数wavread 对语音信号进行采样，提供sound 函数对语音信号进行回放。 [y,fs,nbits]=wavread(file), 采样值放在向量y 中，fs 表示采样频率nbits 表示采样位数。Wavread 的更多用法请使用help 命令自行查询。 2）语音信号的频谱分析 利用fft 函数对信号进行频谱分析 3）受白噪声干扰的语音信号的产生与频谱分析 ①白噪声的产生： N1=sqrt （方差值）×randn(语音数据长度，2)（其中2表示2列，是由于双声道的原因） 然后根据语音信号的频谱范围让白噪声信号通过一个带通滤波器得到一个带限的白噪声信号 N2； 带通滤波器的冲激响应为： h B （n ）= ))((sin ))((sin 1122απ ωπωαπωπω---n c n c c c c c
一．填空题(本大题共10空，每空2分，共20分。) 1．()*(2)k k εδ-= . 2． sin()()2 t d π τδττ-∞ + =? . 3. 已知信号的拉普拉斯变换为 1 s a -，若实数a ，则信号的傅里叶变换不存在. 4. ()()()t h t f t y *=，则()=t y 2 . 5. 根据Parseval 能量守恒定律，计算?∞ ∞-=dt t t 2 )sin ( . 6. 若)(t f 最高角频率为m ω，则对 )2()4()(t f t f t y =取样，其频谱不混迭的最大间隔是 . 7. 某因果线性非时变（LTI ）系统，输入)()(t t f ε=时，输出为： )1()()(t t e t y t --+=-εε；则) 2()1()(---=t t t f εε时，输出)(t y f ＝ . 8. 已知某因果连续LTI 系统)(s H 全部极点均位于s 左半平面，则 ∞→t t h )(的值为 . 9. 若)()(ωj F t f ?，已知)2cos()(ωω=j F ，试求信号)(t f 为 . 10.已知某离散信号的单边z 变换为) 3(,)3)(2(2)(2>+-+=z z z z z z F ，试求其反变换 )(k f = . 二．选择题(本大题共5小题，每题4分，共20分。) 1.下列信号的分类方法不正确的是 ： A 、数字信号和离散信号 B 、确定信号和随机信号 C 、周期信号和非周期信号 D 、因果信号与反因果信号 2. )]2()()[2()]()2([2)(1--++-+=t t t t t t f εεεε，则)] 1()2 1()[21()(--+-=t t t f t f εε
数字信号处理模拟试题一 一、单项选择题(本大题共10小题，每小题2分，共20分) 在每小题列出的四个备选项中只有一个是符合题目要求的，请将其代码填写在题后的括号内。错选、多选或未选均无分。 1.在对连续信号均匀采样时，要从离散采样值不失真恢复原信号，则采样角频率Ωs与信号最高截止频率Ωc应满足关系（A ） A.Ωs>2Ωc B.Ωs>Ωc C.Ωs<Ωc D.Ωs<2Ωc 2.下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？（D） A.y(n)=y(n-1)x(n) B.y(n)=x(n)/x(n+1) C.y(n)=x(n)+1 D.y(n)=x(n)-x(n-1) 3.已知某序列Z变换的收敛域为5>|z|>3，则该序列为（D ） A.有限长序列 B.右边序列 C.左边序列 D.双边序列 4.实偶序列傅里叶变换是（A ） A.实偶序列 B.实奇序列 C.虚偶序列 D.虚奇序列 5.已知x(n)=δ(n)，其N点的DFT［x(n)］=X(k)，则X(N-1)=（B） A.N-1 B.1 C.0 D.-N+1 6.设两有限长序列的长度分别是M与N，欲通过计算两者的圆周卷积来得到两者的线性卷积，则圆周卷积的点数至少应取（B ） A.M+N B.M+N-1 C.M+N+1 D.2(M+N) 7.下面说法中正确的是（C） A.连续非周期信号的频谱为周期连续函数 B.连续周期信号的频谱为周期连续函数 C.离散非周期信号的频谱为周期连续函数 D.离散周期信号的频谱为周期连续函数 8.下列各种滤波器的结构中哪种不是IIR滤波器的基本结构？（C ） A.直接型 B.级联型 C.频率抽样型 D.并联型 9.下列关于FIR滤波器的说法中正确的是（C） A.FIR滤波器容易设计成线性相位特性
6.4 根据下列象函数及所标注的收敛域，求其所对应的原序列。 （1）1)(=z F ，全z 平面 （2）∞<=z z z F ,)(3 （3）0,)(1>=-z z z F （4）∞<<-+=-z z z z F 0,12)(2 （5）a z az z F >-= -,11 )(1 （6）a z az z F <-=-,11 )(1 6.5 已知1)(?k δ，a z z k a k -? )(ε，2 )1()(-?z z k k ε，试利用z 变换的性质求下列序
列的z 变换并注明收敛域。 （1）)(])1(1[2 1k k ε-+ （3）)()1(k k k ε- （5）)1()1(--k k k ε （7）)]4()([--k k k εε （9）)()2 cos( )2 1(k k k επ
6.8 若因果序列的z 变换)(z F 如下，能否应用终值定理？如果能，求出)(lim k f k ∞ →。 （1）)3 1 )(21(1 )(2+-+=z z z z F （3）)2)(1()(2--=z z z z F 6.10 求下列象函数的双边逆z 变换。 （1）31 ,)31)(21(1)(2<--+= z z z z z F （2）21 ,)3 1)(21()(2>--= z z z z z F （3）2 1,) 1()2 1 ()(23 < --= z z z z z F
（4） 2 1 3 1 , )1 ( ) 2 1 ( ) ( 2 3 < < - - =z z z z z F
6.11 求下列象函数的逆z 变换。 （1）1,1 1 )(2>+= z z z F （2）1,) 1)(1()(2 2>+--+=z z z z z z z F （5）1,) 1)(1()(2>--= z z z z z F （6）a z a z az z z F >-+=,) ()(3 2
《信号与系统》期末试卷与答案
第 2 页 共 14 页 《信号与系统》期末试卷A 卷 班级： 学号：__________ 姓名： ________ _ 成绩：_____________ 一． 选择题（共10题，20分） 1、n j n j e e n x )3 4( )3 2(][ππ +=，该序列是 。 A.非周期序列 B.周期3=N C.周期 8 /3=N D. 周期24=N 2、一连续时间系统y(t)= x(sint)，该系统是 。 A.因果时不变 B.因果时变 C.非因果时不变 D. 非因果时变 3、一连续时间LTI 系统的单位冲激响应 ) 2()(4-=-t u e t h t ，该系统是 。 A.因果稳定 B.因果不稳定 C.非因果稳定 D. 非因果不稳定
第 3 页 共 14 页 4、若周期信号x[n]是实信号和奇信号，则其傅立叶级数系数a k 是 。 A.实且偶 B.实且为奇 C.纯虚且偶 D. 纯虚且奇 5、一信号x(t)的傅立叶变换? ? ?><=2 ||02 ||1)(ωωω，，j X ，则x(t)为 。 A. t t 22sin B. t t π2sin C. t t 44sin D. t t π4sin 6、一周期信号∑∞-∞ =-=n n t t x )5()(δ，其傅立叶变换)(ωj X 为 。 A. ∑∞ -∞ =-k k )52(5 2πωδπ B. ∑∞ -∞=-k k )5 2(25 πωδπ C. ∑∞ -∞ =-k k )10(10πωδπ D. ∑∞ -∞ =- k k ) 10 (101πωδπ 7、一实信号x[n]的傅立叶变换为)(ω j e X ，则x[n]奇部的傅立叶变换为 。
语音信号处理试验 实验一：语音信号时域分析 实验目的： (1)录制两段语音信号，内容是“语音信号处理”，分男女声。 (2)对语音信号进行采样，观察采样后语音信号的时域波形。 实验步骤： 1、使用window自带录音工具录制声音片段 使用windows自带录音机录制语音文件，进行数字信号的采集。启动录音机。录制一段录音，录音停止后，文件存储器的后缀默认为.Wav。将录制好文件保存，记录保存路径。男生女生各录一段保存为test1.wav和test2.wav。 图1基于PC机语音信号采集过程。 2、读取语音信号 在MATLAB软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数。通过使用wavread函数，理解采样、采样频率、采样位数等概念！ Wavread函数调用格式： y=wavread（file），读取file所规定的wav文件，返回采样值放在向量y中。
[y，fs，nbits]=wavread（file），采样值放在向量y中，fs表示采样频率（hz），nbits表示采样位数。 y=wavread（file，N），读取前N点的采样值放在向量y中。 y=wavread（file，[N1，N2]）,读取从N1到N2点的采样值放在向量y中。 3、编程获取语音信号的抽样频率和采样位数。 语音信号为test1.wav和test2.wav，内容为“语音信号处理”，两端语音保存到工作空间work文件夹下。在M文件中分别输入以下程序，可以分两次输入便于观察。 [y1,fs1,nbits1]=wavread('test1.wav') [y2,fs2,nbits2]=wavread('test2.wav') 结果如下图所示 根据结果可知：两端语音信号的采样频率为44100HZ，采样位数为16。 4、语音信号的时域分析 语音信号的时域分析就是分析和提取语音信号的时域参数。进行语音分析时，最先接触到并且夜市最直观的是它的时域波形。语音信
摘要 (2) 1 设计目的与要求 (3) 2 设计步骤 (4) 3 设计原理及内容 (5) 3.1 理论依据 (5) 3.2 信号采集 (6) 3.3 构造受干扰信号并对其FFT频谱分析 (8) 3.4 数字滤波器设计 (9) 3.5 信号处理 (10) 总结 (12) 致谢 (13) 参考文献 (14)
用MATLAB对语音信号进行分析与处理，采集语音信号后，在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。 数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测和参数估计等处理。IIR数字滤波器最大的优点是给定一组指标时，它的阶数要比相同组的FIR滤波器的低的多。信号处理中和频谱分析最为密切的理论基础是傅立叶变换（FT）。离散傅立叶变换（DFT）和数字滤波是数字信号处理的最基本内容。 关键词：MATLAB;语音信号；加入噪声；滤波器；滤波
1. 设计目的与要求 （1）待处理的语音信号是一个在20Hz~20kHz频段的低频信号 （2）要求MATLAB对语音信号进行分析和处理，采集语音信号后，在MATLAB平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器进行滤除噪声，恢复原信号。
2. 设计步骤 （1）选择一个语音信号或者自己录制一段语音文件作为分析对象； （2）对语音信号进行采样，并对语音信号进行FFT频谱分析，画出信号的时域波形图和频谱图； （3）利用MATLAB自带的随机函数产生噪声加入到语音信号中，对语音信号进行回放，对其进行FFT频谱分析； （4）设计合适滤波器，对带有噪声的语音信号进行滤波，画出滤波前后的时域波形图和频谱图，比较加噪前后的语音信号，分析发生的变化； （5）对语音信号进行回放，感觉声音变化。
信号与线性系统课后答案 第一章 信号与系统（一） 1-1画出下列各信号的波形【式中)() (t t t r ε=】为斜升函数。 （2）∞<<-∞=-t e t f t ,)( （3）)()sin()(t t t f επ= （4）)(sin )(t t f ε= （5）)(sin )(t r t f = （7）)(2)(k t f k ε= （10）)(])1(1[)(k k f k ε-+= 解：各信号波形为 （2）∞<<-∞=-t e t f t ,)( （3）)()sin()(t t t f επ=
（4）) t fε= (sin )(t （5）) t f= r )(t (sin
（7）)( t f kε 2 )(k = （10）)(])1 k (k f kε ( ) 1[ = - +
1-2 画出下列各信号的波形[式中)()(t t t r ε=为斜升函数]。 （1）)2()1(3)1(2)(-+--+=t t t t f εεε （2）)2()1(2)()(-+--=t r t r t r t f （5） )2()2()(t t r t f -=ε （8）)]5()([)(--=k k k k f εε （11）)]7()()[6 sin()(--=k k k k f εεπ （12） )]()3([2)(k k k f k ---=εε 解：各信号波形为 （1） )2()1(3)1(2)(-+--+=t t t t f εεε
（2） )2()1(2)()(-+--=t r t r t r t f （5） )2()2()(t t r t f -=ε
重庆大学信号与线性系统期末考试试题 一、填空题：（30分，每小题3分） 1. =-? ∞ ∞ -dt t t )()5cos 2(δ 。 2. ()dt t e t 12-?+∞ ∞ --δ= 。 3. 已知 f (t )的傅里叶变换为F (j ω), 则f (2t -3)的傅里叶变换为 。 4. 已知 6 51 )(2 +++= s s s s F ，则=+)0(f ; =∞)(f 。 5. 已知 ω ωπδεj t FT 1 )()]([+=，则=)]([t t FT ε 。 6. 已知周期信号 )4sin()2cos()(t t t f +=，其基波频率为 rad/s ； 周期为 s 。 7. 已知 )5(2)2(3)(-+-=n n k f δδ，其Z 变换 =)(Z F ；收敛域为 。 8. 已知连续系统函数1342 3)(2 3+--+= s s s s s H ，试判断系统的稳定性： 。 9．已知离散系统函数1.07.02 )(2+-+=z z z z H ，试判断系统的稳定性： 。 10．如图所示是离散系统的Z 域框图，该系统的系统函数H(z)= 。 二．(15分)如下方程和非零起始条件表示的连续时间因果LTI 系统，
?????==+=++-- 5 )0(',2)0() (52)(4522y y t f dt df t y dt dy dt y d 已知输入 )()(2t e t f t ε-=时，试用拉普拉斯变换的方法求系统的零状态响应 )(t y zs 和零输入响应)(t y zi ，0≥t 以及系统的全响应),(t y 0≥t 。 三．（14分） ① 已知2 36 62)(22++++=s s s s s F ,2]Re[->s ,试求其拉氏逆变换f (t )； ② 已知) 2(2 35)(2>+-=z z z z z X ，试求其逆Z 变换)(n x 。 四 （10分）计算下列卷积： 1. }1,0,6,4,3{}4,1,2,1{)()(21--*=*k f k f ； 2． )(3)(23t e t e t t εε--* 。
信号与线性系统分析 离散信号部分 1. 用MATLAB画出正弦离散序列的时域波形。 N=100; n=-N:N; w0=0.2; f1=cos((pi*n*w0)/8); f2=cos(2*n*w0); subplot(211); stem(n,f1); grid on; title('f1=cos((pi*n*w0)/8)'); xlabel('n'); ylabel('f1(n)'); subplot(212); stem(n,f2); grid on; title('f2=cos(2*n)');
xlabel('n'); ylabel('f2(n)'); 信号运算部分 2.已知信号 ，画出 的波形； t=-20:0.01:20; f1=0.25*(t+1).*(t>-4&t<0)+1.*(t>0&t<2)+0.*(t>=2&t<=-4&t==0); subplot(211); plot(t,f1); grid on; title('f1=(t+1)/4.*(t>-4&t<0)+1.*(t>0&t<2)+0.*(t>=2&t<=-4)'); xlabel('t'); ylabel('f(t)'); %f2=0.25*((-2)*t+5).*(t>4&t<12)+1.*(t>0&t<4)+0.*(t>=12&t<=0&t== 4); f2=-0.25*(t+1).*(t>2&t<4)+1.*(t>1&t<2)+0.*(t>=4&t<=1&t==2); subplot(212); plot(t,f2); grid on;
title('f2=0.25*(-2*t+5).*(t>-4&t<0)+1.*(t>0&t<2)+0.*(t>=2&t<=-4&t= =0)'); xlabel('t'); ylabel('f(-2t+4)'); 系统响应运算 3、已知描述系统的微分方程和激励信号e(t) 分别如下，试用解析方法求系统的单位冲激响应h(t)和零状态响应r(t)，并用MATLAB绘出系统单位冲激响应和系统零状态响应的波形。 ； a=[1 4 4]; b=[1 3]; subplot(211) impulse(b,a,4); %冲激响应函数 title('?μí3μ￥??3??¤?ìó|'); c=[1 4 4]; d=[1 3]; p1=0.001; t=0:p1:10;
信号与系统期末考试试题6 课程名称： 信号与系统 一、选择题（共10题，每题3分 ，共30分，每题给出四个答案，其中只有一个正确的） 1、 卷积f 1(k+5)*f 2(k-3) 等于 。 （A ）f 1(k)*f 2(k) （B ）f 1(k)*f 2(k-8)（C ）f 1(k)*f 2(k+8)（D ）f 1(k+3)*f 2(k-3) 2、 积分dt t t ?∞ ∞--+)21()2(δ等于 。 （A ）1.25（B ）2.5（C ）3（D ）5 3、 序列f(k)=-u(-k)的z 变换等于 。 （A ） 1 -z z （B ）- 1 -z z （C ） 1 1-z （D ） 1 1--z 4、 若y(t)=f(t)*h(t),则f(2t)*h(2t)等于 。 （A ） )2(4 1t y （B ） )2(2 1t y （C ） )4(4 1t y （D ） )4(21t y 5、 已知一个线性时不变系统的阶跃相应g(t)=2e -2t u(t)+)(t δ,当输入f(t)=3e —t u(t)时，系 统的零状态响应y f (t)等于 （A ）(-9e -t +12e -2t )u(t) （B ）(3-9e -t +12e -2t )u(t) （C ）)(t δ+(-6e -t +8e -2t )u(t) （D ）3)(t δ +(-9e -t +12e -2t )u(t) 6、 连续周期信号的频谱具有 （A ） 连续性、周期性 （B ）连续性、收敛性 （C ）离散性、周期性 （D ）离散性、收敛性 7、 周期序列2)455.1(0 +k COS π的 周期N 等于 （A ） 1（B ）2（C ）3（D ）4 8、序列和 ()∑∞ -∞ =-k k 1δ等于 （A ）1 (B) ∞ (C) ()1-k u (D) ()1-k ku 9、单边拉普拉斯变换()s e s s s F 22 12-+= 的愿函数等于 ()()t tu A ()()2-t tu B ()()()t u t C 2- ()()()22--t u t D 10、信号()()23-=-t u te t f t 的单边拉氏变换()s F 等于 ()A ()()()232372+++-s e s s ()()2 23+-s e B s