什么是频率响应函数
动态信号分析仪的一个常见应用是测量机械系统的频率响应函数(FRF)。这也称为网络分析，系统的输入和输出同时测量。通过这些多通道测量，分析仪可以测量系统如何“改变”输入。一个常见的假设是，如果系统是线性的，那么这个“变化”被频率响应函数(FRF)充分描述。事实上，对于线性和稳定的系统，只要知道频率响应函数，就可以预测系统对任何输入的响应。
宽带随机、正弦、阶跃或瞬态信号在测试和测量应用中被广泛地用作激励信号。图1说明了一个激励信号x，可以应用于一个UUT(测试单元)，并生成一个或多个由y表示的响应，输入和输出之间的关系称为传递函数或频率响应函数，由H(y,x)表示。一般来说，传递函数是一个复杂的函数，描述系统如何将输入信号的大小和相位作为激励频率的函数。
在各种激励条件下，对UUT系统的特性进行了实验测量。这些特征包括:频率响应函数(FRF)，通过以下参量描述: 增益频率函数。相位频率函数。共振频率，阻尼因素，总谐波失真，非线性。
利用宽带随机激励的FFT、交叉功率谱法测量频率响应。宽带激励可以是高斯分布的真随机噪声信号，也可以是一个伪随机信号，其振幅分布可以由用户来定义。宽带这一术语可能具有误导性，因为一个好的实现的随机激励信号应该是频带有限的，并由分析频率范围的上限控制。也就是说，激励不应该激发高于测
量仪器所能测量的频率。随机发生器只产生频宽在分析频率范围内随机信号。这也将把激发能量集中在有用的频率范围，以提高测试动态范围。
宽带随机激励的优点是它能在短时间内激发宽频段，因此总测试时间较短。宽带激励的缺点是其频率能量在短时间内广泛传播。每个频率点激发的能量贡献远小于总信号能量(大概是-30到-50dB小于总数)。即使对于频率响应函数(FRF)估计有一个大的平均数字，宽带信号也不能有效地测量UUT的极端动态特性。
扫频正弦测量，优化了每个频率点的测量值。由于激励信号是一个正弦波，在某一时刻其所有的能量都集中在一个频率上，改进了宽带激励中的动态范围不足的缺点。此外，如果频率响应幅值大小下降，响应的跟踪滤波器可以帮助接收到非常小的正弦信号。只要优化每个频率的输入范围，就可以将测量的动态范围扩展到150分贝以上。
频率响应函数的应用很广，其中测试试件的固有频率是基础应用，可以有效的避免共振频率。试件由于材质、材料属性、形状的不同会影响自身刚度和质量。它的固有频率只受刚度分布和质量分布的影响，阻尼对固有频率的影响有限。质量增大固有频率必然降低，刚度增大固有频率必然增大。
理论上讲，试件有多阶固有频率。在二维频谱图中，并不是所有的峰值对应的都是固有频率，因为有可能是激励频率或是它的倍频。因此通常通过测量频响函数的方式来测量固有频率，频响函数对应的峰值都是系统的固有频率。多数情况下，我们只关心低阶或特定阶固有频率。
常用两种方法测试频率响应函数，锤击法和正弦扫频法。
CoCo-80X 动态信号分析仪利用锤击法测试
Spider-81B 振动控制器+振动台利用正弦扫频法测试
杭州锐达数字技术有限公司是美国晶钻仪器公司中国总代理，负责产品销售、技术支持与产品维护，是机械状态监测、振动噪声测试、动态信号分析、动态数据采集、应力应变测试等领域的供应商，提供手持一体化动态信号分析系统、多通道动态数据采集系统、振动控制系统、多轴振动控制系统、三综合试验系统和远程状态监测系统等。
二阶系统的阶跃响应及频率特性
实验二二阶系统的阶跃响应及频率特性 实验简介：通过本实验学生能够学习二阶系统的频率响应和幅频特性的测试方法，对实验装置和仪器的调试操作，具备对实验数据、结果的 处理及其与理论计算分析比较的能力。 适用课程：控制工程基础 实验目的：A 学习运算放大器在控制工程中的应用及传递函数的求取。 B 学习二阶系统阶跃响应曲线的实验测试方法。 C 研究二阶系统的两个重要参数ζ、ω n 对阶跃瞬态响应 指标的影响。 D 学习频率特性的实验测试方法。 E 掌握根据频率响应实验结果绘制Bode图的方法。 F 根据实验结果所绘制的Bode图，分析二阶系统的主要 动态特性（M P ，t s ）。 面向专业：机械类 实验性质：综合性/必做 知 识 点：A《模拟电子技术》课程中运算放大器的相关知识； B《数字电子技术》课程中采样及采样定理的相关知识； C《机械工程控制基础》课程中，传递函数，时域响应， 频率响应三章的内容。 学 时 数：2 设备仪器：XMN-2自动控制原理学习机，CAE-98型微机接口卡，计算机辅助实验系统2.0软件，万用表。 材料消耗：运算放大器，电阻，电容，插接线。 要 求：实验前认真预习实验指导书的实验内容，完成下述项目， 做实验时交于指导教师检查并与实验报告一起记入实验成绩。 B推导图2所示积分放大器的输出输入时域关系和传递函数。
C 推导图3所示加法和积分放大器的输出输入时域关系（两输入单输出） 和S <1>.写出op1,op2,op9,0p6对应的微分方程组（4个方程）。 <2>.画出系统方框图。 <3>.用方框图化简或方程组联立消元的方法求取实验电路所示系统的 传递函数，写出求解过程。 和ζ。 <4>.求取该系统的ω n 实验地点：教一楼327室 实验照片：实验装置及仪器
各种谱计算,频响函数,传递率
各种谱计算，频响函数，传递率 阅读：22802006-05-25 22:01 A．信号与谱的分类 由于时域信号有不同的分类, 变换后对应的频域也有不同的谱 信号可分为模拟(连续)信号和数字(离散)信号, 连续信号变换后称为谱密度, 离散信号变换 后称为谱. 连续信号又可分为绝对可积,平方可积(能量有限),均方可积(功率有限) 绝对可积信号有傅里叶谱(线性谱)和傅里叶谱密度(线性谱密度),如时域信号单位为电压V, 则前者单位为V,后者单位为V/Hz. 均方可积信号有功率谱PS(单位为V2)和功率谱密度PSD(单位为V2/ Hz.). 平方可积信号有能量谱密度ESD(单位为V2 s / Hz.). 注1 平方量称为功率,平方量乘秒称为能量,谱分量除以频率称为谱密度. 注2 功率谱密度另一定义(离散信号的功率谱密度)见下述, 连续信号的功率谱密度. 为连续(光滑)曲线, 离散信号的功率谱密度为不连续的阶梯形.. 注3 随机信号求功率谱密度时为减少方差,可采用平均,重叠和加窗处理(Welch法). 数字信号又可分为绝对可和,平方可和,均方可和. B.各种谱计算 1. 线性谱Linear Spectrum: 对时域离散信号作DFT(离散傅里叶变换)得到, 采用方法为FFT(快速傅里叶变换)法.X(f)=FFT(x(t)) 2. 自功率谱APS=Auto Power Spectrum: 离散信号的线性谱乘其共轭线性谱 APS(f)=X(f)*conj(X(f)), conj=conjugate共轭(实部不变,虚部变符号). 3. 互功率谱CPS=Cross Power Spectrum::x(t)的线性谱乘y(t)的共轭线性谱 互功率谱是复数,可表示为幅值和相位或实部和虚部等. CPS(f)=X(f) *conj(Y(f)) Y(f)=FFT(y(t)) 4. (自)功率谱密度PSD(=Power Spectrum Density): PSD(f)=APS(f)/Δf Δf—频率分辨率(Hz), 自功率谱密度与自相关函数成傅立叶对应关系 故功率谱密度也称为规一化的功率谱. 5. 互功率谱密度CSD=CPS(f)/Δf A.频响函数FRF, 传递率 A1.频响函数.FRF为响应的傅里叶变换与力的傅里叶变换之比或力和响应的互谱与力的自谱之比后者可通过平均减少噪声,故较常用. H(f)=X(f ) / F(f)=X(f)*conj(F(f)) / F(f)*conj(F(f))=CPS / APS. A2. 频响函数有三种表达形式 频响函数表达成分子多项式与分母多项式(特征多项式)之比,也称有理分式. (两多项式求根后) 频响函数表达成极点,零点和增益ZPK形式. 频响函数表达成部分分式,也称极点留数形式,( 部分分式的分子项称为留数.), 例如:最常见的单自由度(位移)频响函数H(ω)=X(ω)/F(ω)
求函数极限的方法
一、求函数极限的方法 1、运用极限的定义 例: 用极限定义证明: 12 23lim 22=-+-→x x x x 证: 由 2 4 4122322-+-= --+-x x x x x x ()2 2 22 -=--= x x x 0>?ε 取εδ= 则当δ <-<20x 时,就有 ε<--+-12 2 32x x x 由函数极限δε -定义有: 12 23lim 22=-+-→x x x x 2、利用极限的四则运算性质 若 A x f x x =→)(lim 0 B x g x x =→)(lim 0 (I)[]=±→)()(lim 0 x g x f x x )(lim x f x x →±B A x g x x ±=→)(lim 0 (II) []B A x g x f x g x f x x x x x x ?=?=?→→→)(lim )(lim )()(lim 0 (III)若 B ≠0 则： B A x g x f x g x f x x x x x x ==→→→)(lim ) (lim )()(lim 0 00 （IV ）cA x f c x f c x x x x =?=?→→)(lim )(lim （c 为常数） 上述性质对于时也同样成立 -∞→+∞→∞→x x x ,,
例：求 4 5 3lim 22+++→x x x x 解: 4 53lim 22+++→x x x x = 25 4252322=++?+ 3、约去零因式（此法适用于型时0 ,0x x →） 例: 求12 16720 16lim 23232+++----→x x x x x x x 解:原式=() () ) 12102(65) 2062(103lim 2 23223 2 +++++--+---→x x x x x x x x x x x =) 65)(2() 103)(2(lim 222+++--+-→x x x x x x x =) 65() 103(lim 222++---→x x x x x =)3)(2()2)(5(lim 2+++--→x x x x x =2 lim -→x 73 5 -=+-x x 4、通分法（适用于∞-∞型） 例: 求 )21 44( lim 2 2 x x x ---→ 解: 原式=) 2()2() 2(4lim 2x x x x -?++-→ =) 2)(2() 2(lim 2x x x x -+-→ =4 1 21lim 2=+→x x 5、利用无穷小量性质法（特别是利用无穷小量与有界量之乘积仍为无穷小量的性质） 设函数f(x)、g(x) 满足：
频响函数用于转子振动信号诊断
A frequency response function-based structural damage identi?cation method Usik Lee *,Jinho Shin Department of Mechanical Engineering,Inha University,253Yonghyun-Dong,Nam-Ku,Incheon 402-751,South Korea Received 9March 2001;accepted 9October 2001 Abstract This paper introduces an frequency response function (FRF)-based structural damage identi?cation method (SDIM)for beam structures.The damages within a beam structure are characterized by introducing a damage distribution function.It is shown that damages may induce the coupling between vibration modes.The e?ects of the damage-induced coupling of vibration modes and the higher vibration modes omitted in the analysis on the accuracy of the predicted vibration characteristics of damaged beams are numerically investigated.In the present SDIM,two feasible strategies are introduced to setup a well-posed damage identi?cation problem.The ?rst strategy is to obtain as many equations as possible from measured FRFs by varying excitation frequency as well as response measurement point.The second strategy is to reduce the domain of problem,which can be realized by the use of reduced-domain method in-troduced in this study.The feasibility of the present SDIM is veri?ed through some numerically simulated damage identi?cation tests.ó2002Elsevier Science Ltd.All rights reserved. Keywords:Structural damage;Damage identi?cation;Beams;Frequency response function;Damage-induced modal coupling;Reduced-domain method 1.Introduction Existence of structural damages within a structure leads to the changes in dynamic characteristics of the structure such as the vibration responses,natural fre-quencies,mode shapes,and the modal dampings.Therefore,the changes in dynamic characteristics of a structure can be used in turn to detect,locate and quantify the structural damages generated within the structure.In the literature,there have been appeared a variety of structural damage identi?cation methods (SDIM),and the extensive reviews on the subject can be found in Refs.[1–3]. The ?nite element model (FEM)update techniques have been proposed in the literature [4–9].As a draw- back of FEM-update techniques,the requirement of reducing FEM degrees of freedom or extending the measured modal parameters may result in the loss of physical interpretability and the errors due to the sti?-ness di?usion that smears the damage-induced localized changes in sti?ness matrix into the entire sti?ness matrix.Thus,various experimental-data-based SDIM have been proposed in the literature as the alternatives to the FEM-update techniques. The experimental-data-based SDIM depends on the type of data used to detect,locate,and/or quantify structural damages.They include the changes in modal data [10–18],the strain energy [19,20],the transfer function parameters [21],the ?exibility matrix [22,23],the residual forces [24,25],the wave characteristics [26],the mechanical impedances [27,28],and the frequency response functions (FRFs)[29–31].Most of existing modal-data-based SDIM have been derived from FEM model-based eigenvalue problems. As discussed by Banks et al.[32],the modal-data-based SDIM have some shortcomings.First,the modal * Corresponding author.Tel.:+82-32-860-7318;fax:+82-32-866-1434. E-mail address:ulee@inha.ac.kr (U.Lee). 0045-7949/02/$-see front matter ó2002Elsevier Science Ltd.All rights reserved.PII:S 0045-7949(01)00170-5
一二阶系统频率特性测试与分析
广西大学实验报告纸 姓名： 指导老师：胡老师 成绩： 学院：电气工程学院 专业：自动化 班级：121 实验内容：零、极点对限性控制系统的影响 2014年 11月 16 日 【实验时间】2014年11月14日 【实验地点】宿舍 【实验目的】 1. 掌握测量典型一阶系统和二阶系统的频率特性曲线的方法； 2. 掌握软件仿真求取一、二阶系统的开环频率特性的方法； 3. 学会用Nyquist 判据判定系统的稳定性。 【实验设备与软件】 1. labACT 实验台与虚拟示波器 2. MATLAB 软件 【实验原理】 1.系统的频率特性测试方法 对于现行定常系统，当输入端加入一个正弦信号)sin()(t X t X m ωω=时，其稳态输出是一个与输入信号频率相同，但幅值和相位都不同的正弦信号 )sin()()sin()(ψωωψω+=+=t j G X t Y s Y m m 。 幅频特性：m m X Y j G /)(=ω，即输入与输出信号的幅度比值，通常转换成)(lg 20ωj G 形式。 相频特性：)(arg )(ωω?j G =，可以直接基于虚拟示波器读取，也可以用“李沙育图行”法得到。 可以将用Bode 图或Nyquist 图表示幅频特性和相频特。 在labACT 试验台采用的测试结构图如下：
被测定稳定系统对于实验就是有源放大电路模拟的一、二阶稳定系统。 2.系统的频率测试硬件原理 1）正弦信号源的产生方法 频率特性测试时，一系列不同频率输入正弦信号可以通过下图示的原理产生。按照某种频率不断变化的数字信号输入到DAC0832，转换成模拟信号，经一级运放将其转换为模拟电压信号，再经过一个运放就可以实现双极性电压输出。 根据数模转换原理，知 R V N V 8 012- = (1) 再根据反相加法器运算方法，得 R R R V N V N V R R V R R V 1281282282201210--=??? ??+-?-=??? ? ??+-= (2) 由表达式可以看出输出时双极性的：当N 大于128时，输出为正；反之则为负；当输入为128时，输出为0. 在labACT 实验箱上使用的参考电压时5V 的，内部程序可以产生频率范围是对一阶系统是0.5 H Z ~64H Z 、对二阶系统是0.5 H Z ~16 H Z 的信号，并由B2单元的OUT2输出。
频谱分析仪的响应函数
什么是频率响应函数 动态信号分析仪的一个常见应用是测量机械系统的频率响应函数(FRF)。这也称为网络分析，系统的输入和输出同时测量。通过这些多通道测量，分析仪可以测量系统如何“改变”输入。一个常见的假设是，如果系统是线性的，那么这个“变化”被频率响应函数(FRF)充分描述。事实上，对于线性和稳定的系统，只要知道频率响应函数，就可以预测系统对任何输入的响应。 宽带随机、正弦、阶跃或瞬态信号在测试和测量应用中被广泛地用作激励信号。图1说明了一个激励信号x，可以应用于一个UUT(测试单元)，并生成一个或多个由y表示的响应，输入和输出之间的关系称为传递函数或频率响应函数，由H(y,x)表示。一般来说，传递函数是一个复杂的函数，描述系统如何将输入信号的大小和相位作为激励频率的函数。 在各种激励条件下，对UUT系统的特性进行了实验测量。这些特征包括:频率响应函数(FRF)，通过以下参量描述: 增益频率函数。相位频率函数。共振频率，阻尼因素，总谐波失真，非线性。 利用宽带随机激励的FFT、交叉功率谱法测量频率响应。宽带激励可以是高斯分布的真随机噪声信号，也可以是一个伪随机信号，其振幅分布可以由用户来
定义。宽带这一术语可能具有误导性，因为一个好的实现的随机激励信号应该是频带有限的，并由分析频率范围的上限控制。也就是说，激励不应该激发高于测量仪器所能测量的频率。随机发生器只产生频宽在分析频率范围内随机信号。这也将把激发能量集中在有用的频率范围，以提高测试动态范围。 宽带随机激励的优点是它能在短时间内激发宽频段，因此总测试时间较短。宽带激励的缺点是其频率能量在短时间内广泛传播。每个频率点激发的能量贡献远小于总信号能量(大概是-30到-50dB小于总数)。即使对于频率响应函数(FRF)估计有一个大的平均数字，宽带信号也不能有效地测量UUT的极端动态特性。 扫频正弦测量，优化了每个频率点的测量值。由于激励信号是一个正弦波，在某一时刻其所有的能量都集中在一个频率上，改进了宽带激励中的动态范围不足的缺点。此外，如果频率响应幅值大小下降，响应的跟踪滤波器可以帮助接收到非常小的正弦信号。只要优化每个频率的输入范围，就可以将测量的动态范围扩展到150分贝以上。 频率响应函数的应用 频率响应函数的应用很广，其中测试试件的固有频率是基础应用，可以有效的避免共振频率。试件由于材质、材料属性、形状的不同会影响自身刚度和质量。它的固有频率只受刚度分布和质量分布的影响，阻尼对固有频率的影响有限。质量增大固有频率必然降低，刚度增大固有频率必然增大。 理论上讲，试件有多阶固有频率。在二维频谱图中，并不是所有的峰值对应的都是固有频率，因为有可能是激励频率或是它的倍频。因此通常通过测量频响
自动控制原理学生实验：二阶开环系统的频率特性曲线
实验三 二阶开环系统的频率特性曲线 一．实验要求 1．研究表征系统稳定程度的相位裕度γ和幅值穿越频率c ω对系统的影响。 2．了解和掌握欠阻尼二阶开环系统中的相位裕度γ和幅值穿越频率c ω的计算。 3．观察和分析欠阻尼二阶开环系统波德图中的相位裕度γ和幅值穿越频率ωc ，与计算值作比对。 二．实验内容及步骤 本实验用于观察和分析二阶开环系统的频率特性曲线。 由于Ⅰ型系统含有一个积分环节，它在开环时响应曲线是发散的，因此欲获得其开环频率特性时，还是需构建成闭环系统，测试其闭环频率特性，然后通过公式换算，获得其开环频率特性。 自然频率：T iT K = n ω 阻尼比：KT Ti 2 1= ξ （3-2-1） 谐振频率： 2 21ξωω-=n r 谐振峰值：2 121lg 20)(ξ ξω-=r L （3-2-2） 计算欠阻尼二阶闭环系统中的幅值穿越频率ωc 、相位裕度γ： 幅值穿越频率： 24241ξξωω-+? =n c （3-2-3） 相位裕度： 4 24122arctan )(180ξξξω?γ++-=+=c （3-2-4） γ值越小，Mp%越大，振荡越厉害；γ值越大，Mp%小，调节时间ts 越长，因此为使 二阶闭环系统不致于振荡太厉害及调节时间太长，一般希望： 30°≤γ≤70° （3-2-5） 本实验所构成的二阶系统符合式（3-2-5）要求。 被测系统模拟电路图的构成如图1所示。 图1 实验电路 本实验将数/模转换器（B2）单元作为信号发生器，自动产生的超低频正弦信号的频率从低到高变化（0.5Hz~16Hz ），OUT2输出施加于被测系统的输入端r (t)，然后分别测量被测系统的输出信号的开环对数幅值和相位，数据经相关运算后在虚拟示波器中显示。 实验步骤： （1）将数/模转换器（B2）输出OUT2作为被测系统的输入。 （2）构造模拟电路：安置短路套及测孔联线表同笫3.2.2 节《二阶闭环系统的频率特性曲线测试》。 （3）运行、观察、记录： ① 将数/模转换器（B2）输出OUT2作为被测系统的输入，运行LABACT 程序，在界面 的自动控制菜单下的线性控制系统的频率响应分析-实验项目，选择二阶系统，就会弹出虚拟示波器的界面，点击开始，实验开始后，实验机将自动产生0.5Hz~16H 等多种频率信号，等待将近十分钟，测试结束后，观察闭环对数幅频、相频曲线和幅相曲线。 ② 待实验机把闭环频率特性测试结束后，再在示波器界面左上角的红色‘开环’或‘闭
复合函数极限条件
书中这样定义： 设函数y = f[g(x)]是由函数u = g(x)与函数y = f(u)复合而成，f[g(x)]在点x0的某去心邻域内有定义，若lim(x->x0)g(x) = u0, lim(u->u0)f(u) = A，且存在δ > 0,当x属于x0的去心δ邻域时，有g(x)不等于u0，则lim(x->x0)f[g(x)] = A u 与u0的接近程度是用0 < |u - u0| < δ描述的，u -> u0的过程中不等于u0 函数在某点的极限值是自变量逼近这一点时函数值无限接近的一个值，这个值与函数在这一点的函数值无关 如果能进一步针对这条举出反例就更好了， g(x)=xsin(1/x) 若u≠0，f(u)=0 若u=0，f(u)=1 在0的去心邻域中，f(g(x))有定义 (*) 对任意的正数δ，在0的去心δ邻域中，都有无数个点使得g(x)=0, 而f(g(x))=f(0)=1 lim{x→0}g(x)=0 lim{u→0}f(u)=0 而根据(*)，lim{x→0}f(g(x))不存在。 可见这个条件确实不能去掉。如果f(u)在u0处连续，那么这个复合函数的极限运算法则仍然是成立的，g(x)是否在其他点取值u0并无影响，因而很多时候在实际应用这条法则时并不去验证这条，因为我们通常面对的是连续函数。确实是这样的，因为g(x)在0的任意去心邻域内总是存在使得g(x)为0的点，而f(0) = 1 =/= lim(u->0)f(u)。所以就不存在0的某个去心邻域使得|f(g(x))-0|能够小于任意ε>0，自然极限也就不存在了。 另一种情况：设lim(u->u0)f(u) = A，且f(u)在u0的某个去心邻域是连续函数，那么就有f(u0) = lim(u->u0)f(u) = A，再设lim(x->x0)g(x) = u0，那这时候就不用考虑在x0的某个去心邻域中，g(x) =/= u0这个条件了，因为g(x) =u0时，|f(g(x)) - A| = 0 < 任意ε>0 。
离散系统的频率响应分析
实验报告 实验课程：数字信号处理 实验内容：实验4 离散系统的频率响应分析和零、极点分布 院（系）：计算机学院 专业：通信工程 班级：111班 2013年6 月7日
一、实验目的：加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验原理： 离散系统的时域方程为 ∑∑==-=-M k k N k k k n x p k n y d ) ()( 其变换域分析方法如下： 时频域变换 )()()(][][][][][ω ωωj j j m e H e X e Y m n h m x n h n x n y =?-= *=∑ ∞ -∞ = 系统的频率响应为 ω ωω ωω ωω jN N j jM M j j j j e d e d d e p e p p e D e p e H ----++++++==......)()()(1010 时域Z 域变换 ) ()()(][][][][][z H z X z Y m n h m x n h n x n y m =?-= *=∑∞ -∞ = 系统的转移函数为 N N M M z d z d d z p z p p z D z p z H ----++++++= =......)()()(110110 分解因式 ∏-∏-=∑∑= =-=-=-=-N i i M i i N i i k M i i k z z K z d z p z H 1111 0)1()1()(λξ 上式中的i ξ和i λ称为零、极点。 在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统转移函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的级联。
函数极限的运算法则
教学目标：掌握函数极限的运算法则，并会求简单的函数的极限 教学重点：运用函数极限的运算法则求极限 教学难点：函数极限法则的运用 教学过程： 一、引入： 一些简单函数可从变化趋势找出它们的极限，如o x x x x x x o ==→∞ →lim ,01lim .若求极限的函 数比较复杂，就要分析已知函数是由哪些简单函数经过怎样的运算结合而成的，已知函数 二 0）. 说明：当三 例1 求)3(lim 2 2 x x x +→ 例2 求1 1 2lim 2 31 ++-→x x x x 例3 求4 16lim 2 4 --→x x x
分析：当4→x 时，分母的极限是0，不能直接运用上面的极限运用法则.注意函数4 16 2 --= x x y 在定义域4≠x 内，可以将分子、分母约去公因式4-x 后变成4+x ，由此即 可求出函数的极限. 例4 求1 33lim 22 ++-∞ →x x x x 分析：当∞→x 时，分子、分母都没有极限，不能直接运用上面的商的极限运算法则.如果分子、2 总结：lim x x o →lim x ∞ →例5 求lim ∞ →x 分析：同例计算了。 四 （1）lim 2 1 → x （3）lim 4 →x 1 432 1 -+→x x x （5）1 1lim 2 1 +--→x x x （6）9 65lim 2 2 3 -+-→x x x x （7）1 3322lim 2 3 2 +--+∞ →x x x x x （8）5 2lim 3 2 --∞ →y y y y
五 小结 1 有限个函数的和（或积）的极限等于这些函数的和（或积）； 2 函数的运算法则成立的前提条件是函数 )(),(x g x f 的极限存在，在进行极限运算时， 要特别注意这一点. 3 两个（或几个）函数的极限至少有一个不存在时，他们的和、差、积、商的极限不一定不存在. 4 在求几个函数的和（或积）的极限时，一般要化简，再求极限. 六 作业（求下列极限） （1） lim -→x 2 （4）lim 0 →x （7）lim 2 →x （10）x → （13）1 3lim 2 4 3 +++∞ →x x x x x （14）2 3 3 2 )2 312( lim -+→x x x （15）3 526113lim 2 2 1 --+-→x x x x x （16） 3 526113lim 22 --+-∞ →x x x x x （17） 3 2 320 3526lim x x x x x x x ----→ （18） 3 2 323526lim x x x x x x x ----∞ →
实验三 二阶系统频率响应
实验三 二阶系统频率响应 一、实验目的 （1）学习系统频率特性响应的实验测试方法。 （2）了解二阶闭环系统中的对数幅频特性和相频特性的计算。 （3）掌握根据频率响应实验结果绘制波特图的方法。 （4）掌握欠阻尼二阶闭环系统中的自然频率、阻尼比对谐振频率、谐振峰值和带宽的影响及对应的计算。 二、实验设备 （1）XMN-2型学习机； （2）CAE-USE 辅助实验系统 （3）万用表 （4）计算机 三、实验内容 本实验用于观察和分析二阶系统瞬态响应的稳定性。 二阶闭环系统模拟电路如图3-1所示，它由两个积分环节（OP1和OP2）及其反馈回路构成。 图3-1 二阶闭环系统模拟电路图 OP1和OP2为两个积分环节，传递函数为s T s G i 1 )(-=（时间常数RC T i =）。二阶闭环系统等效结构图如图3-2所示。 图3-2 二阶闭环系统等效结构图 该二阶系统的自然振荡角频率为RC T n 11==ω，阻尼为i f R R K 22= =ζ。 四、实验步骤 （1）调整Rf=40K ，使K=0.4（即ζ=0.2）；取R=1M ，C=1μ，使T=1秒（ωn=1/1）。 （2）输入信号位)sin(t X ω=，改变角频率使ω分别为 0.2,0.6,0.8,0.9,1.0,1.2,1.6,2.0,3.0rad/s 。稳态时，记录下输出响应)sin(φω+=t Y y 五、数据采集及处理 输出信号幅值Y 输出信号初相φ L(ω) φ(ω) ω(rad/s) T 0.2 0.6 0.8 0.9 1.0 1.2
1.6 2.0 3.0 六、实验报告 1、绘制系统结构图，并求出系统传递函数，写出其频率特性表达式。 2、用坐标纸画出二阶闭环系统的对数幅频、相频曲线（波特图）。 3、其波特图上分别标示出谐振峰值（Mr）、谐振频率（ωr）和带宽频率（ωb）。 4、观察和分析曲线中的谐振频率（ωr）、谐振峰值（Mr）和带宽（ωb），并与理论计算值作对比。
函数极限的十种求法
函数极限的十种求法 信科2班江星雨20140202250 函数极限可以分成而运用ε-δ定义更多的见诸于已知极限值的证明题中。掌握这类证明对初学者深刻理解运用极限定义大有裨益。以的极限为例，f(x) 在点以A为极限的定义是：对于任意给定的正数ε（无论它多么小），总存在正数，使得当x满足不等式时，对应的f(x)函数值都满足不等式：，那么常数A就叫做函数f(x)当x→x。时的极限。 1.利用极限的四则运算法则： 极限四则运算法则的条件是充分而非必要的，因此，利用极限四则运算法则求函数极限时，必须对所给的函数逐一进行验证它是否满足极限四则运算法则条件，满足条件者。方能利用极限四则运算法则进行求之。不满足条件者，不能直接利用极限四则运算法则求之。但是，井非不满足极限四则运算法则条件的函数就没有极限，而是需将函数进行恒等变形，使其符合条件后，再利用极限四则运算法则求之。而对函数进行恒等变形时，通常运用一些技巧如拆项、分子分母同时约去零因子、分子分母有理化、通分、变量替换等等。例 1 求lim( x 2 ? 3x + 5). x→ 2 解：lim( x 2 ? 3x + 5) = lim x 2 ? lim 3x + lim 5 = (lim x) 2 ? 3 lim x + lim 5 = 2 2 ? 3 ? 2 + 5 = 3. x→2 x →2 x →2 x →2 x →2 x →2 x →2 2.利用洛必达法则 洛必达(L 'Hopital)法则是在一定条件下通过分子分母分别求导再求极限来确定未定式值的方法.简单讲就是，在求一个含分式的函数的极限时，分别对分子和分母求导，在求极限，和原函数的极限是一样的。一般用在求导后为零比零或无穷比无穷的类型。 利用洛必达求极限应注意以下几点： 设函数f(x)和F(x)满足下列条件： （1）x→a时,lim f(x)=0,lim F(x)=0; （2）在点a的某去心邻域内f(x)与F(x)都可导,且F(x)的导数不等于0; （3）x→a时,lim(f'(x)/F'(x))存在或为无穷大 则x→a时,lim(f(x)/F(x))=lim(f'(x)/F'(x)) 例1： 1-cosx = 1-{1-2[sin(x/2)]^2} = 2[sin(x/2)]^2 xsinx = 2xsin(x/2)cos(x/2) 原式= lim 2[sin(x/2)]^2 / [2xsin(x/2)cos(x/2)] = tgx / x 对分子分母同时求导（洛必达法则） (tgx)' = 1 / (cosx)^2 (x)' = 1 原式= lim 1/(cosx)^2 当x --> 0 时,cosx ---> 1 原式= 1 3.利用两个重要极限： 应用第一重要极限时，必须同时满足两个条件： ①分子、分母为无穷小，即极限为0 ； ②分子上取正弦的角必须与分母一样。 应用第二重要极限时，必须同时满足四个条件：
线性控制系统的频率响应分析
一．实验目的 1．了解和掌握对数幅频曲线和相频曲线（波德图）、幅相曲线（奈奎斯特图）的构造及绘制方法。 2．二阶开环系统中的相位裕度和幅值穿越频率的计算。 二．实验内容及要求 1．一阶惯性环节的频率特性曲线测试。 2．二阶开环系统的频率特性测试，研究表征系统稳定程度的相位裕度和 幅值穿越频率对系统的影响。 三、实验主要仪器设备和材料 1．labACT自控/计控原理实验机一台 2．数字存储示波器一台 四、实验方法、步骤及结果测试 1．一阶惯性环节的频率特性曲线 惯性环节的频率特性测试模拟电路见图4-1。 图4-1 惯性环节的频率特性测试模拟电路 实验步骤：注：‘S ST'不能用“短路套”短接！ （1）将数/模转换器（B2）输出OUT2作为被测系统的输入。 （2）按图4-1安置短路套及测孔联线。 （3）运行、观察、记录： ①运行LABACT程序，选择自动控制菜单下的线性控制系统的频率响应分析-实验项目，选择一阶系统，再选择开始实验，点击开始，实验机将自动产生0.5Hz~64Hz多个频率信号，测试被测系统的频率特性，等待将近十分钟，测试结束。 ②测试结束后，可点击界面下方的“频率特性”选择框中的任意一项进行切换，将显示被测系统的对数幅频、相频特性曲线（伯德图）和幅相曲线（奈 奎斯特图），同时在界面上方将显示点取的频率点的L、、Im、Re等相关数
据。如点击停止，将停止示波器运行，不能再测量数据。 ③分别改变惯性环节开环增益与时间常数，观察被测系统的开环对数幅频曲线、相频曲线及幅相曲线，在幅频曲线或相频曲线上点取相同的频率点，测量、记录数据于实验数据表中。 实验数据表1：改变惯性环节开环增益，（T=0.05，C=1u，R2=50K） 实验数据表2： 改变惯性环节时间常数， K=1（R1=50K、R2=50K） 2．二阶开环系统的频率特性曲线 二阶系统模拟电路图的构成如图4-2所示。
一、二阶系统频率特性测试与分析
自动频率序列正弦信号产生单元 广西大学实验报告纸 姓名： 指导老师： 成绩： 学院：电气工程学院 专业：自动化 班级： 实验内容：一、二阶系统频率特性测试与分析 年月日 其他组员及各自发挥的作用： 【实验时间】 【实验地点】 【实验目的】 1. 掌握测量典型一阶系统和二阶系统的频率特性曲线的方法； 2. 掌握软件仿真求取一、二阶系统的开环频率特性的方法； 3. 学会用Nyquist 判据判定系统的稳定性。 【实验设备与软件】 1. labACT 实验台与虚拟示波器 2. MATLAB 软件 【实验原理】 1.系统的频率特性测试方法 对于现行定常系统，当输入端加入一个正弦信号)sin()(t X t X m ωω=时，其稳态输出是一个与输入信号频率相同，但幅值和相位都不同的正弦信号)sin()()sin()(ψωωψω+=+=t j G X t Y s Y m m 。 幅频特性：m m X Y j G /)(=ω，即输入与输出信号的幅度比值，通常转换成)(lg 20ωj G 形式。 相频特性：)(arg )(ωω?j G =，可以直接基于虚拟示波器读取，也可以用“李沙育图行”法得到。 可以将用Bode 图或Nyquist 图表示幅频特性和相频特。 在labACT 试验台采用的测试结构图如下。 图1 频率特性测试结构图 被测定稳定系统对于实验就是有源放大电路模拟的一、二阶稳定系统。 2.系统的频率测试硬件原理 被测稳定系统 频率特性测试 模块单元
1）正弦信号源的产生方法 频率特性测试时，一系列不同频率输入正弦信号可以通过下图示的原理产生。按照某种频率不断变化的数字信号输入到DAC0832，转换成模拟信号，经一级运放将其转换为模拟电压信号，再经过一个运放就可以实现双极性电压输出。 根据数模转换原理，知 R V N V 8012 - = (1) 再根据反相加法器运算方法，得 R R R V N V N V R R V R R V 1281282282201210--=??? ??+-?-=??? ? ??+-= (2) 由表达式可以看出输出时双极性的：当N 大于128时，输出为正；反之则为负；当输入为128时，输出为0. 在labACT 实验箱上使用的参考电压时5V 的，内部程序可以产生频率范围是对一阶系统是0.5 H Z ~64H Z 、对二阶系统是0.5 H Z ~16 H Z 的信号，并由B2单元的OUT2输出。 2）被测对象输出信号的采样方法 对被测对象的输出信号夏阳，首先将其通过LM324与基准电压进行比较嵌位，再通过CD14538进行脉冲整形，一保证有足够的IRQ 采样时间，最后将信号送到处理器的IRQ6脚，向处理器申请中断，在中断中对模拟量V y 进行采样并模数转换，进而进行处理与计算幅值与相位。途中采用ADC089采集模拟量，以单极性方式使用，所以在出现振荡的情况下需要加入一个二极管，将V y 出现负值时将其直接拉倒0。
什么是频率响应函数
动态信号分析仪的一个常见应用是测量机械系统的频率响应函数(FRF)。这也称为网络分析，系统的输入和输出同时测量。通过这些多通道测量，分析仪可以测量系统如何“改变”输入。一个常见的假设是，如果系统是线性的，那么这个“变化”被频率响应函数(FRF)充分描述。事实上，对于线性和稳定的系统，只要知道频率响应函数，就可以预测系统对任何输入的响应。 宽带随机、正弦、阶跃或瞬态信号在测试和测量应用中被广泛地用作激励信号。图1说明了一个激励信号x，可以应用于一个UUT(测试单元)，并生成一个或多个由y表示的响应，输入和输出之间的关系称为传递函数或频率响应函数，由H(y,x)表示。一般来说，传递函数是一个复杂的函数，描述系统如何将输入信号的大小和相位作为激励频率的函数。 在各种激励条件下，对UUT系统的特性进行了实验测量。这些特征包括:频率响应函数(FRF)，通过以下参量描述: 增益频率函数。相位频率函数。共振频率，阻尼因素，总谐波失真，非线性。 利用宽带随机激励的FFT、交叉功率谱法测量频率响应。宽带激励可以是高斯分布的真随机噪声信号，也可以是一个伪随机信号，其振幅分布可以由用户来定义。宽带这一术语可能具有误导性，因为一个好的实现的随机激励信号应该是频带有限的，并由分析频率范围的上限控制。也就是说，激励不应该激发高于测
量仪器所能测量的频率。随机发生器只产生频宽在分析频率范围内随机信号。这也将把激发能量集中在有用的频率范围，以提高测试动态范围。 宽带随机激励的优点是它能在短时间内激发宽频段，因此总测试时间较短。宽带激励的缺点是其频率能量在短时间内广泛传播。每个频率点激发的能量贡献远小于总信号能量(大概是-30到-50dB小于总数)。即使对于频率响应函数(FRF)估计有一个大的平均数字，宽带信号也不能有效地测量UUT的极端动态特性。 扫频正弦测量，优化了每个频率点的测量值。由于激励信号是一个正弦波，在某一时刻其所有的能量都集中在一个频率上，改进了宽带激励中的动态范围不足的缺点。此外，如果频率响应幅值大小下降，响应的跟踪滤波器可以帮助接收到非常小的正弦信号。只要优化每个频率的输入范围，就可以将测量的动态范围扩展到150分贝以上。 频率响应函数的应用很广，其中测试试件的固有频率是基础应用，可以有效的避免共振频率。试件由于材质、材料属性、形状的不同会影响自身刚度和质量。它的固有频率只受刚度分布和质量分布的影响，阻尼对固有频率的影响有限。质量增大固有频率必然降低，刚度增大固有频率必然增大。 理论上讲，试件有多阶固有频率。在二维频谱图中，并不是所有的峰值对应的都是固有频率，因为有可能是激励频率或是它的倍频。因此通常通过测量频响函数的方式来测量固有频率，频响函数对应的峰值都是系统的固有频率。多数情况下，我们只关心低阶或特定阶固有频率。 常用两种方法测试频率响应函数，锤击法和正弦扫频法。
深圳大学一阶、二阶系统的幅频特性测试实验
深圳大学实验报告 课程名称：信号与系统 实验项目名称：一阶、二阶系统的幅频特性测试实验学院：信息工程 专业：通信工程 指导教师： 报告人：学号：班级： 实验时间： 实验报告提交时间： 2015.6.23 教务部制
一、实验目的与要求： 1、学会利用基本的运算电路单元，搭建一些简单的实验系统。 2、学会测试系统的频率响应的方法。 3、了解一阶、二阶系统的阶跃响应特性。 二、实验仪器 1、信号与系统实验箱一台（主板）。 2、线性系统综合设计性模块一块。 3、20M双踪示波器一台。 三、实验原理 1、基本运算单元 （1）比例放大 1）反相数乘器 由: 2 2 1 1 R U R U - =则有： 1 1 2 2R U R U- = 2）同相数乘器 由： 5 4 4 4 3 R R U R U + =则有： 4 5 4 3 4 ) ( R R R U U + = （2）积分微分器 1）积分器：由：12 1 2 1 1 // U U R R sC =-则有：2 21 121 (1) R U U R sR C =- +
2）微分器：由：34 1 1 1 U U R sC =-则有： 4311 U U R C s =- （3）加法器 1）反向加法器 有：) ( 2 1 1 3 2R U R U R U+ - = 2）正向加法器 由： ? ? ? ?? ? ? + = - = + - + 8 7 5 7 6 4 4 3 3 R R U R U R U R U R U 则有) ( ) (* 4 4 3 3 7 8 7 6 5R U R U R R R R U+ + = 2、N阶系统系统 1 011 1 1 011 1 ()()()() ()()()() n n n n n n m m m m m m d d d C y t C y t C y t C y t dt dt dt d d d E x t E x t E x t E x t dt dt dt - - - - - - ++++= ++++ 根据零状态响应（起始状态为零），则对其进行拉氏变换有： 1 011 1 011 ()()()() ()()()() n n n n m m m m C s Y s C s Y s C sY s C Y s E s X s E s X s E sX s E X s - - - - ++++= ++++ 则其传函数可表达为： -1 01-1 -1 01-1 s s s (s) (s) (s)s s s m m m m n n n n E E E E Y H X C C C C ++++ == ++++ 3、作为一阶系统，一般可表达为： 01 01 () E s E H s C s C + = + 一阶系统是构成复杂系统的基本单元，学习一阶的特点有助于对一般系统特性的了