频域分析频谱分析：
上式表明
为基频的整数倍的谐波分量之和。
幅值谱
相位谱
例：
傅氏级数也可以写成复指数函数形式。根据欧拉公式
其中，
非周期信号
第5章频域分析法 学习要点 1 频率特性的概念，常用数学描述与图形表示方法； 2 典型环节的幅相频率特性与对数频率特性表示及特点； 3 系统开环幅相频率特性与对数频率特性的图示要点； 4 应用乃奎斯特判据判断控制系统的稳定性方法； 5 对数频率特性三频段与系统性能的关系； 6 计算频域参数与性能指标； 思考与习题祥解 题判断下列概念的正确性 ω的正弦信号加入线性系统，这个系统的稳态输出也将是同 (1) 将频率为 一频率的。 M仅与阻尼比ξ有关。 (2) 对于典型二阶系统，谐振峰值 p (3) 在开环传递函数中增加零点总是增加闭环系统的带宽。 (4) 在开环传递函数中增加极点通常将减少闭环系统的带宽并同时降低稳定性。 (5) 对于最小相位系统，如果相位裕量是负值，闭环系统总是不稳定的。 (6) 对于最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (7) 对于最小相位系统，如果幅值裕量是负分贝值，闭环系统总是不稳定的。 (8) 对于非最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (9) 对于非最小相位系统，须幅值裕量大于1且相位裕量大于0，闭环系统才是稳定的。 (10) 相位穿越频率是在这一频率处的相位为0。 (11) 幅值穿越频率是在这一频率处的幅值为0dB。 (12) 幅值裕量在相位穿越频率处测量。 (13) 相位裕量在幅值穿越频率处测量。 (14) 某系统稳定的开环放大系数25 K<，这是一个条件稳定系统。 (15) 对于（-2/ -1/ -2)特性的对称最佳系统，具有最大相位裕量。 (16) 对于（-2/ -1/ -3)特性的系统，存在一个对应最大相位裕量的开环放大系数值。 (17) 开环中具有纯时滞的闭环系统通常比没有时滞的系统稳定性低些。 (18) 开环对数幅频特性过0分贝线的渐近线斜率通常表明了闭环系统的相对稳定性。 M和频带宽BW (19) Nichols图可以用于找到一个闭环系统的谐振峰值 p 的信息。
实验二连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义； 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质； 5、学习掌握利用Matlab语言编写计算CTFS、CTFT和DTFT的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT、DTFT的若干重要性质。 基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用Matlab编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS分析 任何一个周期为T 1 的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。 三角傅里叶级数为： ∑∞=+ + = 1 0 0 )] sin( ) cos( [ )( k k k t k b t k a a t xω ω 2.1
或： ∑∞ =++=100)cos()(k k k t k c a t x ?ω 2.2 其中1 02T πω=，称为信号的基本频率（Fundamental frequency ），k k b a a ，和，0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度，k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率0ωk 的函数，绘制出它们与0ωk 之间的图像，称为信号的频谱图（简称“频谱”），k c －0ωk 图像为幅度谱，k ?－0ωk 图像为相位谱。 三角形式傅里叶级数表明，如果一个周期信号x(t)，满足狄里克利条件，就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的正弦信号所组成，其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component)，其幅度（amplitude ）为k c 。也可以反过来理解三角傅里叶级数：用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。 指数形式的傅里叶级数为： ∑∞-∞== k t jk k e a t x 0)(ω 2.3 其中，k a 为指数形式的傅里叶级数的系数，按如下公式计算： ? --=2/2/1110)(1T T t jk k dt e t x T a ω 2.4 指数形式的傅里叶级数告诉我们，如果一个周期信号x(t)，满足狄里克利条件，那么，它就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的周期复指数信号所组成，其中每一个不同频率的周期复指数信号称为基本频率分量，其复幅度
实 验 报 告 学生姓名： 学 号： 指导教师： 一、实验室名称：数字信号处理实验室 二、实验项目名称：采样的时域及频域分析 三、实验原理： 1、采样的概念：采样是将连续信号变化为离散信号的过程。 1. A 、理想采样：即将被采样信号与周期脉冲信号相乘 B 、实际采样：将被采样信号与周期门信号相乘，当周期门信号的宽度很小，可近似为周期脉冲串。 根据傅里叶变换性质 00 0()() ()() ??()()()()()()(()) FT FT a a T n n FT a a T a T a a n n x t X j T j x t x t T x nT t nT X j X j n ωδωδδδω=+∞=+∞=-∞ =-∞ ←?→Ω←?→Ω==-←?→Ω=Ω-Ω∑ ∑式中T 代表采样间隔，01 T Ω= 由上式可知：采样后信号的频谱是原信号频谱以0Ω为周期的搬移叠加 结论：时域离散化，频域周期化；频谱周期化可能造成频谱混迭。 ) (t T δ^ T ^)t
C 、低通采样和Nyquist 采样定理 设()()a a x t X j ?Ω且()0,2a M M X j f πΩ=Ω>Ω=当， 即为带限信号。则当采样频率满足2/22s M M f f π≥Ω=时,可以从采样后的 ^ ()()()a a s s n x t x nT t nT δ∞ =-∞ = -∑信号无失真地恢复()a x t 。称2M f 为奈奎斯特频率， 1 2 N M T f = 为奈奎斯特间隔。 注意： 实际应用中，被采信号的频谱是未知的，可以在ADC 前加一个滤波器（防混迭滤波器）。 2、低通采样中的临界采样、欠采样、过采样的时域及频域变化情况。 低通采样中的临界采样是指在低通采样时采样频率2s M f f = 低通采样中的欠采样是指在低通采样时采样频率2s M f f ≤ 低通采样中的欠采样是指在低通采样时采样频率2s M f f ≥ 设一带限信号的频谱如下： ) () a G j Ω0 m -ΩΩ m Ω0 T T
1.周期信号的频谱 周期信号在满足一定条件时，可以分解为无数三角信号或指数之和。这就是周期信号的傅里叶级数展开。在三角形式傅里叶级数中，各谐波分量的形式为()1cos n n A n t ω?+；在指数形式傅里叶级数中，分量的形式必定为1j n t n F e ω 与1-j -n t n F e ω 成对出现。为了把周期信号所具有的各 次谐波分量以及各谐波分量的特征（如模、相角等）形象地表示出来，通常直接画出各次谐波的组成情况，因而它属于信号的频域描述。 以周期矩形脉冲信号为lifenxi 周期信号频谱的特点。周期矩形信号在一个周期（-T/2,T/2）内的时域表达式为 ,2 0,>2 ()A t T t f t ττ ≤?=?? （2-6） 其傅里叶复数系数为 12 n n A F Sa T ωττ?? = ??? （2-7） 由于傅里叶复系数为实数，因而各谐波分量的相位为零（n F 为正）或为π±（n F 为负），因此不需要分别画出幅度频谱n F 与相位频谱n φ。可以直接画出傅里叶系数n F 的分布图。 如图2.4.1所示。该图显示了周期性矩形脉冲信号()T f t 频谱的一些性质，实际上那个也是周期性信号频谱的普遍特性： ① 离散状频谱。即谱线只画出现在1ω的整数倍频率上，两条谱线的间隔为1ω（等于2π/t ）。 ② 谱线宽度的包络线按采样函数()1/2a S n ωτ的规律变化。如图2.4.2所示。但1ω 为 2π τ 时，即( )2m π ωτ =（m=1，2，……）时，包络线经过零点。在两相邻 零点之间，包络线有极值点，极值的大小分别为-0.212()2A T τ，
信号时域与频域分析 实验报告 姓名：杨 班级：机械 学号： 213
实验数据中，电机转速为1200r/min，采样频率为1280Hz。Hz3为X位移振幅数据，Hz4为Y位移振幅数据，Hz5为速度振幅数据。 Matlab中信号特征对应函数编程 ma = max(Hz) %最大值 mi = min(Hz) %最小值 me = mean(Hz) %平均值 pk = ma-mi %峰-峰值 va = var(Hz); %方差 st = std(Hz); %标准差 ku = kurtosis(Hz); %峭度 rm = rms(Hz); %均方根 一、X轴位移测量分析 plot(Fs3,Hz3)时域图： ma =52.0261 mi =56.7010 me =1.8200 pk =108.7271 va =1.3870e+03 st =37.2431 ku =1.5462 rm =37.2693 频域图： fs=1280; x=Hz3; N=length(Hz3); df=fs/N; f=0:df:N*df-df; y=fft(x); y=abs(y)*2/N; figure(1); plot(f,y); xlabel('频率/Hz') ylabel('幅值') 频谱幅值取得最大值51.9847um，频率为20Hz，与电机转速对应频率一致，应为电机轴未动平衡所致；二倍频处有较大振幅，可能为轴承间隙过大所致。
二、Y轴位移测量分析 plot(Fs4,Hz4)时域图： ma =61.3987 mi =-74.6488 me =-1.1948 pk =136.0475 av =42.6109 va =2.2428e+03 st =47.3582 ku =1.5135 rm =47.3501 频域图： fs=1280; x=Hz4; N=length(Hz4); df=fs/N; f=0:df:N*df-df; y=fft(x); y=abs(y)*2/N; figure(1); plot(f,y); xlabel('频率/Hz') ylabel('幅值') 频谱幅值取得最大值66.6319um，频率为20Hz，与电机转速对应频率一致，应为电机轴未动平衡所致；二倍频处有较大振幅，可能为轴承间隙过大所致。
周期信号的时域及其频域分析 姓名：张敏靓学号：1007433014 一、实验目的 1.掌握Multisim软件的应用及用虚拟仪器对周期信号的频谱测量 2.掌握选频电平表的使用，对信号发生器输出信号（方波、矩形波、 三角波等）频谱的测量 二、实验原理 周期信号的傅里叶级数分析法，可以把周期信号表示为三角傅里叶级数或指数傅里叶级数，其中周期信号满足。 1. 周期信号表示为三角傅里叶级数 2. 周期信号表示为指数傅里叶级数 其中， 周期矩形信号的频谱
三、实验内容 1.在Multisim上实现周期信号的时域、频域测量及分析 （1）绘制测量电路 （2）周期信号时域、频域（幅度频谱）的仿真测量 虚拟信号发生器分别设置如下参数： 周期方波信号：周期T=100μs，脉冲宽度τ=50μs，脉冲幅度 V P=5V； 周期矩形信号：周期T=100μs，脉冲宽度τ=20μs，脉冲幅度 V P=5V； 周期三角波信号：周期T=200μs，脉冲幅度V P=5V； 采用虚拟示波器及虚拟频谱仪分别测量上述信号的时域、频域波形并保存测试波形及数据。
2.周期信号时域、频域（幅度频谱）的测量 信号发生器、示波器、选频电平表的连线如上图所示。信号发生器的输出信号分别为周期分别信号、周期矩形信号、周期三角波信号，参数设置同仿真测量。采用示波器及选频电平表对信号发生器的输出信号分别测量，并将测量数据记录下表中。
四、实验总结 1.在周期矩形信号的实验中，信号频率减小，频谱减小；信号占空 比减小，频谱减小；幅度值减小，频谱减小。 2.未安装Origin绘图软件，Excel绘图未能达到理想效果。
第五章 频域分析法 时域分析法具有直观、准确的优点。如果描述系统的微分方程是一阶或二阶的，求解后可利用时域指标直接评估系统的性能。然而实际系统往往都是高阶的，要建立和求解高阶系统的微分方程比较困难。而且，按照给定的时域指标设计高阶系统也不是一件容易的事。 本章介绍的频域分析法，可以弥补时域分析法的不足。因为频域法是基于频率特性或频率响应对系统进行分析和设计的一种图解方法，故其与时域分析法相比有较多的优点。首先，只要求出系统的开环频率特性，就可以判断闭环系统是否稳定。其次，由系统的频率特性所确定的频域指标与系统的时域指标之间存在着一定的对应关系，而系统的频率特性又很容易和它的结构、参数联系起来。因而可以根据频率特性曲线的形状去选择系统的结构和参数，使之满足时域指标的要求。此外，频率特性不但可由微分方程或传递函数求得，而且还可以用实验方法求得。对于某些难以用机理分析方法建立微分方程或传递函数的元件(或系统)来说，具有重要的意义。因此，频率法得到了广泛的应用，它也是经典控制理论中的重点内容。 5.1 频率特性 对于线性定常系统，若输入端作用一个正弦信号 t U t u ωsin )(= (5—1) 则系统的稳态输出y(t)也为正弦信号，且频率与输人信号的频率相同，即 ) t Y t y ?ω+=sin()( (5—2) u(t)和y(t)虽然频率相同，但幅值和相位不同，并且随着输入信号的角频率ω的改变，两者之间的振幅与相位关系也随之改变。这种基于频率ω的系统输入和输出之间的关系称之为系统的频率特性。 不失一般性，设线性定常系统的传递函数G(s)可以写成如下形式 ) () () () () ())(() ()()()(1 21s A s B p s s B p s p s p s s B s U s Y s G n j j n = +=+++== ∏=Λ (5—3) 式中B(s)——传递函数G(s)的m 阶分子多项式，s 为复变量； A(s)——传递函数G(s)的n 阶分母多项式 (n ≥m)； n p p p ---,,,21Λ—传递函数G(s)的极点，这些极点可能是实数，也可能是复数，对稳定的系统采说，它们都应该有负的实部。 由式(5—1)，正弦输入信号u(t)的拉氏变换为(查拉氏变换表) ) )(()(22ωωω ωωj s j s U s U s U -+=+= （5—4）
实验二 连续时间信号的频域分析 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义； 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质； 5、学习掌握利用Matlab 语言编写计算CTFS 、CTFT 和DTFT 的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT 、DTFT 的若干重要性质。 基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用Matlab 编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS 分析 任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。 三角傅里叶级数为： ∑∞ =++=1 000)]sin()cos([)(k k k t k b t k a a t x ωω 2.1 或： ∑∞=++=1 00)cos()(k k k t k c a t x ?ω 2.2 其中1 02T πω=，称为信号的基本频率（Fundamental frequency ），k k b a a ，和，0分别是信号)(t x 的直流分量、 余弦分量幅度和正弦分量幅度，k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率0ωk 的函数，绘制出它们与0ωk 之间的图像，称为信号的频谱图（简称“频谱”），k c －0ωk 图像为幅度谱，k ?－0ωk 图像为相位谱。 三角形式傅里叶级数表明，如果一个周期信号x(t)，满足狄里克利条件，就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的正弦信号所组成，其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component)，其幅度（amplitude ）为k c 。也可以反过来理解三角傅里叶级数：用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。 指数形式的傅里叶级数为：
课程设计任务书 题目 专业、班级电信1班学号姓名 主要内容、基本要求、主要参考资料等： 基于钟表设计的常识，给出时、分、秒的设计思路，并利用硬件编程语言VHDL或者Verilog-HDL来实 现。要求具有基本功能如调整时间对表、闹铃、计时器等，给出完成控制电路所需要的设计模块；给出硬 件编程语言的实现，并进行仿真；给出下载电路的设计，设计为2种下载方法，其中一种必须为JTAG；同 时设计者报告不允许雷同。 参考资料： 1、潘松、黄继业《EDA技术及其应用》（第四版）科学出版社 2009 2、樊昌信《通信原理》电子出版社 完成期限： 指导教师签名： 课程负责人签名： 年月日
目录 摘要…………………………………………………………………………………II
ABSTRACT……………………………………………………………………………III 绪论…………………………………………………………………………………III 1傅里叶变换原理概述 (1) 1.1 傅里叶变换及逆变换的MATLAB实现 (2) 2 用MATLAB实现典型非周期信号的频域分析 (3) 2.1 单边指数信号时域波形图、频域图 (3) 2.2 偶双边指数信号时域波形图、频域图 (4) 2.3 奇双边指数信号时域波形图、频域图 (4) 2.4 直流信号时域波形图、频域图 (5) 2.5 符号函数信号时域波形图、频域图 (5) 2.6 单位阶跃信号时域波形图、频域图 (6) 2.7 单位冲激信号时域波形图、频域图 (6) 2.8 门函数信号时域波形图、频域图 (7) 3 用MATLAB实现信号的幅度调制 (8) 3.1 实例1 (8) 3.2 实例2 (10) 4 实现傅里叶变换性质的波形仿真 (11) 4.1 尺度变换特性 (11) 4.2 时移特性 (14) 4.3 频移特性 (16) 4.4 时域卷积定理 (18) 4.5 对称性质 (20) 4.6 微分特性 (22) 心得体会 (25) 参考文献 (26) 附录 (27)
第3章连续时间信号与系统的频域分析3.1 学习要求 1、掌握周期信号的频谱及其特点； 2、了解周期信号的响应问题； 3、掌握非周期信号的频域描述——傅立叶变换； 4、熟练掌握傅立叶变换的性质与应用； 5、掌握系统的频域特性及响应问题； 6、了解系统的无失真传输和理想滤波。 3.2 本章重点 1、频谱的概念及其特性； 2、傅里叶变换及其基本性质； 3、响应的频域分析方法； 4、系统频率响应的概念。 3.3 知识结构
3.4内容摘要 3.4.1信号的正交分解 两个矢量1V 和2V 正交的条件是这两个矢量的点乘为零，即： o 1212cos900?=?=V V V V 若有一个定义在区间()12,t t 的实函数集{}()(1,2,,)i g t i n =L ，在该集合中所有的函数满足 ?????=≠===??2 1 21,,2,1,0)()(,,2,1)(2t t j i t t i i n j j i dt t g t g n i k dt t g ΛΛ 则称这个函数集为区间()12,t t 上的正交函数集。式中i k 为常数，当1i k =时，称此函数集为归一化正交函数集。 若实函数集{}(),1,2,,i g t i n =L 是区间()12,t t 内的正交函数集，且除()i g t 之外 {}(),1,2,,i g t i n =L 中不存在()x t 满足下式 2 1 20()t t x t dt <<∞?且2 1 ()()0t i t x t g t dt =? 则称函数集{}(),1,2,,i g t i n =L 为完备正交函数集。 若在区间()12,t t 上找到了一个完备正交函数集{}(),1,2,,i g t i n =L ，那么，在此区间的信号()x t 可以精确地用它们的线性组合来表示 11221 ()()()()()n n i i i x t C g t C g t C g t C g t ∞ ==++++=∑L L 各分量的标量系数为 2 1 21 2 ()()d ()d t i t i t i t x t g t t C g t t = ?? 系数i C 只与()x t 和()i g t 有关，而且可以互相独立求取。 3.4.2周期信号的傅里叶级数 1、三角形式的傅里叶级数 0001 ()(cos sin )n n n x t a a n t b n t ωω∞ ===++∑
信号与系统 实验报告 实验三周期信号的频谱分析 实验报告评分：_______ 实验三周期信号的频谱分析 实验目的： 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；
3、掌握各种典型的连续时间非周期信号的频谱特征。 实验内容： （1）Q3-1 编写程序Q3_1，绘制下面的信号的波形图： 其中，0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos( 0t)、cos(3 0t)、cos(5 0t)和x(t) 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。 程序如下: clear,%Clear all variables close all,%Close all figure windows dt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of time w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t); N=input('Type in the number of the harmonic components N='); x=0; for q=1:N; x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; end subplot(221) plot(t,x1)%Plot x1 axis([-2 4 -2 2]); grid on, title('signal cos(w0.*t)') subplot(222) plot(t,x2)%Plot x2 axis([-2 4 -2 2]); grid on, title('signal cos(3*w0.*t))') subplot(223) plot(t,x3)%Plot x3 axis([-2 4 -2 2])
第5章频域分析法 5.1 学习要点 1 频率特性的概念，常用数学描述与图形表示方法； 2 典型环节的幅相频率特性与对数频率特性表示及特点； 3 系统开环幅相频率特性与对数频率特性的图示要点； 4 应用乃奎斯特判据判断控制系统的稳定性方法； 5 对数频率特性三频段与系统性能的关系； 6 计算频域参数与性能指标； 5.2 思考与习题祥解 题5.1 判断下列概念的正确性 ω的正弦信号加入线性系统，这个系统的稳态输出也将是同 (1) 将频率为 一频率的。 M仅与阻尼比ξ有关。 (2) 对于典型二阶系统，谐振峰值 p (3) 在开环传递函数中增加零点总是增加闭环系统的带宽。 (4) 在开环传递函数中增加极点通常将减少闭环系统的带宽并同时降低稳定性。 (5) 对于最小相位系统，如果相位裕量是负值，闭环系统总是不稳定的。 (6) 对于最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (7) 对于最小相位系统，如果幅值裕量是负分贝值，闭环系统总是不稳定的。 (8) 对于非最小相位系统，如果幅值裕量大于1，闭环系统总是稳定的。 (9) 对于非最小相位系统，须幅值裕量大于1且相位裕量大于0，闭环系统才是稳定的。 (10) 相位穿越频率是在这一频率处的相位为0。 (11) 幅值穿越频率是在这一频率处的幅值为0dB。 (12) 幅值裕量在相位穿越频率处测量。 (13) 相位裕量在幅值穿越频率处测量。 (14) 某系统稳定的开环放大系数25 K<，这是一个条件稳定系统。 (15) 对于（-2/ -1/ -2)特性的对称最佳系统，具有最大相位裕量。 (16) 对于（-2/ -1/ -3)特性的系统，存在一个对应最大相位裕量的开环放大系数值。 (17) 开环中具有纯时滞的闭环系统通常比没有时滞的系统稳定性低些。 (18) 开环对数幅频特性过0分贝线的渐近线斜率通常表明了闭环系统的相对稳定性。 M和频带宽BW (19) Nichols图可以用于找到一个闭环系统的谐振峰值 p 的信息。
第一章引言 语音信号是一种非平稳的时变信号，它携带着各种信息。在语音编码、语音合成、语音识别和语音增强等语音处理中无一例外需要提取语音中包含的各种信息。语音信号分析的目的就在与方便有效的提取并表示语音信号所携带的信息。语音信号分析可以分为时域和频域等处理方法。语音信号可以认为在短时间内(一般认为在 10~30ms 的短时间内)近似不变，因而可以将其看作是一个准稳态过程, 即语音信号具有短时平稳性。任何语音信号的分析和处理必须建立在“短时”的基础上, 即进行“短时分析”。 时域分析:直接对语音信号的时域波形进行分析，提取的特征参数有短时能量，短时平均过零率，短时自相关函数等。 频域分析:对语音信号采样，并进行傅里叶变换来进行频域分析。主要分析的特征参数：短时谱、倒谱、语谱图等。 本文采集作者的声音信号为基本的原始信号。对语音信号进行时频域分析后，进行加白噪声处理并进行了相关分析，设计滤波器并运用所设计的滤波器对加噪信号进行滤波, 绘制滤波后信号的时域波形和频谱。整体设计框图如下图所示： 图1.1时频域分析设计图 图1.2加噪滤波分析流程图
第二章 语音信号时域分析 语音信号的时域分析可直接对语音信号进行时域波形分析，在此只只针对语音信号的短时能量、短时平均过零率、短时自相关函数进行讨论。 2.1窗口选择 由人类的发生机理可知，语音信号具有短时平稳性，因此在分析讨论中需要对语音信号进行加窗处理进而保证每个短时语音长度为10~30ms 。通常选择矩形窗和哈明窗能得到较理想的“短时分析”设计要求。两种窗函数的时域波形如下图2.1所示： sample w （n ） sample w （n ） 图2.1 矩形窗和Hamming 窗的时域波形 矩形窗的定义：一个N 点的矩形窗函数定义为如下 {1,00,()n N w n ≤<=其他 （2.1） 哈明窗的定义：一个N 点的哈明窗函数定义为如下 0.540.46cos(2),010,()n n N N w n π-≤<-??? 其他 = （2.2） 这两种窗函数都有低通特性，通过分析这两种窗的频率响应幅度特性可以发现（如图2.2）：矩形窗的主瓣宽度小（4*pi/N ），具有较高的频率分辨率，旁瓣峰值大（-13.3dB ），会导致泄漏现象；哈明窗的主瓣宽8*pi/N ，旁瓣峰值低（-42.7dB ），可以有效的克服泄漏现象，具有更平滑的低通特性。因此在语音频谱分析时常使用哈明窗，在计算短时能量和平均幅度时通常用矩形窗。表2.1对比了这两种窗函数的主瓣宽度和旁瓣峰值。
课程实验报告 题目：用Matlab进行 信号与系统的时、频域分析 学院 学生姓名 班级学号 指导教师 开课学院 日期 用Matlab进行信号与系统的时、频域分析 一、实验目的 进一步了解并掌握Matlab软件的程序编写及运行； 掌握一些信号与系统的时、频域分析实例； 了解不同的实例分析方法，如：数值计算法、符号计算法； 通过使用不同的分析方法编写相应的Matlab程序； 通过上机，加深对信号与系统中的基本概念、基本理论和基本分析方法的理解。 二、实验任务 了解数值计算法编写程序，解决实例； 在Matlab上输入三道例题的程序代码，观察波形图； 通过上机实验，完成思考题； 完成实验报告。 三、主要仪器设备
硬件：微型计算机 软件：Matlab 四、 实验内容 （1） 连续时间信号的卷积 已知两个信号)2()1()(1---=t t t x εε和)1()()(2--=t t t x εε，试分别画出)(),(21t x t x 和卷积)()()(21t x t x t y *=的波形。 程序代码： T=0.01; t1=1;t2=2; t3=0;t4=1; t=0:T:t2+t4; x1=ones(size(t)).*((t>t1)-(t>t2)); x2=ones(size(t)).*((t>t3)-(t>t4)); y=conv(x1,x2)*T; subplot(3,1,1),plot(t,x1); ylabel('x1(t)'); subplot(3,1,2),plot(t,x2); ylabel('x2(t)'); subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1)); ylabel('y(t)=x1*x2'); xlabel('----t/s'); （2）已知两个信号)()(t e t x t ε-=和)()(2/t te t h t ε-=，试用数值计算法求卷积，并分别画出)(),(t h t x 和卷积)()()(t h t x t y *=的波形。 程序代码： t2=3;t4=11; T=0.01; t=0:T:t2+t4; x=exp(-t).*((t>0)-(t>t2)); h=t.*exp(-t/2).*((t>0)-(t>t4)); y=conv(x,h)*T; yt=4*exp(-t)+2*t.*exp(-1/2*t)-4*exp(-1/2*t); subplot(3,1,1),plot(t,x); ylabel('x(t)'); subplot(3,1,2),plot(t,h); ylabel('h(t)'); subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1),t,yt,'--r'); legend('by numberical','Theoretical'); ylabel('y=x*h'); xlabel('----t/s'); （3）求周期矩形脉冲信号的频谱图，已知s T s A 5.0,1.0,1===τ
第四章 连续信号的频域分析 将信号分解为若干不同频率的正弦信号或虚指数信号，实质上是将信号在频率域上进行分解，因此根据这种基本思想对信号和系统的分析称为频域分析。这种分解过程是通过傅里叶级数和傅里叶变换这一数学工具来实现的。 本章首先介绍连续信号的傅里叶级数和傅里叶变换，熟悉信号频谱的概念。 4.1 基本要求 1．基本要求 ? 了解傅里叶级数和傅里叶变换的定义及其物理含义； ? 掌握信号频谱和频谱密度的概念； ? 了解连续谱和离散谱的特点和区别； ? 掌握傅里叶变换的常用性质； ? 掌握周期信号傅里叶变换的求解方法。 2．重点和难点 ? 傅里叶变换的性质及其应用 4.2 知识要点 1．周期信号的傅里叶级数 （1）傅里叶级数展开式 三角形式：∑∑∞ =∞=+Ω+=Ω+Ω+=1010)cos(2)]sin()cos([2)(n n n n n n t n A A t n b t n a a t f ?（4-1） 指数形式： ∑∑∞ -∞ =+Ω∞ -∞ =Ω= =n t n n n t n n n F F t f )j(j e e )(? （4-2） 其中 ? +Ω= T t t n t t n t f T a 00 d cos )(2 ，n =0，1，2，? （4-3） ? +Ω= T t t n t t n t f T b 00 d sin )(2，n =1，2，? （4-4） 且
n n n n n n a b b a A a A arctg , ,2 200-=+==? （4-5） ?+Ω-= T t t t n n t t f T F 00 d e )(1j （4-6） （2）两种形式之间的转换关系 0)( e 2 1 j ≥=n A F n n n ? （4-7） 并且|F n |为偶函数，?n 为奇函数，即 ||||n n F F -=，||||n n -=?? （4-8） （3）傅里叶级数的物理含义 通过傅里叶级数可以将任意周期信号f (t )分解为若干个正弦信号（三角形式）或复简谐信号（指数形式）的叠加。每个正弦信号分量的频率为周期信号基波频率的n 倍（n ?0），即n ?，而幅度为A n 或者2|F n |，相位为?n ，将其称作第n 次谐波分量。特别地，将频率为0（即n =0）的分量称为直流分量，幅度为A 0/2或者F 0；频率等于基波频率?（即n =1）的分量称为基波分量。 2．周期信号的频谱 通过傅里叶级数可以将时域中的周期信号分解为直流分量、基波分量和各次谐波分量之和，傅里叶级数展开式中的A n 、?n 或傅里叶系数F n 分别代表了各分量的幅度和相位随谐波次数n （从而频率n ?）的变化关系，称为周期信号的频谱，其中A n 或|F n |称为幅度谱，?n 称为相位谱。 A n 或|F n |、?n 都是关于整型变量n 的实函数，分别以其为纵轴，以n （或者n ?）为横轴，得到的图形称为周期信号的幅度谱图和相位谱图，合称为周期信号的频谱图。 但是，在三角形式的傅里叶级数中，A n 和?n 的自变量n 只能取非负的整数，因此称为单边频谱，而在F n 中，n 可以为任意的整数，相应地将F n 称为双边频谱。对同一个周期信号，其单边和双边频谱可以通过式（4-7）进行相互转换。 所有周期信号的频谱都具有离散性，因此称为离散谱。 3．非周期信号的傅里叶变换及其频谱密度 非周期信号的傅里叶变换及傅里叶反变换的定义为 ?∞ ∞--=t t f F t d e )()j (j ωω （4-9） ?∞ ∞ -= ωωωd )e (j 2π1)(j t F t f （4-10） 其中正变换用于根据信号的时域表达式求其频谱表达式，反变换用于根据其频谱表达式求时域表达式。 通过傅里叶变换可以将信号分解为不同频率的复简谐信号的叠加，而信号的傅里叶变换F (j ?)反映了信号中各分量的幅度和相位随其频率? 的变化关系，称为信号的频谱密度，又称为频谱密度函数或频谱函数。 教材表4-1中列出了一些基本信号的傅里叶变换，在求解复杂信号的傅里叶变换和频谱密度时经常用到。 4．傅里叶变换的性质
一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义； 4、学习利用MATLAB 语言编写计算CTFS 和CTFT 的仿真程序。 基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用MATLAB 编程完成相关的傅里叶变换的计算。 二、实验原理及方法 1、连续时间周期信号的傅里叶级数CTFS 分析 任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。 其中三角傅里叶级数为： ∑∞ =++=1000)]sin()cos([)(k k k t k b t k a a t x ωω 9.1 或： ∑∞ =++ =1 00)cos()(k k k t k c a t x ?ω 9.2 其中1 02T π ω= ，称为信号的基本频率，k k b a a ，和，0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度，k k c ?、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率0ωk 的函数，绘制出它们与0ωk 之间的图像，称为信号的频谱图（简称“频谱”）， k c －0ωk 图像为幅度谱，k ?－0ωk 图像为相位谱。 指数形式的傅里叶级数为： ∑∞ -∞ == k t jk k e a t x 0)(ω 9.3 其中，k a 为指数形式的傅里叶级数的系数，按如下公式计算：
? --= 2 /2 /1 110)(1 T T t jk k dt e t x T a ω 9.4 假设谐波项数为N ，则上面的和成式为： ∑-== N N k t jk k e a t x 0)(ω 9.5 显然，N 越大，所选项数越多，有限项级数合成的结果越逼近原信号x(t)。 2、连续时间信号傅里叶变换----CTFT 傅里叶变换在信号分析中具有非常重要的意义，它主要是用来进行信号的频谱分析的。傅里叶变换和其逆变换定义如下： ?∞ ∞ --= dt e t x j X t j ωω)()( 9.6 ? ∞ ∞ -= ωωπ ωd e j X t x t j )(21 )( 9.7 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。任意非周期信号，如果满足狄里克利条件，那么，它可以被看作是由无穷多个不同频率（这些频率都是非常的接近）的周期复指数信号e j ωt 的线性组合构成的，每个频率所对应的周期复指数信号e j ωt 称为频率分量其相对幅度为对应频率的|X(j ω)|之值，其相位为对应频率的X(j ω)的相位。 给定一个连续时间非周期信号x(t)，它的频谱也是连续且非周期的。对于连续时间周期信号，也可以用傅里变换来表示其频谱，其特点是，连续时间周期信号的傅里叶变换时有冲激序列构成的，是离散的——这是连续时间周期信号的傅里叶变换的基本特征。 3、连续周期信号的傅里叶级数CTFS 的MATLAB 实现 3.1 傅里叶级数的MATLAB 计算 设周期信号x(t)的基本周期为T 1，且满足狄里克利条件，则其傅里叶级数的系数可由式9.4计算得到。式9.4重写如下： ?--= 2 /2 /1 110)(1T T t jk k dt e t x T a ω 基本频率为： 1 02T πω= 对周期信号进行分析时，我们往往只需对其在一个周期内进行分析即可，通常选择主周期。假定x 1(t)是x(t)中的主周期，则
福建农林大学计算机与信息学院 信息工程类 课程设计报告 课程名称：信号与系统 课程设计题目：连续时间信号的频域分析 姓名： 系：电子信息工程 专业：电子信息工程 年级：2008 学号： 指导教师： 职称： 2011 年 1 月10 日
福建农林大学计算机与信息学院信息工程类 课程设计结果评定
目录 1课程设计的目的 (1) 2课程设计的要求 (1) 3课程设计报告内容.....................................................................1-13 3.1连续信号的设计..................................................................1-11 3.2验证傅里叶变换的调制定理 (11) 3.3周期信号及其频谱 (12) 4总结 (13) 参考文献 (14)
连续时间信号的频域分析 1.课程设计的目的 （1）熟悉MATLAB语言的编程方法及MATLAB指令； （2）掌握连续时间信号的基本概念； （3）掌握门函数、指数信号和抽样信号的表达式和波形； （4）掌握连续时间信号的傅里叶变换及其性质； （5）掌握连续时间信号频谱的概念以及幅度谱、相位谱的表示； （6）掌握利用MATLAB进行信号的傅里叶变换以及时域波形和频谱的表示；（7）通过连续时间信号的频域分析，更深刻地理解了连续时间信号的时域和频域间的关系，加深了对连续时间信号的理解。 2.课程设计的要求 （1）自行设计以下连续信号：门函数、指数信号和抽样信号。要求：（a）画出以上信号的时域波形图； （b）实现以上信号的傅里叶变换，画出以上信号的幅度谱及相位谱，并对相关结果予以理论分析； （c）对其中一个信号进行时移和尺度变换，分别求变换后信号的傅里叶变换，验证傅里叶变换的时移和尺度变换性质。 （2）自行设计信号，验证傅里叶变换的调制定理。 （3）自行设计一个周期信号，绘出该信号的频谱，并观察周期信号频谱的特点。 3.课程设计报告内容 3.1（a）①门函数(矩形脉冲)： MATLAB中矩形脉冲信号用rectpuls函数表示： y=rectpuls (t,width) %width缺省值为1 >> t=-2:0.001:2; T=2; yt=rectpuls (t,T); plot(t,yt); axis([-2,2,0,1.5]); grid on; %显示格线
一，实验目的四，心得体会 了解信号频谱和信号频域，掌握其特性。 一，实验原理 实验主要分为四个部分，分别分析了连续和离散信号的周期、非周期情况下特性。 1.连续周期信号的频谱分析 首先手算出信号的傅里叶级数，得出信号波形，然后通过代码画出信号波形图。 2.连续非周期信号的频谱分析 先由非周期信号的时域信号得到它的频谱X（w），再通过MATLAB 求出其傅里叶变换并绘出图形。 X=fourier(x) x=ifourier(x) ①符号运算法 syms t ②数值积分法 quad（fun，a，b） ③数值近似法 3.离散周期信号的频谱分析 X=fft（x） 4.离散非周期信号的频谱分析 可以化为两个相乘的矩阵，从而由MATLAB实现。
三，实验内容 （1）已知x（t）是如图周期矩形脉冲信号。 1）.计算该信号的傅里叶级数。 2）.利用MATLAB绘出由前N次谐波合成的信号波形，观察随着N的变化合成信号波形的变化规律。 3）.利用MATLAB绘出周期矩形脉冲信号的频谱，观察参数T和τ变化时对频谱波形的影响。 思考下列问题： ①什么是吉伯斯现象？产生吉伯斯现象的原因是什么？ ②以周期矩形脉冲信号为例，说明周期信号的频谱有什么特点。 ③周期矩形脉冲信号参数τ/T的变化，其频谱结构（如频谱包络形状、过零点、频谱间隔等）如何变化？ （2）已知x（t）是如图所示矩形脉冲信号。 1）.求该信号的傅里叶变幻。 2）. 利用MATLAB绘出周期矩形脉冲信号的频谱，观察参数T和τ变化时对频谱波形的影响。 3）. 让矩形脉冲宽度始终等于一，改变矩形脉冲宽度，观察矩形脉冲信号时域波形和频谱随矩形脉冲宽度的变化趋势。 ①比较矩形脉冲信号和周期矩形脉冲信号的频谱，两者之间有何异同。 ②让矩形脉冲的面积始终等于一，改变矩形脉冲的宽度，观察矩形脉冲信号时域波形和频谱波形随矩形脉冲宽度的变化趋势。