控制系统仿真实验指导书.doc

精品文档控制系统仿真实验指导书 梅秀庄 主编内蒙古工业大学机械学院2009年12月2欢迎下载2欢迎下载。精品文档目 录实验一 MATLAB基本操作与基本运算1实验二 MATLAB 语言的程序设计3实验三 MATLAB 的图形绘制4实验四 函数文件设计和控制系统模型的描述6实验五 控制系统的频域与时域分析7实验六 控制系统PID 校正器设计法8附 录11精品文档实验一 MATLAB基本操作与基本运算 一、实验目的及要求： 1熟悉MATLAB7.0的开发环境； 2掌握MATLAB7.0的一些常用命令； 3掌握矩阵、变量、表达式的输入方法及各种基本运算。 二、实验内容： 1.熟悉MATLAB7.0的开发环境： MATLAB的各种窗口：命令窗口、命令历史窗口、启动平台窗口、工作空间窗口、当前路径窗口。 图1 MATLAB界面窗口路径的设置： 建立自己的文件夹，加入到MATLAB路径中，并保存。 设置当前路径，以方便文件管理。改变命令窗口数据的显示格式 format short format long然后键入特殊变量：pi （圆周率），比较显示结果。2.掌握MATLAB常用命令 who %列出工作空间中变量 whos %列出工作空间中变量，同时包括变量详细信息 save test %将工作空间中变量存储到test.mat文件中 load test %从test.mat文件中读取变量到工作空间中 clear %清除工作空间中变量 help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 lookfor %查找具有某种功能的函数但却不知道该函数的准确名称 如： lookfor Lyapunov 可列出与Lyapunov有关的所有函数。 3. 在MATLAB的命令窗口计算:1) 2) 4. 设计M文件计算: 已知求出： 1) a+b a-b a+b*5 a-b+I（单位阵） 2) a*b a.*b a/b3) a2 a.2 注意：点运算的功能，比较结果。5. 设计M文件计算: x=0:0.1:10当sum1000时停止运算，并显示求和结果及计算次数。三、实验报告要求： 1）体会1、2的用法；2）对3、4、5写出程序及上机的结果。 实验二 MATLAB 语言的程序设计一、实验目的1、熟悉 MATLAB 程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计方法3、函数文件的编写和设计4、了解和熟悉跨空间变量传递和赋值二、实验内容1、用 for 循环语句实现求1100 的和2、用 for 循环语句实现编写一个求n 阶乘的函数文件3、找到一个 n! 10100 的值（利用上题的n阶乘函数文件）4、编写一个乘法表（99）5、已知一维数组 A = 2,4,5,8,10、B =4,9,6,7,4，用for 循环语句实现，求和函数可用sum()6、编写验证魔方矩阵的函数文件，输出要求如下：（1） 如果输入矩阵的维数小于3，输出显示error（2） 如果输入矩阵的不是方阵，输出显示the size of matrix X mustbe N-by-N matrix（3） 显示行、列和及其对角线求和后的值，并判断其和是否相同。若不同，显示No,相同显示Yes。三、实验报告要求： 在M 文件编辑器中，编写程序代码并调试实验三 MATLAB 的图形绘制一、实验目的及要求：1.掌握MATLAB 绘图的基本方法，熟悉各种绘图函数的使用；2.掌握图形的修饰方法和标注方法；3.了解MATLAB 中图形窗口的操作。二、实验内容： x=2，2，y1=sinx、y2=cosx、y3=sin2x、y4=cos 2x用MATLAB语言分四个区域分别绘制的曲线，并且对图形标题及横纵坐标轴进行标注（如下图所示）。图2 四分区绘制曲线另建一个窗口，不分区，用不同颜色、线型绘出四条曲线，并标注图例注解。 图3 同一窗口绘制多条曲线(选做)在图形窗口可利用Figure窗口菜单提供的功能进行操作，并保存成* .fig文件。三、实验报告要求： 写出相应的的程序及上机结果。实验四 函数文件设计和控制系统模型的描述一、实验目的及要求： 1. 掌握函数文件的设计方法；2. 掌握控制系统数学模型的基本描述方法。二、实验内容： 1、函数文件设计： 设计一个函数文件实现一个阶乘运算 ，并设计程序调用该函数。 为保证函数的通用性，当输入负数或小数时，显示出错提示： disp(Input parameter must be a positive integer!) 提示：fix(x) 对零方向取整数 ceil(x) 对+方向取整数 round(x) 四舍五入取整数2、已知两个传递函数分别为：在MATLAB中分别用传递函数、零极点表示；在MATLAB中分别求出通过反馈、串联、并联后得到的系统模型；三、实验报告要求： 写出程序及上机的结果。实验五 控制系统的频域与时域分析一、实验目的：1、掌握控制系统数学模型的基本描述方法和相互转化2、了解控制系统的稳定性分析方法3、掌握控制系统频域与时域分析基本方法二、实验内容1、表示下列传递函数模型，并转化成其它的数学模型2、求下面系统在阶跃信号时系统的响应并求系统性能指标：稳态值、上升时间、调节时间、超调量3、已知两个单位负反馈系统中的开环传递函数分别为： ，分别用传递函数、零极点法表示开环模型；分别绘制开环系统的bode图和nyquist曲线；分别绘制闭环系统的阶跃响应曲线，分析静态误差；分别绘制闭环系统的单位斜坡输入响应曲线，分析静态误差。三、实验报告要求： 写出程序及上机的结果，并对结果进行分析。实验六 控制系统PID 校正器设计法一、实验目的1、熟悉常规 PID 控制器的设计方法2、掌握 PID 参数的调节规律3、学习编写程序求系统的动态性能指标二、 相关知识临界比例度法（边界稳定法）用系统的等幅振荡曲线来整定控制器的参数。先测出系统处于闭环状态下对象的等幅振荡曲线，根据等幅振荡曲线定出一些能反映控制对象动态特性的参数，具体做法是将比例增益K(或比例度=1/K)调在比较小的位置上（对应为比较大位置上），逐渐增大K值（或逐渐减小），直到出现等幅振荡曲线，此时的比例增益为Km，称为临界比例增益， 称为临界比例度。从振荡曲线上读出临界周期Tm。根据得到的Km（或）、Tm两个参数，利用下表来计算控制器的控制参数。控制器类型控制器的控制参数比例增益Kp比例度积分时间Ti微分时间TdP0.5Km20PI0.45Km2.2Tm/1.20PID0.6Km1.70.5Tm0.125Tm三、实验内容1、在SIMULINK 窗口建立如下页模型。2、设计PID控制器，传递函数模型如下 3、修改PID参数讨论参数对系统的影响4、利用临界比例度法(即：稳定边界法)对PID参数校正设计。参 数KmTmKpTiTd数 值5、根据PID参数对系统的影响，调节PID参数实现系统的超调量小于10%。参数KpTiTd数值6、通过程序求得系统的超调量、上升时间和调节时间(误差带选为5%)四、实验报告要求： 写出程序及上机的结果，并对结果进行分析。附 录1、MATLAB基本函数和语言 MATLAB的原意为Matrix Laboratory，即矩阵实验室，是由在数值线性代数领域颇有影响的Cleve Moler博士首创的。后来由Moler博士等一批数学家和软件专家组建了MathWorks软件公司，专门从事MATLAB的扩展与改进。 MATLAB不仅具有强大的数值计算能力，而且具有数据图示功能和符号运算功能。特别是大量的工具箱，扩展了应用领域，是高校学生、教师、科研人员和工程计算人员的最好选择，是数学建模必不可少的工具。一、变量及其赋值（1）MATLAB中的变量不需要事先定义，在遇到新的变量名时，MATLAB会自动建立该变量并分配存储空间。对变量赋值可采用赋值语句：变量名表达式（或常数）例如：(2)利用冒号表达式生成向量，基本形式为：x=x0:step:xn如：a=1:2:12则：a = 1 3 5 7 9 11(3)基本矩阵：zeros 全零矩阵(mn阶)，ones 全一矩阵(mn阶)，eye(n) 单位矩阵(方阶)，length 一维矩阵的长度，i、j虚数单位。MATLAB中提供了一些用户不能清除的固定变量，最常用的固定变量为pi：即3.1415二、基本运算1、算术运算：+ 加 - 减 * 乘 / 除 左除 乘方 2、矩阵运算：(1) 矩阵的转置 A 求矩阵A的转置。(2) 矩阵的逆运算 inv(A) 求方阵A的逆矩阵。(3) 矩阵的行列式运算 det(A) 求方阵A 的行列式的值。(4) 矩阵的指数运算，即eA 矩阵的指数运算的最常用的命令为expm（A）(5) 矩阵的对数运算 logm(A) 求方阵A的对数。(6) 矩阵的开方运算，即 sqrtm(A) 求矩阵A的算术平方根。3、 数组运算： 4、常用函数：函数名功能函数名功能sin正弦函数expe的指数函数asin反正弦函数log自然对数函数cos余弦函数log10常用对数函数acos反余弦函数log2以2为底的对数函数tan正切函数pow2以2为底的幂函数atan反正切函数sqrt平方根函数abs模函数imag复矩阵虚部函数angle相角函数real复矩阵实部函数conj复共轭函数三、逻辑判断逻辑运算是数组运算所特有的一种运算形式，包括逻辑关系运算和逻辑函数运算。结果为真用1表示，结果为假用0表示。基本逻辑运算表符号运算符功能函数名=等于eq=不等于ne大于gt=大于等于ge&逻辑与and|逻辑或or逻辑非not四、基本绘图方法plot(t,y) plot(t,y,*) title(图1) xlabel(t(s) xlabel(V(mv) grid on grid offhold hold off subplot(n,m,p) gtext() stem(t,y) stairs(t,y) bar(t,y)五、流程控制语句(1) 选择语句if语句if表达式 语句组 A，endif 表达式 语句组 A，else 语句组B，end多分支(switch)语句switch 表达式 case数值1, 语句组A, otherwise语句组B, end(2)循环语句 for循环for k= 初值:增量:终值 语句组A，end while循环while 表达式 语句组A，end2、常用矩阵函数命令说明d=eig(A)矩阵特征值v,d=eig(A)矩阵特征值和特征向量det(A)行列式计算inv(A)求逆poly(A)特征多项式trace(A)对角元素之和orth(A)正交化3、常用特殊矩阵命令说明A= 空矩阵A=eye( n )N维单位矩阵A=ones(n,m)全部元素都为1的矩阵A=zeros(n,m)全部元素都为0的矩阵A=rand(n,m)元素为0到1之间均匀分布的随机矩阵A=randn(n,m)元素为零均值单位方差正态分布的随机矩阵4、图形绘制和修饰（1）坐标轴调整axis(xmin xmax ymin ymax) 图形的x轴范围限定在xmin xmax y轴范围限定在ymin ymaxaxis(控制字符串） 根据字符串控制图形set(gca,xtick,标示向量) 按照向量设置X坐标的刻度标示set(gca,ytick,标示向量) 按照向量设置Y坐标的刻度标示set(gca,xticklabel，字符串|字符串) 按字符串设置X坐标的刻度标示set(gca,yticklabel，字符串|字符串) 按字符串设置Y坐标的刻度标示 （2）文字标示title(字符串) 图形标题xlabel(字符串) x轴标注ylabel(字符串) y轴标注 text(x,y,字符串) 在坐标处标注说明文字gtext(字符串) 用鼠标在特定处标注说明文字（3）网格控制grid on在所画的图形中添加网格线grid off在所画的图形中去掉网格线（4）图例注解legend(字符串1,字符串2,参数) 为区分各条曲线以字符串进行注解说明，参数确定注解在图形中的位置。（5）图形的保持hold on保持当前图形及轴系的所有特性hold off解除命令（6）图形窗口的分割 可把一个绘图窗口分割成几个区域，分别绘图。subplot(m，n，p) 将当前绘图窗口分割成m行n列区域，指定第p个编号区域为当前绘图区。区域编号原则：“先上后下，先左后右”。5、线性系统分析设计 （1）模型描述状态空间：写为A、B、C、D四个矩阵的形式；也可用SYS = SS(A,B,C,D) 传递函数：用分子、分母的多项式表示；可用SYS = TF(NUM,DEN)建立tf模型零极点：用增益k、分子零点向量z、分母极点向量p表示，同样可用 SYS = ZPK(Z,P,K)建立zpk模型部分分式：传递函数没有相同极点时，可与部分分式相互转换： r,p,k=residue(num,den) 传递函数变换为部分分式形式 num,den=residue(r,p,k) 部分分式变换为传递函数形式闭环系统feedback 反馈连接 SYS = feedback(SYS1,SYS2,sign) A,B,C,D=feedback(A1,B1,C1,D1,A2,B2,C2,D2,sign)num,den=feedback(num1,den1,num2,den2,sign) sign=+(-)1反馈极性,缺省-1为负反馈类似地，有 parallel系统并联series系统串联（2）模型转换ss状态空间、tf传递函数、zp零极点 num,den=ss2tf(a,b,c,d,iu) 状态空间到传递函数 z,p,k=ss2zp(a,b,c,d,iu) 状态空间到零极点 a,b,c,d=tf2ss(num,den) 传递函数到状态空间 z,p,k=tf2zp(num,den) 传递函数到零极点 a,b,c,d=zp2ss(z,p,k) 零极点到状态空间 num,den=zp2tf(z,p,k) 零极点到传递函数 r,p,k=residue(