MATLAB自动控制原理相关编程

　　传递函数表示
　　显示多项式传递函数
　　matlab语句：分子分母的各项式系数按照降次排序分别计入num〔〕，den〔〕中，缺项系数补零。传递函数systf（num，den）
　　程序clc；clearall；closeall；num〔13406〕；den〔531617〕；systf（num，den）Gtf（num，den）
　　结果
　　因子形式的传递函数
　　matlab语句：利用conv（）实现两个多项式降次系数乘积运算。
　　程序clc；clearall；closeall；numconv（〔13〕，〔176〕）；denconv（〔11〕，〔1305〕）；systf（num，den）Gtf（num，den）
　　结果
　　显示零极点式传递函数
　　matlab语句：zpk（）函数用法。z〔〕输入零点，p〔〕输入极点；k〔〕输入增益。syszpk（z，p，k）；
　　程序clc；clearall；closeall；z〔3〕；p〔245〕；k7；syszpk（z，p，k）
　　结果
　　直接传递函数输入matlab语句：stf（‘s’）；
　　程序clc；clearall；closeall；stf（s）；sys1（s（s1）（s2））
　　结果
　　传递函数多项式形式和零极点形式的转换
　　1零极点形式转换为多项式形式
　　zp2tf函数
　　程序clc；clearall；closeall；z〔7，2〕；p〔3，5，9〕；k4；spkzpk（z，p，k）〔num，den〕zp2tf（z，p，k）stftf（num，den）
　　2多项式形式转换为零极点形式
　　tf2zp（）函数
　　程序clc；clearall；closeall；num〔4，36，56〕；den〔1，17，87，135〕；stftf（num，den）〔z，p，k〕tf2zp（num，den）；spkzpk（z，p，k）
　　结果
　　求系统的单位冲击响应和单位阶跃响应
　　1单位冲激响应
　　impulse函数
　　impulse函数可以求得系统的单位冲击响应，参数为sys和t，其中sys为系统对应的微分方程，t为持续时间。
　　sys变量由tf函数生成，其参数为输入部分的方程系数矩阵和响应部分的方程系数矩阵。使用格式：impulse（sys）；impulse（sys，t）；
　　程序clc；clearall；closeall；systf（〔2，1〕，〔1，4，3〕）；t0：0。1：10；yimpulse（sys，t）；plot（t，y）；figure；impulse（sys）
　　结果
　　2单位阶跃响应
　　step函数可以求得系统的单位阶跃响应，其用法与impulse函数类似。
　　step：动态系统的阶跃响应绘图。具体描述：1。step：计算一个动态系统的阶跃响应。在状态空间的情况下，假定初始状态为零。当它没有输出参数时，调用这个函数在屏幕上的画出阶跃响应。2。step（sys）：画出任意一个动态系统模型sys的阶跃响应。这个模型可以是连续的或离散的，和单输入单输出或多输入多输出。多端输入系统的阶跃响应对于每一个输入方式来说都是阶跃响应的集合。模拟持续的时间根据系统的极点和零点自动地被确定。3。step（sys，Tfinal）：模拟了系统sys从时间t0到tTfinal的阶跃响应。在系统时间单位中，表达式Tfinal在sys的时间单位属性中是被指定的。对于未指定采样时间（Ts1）的离散时间系统，阶跃将Tfinal作为采样周期的数量来模拟。4。step（sys，t）：使用用户提供的时间矢量t来模拟。在系统时间单位中，表达式t在sys的时间单位属性中是指定的。对于离散时间模型，t应该是Ti：Ts：Tf，这里Ts是采样时间。对于持续的时间模型，t应该是Ti：dt：Tf，这里dt变成近似于持续系统的一个离散采样。阶跃命令总是适用于在t0时的阶跃输入，而忽视Ti。5。在一个图上画出一系列模型sys1，sys2，。。。sysN的阶跃响应：step（sys1，sys2，。。。，sysN）；step（sys1，sys2，。。。，sysN，Tfinal）；step（sys1，sys2，。。。，sysN，t）；所有将被画在一个图上的系统必须有相同数量的输入和输出。当然也可以在一个单一的图上画连续时间系统和离散时间系统的混合。这个语法对于比较多个系统的阶跃响应是有用的。6。你也可以给每个系统指定一个独特的颜色，线型，标记，或者三者都有。例如：step（sys1，y：，sys2，g）；当调用输出参数时：ystep（sys，t）；〔y，t〕step（sys）；〔y，t〕step（sys，Tfinal）；〔y，t，x〕step（sys）；step返回输出响应y，模拟时间向量t（如果没有作为一个输入参数被提供）和状态轨迹x（仅为状态空间模型）。在屏幕上没有绘图生成。对于单输入系统来说，y有和采样时间一样多的行，和输出一样多的列。在多端输入的情况下，每个输入方式的阶跃响应沿着输出y的第三个维度被累积起来。
　　程序clc；clearall；closeall；systf（〔2，1〕，〔1，4，3〕）；t0：0。1：10；ystep（sys，t）；plot（t，y）；figure；step（sys）
　　结果
　　3求系统的阶跃响应动态指标
　　实例程序clc；清除命令行空间clearall；清除变量closeall；关闭图像num6；传递函数分子denconv（conv（〔10〕，〔0。11〕），〔0。31〕）；systf（num，den）开环传递函数syscfeedback（sys，1，1）系统闭环传递函数figure；新建一个窗口bode（sys）；画Bode图margin（sys）画出bode所有参数figure；新建一个窗口step（sysc）闭环系统阶跃响应曲线〔y，t，x〕step（sysc）；mpmax（y）；上升时间tr0。518峰值时间tpspline（y，t，mp）cslength（t）；稳态值yssy（cs）稳态误差值ess1yss超调量ct（mpyss）yss
　　结果
　　零极点分布图、根轨迹图、奈奎斯特曲线图、bode图
　　1、零极点分布图和根轨迹图
　　零极点分布图pzmap（）函数
　　根轨迹图rlocus（）函数
　　结果
　　2奈奎斯特曲线图
　　奈奎斯特曲线图nyquist（）函数
　　3伯德图伯德图bode（）和margin（）函数〔Lg，Y，Wg，Wc〕margin（G）四者分别为幅值裕度、相位裕度、相位穿越频率、截止频率
　　程序clc；clearall；closeall；num〔13406〕；den〔531617〕；systf（num，den）；Gtf（num，den）figure；pzmap（G，r）画图figure；rlocus（G）figure；nyquist（G）figure；bode（G）margin（G）；〔Lg，Y，Wg，Wc〕margin（G）四者分别为幅值裕度、相位裕度、相位穿越频率、截止频率相位裕度：对于闭环稳定系统，如果系统开环相频特性再滞后Y度，则系统将处于临界稳定状态相位穿越频率：系统开环频率特性上相位等于180度时所对应的角频率称为相位穿越频率截止频率：系统开环频率特性上幅值为1时所对应的角频率称为幅值穿越频率或截止频率
　　结果如上四个图
　　本文内容来源于网络，仅供参考学习，如内容、图片有任何版权问题，请联系处理，24小时内删除。
　　作者郭志龙
　　编辑郭志龙
　　校对郭志龙