MATLAB教程及实训曹弋课件8章节

《MATLAB教程及实训曹弋课件8章节》由会员分享，可在线阅读，更多相关《MATLAB教程及实训曹弋课件8章节（31页珍藏版）》请在金锄头文库上搜索。

1、第8章 线性控制系统的分析,8.1 控制系统的数学模型 8.2 时域分析的MATLAB实现 8.3 频域分析的MATLAB实现 8.4 根轨迹分析的MATLAB实现 8.5 稳定性分析的MATLAB实现 8.6 稳态误差分析的MATLAB实现 8.7 线性定常系统分析与设计的图形工具 8.8 综合实例介绍,8.1 控制系统的数学模型 8.1.1 创建系统的模型并相互转换,1. 传递函数模型 sys=tf(num,den,Ts) %由分子分母得出传递函数 sys=tf(num,den,Ts,Property1,v1,Porperty2,v2,) %创建传递函数并设置属性 2. 零极点增益模型 G

2、=zpk(z,p,k,Ts) %由零点、极点和增益创建模型 G=zpk(z,p,k,Property1,v1,Porperty2,v2,) %创建模型并设置属性,3. 状态方程模型 G=ss(a,b,c,d,Ts) %由a、b、c、d参数创建模型 G=ss(a,b,c,d,Ts,Property1,v1,Porperty2,v2,) %创建模型并设置属性,4. 模型的转换 5. 连续系统与离散系统模型的转换 MATLAB控制工具箱提供了c2d、d2c和d2d函数实现连续系统和离散系统的相互转换。,8.1.2系统的模型参数,2. 获取模型的参数 （1）获取模型参数的函数 tfdata、zpkda

3、ta和ssdata、dssdata函数分别用来获取传递函数模型、零极点增益模型和状态方程模型的参数，这些函数名都是在创建模型的函数名后面加“data”。 （2）获取模型尺寸的函数 n=ndims(sys) %获取模型的维数 d=size(sys,n) %获取模型的参数 d=size(sys,order) %获取模型的阶数,3. 使用get和set函数 （1）set函数 set(sys,property1,value1,property2,value2,) %设置系统属性 （2）get函数 value=get(sys,property) %获取当前系统的属性,8.1.3系统模型的连接和简化,1.

4、 串联环节 G=series(G1,G2,outputs1,inputs1) %计算串联模型 2. 并联环节 G=parallel(G1,G2,in1,in2,out1,out2) %计算并联模型,u(t),3. 反馈环节 G=feedback(G1,G2,feedin,feedout,sign) %计算反馈模型 说明：G1和G2模型必须都是连续系统或采样周期相同的离散系统；sign表示反馈符号，当sign省略或=-1时为负反馈；feedin和feedout分别是G2的输入端口和G1的输出端口，可省略，当省略时G1与G2端口正好对应连接。,4. 复杂模型的连接 （1）对框图中的每个环节进行编号

5、并建立它们的对象模型，环节是指一条单独的通路； （2）建立无连接的状态空间模型，使用append函数实现，append的命令格式如下： G=append(G1,G2,G3,) （3）写出系统的联接矩阵Q Q是联接矩阵，第一列是各环节的编号，其后各列是与该环节连接的输入通路编号，如果是负连接则加负号； （4）列出系统总的输入和输出端的编号，使用inputs列出输入端编号，outputs列出输出端的编号； （5）使用connect函数生成组合后系统的状态空间模型，connect函数的命令格式为： Sys=connect(G,Q,inputs,outputs),例8-7 根据图8-6所示的模型结构框

6、图计算模型的总传递函数，其中R1=1，R2=2，C1=3，C2=4。,练习,8.1.4 将Simulink模型结构图转化为系统模型,MATLAB提供了linmod和linmod2函数命令将Simulink模型转换为数学模型。 例8-8 根据图8-6的模型结构框图在Simulink环境中创建系统模型，使用函数命令转化为传递函数。,8.2 时域分析的MATLAB实现 8.2.1 使用拉氏变换和逆变换计算时域响应,例8-9 使用拉氏变换和逆变换计算输入信号为阶跃信号和脉冲信号的系统输出响应，已知系统的传递函数为 。 syms s t; rs=laplace(heaviside(t); %阶跃信号 f

7、s=1/(0.5*s2+s+1); cs=rs*fs; ct=ilaplace(cs) %阶跃响应表达式 ct = -exp(-t)*cos(t)-exp(-t)*sin(t)+1,8.2.2线性系统的时域分析,1. 阶跃响应 step(G, T) %绘制系统G的阶跃响应曲线 y,t,x=step(G, T) %得出系统G的阶跃响应数据 2. 脉冲响应 系统的脉冲响应使用impulse函数命令绘制，命令格式与step函数相同。 3. 斜坡响应和加速度响应 斜坡响应=阶跃响应*s 加速度响应=阶跃响应*s2,4. 任意输入响应 lsim(G,U,T) %绘制系统G的任意响应曲线 y,t,x=ls

8、im(G,U,T) %得出系统G的任意响应数据 5. 零输入响应 initial(G,x0, T) %绘制系统G的零输入响应曲线 y,t,x=initial(G,x0, T) %得出系统G的零输入响应的数据 6. 离散系统响应,例8-12 使用lsim函数绘制正弦响应曲线，系统的传递函数为，绘制的曲线如图8-11所示。 t=0:0.1:10; u=sin(t); G=tf(1,2 1) Transfer function: 1 - 2 s + 1 lsim(G,u,t) %正弦响应,8.2.3 线性系统的结构参数与时域性能指标,1. 线性系统的结构参数 （1）pole和zero p=pole(

9、G) %获得系统G的极点 z=zero(G) %得出系统G的零点 （2）pzmap pzmap(G) %绘制系统的零极点分布图 p,z=pzmap(G) %获得系统的零极点值 （3）damp wn,zeta,p=damp(G) %获得G的阻尼系数、固有频率和极点 （4）sgrid sgrid(z,wn) %绘制s平面网格并指定z和wn,2. 时域分析的性能指标 在自动控制原理中，时域分析常用的系统性能指标有超调量p 、上升时间tr、峰值时间tp和过渡时间ts，通过性能指标来分析系统暂态性能的稳定性。 例： detap=exp(-pi*zeta/sqrt(1-zeta2)*100 %计算超调量

10、tr=(pi-acos(zeta)/(wn*sqrt(1-zeta2) %计算上升时间,8.3 频域分析的MATLAB实现 8.3.1 线性系统的频域分析,1. 线性系统的幅频特性和相频特性 Gw=freqresp(G,w) %计算w处系统G的频率特性 2. nyquist曲线 nyquist(G,w) %绘制系统G的nyquist曲线 3. bode图 bode(G,w) %绘制系统G的bode图 mag,pha=bode(G,w) %得出w处的幅值和相角 4. nichols图,例8-19 根据传递函数绘制系统的bode图和对数幅频特性曲线，已知系统的传递函数分别为 。 G=tf(10,c

11、onv(1 1,5 1 0); subplot(1,2,1); bode(G) %绘制系统的bode图 subplot(1,2,2); bodemag(G) %绘制对数幅频特性曲线,8.3.2 频域分析性能指标,1. 开环频率特性的相角域度和幅值域度 margin(G) %绘制bode图并标出幅值裕度和相角裕度 Gm,Pm,Wcg,Wcp=margin(G) %得出幅值裕度和相角裕度和相应的频率 例： Gm,Pm,Wcg,Wcp=margin(G) Gm = 0.1200 Pm = -38.7743 Wcg = 0.4472 Wcp = 1.1395,8.4 根轨迹分析的MATLAB实现 8.

12、4.1 线性系统的根轨迹分析,1. 绘制根轨迹 rlocus(G) %绘制根轨迹 r,k=rlocus(G) %得出系统G的闭环极点r和增益k 2. 得出给定根的根轨迹增益 rlocfind函数可以获得根轨迹上给定根的增益和闭环根。 k,p=rlocfind(G) %得出根轨迹上某点的闭环极点和k,8.4.2 根轨迹设计工具,MATLAB控制工具箱还提供了根轨迹设计器，根轨迹设计器是一个分析系统根轨迹的图形界面。 rltool(G) %打开系统G的根轨迹设计器,8.5 稳定性分析的MATLAB实现 8.5.1 根据闭环特征方程判定系统稳定性,1. 使用特征根判定系统的稳定性 可以通过roots

13、函数计算闭环传递函数分母的根，也可以使用pzmap显示系统的闭环零极点分布图或获取系统的闭环零极点，查看是否在s的左半平面，从而判断系统的稳定性。 2. 使用劳斯判据判定系统的稳定性 对于闭环特征方程，劳斯判据是根据闭环特征方程的系数a0、a1的关系来判定系统稳定性的，对于二阶系统a0、a1、a2必须大于0；对于三阶系统首先满足系数都大于0，然后满足a1*a2a0*a3。,8.5.2 用频率特性法判定系统稳定性,1. 使用bode图判定系统稳定性 在bode图上判定系统稳定性可以针对开环SISO的连续或离散系统进行，稳定判据是在bode图上查看在Lg()0的范围内相频特性没有穿越-180度线，

14、即0时系统稳定。 例： Gm,Pm,Wcg,Wcp=margin(G) Gm = Inf Pm = 101.4434 Wcg = Inf Wcp = 0.0102,8.5.3 用根轨迹法判定系统稳定性,根轨迹是开环系统中某个参数从0变化时，闭环根的移动轨迹，根轨迹在s的左半平面时系统稳定。判定稳定的方法是先用rlocus函数绘制系统的闭环根轨迹，然后使用rlocfind函数找出临界稳定点时的对应参数值。,8.6 稳态误差分析的MATLAB实现,稳态误差是系统稳定误差的终值，稳态误差可表示为，在自动控制原理中使用位置误差系数、速度误差系数和加速度误差系数来计算稳态误差， ， ， 。 MATLAB

15、没有提供专门的计算函数，可以使用求极限的limit函数来计算稳态误差。,练习,8.7 线性定常系统分析与设计的图形工具 8.7.1 线性定常系统仿真图形工具LTI Viewer,LTI Viewer是线性定常系统仿真的图形工具，可以方便地获得系统的各种时域响应和频率特性等曲线，并可以得到系统的性能指标。,8.7.2 SISO设计工具sisotool,sisotool(views,G,C,H,F) %打开SISO设计工具 说明： views是指定SISO设计工具窗口的初始显示图形，可以是一个或多个图形，rlocus是根轨迹图，bode是开环系统的bode图，nichols是nichols图，filter是F的bode图和从F输入到G输出的闭环响应，views参数可以省略，省略时指所有图形都显示； G是前向通道的系统模型，可以是传递函数、零极点形式和状态空间模型； C是前向通道中串联的补偿器模型； H是反馈通道的模型； F是预滤器的模型；这些参数都可以省略，省略时显示空白界面。,8.8 综合实例介绍,例8-30 使用超前校正环节来校正系统，已知系统的开环传递函数为，要求校正后系统的速度误差系数小于10，相角域度为450。,根据系统框图得出系统总传递函数，系统框图如图8-35所示，其中 ， ， ， ， 。,，,，,，,。,G1=tf(1,1 2 5); G2=t