- 1、本文档共45页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第六章Matlab在自动控制理论中应用
Matlab在自动控制理论中的应用;MATLAB语言基础 ;1.矩阵运算功能;3.丰富的绘图功能与计算结果的可视化;5.丰富的MATLAB工具箱;6.MATLAB的开放式可扩充结构 ;应用MATLAB进行模型处理;1、拉氏变换与反变换 ;拉氏反变换
“iLaplace”的调用格式如下:
F=iLaplace(L):是缺省独立变量s的关于符号向量L的拉氏
反变换,缺省返回关于t的函数。
F=iLaplace(L,y):是一个关于y代替缺省t项的拉氏变换。
F=iLaplace(L,y,x):是一个关于x代替缺省t项的拉氏变换。;2 、传递函数 ;例 已知某一系统的微分方程如下,试求其传递函数。;例 将传递函数; 零极点模型
传递函数的分子和分母均为因式的形式称为零极点模型,如下式所示。
;例 将传递函数;3 、控制系统的特征根及零极点图 ;零极点图 ;例 某系统开环传递函数如下,试求其零、极点图。;4 、 控制系统的方框图模型 ;用MATLAB进行动态响应分析;例2-1:用Matlab绘制系统;例2-2:求系统;2、 阶跃响应性能分析 ;例2-3: 用Matlab求系统;已知单位负反馈控制系统的开环传递函数为; 例3-4:控制系统结构图如下图所示,试在Simulink环境下构建系统方框图,并对系统的阶跃响应进行仿真。;用MATLAB绘制系统的根轨迹; 如果引入左端变量,即
[r, K]=rlocus (num, den, K)
则r向量中元素为指定K值闭环极点的位置。如果不指定K,即
[r, K]=rlocus (num, den)
则r矩阵的行数与K向量的元素相同,列数与极点的个数相同。r中元素分别为各对应K值的闭环极点。这时,使用绘图命令
plot(r, ‘’)也可绘制系统根轨迹。
如果希望在绘制根轨迹时标上符号“o”或“x”,则需要采用下列绘图命令
plot(r, o) 或 plot(r, x)
这时,每一个计算的闭环极点都被图解表示出来。 ;例3-1 已知系统的开环传递函数为;例3-2 已知系统的开环传递函数;a=[1 4 0];
b=[1 2 2];
den=conv(a,b);
num=[1 4 16];
K1=0:0.2:20;
K2=20:0.1:30;
K3=30:5:1000;
K=[K1 K2 K3];
r=rlocus(num,den,K)
plot(r,o)
v=[-4 0.5 -15 15];
axis(v);
grid
title(‘Root Locus Plot of G(s)H(s)=K(s^2+4s+16)/[s(s+4)(s^2+2s+2)]);MATLAB频率特性分析;MATLAB中绘制Bode图的函数是bode( ),调用格式为
bode (num,den)
bode (num,den, w )
[mag, phase, w]=bode (num,den)
[ mag, phase ]=bode (num,den, w)
Bode(num,den)可以绘制传递函数为; Bode函数具有自动频率选择功能,函数的输入变量部分未给出频率的范围,则该函数能根据系统模型的特性自地选择频率的变化范围。若需要人为地指定频率范围或频率点 w,可以在函数的输入变量部分包含所定义 w。 w的定义可以采用logspace函数,其格式为W=logspace(a,b,n) 其中,a表示最小频率10a;b表示最大频率10b;n表示10a~10b之间的频率点数。; 在MATLAB中可采用函数 margin( )来求取相对稳定性,调用格式为:
[Gm,Pm,wcg,wcp]=(num,den)
式中,Gm为幅值裕度;Pm为相角裕度;wcg为穿越 线 ;2、 Nyquist图;例如;3、 Nichols 图;MATLAB 离散控制系统分析;R(s);2、采样系统的响应 ;R(s);系统在采样时刻的单位脉冲响应为:
num=[0.368 0.264];
den=[1 -1 0.632];
dimpulse(num,den);3、Z变换和Z反变换;例:计算下列函数的Z变换。
(1) (2) (3)
(4) (5);例:已知连续函数x(t)拉氏变换为 ,
您可能关注的文档
- 第03章 春秋战国时代城市.ppt
- 第9课 北美大陆上新体制(公开课).ppt
- 在EN15085标准下焊接供应商管理体系建立.pdf
- 第一章第三节 饭店分类与等级.ppt
- 第一讲 音与音名.ppt
- 钢筋溷凝土梁建模与分析过程.doc
- 第七章 重大错报风险评估过程.ppt
- 2017届初中学业水平测试模拟试卷2.doc
- 金蝶认证与部分答案.pdf
- 高成形密度粉末冶金零件用润滑剂开发.pdf
- 金融产品2024年投资策略报告:积极适应市场风格,行为金融+机器学习新发现.pdf
- 交运物流2024年度投资策略:转型十字路,峰回路又转(2023120317).pdf
- 建材行业2024年投资策略报告:板块持续磨底,重点关注需求侧复苏.pdf
- 宏观2024年投资策略报告:复苏之路.pdf
- 光储氢2024年投资策略报告:复苏在春季,需求的非线性增长曙光初现.pdf
- 公用环保2024年投资策略报告:电改持续推进,火电盈利稳定性有望进一步提升.pdf
- 房地产2024年投资策略报告:聚焦三大工程,静待需求修复.pdf
- 保险2024年投资策略报告:资产负债匹配穿越利率周期.pdf
- 政策研究2024年宏观政策与经济形势展望:共识与分歧.pdf
- 有色金属行业2024年投资策略报告:新旧需求共振&工业原料受限,构筑有色大海星辰.pdf
文档评论(0)