线性系统的频域分析
《自动控制原理》实验3.线性系统的频域分析
《自动控制原理》实验3.线性系统的频域分析实验三线性系统的频域分析一、实验目的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) 返回奈氏曲线的实部和虚部向量，不作图2s?6例4-1：已知系统的开环传递函数为G(s)?3，试绘制Nyquists?2s2?5s?2图，并判断系统的稳定性。
num=[2 6]; den=[1 2 5 2]; nyquist(num,den)极点的显示结果及绘制的Nyquist图如图4-1所示。
由于系统的开环右根数P=0，系统的Nyquist曲线没有逆时针包围（-1，j0）点，所以闭环系统稳定。
p =-0.7666 + 1.9227i -0.7666 - 1.9227i -0.4668图4-1 开环极点的显示结果及Nyquist图若上例要求绘制??(10?2,103)间的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图有两张图，分别绘制开环频率特性的幅值和相位与角频率?的关系曲线，称为对数幅频特性曲线和对数相频特性曲线。
四、线性系统的频域分析法
其中： A()Ac (j) 幅频特性
A
() (j) 相频特性
RC网络频率特性的物理意义：
1 A()
0.707
频带宽度
b
01 2 3 4 5
TTTT T
() 0
相角迟后
90
01 2 3 4 5
TTTT T
对稳定的线性系统，其频率特性如下：
设： (s)C R ((s s))b a 0 0 ssm n b a 1 1 s sm n 1 1 .... a .b .m n 1 1 s s a b n m
微分环节： s 惯性环节： 1/(Ts1) 一阶微分环节： Ts1
振荡环节： 1 /s (2/ n 2 2s/ n 1 )0 , 1
二阶微分环节： s2/n22 s/n 1 ,01
例如：G(s)s(0.5s K 1()ss( 21 )2s5) 由上述的5个环节组成。
A()1/ ()900
db 60 40 20 0 900
[20]
0.1
1
j
0
幅相曲线
对数频率特性曲线
L()2l0g A()
20lg () 900
10
3）微分环节： s 由 G(s)s
A() ()900
db 60 40 20 0 90 0 00
uc
ur
ur Asi nt c u c
设初值为0， 对上式拉氏变换，设A=1，得：
Uc(s)RC 1s1Ur(s) s1/1T/Ts2 2
RC网络
TRC
s1x/Tsy2sz2 (xy)s2( s (z1 /T y)/T s(2) s x 2 )2z/T
线性系统的频域分析
第五章 线性系统的频域分析频域分析法是应用频率特性研究线性系统的一种经典方法。
它以控制系统的频率特性作为数学模型，以伯德图或其他图表作为分析工具，来研究、分析控制系统的动态性能与稳态性能。
频域分析法由于使用方便，对问题的分析明确，便于掌握，因此和时域分析法一样，在自动控制系统的分析与综合中，获得了广泛的应用。
本章研究频率特性的基本概念、典型环节和控制系统的频率特性曲线、奈奎斯特稳定判据以及开环频域性能分析等内容。
§5-1 频率特性的基本概念一、频率特性的基本概念频率特性又称频率响应，它是系统（或元件）对不同频率正弦输入信号的响应特性，对于线性系统，若其输入信号为正弦量，则其稳态输出信号也将是同频率的正弦量，但其幅值和相位都不同与输入量。
下面以RC 电路为例，说明频率特性的基本概念。
图5-1所示的RC 电路，)(t u i 和)(0t u 分别为电路的输入电压和输出电压，电路的微分方程为：)()()(00t u t u dtt du Ti =+ 式中T=RC 为电路的时间常数。
RC 电路的传递函数为11)()(0+=Ts s U s U i （5-1） Rui )t图 5-1 RC 电路当输入电压为正弦函数t U t u i i ωsin )(=，则由式（5-1）可得22011)(11)(ωω+⋅+=+=s U Ts s U Ts s U i i 经拉氏反变换得电容两端的输出电压)sin(11)(122/220T tg t T U e T T U t u iT t i ωωωωω---+++=式中，第一项为输出电压的暂态分量，第二项为稳态分量，当∞→t 时，第一项趋于零，于是)sin(1|)(1220T tg t T U t u i t ωωω-∞→-+=)](sin[)(ωϕωω+=t A U i （5-2）式中：2211)(TA ωω+=，T tgωωϕ1)(--=，分别反映RC 网络在正弦信号作用下，输出稳态分量的幅值和相位的变化，二者皆是输入正弦信号频率ω的函数。
自动控制原理课件：线性系统的频域分析
包围坐标原点 − 周。
m
F (s)
K1 ( s z j )
j 1
n
i 1
( s pi )
24
• 02
基本概念
m
1 G ( s) H ( s) F ( s)
K1 ( s z j )
j 1
在 平面上的映射曲线 F 1 G ( j ) H ( j )将按逆时针方向
围绕坐标原点旋转 = − 周。
如果在s平面上，s沿着奈奎斯特回线顺时针方向移动一周时，
在 平面上的映射曲线围绕坐标原点按逆时针方向旋转 =
周，则系统为稳定的。
26
根据
( 1, j 0)
L( ) 20 lg K 20 lg 1 12 2 20 lg 1 22 2
( ) arctg 1 arctg 2
τ2
20dB / dec 1
2
L3 ( )
L2 ( )
40dB / dec
( )
0
L( )
90
A( ) 1, ( )
L ( ) 20 lg A( ) 0
L( )
jQ( )
L( ) 0
0
( )
1
0
1
P( )
1
0
30
60
16
5.3
系统开环频率特性图
设开环系统由n个典型环节串联组成
G(s ) G 1(s )G 2(s ) G n(s )
这意味着 的映射曲线 F 围绕原点运动的情况，相当于
自动控制理论 线性系统的频域分析法
() tg 1 Q() P( )
线性系统的频域分析法>>线性系统的频域特性
频率特性与传递函数的关系为：
G( j ) G(s) |s j
由于这种简单关系的存在，频率响应法和利用传递函数的时域 法在数学上是等价的。
[结论]：当传递函数中的复变量s用 j代替时，传递函数就转变
第六章 线性系统的频域分析法
1 线性系统的频率特性及图示 2 开环系统的典型环节 3 频率域稳定判据 4 稳定裕度 5 闭环系统的频域特性
线性系统的频域分析法>>线性系统的频域特性
6.1 频率特性的基本概念
考察一个系统的好坏，通常用阶跃输入下系统的阶跃响应 来分析系统的动态性能和稳态性能。
有时也用正弦波输入时系统的响应来分析，但这种响应并 不是单看某一个频率正弦波输入时的瞬态响应，而是考察频率 由低到高无数个正弦波输入下所对应的每个输出的稳态响应。 因此，这种响应也叫频率响应。
N (s)
Rm
(s p1)(s p2 )...(s pn ) (s p1)(s p2 )...(s pn ) (s j )(s j )
k1 k2 ... kn kc1 kc2
s p1 s p2
s pn s j s j
拉氏反变换为：
c(t) k1e p1t k2e p2t ... kne pnt kc1e jt kc2e jt
频率响应法的优点之二在于它可以用图来表示，这在控制 系统的分析和设计中有非常重要的作用。
由实验方法求频率特性
正弦信号 发生器
实验装置 （系统或元件）
双踪 示波器
图 求频率特性的实验方法
系统的幅频特性： | G( j) | Y
第五章线性系统的频域分析法
对 A(ω ) 求导并令等于零，可解得 A(ω ) 的极值对应的频率 ω r 。
ω r = ω n 1 2ζ 2
该频率称为谐振峰值频率。可见，当 ζ = 当ζ
> 1 2
s = jω
G( jω) =| G( jω) | e
j∠G( jω)
= A(ω)e
j (ω)
G( jω) = G(s) |s= jω
G( jω) = G(s)|s= jω =| G( jω)| e j∠G( jω) = A(ω)e j(ω)
A A j (ω ) k1 = G( jω ) e k2 = G( jω ) e j (ω ) 2j 2j
可以作为系统模型
G( jω) = G(s) |s= jω = G( jω) e j(ω)
定义 幅频特性
A(ω ) =| G( jω ) |
(ω ) = ∠G ( jω )
它描述系统对不同频率输入信号在稳态时的放大特性； 它描述系统对不同频率输入信号在稳态时的放大特性； 相频特性
它描述系统的稳态响应对不同频率输入信号的相位移特性； 它描述系统的稳态响应对不同频率输入信号的相位移特性； 幅频特性和相频特性可在复平面上构成一个完整的向量 G ( jω ), 频率特性。 频率特性 G ( jω ) = A(ω )e j (ω ) ，它也是 ω 的函数。G( jω) 称为频率特性 还可将 G ( jω ) 写成复数形式，即
A(ω ) = 1 1 + T 2ω 2 ,
G (s) =
1 Ts + 1
G ( jω ) =
1 jT ω + 1
(ω ) = tg 1T ω
幅频特性 L(ω) = 20log A(ω) = 20log K 20log 1+ T 2ω2 低频段：当Tω << 1时，ω 高频段：当 Tω >> 1时， ω
线性系统的频域分析法
5.1 频率特性
lg
1 0
2
0.301
3
0.477
4
0.602
5
0.699
6
0.778
7
0.845
8
0.903
9
0.954
10
1
※※
( )
40
20 0dB -20 -40
2、对数频率特性曲线 [ 伯德（Bode）图 ]
L ( ) 20 lg A( ) 20 lg G ( j ) ( dB )
L ( ) 20 lg (T ) 1 20 lg T
2
当 T 即 T 1 时
L(ω)dB 40 20 0dB -20 - 40
1
T
1 T
当
1 T
时 时
20 lg T 0
20 lg T 20
dB
dB
10 T
频 率 特 性 ： G ( j ) 1 j T 1
( ) tg T
1
A ( )
1 T 1
2 2
ω 1/10T φ (ω )(度) -5.7 L(ω )（dB）
从到值 取 代入计算,得
对数幅频特性曲线 Bode图如右
1/5T -11.3
1/2T -26.6
2.频域法的基本思想：利用系统的开环频率特 性来分析闭环响应。对系统进行定性分析和定量 计算。
3.频率特性的性质 考察一个系统的好坏，通常用阶跃输入下系统的阶跃响应 来分析系统的动态性能和稳态性能。
有时也用正弦波输入时系统的响应来分析，但这种响应并 不是单看某一个频率正弦波输入时的瞬态响应，而是考察频率 由低到高无数个正弦波输入下所对应的每个输出的稳态响应。 因此，这种响应也叫频率响应。
自动控制原理第五章线性系统的频域分析法
自动控制原理第五章线性系统的频域分析法1、基本内容和要点（l）频率特性系统的稳态频率响应，频率响应的物理概念及数学定义；求取频率特性的分析法和实验法。
（2）典型环节的频率特性比例、惯性、积分、微分、振荡、延迟环节的频率特性和对数频率特性。
非最小相位环节的频率特性。
（3）反馈控制系统的开环频率特性研究系统开环频率特性的意义。
单环系统开环对数频率持性的求取与绘制。
最小相位系统开环对数幅频特性与相频特性间的对应关系。
（4）奈奎斯特稳定判据幅角定理。
S平面与F平面的映射关系。
根据开环频率特性判别闭环系统稳定性的奈氏判据。
奈氏判据在多环系统中的应用和推广。
系统的相对稳定性。
相角与增益稳定裕量。
（5）二阶和高阶系统的频率域性能指标与时域性指标。
系统频率域性能指标。
二阶和高阶系统暂态响应性能指标与频率域性能指标间的解析关系及近似关系。
（6）系统的闭环频率特性开环频率特性与闭环频率特性间的解析关系。
用等M圆线从开环频率特性求取闭环频率特性。
用尼氏图线从开环对数频率特性求取闭环频率特性。
2、重点（l）系统稳态频率响应和暂态时域响应的关系。
（2）系统开环频率特性的绘制，最小相位系统开环频率特性的特点。
（3）奈奎斯特稳定判据和稳定裕量。
5-1引言第三章，时域分析，分析系统零、极点与系统时域指标的关系；典型二阶系统极点或和n与时域指标tp、和t、tr及稳态误差等的关系，及高阶系统的近似指标计算；第四章，根轨迹分析，研究系统某一个参数变化对系统闭环极点的影响；本章讨论系统零、极点对系统频率域指标的关系，频域指标又分开环频域指标和闭环频域指标，它们都是在频域上评价系统性能的参数。
频域分析是控制理论的一个重要分析方法。
5-2频率特性1．频率特性的基本概念理论依据定理：设线性定常系统G()的输入信号是正弦信号某(t)某int，在过度过程结束后，系统的稳态输出是与输入同频率的正弦信号，其幅值和相角都是频率的函数，即为c(t)Y()in[t()]。
实验四 线性系统的频域分析
实验四线性系统的频域分析
线性系统的频域分析是一种利用线性系统的响应特性来提高系统性能的有效手段，它
在系统设计中起着重要的作用。
其主要思想是将系统的响应特性根据其与频率之间的关系
进行分割，从而更好地理解该响应的物理规律。
本文的目的是介绍线性系统的频域分析方法。
线性系统的频域分析分为时域分析和频域分析两种技术。
时域分析是检测一个系统在
其他变量没有变化时，系统输出信号形状及其随时间变化趋势的一种分析方法。
时域分析中，将系统的输入和输出逐样本放入示波器进行分析及测试。
频域分析是通过将系统的输
入和输出信号进行频谱分析，将它们映射到频率轴上进行分析的一种方法。
在频域分析中，我们可以通过频谱分析仪、傅里叶变换、系统增益、阶跃响应等技术来检测系统响应的特性，得出系统的频率响应函数，从而研究系统是否属于线性系统。
线性系统的频域分析一般步骤如下：
1、定义时域函数并将其傅里叶变换，从而得到其频域函数；
2、计算系统的增益及其全频响应曲线，以便了解频率和增益之间的关系；
3、根据阶跃响应的拟合结果，利用积分和微分的技巧，确定系统的阶跃函数；
4、选择优化算法，进行系统参数优化调整，使系统达到所需要的设计目标。
以上就是线性系统的频域分析方法介绍，从分析输入输出信号，到频域拟合分析，再
到进行参数优化调整，这一系列的步骤可以帮助我们更好的理解系统的物理机理，实现系
统的最佳设计性能。
自动控制原理第五章 线性系统的频域分析法-5-1
如同收音机、电视机一样，任一系统的频率响应反映系统的频率特性，体现系统的控制性能。
系统频率特性物理意义明确。应用频率特性分析研究系统性能的方法称为频域分析法。
控制系统的频域分析法兼顾动态响应和噪声抑制的要求，可以拓展应用于非线性系统。
频率特性定义
分别称为系统的幅频特性和相频特性。
系统数学模型间的关系
控 制 系 统
傅氏变换
拉氏变换
g(t)
数学建模
例5.1-1
图示系统，设输入为r(t)=sin(5t)，计算系统的频率响应和稳态误差。
当
1
2
3
4
5
6
7
8
9
10
20
100
1
2
3
4
5
6
7
8
9
10
0
0.301
0.477
0.602
0.699
0.788
0.845
0.903
0.954
1
十倍频程
两倍频程
0.1
0.2
200
十倍频程
十倍频程
对数坐标的单位长度
⑶ 对数频率特性曲线
对数幅频特性曲线 纵坐标： ，线性刻度，单位为分贝(dB) 横坐标：ω ，对数刻度，单位为弧度/秒(rad/s)
绘制一阶系统幅相频率特性曲线
解：系统频率特性为
且有
即
复平面上位于第Ⅳ象限圆心为(1/2,j0)，半径为1的半圆。
箭头表示随ω增加,曲线的运动方向
2. 对数频率特性曲线（对数坐标图、伯德(Bode)图）
⑴ 频率特性的常用对数函数
线性系统的频域分析_自动控制原理
X G(-j )X d 1 G(s) 2 (s j ) S -j 2 2j s X G(j )X d 2 G(s) 2 (s - j ) S j 2 2j s G(j ) | G(j ) | e j G(-j ) | G(-j ) | e - j | G(j ) || G(-j ) |
第五章 线性系统的频域分析 §1 频率响应及其描述
一．频率特性 1.频率特性的基本概念 a.RC网络
右图所示的RC 网络的微分方程为
0 T dU dt U 0 U i
R UI C U0
式中
T RC 则
U 0 (S) U i (S)
1 TS 1
设 U i Asin t U 0 (S)
说明: 1.在稳态求出的输出信号 与输入信号的幅值比是 的非 线性函数, 称为幅频特性 Y/X | j ) | 2.输出信号与输入信号的 相位差是的非线性函数 ,称 为相频特性 .它描述在稳态情况下 ,当系统输入不同频率 的谐波信号时 , 其相位产生超前 ( 0 )或滞后( 0 )的 特性. 3.幅频特性和相频特性总 称为频率特性 , 记为 G(j ) G(j ) e jG(j ) 4.频率特性的求取 G(j ) G(s) s j
X(t) xsint Y(S)
B( s ) x ( s s1 )( s s2 ) ( s sn ) (s j )(s - j ) d1 d2 c1 cn s j s j s s1 s sn
y(t) d1e - jt d 2e jt c1e s1t c n e sn t 对于稳定系统 ,由于极点S1 , S2 , , Sn都有负实部 , 所以当t 时 y ss ( t ) d1e jt d 2e jt
第4章 线性系统的频域分析
第4章线频域分析法频域分析方法是根据系统的频率特性来分析系统的性能，也常称为频率特性法或频率法。
频域分析法有以下特点，首先是频率特性有明确的物理意义。
系统的频率响应可以用数学模型算出，也可以通过实际的频率特性实验测出。
这一点在工程实践上价值很大，特别是对结构复杂或机理不明确的对象，频率分析法提供了一个处理这类问题的有效方法。
频率法计算简单，只用很小的计算量和很简单的运算方法，再辅以作图，便可以完成分析与综合的工作。
当前已有一套完整便捷的基于频率法的计算机辅助设计软件，可以代替人工完成绝大部分的设计工作。
频率法也有其缺点和局限性。
频率法只适合用于线性定常系统。
从原理上讲频率法不能用于非线性系统或时变系统。
虽然在研究非线性系统时也借用了频率法的一些思想，但只能在特定的条件下解决一些很有局限性的问题。
本章研究频率特性的基本概念、图示方法、控制系统的稳定性判据、系统性能的频域分析方法。
4.1 频率特性系统的频率特性描述了线性系统在正弦信号输入下其稳态输出和输入的关系。
为了说明频率特性的概念，下面分析线性系统在正弦输入信号的作用下，其输出信号和输入信号间的关系。
设线性定常系统输入信号为()r t ，输出信号为()c t ，如图4-1所示。
图中G(s)为系统的传递函数。
即 1011111()()()mm m m n n n nb s b s b s b C s G s R s s a s a s a ----++⋅⋅⋅++==++⋅⋅⋅++ （n m ≥） （4-1）若在系统输入端作用一个时间的谐波函数，即0()s i n ()r t r t ωϕ=⋅+ ，式中，0r 是振幅；ω是频率；ϕ是相角。
为简便起见，假设0ϕ=，则0()sin r t r t ω=⋅ 图4-1 一般线性定常系统由于0022()()()r r R s s s j s j ωωωωω==++- （4-2）系统输出()C s 为10110111()()()()()m m m m n n n n b s b s b s b r C s G s R s s a s a s a s j s j ωωω----++⋅⋅⋅++==⋅++⋅⋅⋅+++-1()ni i i C B Ds s s j s j ωω==++-+-∑（4-3）式中，i s 为系统特征根，即极点（设为互异）；C i ，B ，D 均为相应极点处留数。
第五章 线性系统的频域分析法
4.还可以推广到研究某些非线性系统。
时域分析法与频域分析法比较：
时域分析法是分析控制系统的直接方法，比较直观、 精确。当往往需要求解复杂的微分方程。 频域分析法是一种图解分析法。它依据系统的又一种 数学模型——频率特性，利用频域指标和时域指标之间的 对应关系，间接地揭示系统的暂态特性和稳态特性，简单 迅速地判断某些环节或者参数对系统的暂态特性和稳态特 性的影响，并能指明改进系统的方向。也是一种工程上常 用的方法。
2 0.707 2
时，谐振峰值 M r 1 。
2 , (0, r ), 0 2 0 2 , ( , ), r 2
4.无谐振时
2 1, (0, ), 2
A( )
1
2 2 2 1 2 4 2 n n 2
参见《信号与系统》
频域分析法的基本介绍 •控制系统的频率特性反映正弦信号作用下系统响应的性能， 是系统的一种数学模型。 •应用频率特性来研究线性系统的经典方法称为频域分析法。 频域分析法具有以下特点：
1.控制系统及其元部件的频率特性可以运用分析法或者实验 法获得，并可用多种形式的曲线来表示，因而系统分析和控 制器设计可以应用图解法进行。
4.系统的开环幅相曲线（Nyquist图）
5.系统的开环对数频率特性曲线（bode图） 6.传递函数的频域实验确定
7.延迟环节和延迟系统
重点掌握最小相位情况的各个知识点，非最小相位情况的考试不考，考研可能考。
1.典型环节
2.最小相位环节的频率特性
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
线性系统的频域分析
1．实验目的
1. 掌握用MATLAB语句绘制各种各样频域曲线。
2. 掌握控制系统的频域分析方法。
二．练习：
1．典型二阶系统
绘制出，，0.3，0.5，0.8，2的bode图，记录并分析对系统bode图的影响。
解：MATLAB编程如下：
>> 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
Current plot held
>> bode(num,den2,w)
>> bode(num,den3,w)
>> bode(num,den4,w)
>> bode(num,den5,w)
（2）系统的开环传递函数为
绘制系统的Nyquist曲线Bode图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。
解：（1）MATLAB如下
>> 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)
（2）MATLAB编程如下：
>> num2=[8,8];den2=conv([1,0],conv([1,0],conv([1,15],[1,6,10]))); >> w=logspace(-1,1,100);
>> nyquist(num2,den2)
（3）MATLAB编程如下
>> num3=[4/3,4];den3=conv([1,0],conv([0.02,1],conv([0.05,1],[0.1,1]))); >> w=logspace(-1,1,100);
>> nyquist(num3,den3)
Bode图：
MATLAB编程如下：
>> num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5])));
>> num2=[8,8];den2=conv([1,0],conv([1,0],conv([1,15],[1,6,10]))); >> num3=[4/3,4];den3=conv([1,0],conv([0.02,1],conv([0.05,1],[0.1,1]))); >> bode(num1,den1)
>> grid
>> hold
Current plot held
>> bode(num2,den2)
>> bode(num3,den3)
阶跃响应曲线：
(1). >> num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5]))); >> step(num1,den1)
>> grid
(2).>> num2=[8,8];den2=conv([1,0],conv([1,0],conv([1,15],[1,6,10]))); >> step(num2,den2)
>> grid
(3)
>> num3=[4/3,4];den3=conv([1,0],conv([0.02,1],conv([0.05,1],[0.1,1]))); >> step(num3,den3)
>> grid
3．已知系统的开环传递函数为。
求系统的开环截止频率、穿越频率、幅值裕度和相位裕度。
应用频率稳定判据判定系统的稳定性。
解：MATLAB编程如下：
>> num=[0 0 1 1];
>> den=[0.1 1 0 0];
>> w=logspace(-1,3,200);
>> bode(num,den,w)
>> [gm,pm,wcg,wcp]=margin(num,den);
>> gm,pm,wcg,wcp
gm =
Inf
pm =
44.4594
wcg =
Inf wcp =
1.2647 >> grid。