时频域分析

时频域分析
时频域分析

时频域分析有关知识

问题

1 现在咱们做的这个数字舵机需要得到的阶跃响应的参数有:超调量(%),时间常数,上升时间,调节时间,稳态误差(°),你给我详细说说那几个量怎么算,

2 数字舵机还要求有频域特性测试功能,跟用户沟通时,客户的意思是做一个扫频功能,发出频率1HZ~60HZ,每个频率求一个最大值,画出频率/幅值曲线,我不知道这个究竟有什么用。

另外幅值裕度,相位裕度如何进行测量?

控制系统性能的评价分为动态性能指标和稳态性能指标。上述两个问题是分析控制系统的动态性能和稳态性能的两类方法,一是时域分析法,二是频域分析法。

1 时域分析法

时域分析法是一种直接在时间域中对系统进行分析的方法。为了求解系统的时间响应,必须了解输入信号的解析表达式。在一般情况下,控制系统的外加信号具有随机性而无法预先确定,需要选择典型信号输入。

工程中常见信号有:单位阶跃函数、单位斜坡函数、单位加速度函数、单位脉冲函数和正弦函数。在实际应用中一般选取最不利的信号作为系统的典型信号。

1.1 动态性能性能指标

通常在阶跃函数作

用下,测定获计算系统

的动态性能。一般认为

阶跃输入对系统来说是

最严峻的工作状态,如

果系统在阶跃函数作用

下的动态性能满足要求,

那么系统在其他形式的

函数作用下,其动态性

能也是令人满意的。

图1单位阶跃响应曲线

以下是稳定系统在单位阶跃函数作用下,动态过程随时间变化状况的动态性能指标,图1是单位阶跃响应曲线。

其动态性能指标通常为:

峰值时间tp,指响应超过其终值到达第一个峰值所需要的时间;

延迟时间td,指响应曲线第一次达到其终值一半所需要的时间;

上升时间tr,指响应从终值10%上升到90%所需要的时间;对于有振荡的系统,亦可定义为响应从0第一次上升到终值所需要的时间。上升时间是系统响应速度的一种度量。上升时间越短,系统响应速度越快。

调节时间ts,指响应到达并保持在终值±5%内所需要的最短时间。有时也用终值的±2%误差范围定义调节时间。

超调量σ%,指响应的最大偏离量与终值之差的百分比,即

σ%=h t p?h(∞)

h(∞)

×100%

若h(t)<h(∞),则响应无超调。超调量亦称为最大超调量,或百分比超调量。在实际应用中,常用的动态性能指标为上升时间、调节时间和超调量。

1.2 稳态误差

稳态误差是描述系统稳态性能的一种指标,在阶跃函数、斜坡函数或加速度函数作用下进行测量或计算。若时间趋于无穷时,系统的输出量不等于输入量的确定函数,则系统存在稳态误差,稳态误差是系统控制精度或扰动能力的一种度量。

稳态误差的计算部分参考《自动控制原理》部分。

1.3 时域动态指标的测量

一般是控制系统发出阶跃信号(一般用脉冲信号代替),同时记录被测产品的输出信号,用波形测量中的“Transition measurement .VI”其帮助如图2所示。

上升时间可以由VI直接测量出来,调节时间和超调量要根据定义进行数据出来,在带入公式计算。

2 频域分析法

频域分析法是应用频率特性研究线性控制系统的另外一种方法,频域分析法的

典型输入信号是扫频信号。扫频信号是指在一定时间内,一规则信号(正弦、锯齿等信号)的频率从低到高地连续变化,其幅值一般保持不变。实际上,就是被测产品的频率响应。

图2 Transition Measurement. VI Help

2.1 频率特性的基本概念

以RC网络为例(如图3)说明频率特性的概念。网络的微分方程为

R

图3 RC网络

T de c

dt

+e c=e r

式中 T=RC。网络的传递函数

E c(s) r =

1

s

1?1

若网络输入为正弦电压为正弦电压,即

e r=Asinωt 则由式1-1得

E c s=1

E r s=

1

?

22

经拉氏变换,得电容两端的输出电压

e c=

AωT

1+ω2T2

e?

1

T+

A

1+ω2T2

?(ωt?arctgωT)

式中,第一项为输出电压的瞬态分量;第二项为稳态分量。随着时间趋于无穷,第一项趋于0,于是

lim t→∞e c=

A

1+ω2T2

?(ωt?arctgωT)

由上式可见,网络的稳定输出仍然是正弦电压,其频率和输入电压的频率相同,幅值是输入的1/2T2倍,相角比输入迟后arctgωT。显然,1/1+ω2T2和-arctgωT皆为输入电压频率ω的函数。前者称为RC网络的频率特性,后者称为RC网络的相频特性。表1是两个特性的计算数据,以ω为横坐标就可画出幅频和相频特性曲线。

1/1+ωT

-arctgωT(°) 0 -26.6 -45 -63.5 -71.5 -76 -78。7 -90

2.2 稳定裕度

奈奎斯特稳定判据是根据开环幅相曲线判断闭环系统稳定性的一种准则,根据奈氏判据知,系统开环幅相曲线临界点附近的形状,对闭环稳定性影响大,曲线越接近临界点,系统的稳定程度就越差。若幅相曲线通过临界点,则稳态后系统的输出端会出现等幅振荡信号。用相角裕度(相位裕度)γ和幅值裕度h表征开环幅相曲线接近临界点的程度,其定义分别如下。

2.2.1 幅值裕度

幅值裕度h定义为幅相曲线上,相角为-180°时对应幅值的倒数,即

h=

1

|G(jωg)H(jωg)|

式中,ω

g

为相角交界频率。

幅值裕度h的含义:如果系统的开环传递函数系统增大到原来的h倍,则系统就处于临界稳定状态。但是,仅仅用幅值裕度不足以表示所有系统的稳定程度,

还要引入相角裕度。

图4 信号传递

2.2.2 相角裕度

相角裕度定义为180°加开环幅相曲线幅值1时的相角,即

γ=180°+∠G jωc H(jωc)

称为系统的截至频率。

式中ω

c

相角裕度的含义:如果系统对频率ω

延迟γ度,则系统将处于临界稳定状

c

态。

2.2.3 带宽频率和带宽

称为闭环幅频特性下降到频率为0时的分贝值以下3dB时,对应的频率ω

b

带宽频率,即

20lg|Φ(jω)|<20lg|Φ(j0)|-3 (ω>ω

b

带宽:系统幅频特性不低于-3dB对应的频率范围

0≤ω≤ω

b

称为系统的带宽。

2.2.4 实际应用

用扫频信号给被测产品一个激励,同时记录被测产品的输出量,一般是输出量随着频率增加而下降。计算出不同频率时输出量相对于某起始频率的输出量之比A,以20lgA为纵坐标,频率为横坐标,绘制幅频特线曲线,其幅值比为0.707时的频率为幅频宽;不同频率相位相对于输入信号的相位差称为相频特性,也绘制在幅频特线曲线上(双坐标)就可以得到相频特性曲线。

信号的相位和频率可以用波形测量.VI测量,不过测量误差比较大。

比较理想的是激励源用信号发生器,波形测量数字化信号仪,这样信号生成和测量比较准确。

大作业1(机电控制系统时域频域分析)

《机电系统控制基础》大作业一 基于MATLAB的机电控制系统响应分析 哈尔滨工业大学 2013年11月4日

1 作业题目 1. 用MATLAB 绘制系统2 ()25()() 425 C s s R s s s Φ== ++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2 ()25 ()()425 C s s R s s s Φ==++的单位阶跃响应性能指标:上升时间、峰值时间、调节时间和超调量。 3. 数控直线运动工作平台位置控制示意图如下: X i 伺服电机原理图如下: L R (1)假定电动机转子轴上的转动惯量为J 1,减速器输出轴上的转动惯量为J 2,减速器减速比为i ,滚珠丝杠的螺距为P ,试计算折算到电机主轴上的总的转动惯量J ; (2)假定工作台质量m ,给定环节的传递函数为K a ,放大环节的传递函数为K b ,包括检测装置在内的反馈环节传递函数为K c ,电动机的反电势常数为K d ,电动机的电磁力矩常数为K m ,试建立该数控直线工作平台的数学模型,画出其控制系统框图; (3)忽略电感L 时,令参数K a =K c =K d =R=J=1,K m =10,P/i =4π,利用MATLAB 分析kb 的取值对于系统的性能的影响。

2 题目1 单位脉冲响应曲线 单位阶跃响应曲线

源代码 t=[0:0.01:1.6]; %仿真时间区段和输入 nC=[25]; dR=[1,4,25]; fi=tf(nC,dR); %求系统模型 [y1,T]=impulse(fi,t); [y2,T]=step(fi,t); %系统响应 plot(T,y1); xlabel('t(sec)'),ylabel('x(t)'); grid on; plot(T,y2); xlabel('t(sec)'),ylabel('x(t)'); grid on; %生成图形 3 题目2 借助Matlab,可得: ans = 0.4330 0.6860 25.3826 1.0000 即

第5章频域分析法习题解答

第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 的信息。

控制系统的频域分析实验报告

实验名称: 控制系统的频域分析 实验类型:________________同组学生姓名:__________ 一、实验目的和要求 用计算机辅助分析的方法,掌握频率分析法的三种方法,即Bode 图、Nyquist 曲线、Nichols 图。 二、实验内容和原理 (一)实验原理 1.Bode(波特)图 设已知系统的传递函数模型: 1 1211121)(+-+-+???+++???++=n n n m m m a s a s a b s b s b s H 则系统的频率响应可直接求出: 1 1211121)()()()()(+-+-+???+++???++=n n n m m m a j a j a b j b j b j H ωωωωω MATLAB 中,可利用bode 和dbode 绘制连续和离散系统的Bode 图。 2.Nyquist(奈奎斯特)曲线 Nyquist 曲线是根据开环频率特性在复平面上绘制幅相轨迹,根据开环的Nyquist 线,可判断闭环系统的稳定性。 反馈控制系统稳定的充要条件是,Nyquist 曲线按逆时针包围临界点(-1,j0)p 圈,为开环传递函数位于右半s 一平面的极点数。在MATLAB 中,可利用函数nyquist 和dnyquist 绘出连续和离散系统的乃氏曲线。 3.Nicho1s(尼柯尔斯)图 根据闭环频率特性的幅值和相位可作出Nichols 图,从而可直接得到闭环系统的频率特性。在 MATLAB 中,可利用函数nichols 和dnichols 绘出连续和离散系统的Nichols 图。 (二)实验内容 1.一系统开环传递函数为 ) 2)(5)(1(50)(-++=s s s s H 绘制系统的bode 图,判断闭环系统的稳定性,并画出闭环系统的单位冲击响应。 2.一多环系统 ) 10625.0)(125.0)(185.0(7.16)(+++=s s s s s G 其结构如图所示 试绘制Nyquist 频率曲线和Nichols 图,并判断稳定性。 (三)实验要求

北京理工大学信号与系统实验实验5连续时间系统地复频域分析报告报告材料

实验5 连续时间系统的复频域分析 一、实验目的 1.掌握拉普拉斯变换及其反变换的定义,并掌握MATLAB 实现方法。 2.学习和掌握连续时间系统系统函数的定义及复频域分析方法。 3.掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号)(t x 的拉普拉斯变换定义为 )1.....(..........)()(dt e t x s X st ? +∞ ∞ --= 拉普拉斯反变换定义为 )2....(..........)(21)(ds e s X j t x j j st ?∞ +∞ -=σσπ 在MATLAB 中,可以采用符号数学工具箱的laplace 函数和ilaplace 函数进行拉氏变换和反拉氏变换。 L=laplace(F)符号表达式F 的拉氏变换,F 中时间变量为t ,返回变量为s 的结果表达式。 L=laplace(F,t)用t 替换结果中的变量s 。 F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换,返回时间变量为t 的结果表达式。 F=ilaplace(L,x)用x 替换结果中的变量t 。 除了上述ilaplace 函数,还可以采用部分分式法,求解拉普拉斯逆变换,具体原理如下: 当 X (s )为有理分式时,它可以表示为两个多项式之比: )3.(..........)()()(0 110 11a s a s a b s b s b s D s N s X N N N N M M M M +?+++?++==---- 式(3)可以用部分分式法展成一下形式 )4.....(.............)(2211N N p s r p s r p s r s X -++-+-= 通过查常用拉普拉斯变换对,可以由式(1-2)求得拉普拉斯逆变换。 利用 MATLAB 的residue 函数可以将 X (s )展成式(1-2)所示的部分分式展开式,该 函数的调用格式为:[r,p,k] = residue(b,a) 其中b 、a 为分子和分母多项式系数向量,r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数

第五章 频域分析法

第五章 频域分析法 时域分析法具有直观、准确的优点。如果描述系统的微分方程是一阶或二阶的,求解后可利用时域指标直接评估系统的性能。然而实际系统往往都是高阶的,要建立和求解高阶系统的微分方程比较困难。而且,按照给定的时域指标设计高阶系统也不是一件容易的事。 本章介绍的频域分析法,可以弥补时域分析法的不足。因为频域法是基于频率特性或频率响应对系统进行分析和设计的一种图解方法,故其与时域分析法相比有较多的优点。首先,只要求出系统的开环频率特性,就可以判断闭环系统是否稳定。其次,由系统的频率特性所确定的频域指标与系统的时域指标之间存在着一定的对应关系,而系统的频率特性又很容易和它的结构、参数联系起来。因而可以根据频率特性曲线的形状去选择系统的结构和参数,使之满足时域指标的要求。此外,频率特性不但可由微分方程或传递函数求得,而且还可以用实验方法求得。对于某些难以用机理分析方法建立微分方程或传递函数的元件(或系统)来说,具有重要的意义。因此,频率法得到了广泛的应用,它也是经典控制理论中的重点内容。 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、掌握测试系统的频率响应的方法。 二、预习内容 1、系统频响的方法。(见第四章波特图的介绍) 三、实验原理 1. N 阶系统系统的传递函数 用微分方程描述的N 阶系统为: 根据零状态响应(起始状态为零),则对其进行拉氏变换有: 则系统传递函数可表达为: 用差分方程描述的N 阶系统为: 根据零状态响应(起始状态为零),则对其进行拉氏变换有: 则系统传递函数可表达为: 2.根据系统传递函数的零极点图分析系统 零点:传递函数分子多项式的根。 极点:传递函数分母多项式的根。 根据零极点图的不同分布分析系统。 3.涉及到的Matlab 函数 (1)freqz 函数:实验六中出现过,可用来求单位圆上的有理z 变换的值。调用格式:同实验六 (2)zplane 函数:得到有理z 变换的零极点图。 调用格式:zplane(num,den)

其中,num和 den是按z ?1 的升幂排列的、z 变换分子分母多项式系数的行向量。 (3)roots 函数:求多项式的根。 调用格式:r=roots(c), c 为多项式系数向量;r 为根向量。 四、实验内容 1.系统零极点的求解 (1)求解系统和的零极点,验 证下面程序的运行结果,根据系统零极点图分析系统性质。 b=[1,0,-1]; a=[1,2,3,2]; zr=roots(b); pr=roots(a); plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2); grid; legend('零点','极点'); figure; zplane(b,a); (2)参考上述程序,绘制系统和 的零极点图,并分析系统性质。与用zplane 函数直接绘制系统零极点图(注:圆心的圆圈并非系统的零点)做比较。

(实验三)连续时间LTI系统的频域分析汇总

实验三 连续时间LTI 系统的频域分析 一、实验目的 1、掌握系统频率响应特性的概念及其物理意义; 2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用; 3、学习和掌握幅度特性、相位特性以及群延时的物理意义; 4、掌握用MA TLAB 语言进行系统频响特性分析的方法。 基本要求:掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法,深刻理解LTI 系统的频率响应特性的物理意义,理解滤波和滤波器的概念,掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。 二、实验原理及方法 1 连续时间LTI 系统的频率响应 所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response ),是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。 上图中x(t)、y(t)分别为系统的时域激励信号和响应信号,h(t)是系统的单位冲激响应,它们三者之间的关系为:)(*)()(t h t x t y =,由傅里叶变换的时域卷积定理可得到: )()()(ωωωj H j X j Y = 3.1 或者: ) () ()(ωωωj X j Y j H = 3.2 )(ωj H 为系统的频域数学模型,它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即 ? ∞ ∞ --= dt e t h j H t j ωω)()( 3.3 由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换,如果h(t)是收敛的,或者说 是绝对可积(Absolutly integrabel )的话,那么H(j ω)一定存在,而且H(j ω)通常是复数,

连续系统的复频域分析

实验四:连续系统的复频域分析 一、实验目的: 1、掌握连续与离散时间系统的正反复频域与Z域变换 2、掌握利用MATLAB进行零极点分析,进一步了解零极点对整个系统的影响 3、掌握simulink环境下系统建模与仿真以及系统求解。 二、实验内容: 1、已知某连续系统的系统函数为: (1)利用[r, p, k]=residue(num, den),求H(s)的极零点以及多项式系数; (2)画出系统的零极点分布图,判断系统得稳定性。 (3)求h(t),判断系统得稳定性。 2、已知某离散系统的系统函数为:, (1)利用[r, p, k]=residuez(num, den)求H(z)的极零点以及多项式系数; (2)画出零极点分布图,判断系统得稳定性。 (3)求单位函数响应用impz(b, a),判断系统是否稳定; 3、已知线性时不变微分方程 在Simulink环境下搭建起系统的仿真模型,并查看仿真结果曲线。(1)写出传递函数H(s),绘出系统模拟框图; (2)当f(t)分别为,,的零状态响应;且当与课本P81的结果进行比较(3)方程的初值为, ,求全响应; 4、已知某信号,n(t)为正态噪声干扰且服从N(0,0.22)分布,对此信号进行采样,采样间隔为0.001s,之后对此信号进行Botterworth低通滤波,从信号中过滤10HZ的输出信号,试对系统进行建模与仿真。 三、实验数据处理与结果分析: 第一题:题1_1:

>> num=[2,5]; den=[1,1,3,2]; [r,p,k]=residue(num,den) r = -0.5750 - 0.7979i -0.5750 + 0.7979i 1.1499 p =-0.1424 + 1.6661i -0.1424 - 1.6661i -0.7152 k =[]

第5章频域分析法习题解答

第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 的信息。

系统频域分析课程设计报告

系统频域分析课程设计 报告 Company number【1089WT-1898YT-1W8CB-9UUT-92108】

《综合仿真》课程设计报告 姓名 学号 同组成员 指导教师 时间 11周至14周

系统的频域分析 【目的】 (1) 加深对系统频域分析基本原理和方法的理解。 (2) 加深对信号幅度调制与解调基本原理和方法的理解。 (3) 锻炼学生综合利用所学理论和技术,分析与解决工程实际 问题的能力。 【研讨内容】 题目1.幅度调制和连续信号的Fourier 变换 本题研究莫尔斯码的幅度调制与解调。本题中信号的形式为 )π2sin()()π2sin()()π2cos()()(132211t f t m t f t m t f t m t x ++= 其中信号x (t )由文件定义,可用命令Load ctftmod 将文件定义的变量装入系统内存。运行命令Load ctftmod 后,装入系统的变量有 af bf dash dot f1 f2 t x 其中 bf af : 定义了一个连续系统H (s )的分子多项式和分母多项式。可利用freqs(bf,af,w)求出该系统的频率响应,也可用sys=tf(bf,af)得到系统的模型,从而用lsim 求出信号通过该系统的响应。 dash dot : 给出了莫尔斯码中的基本信号dash 和dot 的波形 f1 f2: 载波频率 t: 信号x (t )的抽样点 x: 信号x (t )的在抽样点上的值 信号x (t )含有一段简单的消息。Agend 007的最后一句话是

The future of technology lies in ··· 还未说出最后一个字,Agend 007就昏倒了。你(Agend 008)目前的任务就是要破解Agend 007的最后一个字。该字的信息包含在信号x (t )中。信号x (t )具有式(1)的形式。式中的调制频率分别由变量f1和f2给出,信号m 1(t ),m 2(t )和m 3(t )对应于字母表中的单个字母,这个字母表已用国际莫尔斯码进行编码,如下表所示: (1)字母B 可用莫尔斯码表示为b=[dash dot dot dot],画出字母B 莫尔 斯码波形; (2) 用freqs(bf,af,w)画出系统的幅度响应; (3) 利用lsim 求出信号dash 通过由sys=tf(bf,af)定义的系统响应,解释你所获得的结果; (4)用解析法推导出下列信号的Fourier 变换 )π2cos()π2cos()(21t f t f t m )π2sin()π2cos()(21t f t f t m

第5章_用MATLAB进行控制系统频域分析

第5章 用MATLAB 进行控制系统频域分析 一、基于MATLAB 的线性系统的频域分析基本知识 (1)频率特性函数)(ωj G 。 设线性系统传递函数为: n n n n m m m m a s a s a s a b s b s b s b s G ++???++++???++=---1101110)( 则频率特性函数为: n n n n m m m m a j a j a j a b j b j b j b jw G ++???++++???++=---)()()()()()()(1101110ωωωωωω 由下面的MATLAB 语句可直接求出G(jw)。 i=sqrt(-1) % 求取-1的平方根 GW=polyval(num ,i*w)./polyval(den ,i*w) 其中(num ,den )为系统的传递函数模型。而w 为频率点构成的向量,点右除(./)运算符表示操作元素点对点的运算。从数值运算的角度来看,上述算法在系统的极点附近精度不会很理想,甚至出现无穷大值,运算结果是一系列复数返回到变量GW 中。 (2)用MATLAB 作奈魁斯特图。 控制系统工具箱中提供了一个MATLAB 函数nyquist( ),该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。当命令中不包含左端返回变量时,nyquist ()函数仅在屏幕上产生奈氏图,命令调用格式为: nyquist(num,den) nyquist(num,den,w) 或者 nyquist(G) nyquist(G,w) 该命令将画出下列开环系统传递函数的奈氏曲线: ) () ()(s den s num s G = 如果用户给出频率向量w,则w 包含了要分析的以弧度/秒表示的诸频率点。在这些频率点上,将对系统的频率响应进行计算,若没有指定的w 向量,则该函数自动选择频率向量进行计算。 w 包含了用户要分析的以弧度/秒表示的诸频率点,MATLAB 会自动计算这些点的频率响应。 当命令中包含了左端的返回变量时,即: [re,im,w]=nyquist(G) 或

信号与系统报告 实验5 连续系统的复频域分析实验

信号与系统 实验报告 实验五连续系统的复频域分析 实验五连续系统的复频域分析 一、实验目的 1. 深刻理解拉普拉斯变换、逆变换的定义,掌握用MATLAB实现拉普拉斯变换、逆变换的方法。 2会求几种基本信号的拉氏变换。 3 掌握用MATLAB绘制连续系统零、极点的方法。 4 求解系统函数H(s)。 二

1已知连续时间信号f(t)=sin(t)u(t)、求出该信号的拉普拉斯变换,并用MATLAB 绘制拉普拉斯变换的曲面图。 syms t; ft=sin(t)*heaviside(t); Fs=Laplace(ft); a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); c=a+i*b; d=ones(size(a)); c=c.*c; c=c+d; c=1./c; c=abs(c); mesh(a,b,c); surf(a,b,c) axis([-0.5,0.5,-2,2,0,10]) colormap(hsv

) 2求[(1-e^(-at))]/t的拉氏变换。 syms t s a f1=(1-exp(-a*t))/t; F=laplace(f1,t,s) F = log(s+a)-log(s) 3求F(s)=-log(s)+ log(s+a)的拉氏逆变换syms t s a F =log(s+a)-log(s); f1=ilaplace(F,s,t) f1 = (1-exp(-a*t))/t

4已知某连续系统的系统函数为: H(s)=(s^2+3s+2)/(8s^4+2s^3+3s^2+5)试用MATLAB求出该系统的零极点,画出零极点分布图。 b=[1 3 2]; a=[8 2 3 0 5]; zs=roots(b); ps=roots(a); hold on plot(real(zs),imag(zs),'o'); plot(real(ps),imag(ps),'x'); grid axis([-2.5,1,-1,1]) 5已知H(s)=(s+1)/(s^2+s+1),绘制阶跃响应图形,冲激响应图形,频率激响应图形。 syms t s H=(s+1)/(s^2+s+1); f1=ilaplace(H,s,t); f2=heaviside(t);

实验三线性系统的频域分析

自动控制理论 上 机 实 验 报 告 学院:机电工程学院 班级:13级电信一班

: 学号: 实验三 线性系统的频域分析 一、实验目的 1.掌握用MATLAB 语句绘制各种频域曲线。 2.掌握控制系统的频域分析方法。 二、基础知识及MATLAB 函数 频域分析法是应用频域特性研究控制系统的一种经典方法。它是通过研究系统对正弦信号下的稳态和动态响应特性来分析系统的。采用这种方法可直观的表达出系统的频率特性,分析方法比较简单,物理概念明确。 1.频率曲线主要包括三种:Nyquist 图、Bode 图和Nichols 图。 1)Nyquist 图的绘制与分析 MATLAB 中绘制系统Nyquist 图的函数调用格式为: nyquist(num,den) 频率响应w 的围由软件自动设定 nyquist(num,den,w) 频率响应w 的围由人工设定 [Re,Im]= nyquist(num,den) 返回奈氏曲线的实部和虚部向量, 不作图 例4-1:已知系统的开环传递函数为2 526 2)(2 3++++=s s s s s G ,试绘制Nyquist 图,并判断系统的稳定性。

num=[2 6]; den=[1 2 5 2]; [z,p,k]=tf2zp(num,den); p nyquist(num,den) 极点的显示结果及绘制的Nyquist 图如图4-1所示。由于系统的开环右根数P=0,系统的Nyquist 曲线没有逆时针包围(-1,j0)点,所以闭环系统稳定。 p = -0.7666 + 1.9227i -0.7666 - 1.9227i -0.4668 若上例要求绘制)10,10(32-∈ω间的Nyquist 图,则对应的MATLAB 语句为: num=[2 6]; den=[1 2 5 2]; w=logspace(-1,1,100); 即在10-1和101之间,产生100个等距 离的点 nyquist(num,den,w) 2)Bode 图的绘制与分析 系统的Bode 图又称为系统频率特性的对数坐标图。Bode 图有两图,分别绘制开环频率特性的幅值和相位与角频率ω的关系曲线,称为对数幅频特性曲线和对数相频特性曲线。 MATLAB 中绘制系统Bode 图的函数调用格式为: bode(num,den) 频率响应w 的围由软件自动设定 bode(num,den,w) 频率响应w 的围由人工设定 图4-1 开环极点的显示结果及Nyquist 图

控制系统时域与频域性能指标的联系

控制系统时域与频域性能指标的联系 经典控制理论中,系统分析与校正方法一般有时域法、复域法、频域法。时域响应法是一种直接法,它以传递函数为系统的数学模型,以拉氏变换为数学工具,直接可以求出变量的解析解。这种方法虽然直观,分析时域性能十分有用,但是方法的应用需要两个前提,一是必须已知控制系统的闭环传递函数,另外系统的阶次不能很高。 如果系统的开环传递函数未知,或者系统的阶次较高,就需采用频域分析法。频域分析法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法,而且当系统的数学模型未知时,还可以通过实验的方法建立。此外,大量丰富的图形方法使得频域分析法分析高阶系统时,分析的复杂性并不随阶次的增加而显著增加。 在进行控制系统分析时,可以根据实际情况,针对不同数学模型选用最简洁、最合适的方法,从而使用相应的分析方法,达到预期的实验目的。 系统的时域性能指标与频域性能指标有着很大的关系,研究其内在联系在工程中有着很大的意义。 一、系统的时域性能指标 延迟时间t d 阶跃响应第一次达到终值h (∞)的50%所需的时间 上升时间 t r 阶跃响应从终值的10%上升到终值的90%所需的时间;对有振荡的系 统,也可定义为从0到第一次达到终值所需的时间 峰值时间t p 阶跃响应越过终值h (∞)达到第一个峰值所需的时间 调节时间 t s 阶跃响应到达并保持在终值h (∞)的±5%误差带内所需的最短时间 超调量%σ 峰值h( t p )超出终值h (∞)的百分比,即 %σ= () ()() ∞∞-h h h t p ?100% 二、系统频率特性的性能指标 采用频域方法进行线性控制系统设计时,时域内采用的诸如超调量,调整时间等描述系统性能的指标不能直接使用,需要在频域内定义频域性能指标。

噪声中正弦信号的经典法频谱分析

实验报告 一、实验名称 噪声中正弦信号的经典法频谱分析 二、实验目的 通过对噪声中正弦信号的经典法频谱分析,来理解和掌握经典谱估计的知识,以及学会应用经典谱估计的方法。 三、基本原理 1.周期图法:又称直接法。把随机信号)(n x 的N 点观察数据)(n x N 视为一能量有限信号,直接取)(n x N 的傅里叶变换,得)(jw N e X ,然后再取其幅值的平方,并除以N ,作为对)(n x 真 实的功率谱)(jw e P 的估计,以)(?jw PER e P 表示用周期图法估计出的功率谱,则2)(1)(?w X N w P n PER =。 2.自相关法:又称为间接法功BT 法。先由)(n x N 估计出自相关函数)(?m r ,然后对)(?m r 求傅里叶变换得到)(n x N 的功率谱,记之为)(?w P BT ,并以此作为对)(w P 的估计,即1,)(?)(?-≤=--=∑N M e m r w P jwm M M m BT 。 3.Bartlett 法:对L 个具有相同的均值μ和方差2σ的独立随机变量1X ,2X ,…,L X ,新随机变量L X X X X L /)(21+++= 的均值也是μ,但方差是L /2σ,减小了L 倍。由此得 到改善)(?w P PER 方差特性的一个有效方法。它将采样数据)(n x N 分成L 段,每段的长度都是M ,即N=LM ,第i 段数据加矩形窗后,变为L i e n x M w x M n jwn i N I PER ≤≤=∑-=-1,)(1)(?2 10 。把)(?w P PER 对应相加,再取平均,得到平均周期图2 1110 )(1)(?1)(∑∑∑==-=-==L i L i M n jwn i N i PER PER e n x ML w P L w P 。 4.Welch 法:它是对Bartlett 法的改进。改进之一是,在对)(n x N 分段时,可允许每一段的数据有部分的交叠。改进之二是,每一段的数据窗口可以不是矩形窗口,例如使用汉宁窗或汉明窗,记之为)(2n d 。这样可以改善由于矩形窗边瓣较大所产生的谱失真。然后按Bartlett

线性系统的频域分析报告

1 γ = 50 20- =s K0

原系统的伯德图: num/den = 1.2347 s + 1 ------------- 0.20154 s + 1 校正之后的系统开环传递函数为: num/den = 6.1734 s + 5 ------------------------------------------- 0.20154 s^4 + 1.6046 s^3 + 3.4031 s^2 + 2 s alpha =6.1261; P h a s e (d e g ) Bode Diagram Gm = Inf dB (at Inf rad/sec) , P m = 9.04 deg (at 3.14 rad/sec) -200204060 80M a g n i t u d e (d B )

[il,ii]=min(abs(mag1-1/sqrt(alpha))); wc=w( ii); T=1/(wc*sqrt(alpha)); numc=[alpha*T,1]; denc=[T,1]; [num,den]=series(num0,den0,numc,denc); [gm,pm,wcg,wcp]=margin(num,den); printsys(numc,denc) disp('D£?y??oóμ??μí3?a?·′?μYoˉêy?a:');printsys(num,den) [mag2,phase2]=bode(numc,denc,w); [mag,phase]=bode(num,den,w); subplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),'--',w,20*log10(mag2),'-.'); grid; ylabel('·ù?μ(db)'); title('--Go,-Gc,GoGc'); subplot(2,1,2); semilogx(w,phase,w,phase1,'--',w,phase2,'-',w,(w-180-w),':'); grid; ylabel('?à??(0)'); xlabel('?μ?ê(rad/sec)'); title(['D£?y?°£o·ù?μ?£á?=',num2str(20*log10(gm1)),'db','?à???£á?=',num2str(pm1),'0'; 'D£?yoó£o·ù?μ?£á?=',num2str(20*log10(gm)),'db','?à???£á?=',num2s tr(pm),'0']); 10-110 10 1 10 2 -60 -40-20020 40幅值(d b ) --Go,-Gc,GoGc 10 -110 10 1 10 2 -300 -200-1000 100相位(0) 频率(rad/sec) 矫正后系统的伯德图

实验八 连续系统的复频域分析

a=-0.5:0.08:0.5; b=-2:0.08:2; [a,b]=meshgrid(a,b); d=ones(size(a)); c=a+i*b; c=c.*c; c=c+d; c=1./c c=abs(c); surf(a,b,c); axis=([-0.5,0.5,-2,2,0.15]); title('单边正弦信号拉氏变换图'); colormap(hsv); 2. a=0:0.5:5; b=-20:0.1:20; [a,b]=meshgrid(a,b); c=a+i*b; c=(1-exp(-2*c)./c); c=abs(c) mesh(a,b,c); sufr(a,b,c); view(-10,20); axis([-0.5:-20,20,0.2]); title('拉氏变换S域像函数'); w=-20:0.1:20; Fw=(2*sin(w).*exp(i*w)); plot(w,abs(Fw)) title=('傅里叶变换'); xlabel('频率w'); 3. a=[1,2,-3,2,1]; b=[1,4]; sjdt(a,b); a=[1,5,16,30]; b=[5,20,25,0]; sjdt(a,b); 4. a=[8,2,3,15]; b=[1,3,2]; [p,q]=sjdt(a,b); 5.

a=[1,0]; b=[1]; impulse(b,a) a=[1,2]; b=[1]; impulse(b,a) a=[1,-2]; b=[1]; impulse(b,a) a=[1,1,16,25]; b=[1]; impulse(b,a) a=[1,0,16]; b=[1]; impulse(b,a) a=[1,-1,16,25]; b=[1]; impulse(b,a) 6. q=[0,0]; p=[-100,100]; f1=0; f2=1000; k=0.1; splxy(f1,f2,k,p,q) q=[0,0]; f1=0; f2=1000; k=0.1; p=[-500,-1000]; splxy(f1,f2,k,p,q) q=[0,0]; f1=0; f2=1000; k=0.1; p=[-2000,-4000]; splxy(f1,f2,k,p,q) 7. 8

实验二连续时间信号的频域分析

实验二连续时间信号的频域分析 令狐采学 一、实验目的 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因; 3、掌握连续时间傅里叶变换的分析方法及其物理意义; 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质; 5、学习掌握利用Matlab语言编写计算CTFS、CTFT和DTFT 的仿真程序,并能利用这些程序对一些典型信号进行频谱分析,验证CTFT、DTFT的若干重要性质。 基本要求:掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用Matlab编程完成相关的傅里叶变换的计算。 二、原理说明 1、连续时间周期信号的傅里叶级数CTFS分析 任何一个周期为T1的正弦周期信号,只要满足狄利克利条

件,就可以展开成傅里叶级数。 三角傅里叶级数为: ∑∞=++=1000)]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 其中102T π ω=,称为信号的基本频率(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 为指数形式的傅里叶级数的系数,按如下公式计算:

自动控制原理线性系统的频域分析实验报告

实验四 专业 自动化 班号 03班 指导教师 陈艳飞 姓名 胡波 实验名称 线性系统的频域分析 实验日期 第 次实验 一、实验目的 1.掌握用MATLAB 语句绘制各种频域曲线。 2.掌握控制系统的频域分析方法。 二、实验内容 1.典型二阶系统 2 2 22)(n n n s s s G ωζωω++= 绘制出6=n ω,1.0=ζ,0.3,0.5,0.8,2的bode 图,记录并分析ζ对系统bode 图的影响。 解: 程序如下: num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36]; den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36]; w=logspace(-2,3,100); bode(num,den1,w) grid hold bode(num,den2,w) bode(num,den3,w) bode(num,den4,w) bode(num,den5,w)

-100-80-60-40-200 20M a g n i t u d e (d B )10 -2 10 -1 10 10 1 10 2 10 3 P h a s e (d e g ) Bode Diagram Frequency (rad/sec) 分析:随着.0=ζ的增大 ,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大. 2.系统的开环传递函数为 ) 5)(15(10 )(2+-= s s s s G ) 106)(15() 1(8)(22++++= s s s s s s G ) 11.0)(105.0)(102.0() 13/(4)(++++= s s s s s s G 绘制系统的Nyquist 曲线、Bode 图和Nichols 图,说明系统的稳定性,并通过绘制阶跃响应曲线验证。 解: 程序如下 奈氏曲线: (1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); w=logspace(-1,1,100); nyquist(num1,den1,w)

系统的复频域分析实验报告

实验六 系统的复频域分析 信号)(t x 的拉普拉斯变换 ?∞ ∞--=dt e t x s X st )()( (6.1) 是连续时间傅立叶变换地推广。连续时间傅立叶变换在研究连续时间信号与系统中是很有用的。然而,许多信号不存在傅立叶变换而存在拉普拉斯变换,这使得拉普拉斯变换成为线性时不变系统分析的一种有用方法。对一大类信号来说,它们的拉普拉斯变换可以表示为s 的多项式之比,即 ) () ()(s D s N s X = 这里)(s N 和)(s D 分别称作分子和分母多项式。能表示成多项式之比的变换称为有理变换,这里作为满足线性常系数微分方程的LTI 系统的系统函数中常常出现。除了一个标量因子外,有理变换是完全由多项式)(s N 和)(s D 的根决定的,这些根分别称为零点和极点。由于这些根在LTI 系统的研究中起着重要的作用,所以它们以零极点图的方式展现出来的是很方便的。这一章将用拉普拉斯变换在复频域研究LTI 系统的一些性质。 §6.1 MATLAB 函数lsim (用于系统函数) 目的 用lsim 仿真由系统函数表征的因果LTI 系统的输出。 相关知识 第二章所讨论的是如何用lsim 命令仿真一个输出满足一个线性常系数微分方程的因果LTI 连续时间系统。因为系统函数唯一地表征了关联系统输入和输出的微分方程。所以由系统函数表征的因果LTI 系统的输出也能够用lsim 仿真。如果系统函数给出如下形式: ) ()1()1()()1()1()(N a s N a s a M b s M b s b s H N M +-+++-++= (6.2) 那么,对输入地系统)(t x 的系统输出就能用lsim(b,a,x,t)仿真,其中MATLAB 向量b 和a 包含了分子分母s 多项式的系数。

相关文档
最新文档