在自动控制原理课程中,利用折线式伯德图计算截止频率是很常见的题型,下面介绍两种做法。

对于以下传递函数:
G(s)=50s2(s2+s+1)(10s+1)=G1(s)G2(s)G3(s)G4(s)G5(s){G(s)=\frac{50}{s^2(s^2+s+1)(10s+1)}=G_1(s)G_2(s)G_3(s)G_4(s)G_5(s)}G(s)=s2(s2+s+1)(10s+1)50​=G1​(s)G2​(s)G3​(s)G4​(s)G5​(s)
其中:
G1(s)=50,比例环节{G_1(s)=50,比例环节}G1​(s)=50,比例环节
G2(s)=1s,积分环节{G_2(s)=\frac{1}{s},积分环节}G2​(s)=s1​,积分环节
G3(s)=1s,积分环节{G_3(s)=\frac{1}{s},积分环节}G3​(s)=s1​,积分环节
G4(s)=1s2+s+1,振荡环节,ω1=1rad/s{G_4(s)=\frac{1}{s^2+s+1},振荡环节,\omega_1=1rad/s}G4​(s)=s2+s+11​,振荡环节,ω1​=1rad/s
G5(s)=110s+1,延时环节,ω2=0.1rad/s{G_5(s)=\frac{1}{10s+1},延时环节,\omega_2=0.1rad/s}G5​(s)=10s+11​,延时环节,ω2​=0.1rad/s
首先画出伯德图的草图:

在卡西欧计算器上,利用未经化简的L(ω){L(\omega)}L(ω)求解的结果是ωc=1.4193rad/s{\omega_c=1.4193rad/s}ωc​=1.4193rad/s,该值作为本文的真值。

方法一:利用L(ω){L(\omega)}L(ω)的分段函数特点

利用[1]中的方法(具体思路不再重复),可以求出转折频率的估计值:
(1)假设ωc<0.1rad/s{\omega_c<0.1rad/s}ωc​<0.1rad/s,只考虑比例和积分环节,则有L(ω)=20lg(50)−40lg(ω){L(\omega)=20lg(50)-40lg(\omega)}L(ω)=20lg(50)−40lg(ω)
令L(ω)=0{L(\omega)=0}L(ω)=0,解得ω=50{\omega=\sqrt{50}}ω=50​,不符合假设,舍去。
(2)假设0.1rad/s<ωc<1rad/s{0.1rad/s<\omega_c<1rad/s}0.1rad/s<ωc​<1rad/s,不考虑振荡环节,则有
L(ω)=20lg(50)−40lg(ω)−20lg(10ω){L(\omega)=20lg(50)-40lg(\omega)-20lg(10\omega)}L(ω)=20lg(50)−40lg(ω)−20lg(10ω)
令L(ω)=0{L(\omega)=0}L(ω)=0,解得ω=53{\omega=\sqrt[3]{5}}ω=35​,不符合假设,舍去。
(3)假设ωc>1rad/s{\omega_c>1rad/s}ωc​>1rad/s,则有
L(ω)=20lg(50)−40lg(ω)−20lg(10ω)−40lg(ω){L(\omega)=20lg(50)-40lg(\omega)-20lg(10\omega)-40lg(\omega)}L(ω)=20lg(50)−40lg(ω)−20lg(10ω)−40lg(ω)
令L(ω)=0{L(\omega)=0}L(ω)=0,解得ω=55=1.3797rad/s{\omega=\sqrt[5]{5}=1.3797rad/s}ω=55​=1.3797rad/s,符合假设,是正确的截止频率。

方法二:利用ω{\omega}ω靠近零处的特点

当ω{\omega}ω很小时,除了积分环节,其他环节都可以忽略,因此:
L(ω)=20lg(50)−40lg(ω){L(\omega)=20lg(50)-40lg(\omega)}L(ω)=20lg(50)−40lg(ω)
取ω=0.01rad/s{\omega=0.01rad/s}ω=0.01rad/s,有L(0.01)=20lg(50)+80{L(0.01)=20lg(50)+80}L(0.01)=20lg(50)+80
从ω=0.01rad/s{\omega=0.01rad/s}ω=0.01rad/s开始,L(ω){L(\omega)}L(ω)不断下降,因此可以把ωc{\omega_c}ωc​看作是下降过程中的一点。从ω=0.01rad/s{\omega=0.01rad/s}ω=0.01rad/s处下降到ωc{\omega_c}ωc​处,L(ω){L(\omega)}L(ω)首先经过了ω=0.1rad/s{\omega=0.1rad/s}ω=0.1rad/s,然后经过了ω=1rad/s{\omega=1rad/s}ω=1rad/s,由此画出下面的阶梯图:

对于图中折线的斜率,例如-40dB/dec,表示ω{\omega}ω每增大10倍,L(ω){L(\omega)}L(ω)就减小40,其余的同理。
因此,可以直接写出L(ωc){L(\omega_c)}L(ωc​)的表达式:
L(ωc)=L(0.01)−40−60−100lg(ωc)−lg(1)lg(10)−lg(1)=0{L(\omega_c)=L(0.01)-40-60-100\frac{lg(\omega_c)-lg(1)}{lg(10)-lg(1)}=0}L(ωc​)=L(0.01)−40−60−100lg(10)−lg(1)lg(ωc​)−lg(1)​=0
前面很好理解,最后一项是什么呢?不妨将-100dB/dec的直线延伸到ω=10rad/s{\omega=10rad/s}ω=10rad/s看看:

从ω=1rad/s{\omega=1rad/s}ω=1rad/s处到ω=10rad/s{\omega=10rad/s}ω=10rad/s处,L(ω){L(\omega)}L(ω)减小了10。而图中的红色和黄色两个三角形是相似的,假设从ω=1rad/s{\omega=1rad/s}ω=1rad/s处到ωc{\omega_c}ωc​处下降了x{x}x,根据比例关系很容易有:
lg(ωc)−lg(1)lg(10)−lg(1)=x100{\frac{lg(\omega_c)-lg(1)}{lg(10)-lg(1)}=\frac{x}{100}}lg(10)−lg(1)lg(ωc​)−lg(1)​=100x​
用这种方法求出
ωc=10lg(5)5=55=1.3797rad/s{\omega_c=10^\frac{lg(5)}{5}=\sqrt[5]{5}=1.3797rad/s}ωc​=105lg(5)​=55​=1.3797rad/s
与方法一的结论是相同的。

补充:
根据评论区carrot@的建议,最后“从ω=1rad/s{\omega=1rad/s}ω=1rad/s处到ω=ωcrad/s{\omega=\omega_c rad/s}ω=ωc​rad/s处”也可以这么做:
假设w=aw=aw=a时,L(a)=xL(a)=xL(a)=x,ωc\omega_cωc​在aaa的右侧。用d表示对数坐标系下横坐标到原点的长度,则有d1=lg(a),d2=lg(ωc)d_1=lg(a), d_2=lg(\omega_c)d1​=lg(a),d2​=lg(ωc​)。以ddd为自变量,可以写出直线lll的方程如图。

如上,求出y=0y=0y=0时的ωc\omega_cωc​表达式,然后代入xxx即可。

总结

从结果来看,两种方法得到的结论是相同的,且都和计算器求解的结果比较接近,但是方法二的求解过程明显更简单,做题更节省时间。

参考文献

[1] https://wenku.baidu.com/view/f6abcd52dd36a32d73758182.html

自动控制原理:如何用折线式伯德图计算截止频率?相关推荐

  1. 自动控制原理3.6---线性系统的稳态误差计算

    参考书籍:<自动控制原理>(第七版).胡寿松主编. <自动控制原理PDF版下载> 1.误差与稳态误差 控制系统的稳态误差,是系统控制准确度(控制精度)的一种度量,通常称为稳态性 ...

  2. 【经典控制理论】| 自动控制原理知识点概要(上)

    参考资料 胡寿松<自动控制原理>第6版 刘豹<现代控制理论>第3版 斐润<自动控制原理> 哈工大控制学科803考研资料 前言 对于控制工程专业的我来说,经典控制理论 ...

  3. c 语言与自动控制原理,自动控制原理C作业(第二章)答案

    <自动控制原理C作业(第二章)答案>由会员分享,可在线阅读,更多相关<自动控制原理C作业(第二章)答案(20页珍藏版)>请在人人文库网上搜索. 1.第二章 控制系统的数学模型2 ...

  4. 【经典控制理论】| 自动控制原理知识点概要(下)

    参考资料 胡寿松<自动控制原理>第6版 刘豹<现代控制理论>第3版 斐润<自动控制原理> 哈工大控制学科803考研资料 前文 [经典控制理论]| 自动控制原理知识点 ...

  5. 控制系统的滞后校正matlab,自动控制原理课程设计__用MATLAB进行控制系统的滞后校正设计...

    自动控制原理课程设计__用MATLAB进行控制系统的滞后校正设计 1滞后校正特性及校正方法 1.1滞后校正特性 滞后校正就是在前向通道中串联传递函数为的校正装置来校正控制系统,的表达式如下所示. (1 ...

  6. 关于汽车悬挂系统减震控制的进一步研究(自动控制原理课程设计小论文)

    摘要 随着社会的进步.科技的发展,人们生活水平不断提高,乘坐车辆的频率增多,对于汽车乘坐舒适性的要求也越来越高.另一方面,由于最近几年国内外汽车制造业的迅猛发展,汽车制造业的竞争也日益激烈,各个汽车制 ...

  7. 自动控制原理中的MATLAB函数以及用法总结

    MATLAB与自动控制原理 简记 最近要进行自控实验考试,在这里记录一下分析自控问题需要用的到一些函数以及用法,以供自己查阅,之后可能会总结一下如何用MATLAB来解自控的题. 一.数学模型的表示 建 ...

  8. 自动控制原理matlab pdf下载,第五章-matlab在自动控制原理中的应用.pdf

    第五章-matlab在自动控制原理中的应用.pdf 第第五五章章 MATLAB MATLAB MATLAB MATLAB在自动控制原理的应用在自动控制原理的应用 5 5 5 5 1 1 1 1 控制系 ...

  9. 自动控制原理7.5---离散系统的稳定性与稳态误差

    参考书籍:<自动控制原理>(第七版).胡寿松主编. <自动控制原理PDF版下载> 5.离散系统的稳定性与稳态误差 5.1 s s s域到 z z z域的映射 在 z z z变换 ...

最新文章

  1. MMX Intrinsics各函数介绍
  2. mysql系列之2.mysql多实例
  3. 全球及中国自媒体行业营销模式及应用规模前景分析报告2021-2027年
  4. 官方一步解决各种Windows更新问题
  5. java map扩容机制_java中ConcurrentHashMap的扩容机制是怎样的?详细解析
  6. 数据治理项目失败,90%都是被这29条骚操作搞垮的
  7. ADB工具 获取ROOT权限及复制文件方法
  8. C语言实现单例模式,以及使用内存屏障的性能优化方案
  9. 2. 创建多网站和商店(Creating Multiple Websites and Stores)
  10. 小D课堂 - 零基础入门SpringBoot2.X到实战_第11节 Logback日志框架介绍和SpringBoot整合实战_44、新日志框架LogBack介绍...
  11. 没有无线网络设备时如何共享无线网络
  12. Teraterm 脚本
  13. python-scrapy模拟登陆网站--登陆青果教务管理系统(二)
  14. linux开启nginx 443端口无效,http – Nginx正在侦听端口80或443但没有响应
  15. 外汇汇率接口 java_基于java的货币汇率接口调用代码实例
  16. CF1720A Burenka Plays with Fractions 题解
  17. Wireshark-----抓包分析
  18. intel固态硬盘误删文件该如何进行恢复
  19. 齐岳|聚乙二醇-四氧化三铁-二氢卟吩纳米复合物Fe3O4-PEG-Ce6|肝靶向功能的四氧化三铁Fe3O4纳米粒子
  20. IIS7用FastCGI运行PHP配置

热门文章

  1. 逻辑与计算机设计基础实验报告,学位论文_逻辑与计算机设计基础课实验报告.doc...
  2. python中查找文件当前位置定位,Python是怎样定位和修改文件读写位置的?
  3. 2022年 AI 技术成熟度曲线
  4. 满减折扣促销功能代码优化实战
  5. 给大家分享一篇 用Python抓取漫画并制作mobi格式电子书
  6. v8引擎和v12引擎_v8和v12引擎的区别是什么?
  7. html5在线裁剪,html5 Canvas裁剪图片
  8. Origami – 简洁轻快的WordPress主题
  9. 轻快步伐远不足以跟上轻快心情
  10. 谷俊丽:基于大数据的深度学习