振动力学实验讲义2013

汽车振动分析试题1

2008年振动力学期末考试试题 第一题（20分） 1、在图示振动系统中，已知：重物C 的质量m 1，匀质杆AB 的质量m 2，长为L ，匀质轮O 的质量m 3，弹簧的刚度系数k 。当AB 杆处于水平时为系统的静平衡位置。试采用能量法求系统微振时的固有频率。 解： 系统可以简化成单自由度振动系统，以重物C 的位移y 作为系统的广义坐标，在静平衡位置时 y ＝0，此时系统的势能为零。 AB 转角：L y /=? 系统动能： m 1动能：2 1121y m T = m 2动能：2222222 22 222)3 1(21))(31(21)31(2121y m L y L m L m J T ====? ω m 3动能：2322 323 33)2 1(21))(21(212 1y m R y R m J T === ω 系统势能： 2 21)21(21)21( y k y g m gy m V + +-= 在理想约束的情况下，系统的主动力为有势力，则系统的机械能守恒，因而有： E y k gy m gy m y m m m V T =+ +-++= +2 212 321) 2 1(2 12 1)2 13 1(2 1 上式求导，得系统的微分方程为： E y m m m k y '=+ + +) 2 131(4321 固有频率和周期为： ) 2 131(43210m m m k + + = ω 2、质量为m 1的匀质圆盘置于粗糙水平面上，轮缘上绕有不可伸长的细绳并通过定滑轮A 连在质量为m 2的物块B 上；轮心C 与刚度系数为k 的水平弹簧相连；不计滑轮A ，绳及弹簧的质量，系统自弹簧原长位置静止释放。试采用能量法求系统的固有频率。 解：系统可以简化成单自由度振动系统，以重物B 的位移x 作为系统的广义坐标，在静平衡位置时 x ＝0，此时系统的势能为零。 物体B 动能：2 212 1x m T = 轮子与地面接触点为速度瞬心，则轮心速度为x v c 2 1= ，角速度为x R 21=ω，转过的角度为x R 21= θ。轮子动能： )83(21)41)(21(21)4 1( 2 12 1212 122 21212 2 12x m x R R m x m J v m T c =+= + = ω 系统势能： x

河流动力学实验

泥沙颗粒分析及沉降速度实验 一、试验目的 1、了解在实验室进行泥沙颗粒分析及沉速实验的一般方法； 2、掌握筛分法和移液管法的适用性及操作过程； 3、掌握泥沙颗粒级配曲线的绘制方法及沙样特征值的确定方法； 4、泥沙沉降现象可结合移液管法的操作过程进行观察。 二、试验方法及适用范围 1、筛分法：适用于泥沙粒径大于0.075mm的颗粒。 2、移液管法：适用于粒径小于0.075mm的颗粒。 3、若沙样中粗细颗粒兼有，则要联合使用筛分法及移液管法。 三、实验方法原理 １、对d>0.1mm的泥沙，应用筛分法测量泥沙颗粒级配。筛分法原理是利用孔径不同、逐级叠置的筛子，通过振动分选，再分别称出各级筛上的沙重，计算绘出沙样的级配曲线； ２、对d<0.1mm的泥沙，应用移液管法测量泥沙颗粒级配。移液管法原理为根据泥沙在沉降筒中沉降快慢的不同，来测定不同粒级的泥沙的数量，通过计算分析绘出沙样的级配曲线。 四、筛分法实验 1、仪器设备：振筛机、烘箱、天平、盛沙杯、沉降筒、温度计、干燥器等。试验筛：粗筛：圆孔孔径为60mm,40mm,20mm,10mm,5mm,2mm;细筛：孔径为2.0,1.0,0.5,0.25,0.1,0.075mm。 天平：称量1000g与称量200g；台秤：称量5kg。 振筛机：符合GB9909-88的技术条件。 其它：盛沙杯、沉降筒、温度计、干燥器等。 五、筛分法实验方法 1、将252.37g沙样放在精密天秤上称重，放入容量瓶的水中称其体积107ml，测出湿密度2.359g/ml。 2、将252.37g沙样放入干燥器中烘干，将烘干后的沙样作为试样放在天平/台

秤上称重为238.4g。（称量准确至0.1g，当沙样质量多于500g时，准确至1g）. 2、将试样倒入依次叠好的最上层筛中，进行筛析。细筛宜放在振筛机上震摇，震摇时间一般为10-15min。 3、由最大孔径筛开始，顺序将各筛取下，在白纸上用手轻叩摇晃，如仍有土粒漏下，应继续轻叩摇晃，至无土粒漏下为止。漏下的土粒应全部放入下级筛内。并将留在各筛上的试样分别称量，准确至0.1g。称量并计算出小于各级筛的泥沙总重量。 4、各细筛上及底盘内土质量总和与筛前所取试样的质量之差0.92%不大于1%。

语音信号处理实验指导书

语音信号处理实验指导书 实验一 语音信号采集与简单处理 一、 实验目的、要求 （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

信号与系统实验指导书

实验一 常用信号分类与观察 一、实验目的 1、了解单片机产生低频信号源； 2、观察常用信号的波形特点及产生方法； 3、学会使用示波器对常用波形参数的测量。 二、实验内容 1、信号的种类相当的多，这里列出了几种典型的信号，便于观察。 2、这些信号可以应用到后面的“基本运算单元”和“无失真传输系统分析”中。 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、抽样信号、钟形信号、脉冲信号等。 1、正弦信号：其表达式为)sin()(θω+=t K t f ，其信号的参数：振幅K 、角频率ω、与初始相位θ。其波形如下图所示： 图 1-5-1 正弦信号 2、指数信号：指数信号可表示为at Ke t f =)(。对于不同的a 取值，其波形表现为不同的形式，如下图所示：

图 1-5-2 指数信号 3、指数衰减正弦信号：其表达式为 ?? ? ??><=-)0()sin()0(0)(t t Ke t t f at ω 其波形如下图： 图 1-5-3 指数衰减正弦信号 4、抽样信号：其表达式为： sin ()t Sa t t = 。)(t Sa 是一个偶函数，t = ±π,±2π,…,±n π时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

河海大学力学08级振动力学结构动力学试卷

一、 1.在单自由度振动系统中，结构振动响应的频率与外加荷载的频率无关（×） 2.在含有阻尼的单自由度振动系统中，结构振动的固有频率与阻尼无关（×） 3.对于图示简支梁，不计梁的质量，分别将物体M 从在距梁中点正上方高H1和H2处 自由释放，H1=2H2,则振动的频率是一样的（√） 二、 1.如图所示，除支撑不同外，其余均相同。（B ） A.图a 振动周期大 B.图b 振动周期大 C.振动周期一样 D.不能判断 2.一物体从高度为h 的地方落下，系统振动频率是（C ） A.h 越大，频率越大 B.h 越大，频率越小 C.与h 无关 D.不能断定 3.对于一个有阻尼的单自由度强迫振动系统来讲，振动响应频率（C ） A.仅由外荷载频率确定 B.仅由系统固有频率确定 C.在系统振动响应一段时间后，仅与外荷载频率有关 D.在系统振动响应一段时间后，仅与系统固有频率有关 4.对于多自由度系统来讲，假设无重频现象，则两个不同的振型φi 和φj 的关系为（C ） A.j T i φφ?一定为零 B.j T i φφ?一定不为零 C.j T i M φφ一定为零 D.j T i K φφ可能不是零 5.对于一个三自由度系统，设某阶段振型为[]T 1,2,1=φ，骑广义质量为4，则其正则振型为（A ）

A.[]T 5,0,1,5.0=φ B.[]T 25.0,5.0,25.0=φ C. []T 1,2,1=φ D.[]T 2,4,2=φ 三、一个单自由度振动系统，自由振动试验测得经过6周后振幅降为原来的1/10，试求阻尼比和在简谐荷载作用下发生共振时的放大系数（15） 解：ξπδm y y m i i y 2ln '==+ m=6 ∴10ln ln 6 =+i i y y ∴0611.06210ln =?=πξ 197.821==ξ μ 四、试写出图示结构的运动方程和位移动力系数（EI 为常数， t F t F θsin )(=） 解：a 12=F a 2 11=F 2____ M 1____M EI a 38311=δ EI a 65312=δ )(16 5)(1112t F t F Fe ==δδ )(16 5t F ky y m =+?? 3 383)2(3a EI a EI k == 383ma EI m k ==ω 211βμ-= EI ma 322θω?β== 32833a m EI EI θμ-= 五、如图所示结构，层间高度均为L ，m1=m2=m ，求系统的固有圆

语音信号处理实验报告

语音信号处理实验 班级： 学号： 姓名：

实验一 基于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 100.5 1 1.5 2 2.5 3x 10 4 024 N=3200.5 1 1.5 2 2.5 3x 10 4 05 N=6400.5 1 1.5 2 2.5 3x 10 4 0510 N=12800.5 1 1.5 2 2.5 3x 10 4 01020 N=2560 0.5 1 1.5 2 2.5 3x 10 4 02040 N=512 （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'); 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 100.5 1 1.5 2 2.5 3x 10 4 012 N=3200.5 1 1.5 2 2.5 3x 10 4 024 N=6400.5 1 1.5 2 2.5 3x 10 4 024 N=12800.5 1 1.5 2 2.5 3x 10 4 0510 N=2560 0.5 1 1.5 2 2.5 3x 10 4 01020 N=512 2） 短时平均过零率 a=wavread('mike.wav'); a=a(:,1); n=length(a); N=320; subplot(3,1,1),plot(a); h=linspace(1,1,N); En=conv(h,a.*a); %求卷积得其短时能量函数En subplot(3,1,2),plot(En); for i=1:n-1 if a(i)>=0 b(i)= 1;

机械行业振动力学期末考试试题(doc-11页)(正式版)

… 2008年振动力学期末考试试题 第一题（20分） 1、在图示振动系统中，已知：重物C 的质量m 1，匀质杆AB 的质量m 2，长为L ，匀质轮O 的质量m 3，弹簧的刚度系数k 。当AB 杆处于水平时为系统的静平衡位置。试采用能量法求系统微振时的固有频率。 解： 系统可以简化成单自由度振动系统，以重物C 的位移y 作为系统的广义坐标，在静平衡位置时 y ＝0，此时系统的势能为零。 AB 转角： 系统动能： % m 1动能： m 2动能： m 3动能： 系统势能： 在理想约束的情况下，系统的主动力为有势力，则系统的机械能守恒，因而有： 上式求导，得系统的微分方程为： E y m m m k y '=+++) 2 1 31(4321 固有频率和周期为： ~ ) 2 131(43210m m m k ++= ω 2、质量为m 1的匀质圆盘置于粗糙水平面上，轮缘上绕有不可伸长的细绳并通过定滑轮A 连在质量为m 2的物块B 上；轮心C 与刚度系数为k 的水平弹簧相连；不计滑轮A ，绳及弹簧的质量，系统自弹簧原长位置静止释放。试采用能量法求系统的固有频率。 解：系统可以简化成单自由度振动系统，以重物B 的位移x 作为系统的广义坐标，在静平衡位置时 x ＝0，此时系统的势能为零。 物体B 动能：2212 1 x m T = 轮子与地面接触点为速度瞬心，则轮心速度为x v c 21=，角速度为x R 21=ω，转过的角度为x R 21 = θ。轮子动能： )83 (21)41)(21(21)41(212121212221212212x m x R R m x m J v m T c =+=+=ω \ x

振动实验报告

振动力学实验报告 学院：___________________ 班级：___________________ 学号：___________________ 姓名：___________________ 山东科技大学

单自由度系统振动实验报告 实验者姓名：________ 院系：_______系_______专业_______班_______组实验日期：________年________月________日 自由振动法测量单自由度系统的参数 一、实验目的 二、实验对象和装置 三、实验步骤 四、实验数据记录和整理 1、无阻尼单自由度自由振动系统实验测量：

计算单自由度振动的振动频率、周期、固有频率、衰减系数、相对阻尼系数周期、频率和阻尼系数： 2、有阻尼单自由度自由振动系统实验测量： 计算单自由度振动的振动频率、周期、固有频率、阻尼系数、相对阻尼系数： 五、简答 1、上述无阻尼自由振动实验中，为什么振动曲线呈现衰减状态？ 2、简述阻尼对于自由振动周期、频率的影响。

用冲击激励法测量系统的频率响应函数 实验者姓名：________ 院系：_______系_______专业_______班_______组实验日期：________年________月________日 一、实验目的 二、实验对象和装置 三．实验步骤

四、实验数据记录和整理 1、无阻尼单自由度自由振动系统实验测量： 2、有阻尼单自由度自由振动系统实验测量： 五、简答 1、力锤施加力的大小是否影响单自由度系统的振动频率和阻尼，为什么？ 2、实验过程中，力锤敲击质量块时应注意什么？

信号与系统实验指导书

信号与系统实验指导手册 沈阳工业大学信息科学与工程学院 2005年10月

前言 “信号与系统”是电子工程、通信工程、信息工程、微电子技术、自动化、计算机等电类相关专业的一门重要的专业基础课，为国内、外各高等院校相关专业的主要课程。由于本课程的理论性、系统性较强，为使抽象的概念和理论形象化、具体化，使学生能够比较深入的理解《信号与系统》课程的基本理论和分析方法，并提高学生分析问题和解决问题的能力。为此，开设了基于本课程的实验。 《信号与系统》实验指导手册，将《信号与系统》课程的理论知识与“信号与系统”的实验系统设备结合，从内容上对教材起到了一定的补充作用，为学生具体实验进行了指导。 鉴于时间仓促，可能会存在一些不足与错误之处，欢迎大家批评指正，使之完善。 编者 2005年10月

目录 实验一系统的特性测试 (1) 实验二信号的采样与恢复 (8) 实验三模拟滤波器分析 (14) 实验四模拟滤波器的设计 (26)

实验一系统的特性测试 一、实验目的 1、学会利用运算单元，搭建一些简单的实验系统。 2、学会测试系统的频率响应的方法。 3、了解二阶系统的阶跃响应特性。 4、学会对其零状态响应和零输入响应进行分析。 二、实验内容 1、根据要求搭建一阶、二阶实验系统。 2、测试一阶、二阶系统的频响特性和阶跃响应。 三、预备知识 学习使用波特图测试系统频响的方法。 四、实验仪器 1、信号与系统实验箱一台（主板）。 2、线性系统综合设计性模块一块。 3、20MHz双踪示波器一台。 五、实验原理 1、基本运算单元 （1）比例放大 1）反相数乘器

由： 2211R U R U -= 则有：1 122R U R U = 2）同相数乘器 由： 544 43R R U R U += 则有：()4 5434R R R U U += （2） 积分微分器 1）积分器： 由：212 11//1R SC U R U -= 则有：()12121 21C SR R R U U +-= 2）微分器 由： 141 31R U SC U -= 则有：S C R U U 1134-= （3） 加法器

振动力学试验

SINOCERA? YE6251振动力学实验系统 实 验 指 导 书 (参考) 江苏联能电子技术有限公司 一、软件安装 运行光盘中YE6251控制软件目录中的setup.exe即可完成控制软件的安装。 用USB接口线连接计算机的USB口与调理器的USB口,打开调理器电源,这时计算机会提示找到新设备,并需要安装驱动程序,这时请指定驱动程序的安装路径为上面安装程序的目录下的 AQU采集器驱动\Win98(或Win2k或WinXP) 二、软件操作 1、登录输入 2、试验项目选择 选择当前的试验项目，包括系统名称、试验项目名称。下面以“简支梁系统”,用“冲击激励法测量模态参数”以例说明具体使用过程。 选择试验项目后，系统会自动显示本实验的实验向导，这样实验时可按照上面的实验向导的

步骤进行实验。 3、通道参数设置 当选择一个试验项目时，系统已经给出一个大致合理的通道设置，当然用户也可进行部分修改。通道参数包括测点号设置，通道是否测量设置,工程单位设置,以及满量程设置(也即通道增益选择)，本试验实际是用第5通道测力锤信号，用第6通道测量加速度信号. “测点号”是描述测点位置的信息。 一般试验与通道号相对应，如通道号1，对应的测点号为1，通道号2对应的测点为2。 对于模态试验时，一般将实验对象划分为若干个测点，如简支梁划分为10等份（两端固定），共需测量9个测点，分别为测点1到测点9（自左向右），测量的方法一般为：如选择测点3位置作为原点，将加速度传感器置于第3点，用力锤依次敲击测点1到测点9，对应的力锤信号命名为f1到f9,对应的加速度信号命名为1到9,这样f1(激励)与1（响应）对应，f9(激励)与9（响应）对应，共测量了9组传递函数. “测量选择”表示左边的对应的通道是否测量。双击可以在测量与不测量之间切换。 “满量程”表示当前通道所通道测量的满度值。 对于位移通道有2档增益,1倍、10倍,对应的位移满量程为5000um、500um。 对于力通道有3档增益,1倍、 10倍、100倍,对于的力的满量程为5000N 、500N、50N。 对于加速度通道4档增益,1倍、 10倍、100倍、1000倍,对于的力的满量程为5000m/s2 、500m/s2、 50m/s2、 5m/s2 注意：当信号调理单元将力通道的增益设置为10倍时，请记住一定要在通道参数栏将该通道的满度值设置为500N（第2档），力通道的增益设置为100倍时，请记住一定要在通道参数栏将该通道的满度值设置为50N（第3档）. 4、系统参数设置

语音信号处理实验报告

通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程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的工作空间。

信号与系统实验指导书版

《信号与系统》实验指导书 中南大学信息科学与工程学院 2014年3月

目录 实验一基本信号的生成 (1) 实验二信号的基本运算 (8) 实验三系统的时域分析 (13) 实验四周期信号的频域分析 (20)

实验一 基本信号的生成 1．实验目的 ● 学会使用MATLAB 产生各种常见的连续时间信号与离散时间信号； ● 通过MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的 理解； ● 熟悉MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠 定基础。 2．实验原理 ⑴ 连续信号的MATLAB 表示 ① 指数信号 指数信号at Ae 在MATLAB 中可用exp 函数表示，其调用形式为 exp()y A a t =** 【例1】 单边衰减指数信号的MATLAB 表示如下：（取1,0.4A a ==-） % program exa_1.m, decaying exponential A=1;a=-0.4; t=0:0.001:10; ft=A*exp(a*t); plot(t,ft) ② 正弦信号 正弦信号0cos()A t ωφ+和0sin()A t ωφ+分别用MATLAB 的内部函数cos 和sin 表示，其调用形式为

cos(0phi) sin(0phi) A w t A w t **+**+ 【例2】 正弦信号的MATLAB 表示如下：（取01,2,/6A ωπφπ===） % program exa_2.m, sinusoidal signal A=1; w0=2*pi; phi=pi/6; t=0:0.001:8; ft=A*sin(w0*t+phi); plot(t,ft) ③ 抽样函数 抽样函数()Sa t 在MATLAB 中用sinc 函数表示，定义为 sinc()sin()/()t t t ππ= 其调用形式为 y=sinc()t 【例3】 抽样函数的MATLAB 表示如下：（取33t ππ=-:） % program exa_3.m, sample function t=-3*pi:pi/100:3*pi; ft=sinc(t/pi); plot(t,ft)

振动力学期末考试试题和答案

振动力学期末考试试题和答案 振动力学(试题) 2008 一、填空(每空2分) 1、设周期振动信号的周期为，则其傅里叶级数的展开的基频为,T ,,, 2、单自由度粘性阻尼系统的阻尼因子与阻尼系数的关系为,,, , 作用下系统响应的稳态振3、单自由度粘性阻尼系统在简谐力ptsin,0 动的幅值为,,, 4、粘性阻尼一周期内所消耗的能量与频率成,,,比。 5、无阻尼多自由度系统的主振型正交关系为,,,,,, 6、写出多自由度系统再频率域的输入与输出之间的关系,,,,, 7、写出瑞利商的表达式,,,,,, r8、多自由度系统中共存在个主固有频率，其相应的主振型,,, 正交。 9、无阻尼多自由度系统，利用里兹法计算出的主振型关于M、K是 否正交,,,,(答是或否) 10、写出如图T-1所示梁的左端边界条件,,,,,,,,,, y L x K 图T-1 二、(20分)系统如图T-2所示，杆AB为刚性、均质，长度为,总L 质量为，弹簧刚度为，阻尼系数为。求系统的固有频率及阻mck

尼因子。 图T-2 三、系统如图T-3所示。求系统的固有频率与主振型。 k k k k k m m m X X X 123 图T-3 四、 五、(20分)简支梁如图T-5所示，弹性模量为E，质量密度为，, 横截面积为A，截面惯性矩为J。求梁在中央受集中弯矩M下的响应。(假设梁的初始状态为零)

图T-5 答案 一、填空(每空2分) 1、周期振动信号的周期为，则其傅里叶级数的展开的基频为 T2/,T 2、单自由度粘性阻尼系统的阻尼因子与阻尼系数的关系为, c ,, 2mk 作用下系统响应的稳态振3、单自由度粘性阻尼系统在简谐力ptsin,0 p10动的幅值为 ,,B222k,,,,，(1)(2) 4、粘性阻尼一周期内所消耗的能量与频率成,正,比。 5、无阻尼多自由度系统的主振型正交关系为加权(M,K)正交: 0()ij,0()ij,,,TTTT ,,,,M,K,,,ijijMij(),Kij(),pipi,, 6、写出多自由度系统在频率域的输入与输出之间的关系 21,其中 xHP()()(),,,,HKMiC()(),,,,,， TXKX7、写出瑞利商的表达式 ()RX,TXMX r8、多自由度系统中共存在个重固有频率，其相应的主振型,,加 权(M,K)正交。 MK9、无阻尼多自由度系统，利用里兹法计算出的主振型关于、是

《土动力学实验》

《土动力学实验》课程教学大纲 课程编号：031157 学分：0 周数：0.75周 大纲执笔人：黄茂松大纲审核人：李镜培 说明：《土动力学实验》是033027《土动力学》及031349《土动力学与基础抗震》课程的实验部分，该课程有课号、有学时、无学分。 一、实验性质与目的 《土动力学实验》是土木工程专业的实践环节课程。 其主要教学目的为：使学生了解土的动力特性，掌握地震区场地评价方法，了解砂土液化的基本概念及评价方法和处理措施，掌握基础振动分析方法并能够进行动力基础的设计，掌握地基基础的抗震强度验算方法以及抗震措施，掌握一些基本的实验方法如：地基土动力参数的测试、基础动力测试、桩基础动力检测等。 二、课程面向专业 土木工程专业 三、实验基本要求 1．通过实验掌握振动理论以及波动理论中的有关知识； 2．通过实验掌握波速测试的方法； 3．通过实验掌握激振法测试地基土动力特性参数的操作与数据处理。 四、实验教学基本内容 1．振动模型实验：对一质量－弹簧－阻尼体系激振，采集振动信号，采用质量－弹簧－阻尼体系振动理论分析振动数据。 2．波动模型实验：对一杆件的杆端进行激振，采集振动信号，采用应力波在一维杆件中的传播方程对数据进行处理。 3．波速测试现场实验：采用跨孔法测试地基土的剪切波波速。 4．激振法现场实验：采用竖向强迫激振法测试地基土的抗压刚度系数。 五、实验内容和主要仪器设备与器材配置

六、实验预习和实验报告的要求、考核方式 学生应当在实验以前掌握与实验有关的各方面的知识，如理论方面、操作方面、数据处理方面、工程应用方面。实验结束后编写实验报告，实验报告中要包括上面提到的四方面的内容。考核方式包括以下几个方面：1、实验操作；2、数据处理；3、报告编写。 七、学时分配 八、教材、实验指导书与主要参考书 教材： 《土动力学》，自编教材。 主要参考书： 《土与基础振动》，F.E Richart,Jr.U .S .A。 《土动力学》，Dao U.S.A浙江大学翻译。 实验指导书名称： 《土动力学基础》

大学本科语音信号处理实验讲义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、辐射模型 从声道模型输出的是速度波，而语音信号是声压波。二者倒比称为辐射阻抗，它表征了

信号与系统实验指导书

《信号与系统》实验指导书 科技学院电气与信息工程学院 2014年3月

目录 实验一连续时间信号的MATLAB表示 实验二连续时间LTI系统的时域分析 实验三连续时间LTI系统的频率特性及频域分析实验四连续时间LTI系统的零极点分析 实验五典型离散信号及其MATLAB实现

实验一连时间信号的MATLAB表示 实验目的 1．掌握MATLAB语言的基本操作，学习基本的编程功能； 2．掌握MATLAB产生常用连续时间信号的编程方法； 3．观察并熟悉常用连续时间信号的波形和特性。 实验原理 1. 连续信号MATLAB实现原理 从严格意义上讲，MATLAB数值计算的方法并不能处理连续时间信号。然而，可用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能够被MATLAB处理，并且能较好地近似表示连续信号。 MATLAB提供了大量生成基本信号的函数。比如常用的指数信号、正余弦信号等都是MATLAB 的部函数。为了表示连续时间信号，需定义某一时间或自变量的围和取样时间间隔，然后调用该函数计算这些点的函数值，最后画出其波形图。 2. 典型信号的MATLAB表示 (1) 实指数信号：()at ，式中K，a为实数。在MATLAB中可以利用exp()函数实现， f t Ke 其语句格式为：y=K*exp(a*t)。

(2) 正弦信号：()sin()f t K t ω?=+或()cos()f t K t ω?=+，其中K 为振幅，ω是角频率；φ是初相位。在MATLAB 中可以利用sin()或cos()函数实现，其语句格式为K *sin(ω*t+phi)，K *cos(ω*t+phi)。 (3) 单位阶跃信号：单位阶跃信号定义为() () 10()0 0t t t ε>??=? =0)”产生。MATLAB 表达式“y = ( t >=0)”的含义就是t ≥ 0时y =1，而当t < 0时y =0。 此外，也可在MATLAB 的work 目录下创建uCT 的M 文件，其MATLAB 源程序为 function f=uCT(t) f=(t>=0); 保存为uCT.m 后，就可以调用该函数。 自定义函数 注意，在此定义的uCT 函数是阶跃信号数值表示方法，因此在数值计算中我们将调用uCT 函数。而在MATLAB 的MAPLE 核中，将Heaviside 函数定义为阶跃信号符号表达式，在符号运算过程中，若要调用它必须用sym 定义后，才能实现。 (4) 复指数信号：()()cos()sin()st j t t t f t Ke Ke Ke t jKe t σωσσωω+===+，其中，s =σ+j ω是复变量；σ、ω为实数。在MATLAB 中可以利用exp()函数表示复指数信号，MATLAB 中默认变量i 为虚部单位， 求取实部用real()函数，求取虚部用imag()函数，求取模用abs()函数，求取相角用angle()函数。 (5) 抽样信号:()sin()/Sa t t t =，在MATLAB 中用与Sa(t )类似的sinc (t )函数表示，定义为sin ()sin()/()c t t t ππ=。 (6) 矩形脉冲信号：()() 1/2()0 /2t g t t τττ???在MATLAB 中可用rectpuls()函数产生，其语句 格式为y=rectpuls(t,width)，该函数用于产生一个幅度为1、宽度为width ，且以t=0为对称轴的矩形脉冲信号，width 的默认值为1。 周期性矩形波信号或方波在MATLAB 中可用square 函数产生，其语句格式为

上海交通大学2008年振动力学期末考试试题

上海交通大学2008年振动力学期末考试试题 第一题（20分） 1、在图示振动系统中，已知：重物C的质量m1，匀质杆AB的质量m2，长为L，匀质轮O的质量m3，弹簧的刚度系数k。当AB杆处于水平时为系统的静平衡位置。试采用能量法求系统微振时的固有频率。 解： 系统可以简化成单自由度振动系统，以重物C的位移y作为系统的广义坐标，在静平衡位置时y＝0，此时系统的势能为零。 AB转角： 系统动能： m1动能： m2动能： m3动能： 系统势能： 在理想约束的情况下，系统的主动力为有势力，则系统的机械能守恒，因而 有： 上式求导，得系统的微分方程为：

固有频率和周期为： 2、质量为m1的匀质圆盘置于粗糙水平面上，轮缘上绕有不可伸长的细绳并通过 定滑轮A连在质量为m2的物块B上；轮心C与刚度系数为k的水平弹簧相连；不计滑轮A，绳及弹簧的质量，系统自弹簧原长位置静止释放。试采用能量法求 系统的固有频率。 解：系统可以简化成单自由度振动系统，以重物B的位移x作为系统的广义坐标，在静平衡位置时x＝0，此时系统的势能为零。 物体B动能： 轮子与地面接触点为速度瞬心，则轮心速度为，角速度为，转过的角度为。轮子动能： 系统势能： 在理想约束的情况下，系统的主动力为有势力，则系统的机械能守恒，有：上式求导得系统的运动微分方程：

固有频率为： 第二题（20分） 1、在图示振动系统中，重物质量为m，外壳质量为2m，每个弹簧的刚度系数均为k。设外壳只能沿铅垂方向运动。采用影响系数方法：（1）以x1和x2为广义坐标，建立系统的微分方程；（2）求系统的固有频率。 解： 系统为二自由度系统。 当x1＝1，x2＝0时，有：k11＝2k，k21＝－2k 当x2＝1，x2＝1时，有：k22＝4k，k12＝－2k 因此系统刚度矩阵为： 系统质量矩阵为： 系统动力学方程为： 频率方程为： 解出系统2个固有频率： ，

振动力学参考答案

请打双面 习题与综合训练第一章 2-1一单层房屋结构可简化为题2-1图所示的模型，房顶质量为m，视为一刚性杆；柱子 高h，视为无质量的弹性杆， 其抗弯刚度为EJ。求该房屋 作水平方向振动时的固有 频率。 解：由于两根杆都是弹性的，可以看作是两根相同的弹簧的并联。 等效弹簧系数为k 则 其中为两根杆的静形变量，由材料力学易知 = 则= 设静平衡位置水平向右为正方向，则有 所以固有频率 2-2一均质等直杆，长为 l，重量为W，用两根长h的相同的铅垂线悬挂成水平位置，如题2-2图所示。试写出此杆绕通过重心的铅垂轴作微摆动的振动微分方程，并求出振动固有周期。 解：给杆一个微转角θ θ＝hα 2F＝mg 由动量矩定理： 其中 2-3求题2-3图中系统的固有频率，悬臂梁端点的刚度分别是和，悬臂梁的质量忽略不计。 解：悬臂梁可看成刚度分 别为k1和k3的弹簧，因此，k1 与k2串联，设总刚度为k1ˊ。 k 1 ˊ与k3并联，设总刚度为k2 ˊ。k2ˊ与k4串联，设总刚度 为k。即为 ，， mg kδ =δ δ 3 24 mgh EJ = k3 24EJ h " m x kx =- 3 n 24 mh EJ p= 2 a a h a mg a mg Fa M ml I M I 8 2 2 cos sin 12 1 2 2 - = - ≈ ? - = == = α θ α θ&& 1 2 cos sin≈ ≈ θ α α h l ga p h a mg ml n2 2 2 2 2 3 4 12 1 = = ? +θ θ&& g h a l ga h l p T n 3 π2 3 π2 π2 2 2 = = = 1 k3k 2 1 2 1 1k k k k k + = ' 2 1 2 1 3 2k k k k k k + + = ' 4 2 4 1 2 1 3 2 3 1 4 2 1 4 3 2 4 2 1 k k k k k k k k k k k k k k k k k k k k + + + + + + = θ F sinα 2 θ α F h mg θ F