基于MATLAB的控制系统单位阶跃响应分析

电子科技大学学院学生实验报告

学院：机电工程学院专业：17自动化课程名称：自动控制原理实验与仿真

性能指标

Pole Damping Frequency Time Constant

(rad/seconds) (seconds)

-1.50e+00 + 9.89e+00i 1.50e-01 1.00e+01 6.67e-01 -1.50e+00 - 9.89e+00i 1.50e-01 1.00e+01 6.67e-01

2.=n n ζωωζ当0，0.25，0.5，0.75，1，1.25时，对应系统的闭环极点和自然振荡频率见表，编程求取对应系统的阶跃响应曲线，并分析一定时，变化对系统性能的影响。。

ζ

闭环极点

/(/)n rad s ω

阶跃响应曲线 =0ζ j ± 10 等幅振荡 =0.25ζ 2.59.68j -±

10 衰减振动 =0.5ζ

58.66j -± 10 衰减振动 =0.75ζ 7.5 6.61j -±

10 衰减振动 =1ζ

两实重根-10 10 单调上升 =1.25ζ

两不等实根 -5和-20

5,20

单调上升

曲线：

结论：可见当/(/)n rad s ω一定时，系统随着阻尼比ξ增大，闭环极点的实部在S 左半平面的位置更加远离原点，虚部减小到0，超调量减小，调节时间缩短，稳定性更好。

0.25,10,30,50n n ζωζ==3.时，对应点的单位阶跃响应曲线并分析不变，对系统性能的影响。

曲线：

结论：可见当ξ一定时，随着/(/)n rad s ω增大，系统响应加速，振荡频率增大，系统调整时间缩短，但是超

结论：闭环零点使得超调量增大，峰值时间前移，且加入的零点越靠近虚轴，影响越明显

1212

225020

2-5,-2,(5)(210)(2)(210)

====++++++）系统分别增加极点，（）

（）p p G s G s s s s s s s 代码及曲线：

代码：sys=tf(10,[1 2 10]);

step(sys) hold on num=50;

den=conv([1,5],[1 2 10]); sysc=tf(num ,den);

step(sysc) hold on num1=20;

den1=conv([1,2],[1 2 10]); sysx=tf(num1,den1);

step(sysx) hold off

title('单位阶跃系统增加极点比较') lab1='原系统G(s)';text(1,1.3,lab1) lab2='增加极点-2';text(0.25,1.1,lab2) lab3='增加极点-5';text(0.95,0.771,lab3)

结论：

闭环极点使得超调量减小，峰值时间后移，且加入的极点越靠近虚轴，影响越明显。

四.实验作业要求

（1）完成实验容中的实验，编写程序，记录相关数据，并分析，得出结论。 （2）总结闭环零极点对系统阶跃响应影响的规律。

五、 拓展思考

2

(1)(5)

()1()*1()()(2)(1)

++=+=

+++当输入信号为时，求系统的输出响应曲线，并测出动态性能指标。s s u t t t t G s s s s 代码、曲线及性能指标：

代码：

控制系统时间响应分析”实验报告

控制系统时间响应分析”实验报告

实验一、“控制系统时间响应分析”实验报告 一、实验类型 验证性实验 二、实验目的 1、 求系统在时间常数τ不同取值时的单位脉冲、单位阶跃响应和任意输入响应，熟悉系统时间响应的定义和图形曲线 2、 求系统的上升时间、峰值时间、最大超调量和调整时间等性能指标，熟悉系统瞬态性能指标的定义。 三、实验仪器与设备（或工具软件） 计算机，MATLAB 软件 四、实验内容、实验方法与步骤 已知系统传递函数 50 )1(05.050)(2+++=s s s G τ 1、求系统在时间常数τ不同取值时的单位脉冲、单位阶跃响应和任意输入 响应。 应用impulse 函数，可以得到τ=0，τ=0.0125、τ=0.025时系统单位脉冲响 应；应用step 函数，同样可以得到τ=0，τ=0.0125、τ=0.025时系统单位阶跃响应。 2、求系统的瞬态性能指标 五、实验结果 1、系统在时间常数τ不同取值时的单位脉冲、单位阶跃响应和任意输入响 t=[0:0.01:0.8];%仿真时间区段 nG=[50]; tao=0; dG=[0.05 1+50*tao 50]; G1=tf(nG,dG); tao=0.0125; dG=[0.05 1+50*tao 50]; G2=tf(nG,dG); tao=0.025; dG=[0.05 1+50*tao 50]; G3=tf(nG,dG);%三种τ值下，系统的传递函数模型 [y1,T]=impulse(G1,t);[y1a,T]=step(G1,t); [y2,T]=impulse(G2,t);[y2a,T]=step(G2,t);

二阶系统的阶跃响应及频率特性

实验二二阶系统的阶跃响应及频率特性 实验简介：通过本实验学生能够学习二阶系统的频率响应和幅频特性的测试方法，对实验装置和仪器的调试操作，具备对实验数据、结果的 处理及其与理论计算分析比较的能力。 适用课程：控制工程基础 实验目的：A 学习运算放大器在控制工程中的应用及传递函数的求取。 B 学习二阶系统阶跃响应曲线的实验测试方法。 C 研究二阶系统的两个重要参数ζ、ω n 对阶跃瞬态响应 指标的影响。 D 学习频率特性的实验测试方法。 E 掌握根据频率响应实验结果绘制Bode图的方法。 F 根据实验结果所绘制的Bode图，分析二阶系统的主要 动态特性（M P ，t s ）。 面向专业：机械类 实验性质：综合性/必做 知 识 点：A《模拟电子技术》课程中运算放大器的相关知识； B《数字电子技术》课程中采样及采样定理的相关知识； C《机械工程控制基础》课程中，传递函数，时域响应， 频率响应三章的内容。 学 时 数：2 设备仪器：XMN-2自动控制原理学习机，CAE-98型微机接口卡，计算机辅助实验系统2.0软件，万用表。 材料消耗：运算放大器，电阻，电容，插接线。 要 求：实验前认真预习实验指导书的实验内容，完成下述项目， 做实验时交于指导教师检查并与实验报告一起记入实验成绩。 B推导图2所示积分放大器的输出输入时域关系和传递函数。

C 推导图3所示加法和积分放大器的输出输入时域关系（两输入单输出） 和S <1>.写出op1,op2,op9,0p6对应的微分方程组（4个方程）。 <2>.画出系统方框图。 <3>.用方框图化简或方程组联立消元的方法求取实验电路所示系统的 传递函数，写出求解过程。 和ζ。 <4>.求取该系统的ω n 实验地点：教一楼327室 实验照片：实验装置及仪器

控制系统时间响应分析”实验报告

实验一、“控制系统时间响应分析”实验报告 一、实验类型 验证性实验 二、实验目的 1、 求系统在时间常数τ不同取值时的单位脉冲、单位阶跃响应和任意输入响应，熟悉系统时间响应的定义和图形曲线 2、 求系统的上升时间、峰值时间、最大超调量和调整时间等性能指标，熟悉系统瞬态性能指标的定义。 三、实验仪器与设备（或工具软件） 计算机，MATLAB 软件 四、实验内容、实验方法与步骤 已知系统传递函数 50 )1(05.050)(2+++=s s s G τ 1、求系统在时间常数τ不同取值时的单位脉冲、单位阶跃响应和任意输入响应。 应用impulse 函数，可以得到τ=0，τ=0.0125、τ=0.025时系统单位脉冲响应；应用step 函数，同样可以得到τ=0，τ=0.0125、τ=0.025时系统单位阶跃响应。 2、求系统的瞬态性能指标 五、实验结果 1、系统在时间常数τ不同取值时的单位脉冲、单位阶跃响应和任意输入响 t=[0:0.01:0.8];%仿真时间区段 nG=[50]; tao=0; dG=[0.05 1+50*tao 50]; G1=tf(nG ,dG); tao=0.0125; dG=[0.05 1+50*tao 50]; G2=tf(nG ,dG); tao=0.025; dG=[0.05 1+50*tao 50]; G3=tf(nG,dG);%三种τ值下，系统的传递函数模型 [y1,T]=impulse(G1,t);[y1a,T]=step(G1,t); [y2,T]=impulse(G2,t);[y2a,T]=step(G2,t); [y3,T]=impulse(G3,t);[y3a,T]=step(G3,t);%系统响应 subplot(131),plot(T,y1,'--',T,y2,'-.',T,y3,'-') legend('tao=0','tao=0.0125','tao=0.025') xlabel('t(sec)'),ylabel('x(t)');grid on; subplot(132),plot(T,y1a,'--',T,y2a,'-.',T,y3a,'-') legend('tao=0','tao=0.0125','tao=0.025') grid on;xlabel('t(sec)'),ylabel('x(t)');%产生图形 t=[0:0.01:1];u=sin(2*pi*t);% 仿真时间区段和输入 Tao=0.025;

二阶系统阶跃响应实验报告

实验一 二阶系统阶跃响应 一、实验目的 （1）研究二阶系统的两个重要参数：阻尼比ξ和无阻尼自振角频率ωn 对系统动 态性能的影响。 （2）学会根据模拟电路，确定系统传递函数。 二、实验内容 二阶系统模拟电路图如图2-1 所示。 系统特征方程为T 2s 2+KTs+1=0，其中T=RC ，K=R0/R1。根据二阶系统的标准 形式可知，ξ=K/2，通过调整K 可使ξ获得期望值。 三、预习要求 （1） 分别计算出T=0.5,ξ= 0.25，0.5，0.75 时，系统阶跃响应的超调量σP 和过渡过 程时间tS 。 ) 1( p 2 e ζζπσ--=， ζ T 3t s ≈

代入公式得： T=0.5,ξ= 0.25，σp=44.43% ，t s=6s； T=0.5,ξ= 0.5，σp=16.3% ，t s=3s； T=0.5,ξ= 0.75，σp=2.84% ，t s=2s； （2）分别计算出ξ= 0.25，T=0.2,0.5，1.0 时，系统阶跃响应的超调量σP 和过渡过程时间tS。 ξ= 0.25，T=0.2，σp=44.43% ，t s=2.4s； ξ= 0.25，T=0.5，σp=44.43% ，t s=6s； ξ= 0.25，T=1.0，σp=44.43% ，t s=12s； 四、实验步骤 （1）通过改变K，使ξ获得0，0.25，0.5，0.75，1.0 等值，在输入端加同样幅值的阶跃信号，观察过渡过程曲线，记下超调量σP 和过渡过程时间tS，将实验值和理论值进行比较。 （2）当ξ=0.25 时，令T=0.2 秒，0.5 秒，1.0 秒（T=RC,改变两个C），分别测出超调量σP 和过渡过程tS，比较三条阶跃响应曲线的异同。 五、实验数据记录与处理： 阶跃响应曲线图见后面附图。 原始数据记录： （1）T=0.5，通过改变R0的大小改变K值

控制系统的校正研究——频率响应法

论文题目:控制系统的校正研究——频率响应法 专业: 电子信息工程专业 姓名：签名：________ 指导老师：签名: ________ 摘要 摘要：近年来，自动控制系统在如今的工业和生活中，起着越来越重要的作用。所以，据用户要求的性能指标进行自动控制系统的串联校正设计有很重要的现实意义。对于给定的线性定常系统，通常通过加入串联超前、滞后或超前滞后综合校正装置，以达到提高系统的精度和稳定性的目的。该文分别给出基于频率特性法串联校正的具体设计方法，应用MATLAB对系统进行通用程序设计，并对实例进行仿真。仿真实例结果表明了此设计方法的有效性和实用性。 【关键词】：自动控制系统；频率响应法；MATLAB；伯德图 【论文类型】：理论研究型

Title:Correction of control system——Frequency response method Major: Electronic & Information Engineering Name:Signature: Supervisor：Signature: In recent years, automatic control systems play an increasingly important role in today's industrial and domestic.Therefore, the performance according to user requirements for the automatic control system series correcting design has a very important practical significance. For a given linear time-invariant systems, usually by joining the series ahead of lag or lead and lag correction device, in order to achieve the purpose to improve the accuracy and stability of the system. This paper gives specific design series based on the frequency characteristics correction, MATLAB system for generic programming, and simulation instance. The simulation results show the effectiveness and practicality of this design method. 【key word】：Automatic control system;Frequency response method，MATLAB;Bode diagram 【Type of Thesis】：Theory research

matlab中的矩阵的基本运算命令范文

1.1 矩阵的表示 1.2 矩阵运算 1.2.14 特殊运算 1．矩阵对角线元素的抽取 函数diag 格式X = diag(v,k) %以向量v的元素作为矩阵X的第k条对角线元素，当k=0时，v为X的主对角线；当k>0时，v为上方第k条对角线；当k<0时，v为下方第k条对角线。 X = diag(v) %以v为主对角线元素，其余元素为0构成X。 v = diag(X,k) %抽取X的第k条对角线元素构成向量v。k=0：抽取主对角线元素；k>0：抽取上方第k条对角线元素；k<0抽取下方第k条对角线元素。 v = diag(X) %抽取主对角线元素构成向量v。 2．上三角阵和下三角阵的抽取 函数tril %取下三角部分 格式L = tril(X) %抽取X的主对角线的下三角部分构成矩阵L L = tril(X,k) %抽取X的第k条对角线的下三角部分；k=0为主对角线；k>0为主对角线以上；k<0为主对角线以下。函数triu %取上三角部分 格式U = triu(X) %抽取X的主对角线的上三角部分构成矩阵U U = triu(X,k) %抽取X的第k条对角线的上三角部分；k=0为主对角线；k>0为主对角线以上；k<0为主对角线以下。3．矩阵的变维 矩阵的变维有两种方法，即用“：”和函数“reshape”，前者主要针对2个已知维数矩阵之间的变维操作；而后者是对于一个矩阵的操作。 （1）“：”变维 （2）Reshape函数变维 格式 B = reshape(A,m,n) %返回以矩阵A的元素构成的m×n矩阵B B = reshape(A,m,n,p,…) %将矩阵A变维为m×n×p×… B = reshape(A,[m n p…]) %同上 B = reshape(A,siz) %由siz决定变维的大小，元素个数与A中元素个数 相同。 （5）复制和平铺矩阵 函数repmat 格式 B = repmat(A,m,n) %将矩阵A复制m×n块，即B由m×n块A平铺而成。 B = repmat(A,[m n]) %与上面一致 B = repmat(A,[m n p…]) %B由m×n×p×…个A块平铺而成 repmat(A,m,n) %当A是一个数a时，该命令产生一个全由a组成的m×n矩阵。 1.3 矩阵分解 1.3.1 Cholesky分解 函数chol 格式R = chol(X) %如果X为n阶对称正定矩阵，则存在一个实的非奇异上三角阵R，满足R'*R = X；若X非正定，则产生错误信息。 [R,p] = chol(X) %不产生任何错误信息，若X为正定阵，则p=0，R与上相同；若X非正定，则p为正整数，R是有序的上三角阵。 1.3.2 LU分解

线性控制系统的频率响应分析

一．实验目的 1．了解和掌握对数幅频曲线和相频曲线（波德图）、幅相曲线（奈奎斯特图）的构造及绘制方法。 2．二阶开环系统中的相位裕度和幅值穿越频率的计算。 二．实验内容及要求 1．一阶惯性环节的频率特性曲线测试。 2．二阶开环系统的频率特性测试，研究表征系统稳定程度的相位裕度和 幅值穿越频率对系统的影响。 三、实验主要仪器设备和材料 1．labACT自控/计控原理实验机一台 2．数字存储示波器一台 四、实验方法、步骤及结果测试 1．一阶惯性环节的频率特性曲线 惯性环节的频率特性测试模拟电路见图4-1。 图4-1 惯性环节的频率特性测试模拟电路 实验步骤：注：‘S ST'不能用“短路套”短接！ （1）将数/模转换器（B2）输出OUT2作为被测系统的输入。 （2）按图4-1安置短路套及测孔联线。 （3）运行、观察、记录： ①运行LABACT程序，选择自动控制菜单下的线性控制系统的频率响应分析-实验项目，选择一阶系统，再选择开始实验，点击开始，实验机将自动产生0.5Hz~64Hz多个频率信号，测试被测系统的频率特性，等待将近十分钟，测试结束。 ②测试结束后，可点击界面下方的“频率特性”选择框中的任意一项进行切换，将显示被测系统的对数幅频、相频特性曲线（伯德图）和幅相曲线（奈 奎斯特图），同时在界面上方将显示点取的频率点的L、、Im、Re等相关数

据。如点击停止，将停止示波器运行，不能再测量数据。 ③分别改变惯性环节开环增益与时间常数，观察被测系统的开环对数幅频曲线、相频曲线及幅相曲线，在幅频曲线或相频曲线上点取相同的频率点，测量、记录数据于实验数据表中。 实验数据表1：改变惯性环节开环增益，（T=0.05，C=1u，R2=50K） 实验数据表2： 改变惯性环节时间常数， K=1（R1=50K、R2=50K） 2．二阶开环系统的频率特性曲线 二阶系统模拟电路图的构成如图4-2所示。

MATLAB下二阶系统的单位阶跃响应

二阶系统在不同参数下对单位阶跃信号的响应 一、二阶系统 所谓二阶系统就是其输入信号、输出信号的关系可用二阶微分方程来表征的系统。比如常见的RLC电路（图a）、单自由度振动系统等。 图a 图b 二阶系统传递函数的标准形式为 2 22 () 2 n n n H s s s ω ξωω = ++ 二、二阶系统的Bode图（nω=1） MATLAB程序为 >> clear >> num=[1]; >> den=[1 0.2 1]; >> bode(num,den); grid on hold on den=[1 0.4 1]; bode(num,den); >> den=[1 0.6 1]; >> bode(num,den); >> den=[1 0.8 1]; >> bode(num,den); >> den=[1 1.4 1]; >> bode(num,den); >> den=[1 2 1]; >> bode(num,den); >> legend('0.1','0.2','0.3','0.4','0.7','1.0')

运行结果为 三、二阶系统对单位阶跃信号的响应（ =1） n MATLAB程序为 >> clear >> num=[1]; >> den=[1 0 1]; >> t=0:0.01:25; >> step(num,den,t) >> grid on >> hold on >> den=[1 0.2 1]; >> step(num,den,t) >> den=[1 0.4 1]; >> step(num,den,t) >> den=[1 0.6 1]; >> step(num,den,t) >> den=[1 0.8 1]; >> step(num,den,t) >> den=[1 1.0 1]; >> step(num,den,t)

一阶系统的单位阶跃响应

图3－5所示系统。其输入－输出关系为 1 1 111)()(+= +=Ts s K s R s C （3-3） 式中K T 1 = ，因为方程(3-3)对应的微分方程的最高阶次是1，故称一阶系统。 实际上，这个系统是一个非周期环节，T 为系统的时间常数。 一、一阶系统的单位阶跃响应 因为单位阶跃函数的拉氏变换为s 1，将s s R 1)(=代入方程(3-3)，得 s Ts s C 1 11)(+= 将)(s C 展开成部分分式，有 11()1C s s s T =- + （3-4） 对方程(3-4)进行拉氏反变换，并用)(t h 表示阶跃响应)(t C ，有 t T e t h 1 1)(--= 0t ≥ (3-5) 由方程(3-5)可以看出，输出量)(t h 的初始值等于零，而最终将趋于1。常数项“1”是由s 1反变换得到的，显然，该分量随时间变化的规律和外作用相似(本例为相同)，由于它在稳态过程中仍起作用，故称为稳态分量 (稳态响应)。方程(3-5)中第二项由1 1/()s T +反变换得到， 它随时间变化的规律取决于传递函数1/(1)Ts +的极点，即系统特 征方程()10D s Ts =+=的根(1/)T -在复平 面中的位置，若根处在复平面的左半平面 如图3－6(a)所示，则随着时间 t 的增加， 它将逐渐衰减， 最后趋于零 (如图3－6(b) 所示)，称为瞬态响应。可见，阶跃响应曲线具有非振荡特性，故也称为非周期响应。 显然，这是一条指数响应曲线，其初始斜率等于1/T ，即 T e T dt dh t t T t 1 |1|01 0===-= (3-6) 这就是说，假如系统始终保持初始响应速度不变，那么当T t =时， 输出量就能达到稳态值。

实验一 MATLAB基本操作及运算

实验一 MATLAB 基本操作及运算 一、 实验目的 二、 实验的设备及条件 三、 实验内容 1、 建立以下标量： 1） a=3 2） ，（j 为虚数单位） 3） c=3/2πj e 2、 建立以下向量： 1） Vb= 2.71382882????????-???? 2） Vc=[4 3.8 … -3.8 -4 ] （向量中的数值从4到-4，步长为-0.2） 3、 建立以下矩阵： 1） 3 333Ma ????=?????? Ma 为一个7×7的矩阵，其元素全为3. 2） 11191212921020100Mb ??????=??????

Mb 为一个10×10的矩阵. 3） 114525173238Mc ????=?????? 4、 使用题1中的变量计算下列等式的x,y,z 的值： 1） ((15)/6)111a x e --=+ 2） 2x π= 3） 3ln([()()]sin(/3))x b c b c a π=+-R ，其中R 表示复数实部。 5、 求解函数值22/(2.25)ct y e -=，其中c 取值见题1，t 的取值范围为题2中行 向量Vc 。 6、 使用题1和题3中所产生的标量和矩阵计算等式 1()()T Mx a Mc Mc Mc -=?? 其中*为矩阵所对应行列式的值，参考det 。 7、 函数的使用和矩阵的访问。 1） 计算矩阵Mb 每一列的和，结果应为行向量形式。 2） 计算整个矩阵Mb 的平均值。 3） 用向量[1 1…1] 替换Mb 的最上一行的值 4） 将矩阵Mb 的第2~5行，第3到9列的元素所构成的矩阵赋值给矩阵SubMb 。 5） 删除矩阵Mb 的第一行； 6） 使用函数rand 产生一个1×10的向量r ，并将r 中值小于0.5的元素设置为0。 8、 已知CellA （1, 1）＝‘中国’，CellA （1，2）＝‘北京’，CellA （2，1）是一个3乘3的单位阵，CellA （2, 2）＝[1 2 3]，试用MATLAB 创建一个2×2的细胞数组CellA 。 9、 已知结构数组student 中信息包含有姓名，学号，性别，年龄和班级，试用MATLAB 创建相应的结构数组student 。该数组包含有从自己学号开始连续5个同学的信息（如果学号在你后面的同学不足5个则往前排序），创建完成后查看自己的信息。

第四章控制系统的频率特性

第四章控制系统的频率特性 本章要点 本章主要介绍自动控制系统频域性能分析方法。内容包括频率特性的基本概念，典型环节及控制系统Bode图的绘制，用频域法对控制系统性能的分析。 用时域分析法分析系统的性能比较直观，便于人们理解和接受。但它必须直接或间接地求解控制系统的微分方程，这对高阶系统来说是相当复杂的。特别是当需要分析某个参数改变对系统性能的影响时，需反复重新计算，而且还无法确切了解参数变化量对系统性能影响的程度。而频率特性不但可以用图解的方法分析系统的各种性能, 而且还能分析有关参数对系统性能的影响，工程上具有很大的实用意义。 第一节频率特性的基本概念 一、频率特性的定义 频率特性是控制系统的又一种数学模型，它是系统(或元件)对不同频率正弦输入信号的响应特性。对线性系统，若输入信号为正弦量，则其稳态输出信号也将是同频率的正弦量，但是输出信号的幅值和相位一般不同于输入量，如图4-1。 若设输入量为r(t)=A「sin( 3 t+ u r) 其输出量为c(t)=A c sin@ t+ u c) 若保持输入信号的幅值A r不变，改变输入信号的角频率3,则输出信号的角频率 也变化，并且输出信号的幅值和相位也随之变化。 图4-1控制系统的频率响应

我们定义系统(或环节)输出量与输入量幅值之比为幅值频率 特性，简称幅频 M( 3 )表示。输出量与输入量的相位差为相位频率特 3变化，常用U (3 )表示。其数学定义为 M "A U ( 3 )= U c - U 幅频特性和相频特性统称为频率特性，用 G(j 3 )表示。由此，幅频特性 M( 3 )又可 表示为|G(j ;i )，相频特性u (3 )又可表示为Z G(j ■)，三者可表示成下面的形式: G(j a )=|G(j m )|Z G(j s ) M (co ) = G(jco) 「()二/G( j ?) 二、频率特性与传递函数的关系 频率特性和传递函数之间存在密切关系：若系统(或元件)的传递函数为 G(s), 则其频率特性为 G(j 3 )。这就是说，只要将传递函数中的复变量 s 用纯虚数j 3代替， 就可以 得到频率特性。即 G(s) > G(j ■) 三、频率特性的表示方法 1 .数学式表示法 频率特性是一个复数，所以它和其他复数一 | 样，可以表示为极坐标式、直角坐标和指数坐标 三种形式。见图 4-2所示。 G(j ?)二 G(j J- G(j ) 二U (■) jVC ) -M ( )e j () 显然， M =|G( j ⑷)| 2 (co )+V 2?) w G(j "arcta 说 例4-1写出惯性环节的幅频特性、相频特性和频率特性。 特性，它随角频率 3变化，常用 性，简称相频特性，它也随角频率 其中 图4-2频率特性的表示方法

2. 实验二 二阶系统阶跃响应

实验二二阶系统阶跃响应 一、实验目的 1. 研究二阶系统的特征参数，阻尼比ζ和无阻尼自然频率ωn对系统动态性能的影响，定量分析ζ和ωn与最大超调量σp和调节时间ts之间的关系。 2. 进一步学习实验系统的使用。 3. 学会根据系统的阶跃响应曲线确定传递函数。 4. 学习用MATLAB仿真软件对实验内容中的电路进行仿真。 二、实验原理 典型二阶闭环系统的单位阶跃响应分为四种情况： 1）欠阻尼二阶系统 如图1所示，由稳态和瞬态两部分组成：稳态部分等于1，瞬态部分是振荡衰减的过程，振荡角频率为阻尼振荡角频率，其值由阻尼比ζ和自然振荡角频率ωn决定。 （1）性能指标： : 单位阶跃响应C(t)进人±5%(有时也取±2%)误差带，并且不再超出该误差带的调节时间t S 最小时间。 超调量σ% ；单位阶跃响应中最大超出量与稳态值之比。 单位阶跃响应C(t)超过稳态值达到第一个峰值所需要的时间。 峰值时间t P ： 结构参数ξ：直接影响单位阶跃响应性能。 （2）平稳性：阻尼比ξ越小，平稳性越差 长，ξ过大时，系统响应迟钝，（3）快速性：ξ过小时因振荡强烈，衰减缓慢，调节时间t S 也长，快速性差。ξ＝0.7调节时间最短，快速性最好。ξ＝0.7时超调量σ%<5％,调节时间t S 平稳性也好，故称ξ＝0.7为最佳阻尼比。 2）临界阻尼二阶系统（即ξ＝1） 系统有两个相同的负实根，临界阻尼二阶系统单位阶跃响应是无超调的，无振荡单调上升的，不存在稳态误差。

3）无阻尼二阶系统（ξ＝0时）此时系统有两个纯虚根。 4）过阻尼二阶系统（ξ>1）时 此时系统有两个不相等的负实根，过阻尼二阶系统的单位阶跃响应无振荡无超调无稳态误差，上升速度由小加大有一拐点。 三、实验内容 1. 搭建模拟电路 典型二阶系统的闭环传递函数为： 其中，ζ 和ωn对系统的动态品质有决定的影响。 搭建典型二阶系统的模拟电路，并测量其阶跃响应： 二阶系统模拟电路图其结构图为： 系统闭环传递函数为： 式中, T=RC，K=R2/R1。 比较上面二式，可得：ωn=1/T=1/RC ζ=K/2=R2/2R1。 2 2 2 2 ) ( ) ( ) ( n n n w s w s w s R s C S + + = = ξ φ

二阶系统阶跃响应实验报告

实验一二阶系统阶跃响应 一、实验目的 (1)研究二阶系统的两个重要参数：阻尼比E和无阻尼自振角频率3 态性能的影 响。 (2)学会根据模拟电路，确定系统传递函数。 二、实验内容 二阶系统模拟电路图如图2-1所示 a 2-i二阶系疣按拟电帘图 系统特征方程为TV+KTS+仁0其中T=RC K=R0/R1根据二阶系统的标准 形式可知，E =K/2，通过调整K可使E获得期望值 三、预习要求 (1) 分别计算出T=0.5,E = 0.25, 0.5, 0.75时，系统阶跃响应的超调量c P和过渡过程时 间ts。 代入公式得: T=0.5, E : =0.25, c P=44.43%，t s=6s； T=0.5, E : =0.5 ， d P=16.3% ，t s=3s； T=0.5, E : =0.75, c p=2.84% ，t s=2s; (2) 分别计算出E = 0.25，T-0.2,0.5，1.0时，系统阶跃响应的超调量c P和过渡 过程时间ts。 E = =0.25,T-0.2, c p-44.43% ,t s- 2.4s； E = =0.25,T-0.5, c P-44.43% ,t s-6s； E = =0.25,T-1.0, c P-44.43% ,t s- 12s； 四、 (1) 实验步骤 通过改变K,使E获得0, 0.25, 0.5, 0.75, 1.0等值，在输入端加同样幅值的阶跃 信号，观察过渡过程曲线，记下超调量b P和过渡过程时间ts,将实验值和理论值 进行比较。 n对系统动 ) 2 t s 3T

(2)当E =0.25时，令T=0.2秒，0.5秒，1.0秒(T=RC改变两个C),分别测出超调量b P和过渡过 程tS,比较三条阶跃响应曲线的异同。 五、实验数据记录与处理： 阶跃响应曲线图见后面附图。 原始数据记录： (1) T=0.5，通过改变R0的大小改变K值 理论值与实际值比较: 对误差比较大，比如T=0.5,E =0.75时，超调量的相对误差为30%左右。造成误差的原因主要有以下几个方面： (1)由于R0是认为调整的阻值，存在测量和调整误差，且不能精确地保证E的大小等于 要求的数值； (2)在预习计算中我们使用了简化的公式，例如过渡时间大约为3~4T/ E，这并不是一个 精确的数值，且为了计算方便取3T/E作统一计算； (3)实际采样点的个数也可能造成一定误差，如果采样点过少，误差相对会大。 六、实验总结 通过本次实验，我们从图形上直观的二阶系统的两个参数对系统动态性能的影响，巩固了理论知识。其次我们了解了一个简单的系统是如何用电路方式实现的，如何根据一个

单位阶跃响应

20 7619 )(2++= Φs s s 单位阶跃响应 num=[19]; den=[6,7,20]; %描述闭环系统传递函数的分子、分母多项式 sys=tf(num,den); %高阶系统建模 step(sys); %绘制高阶系统的单位阶跃响应曲线 grid; %添加栅格 title('单位阶跃响应'); %标注标题 xlabel('t'); ylabel('c(t)'); %标注横、纵坐 标轴 12345678910 00.20.40.60.8 1 1.2 1.4 单位阶跃响应 t (sec) c (t ) 单位斜坡响应 sys=tf([19],[6,7,20]); %系统建模 t=0:0.01:10; %响应时间 u=t; %单位斜坡输入 lsim(sys,u,t) %单位斜坡响应 grid xlabel(‘t ’); ylabel(‘c(t)’) %标注横、纵坐标轴 title(‘单位斜坡响应’); %标注标题

012345678910 123456 789 10单位斜坡响应 t (sec) c (t ) 单位加速度响应 num=[19]; den=[6,7,20]; sys=tf(num,den); %系统建模 t=0:0.01:10; %响应时间序列 u=0.5*t.^2; %单位加速度输入 lsim(sys,u,t) %绘制单位加速度响应曲线 grid xlabel('t'); ylabel('c(t)'); title('单位加速度响应');

012345678910 51015202530 354045 50单位加速度响应 t (sec) c (t ) 1 4014121 38)(2 3++++= Φs s s s s 单位阶跃响应 num=[38,1]; den=[13,14,40,1]; %描述闭环系统传递函数的分子、分母多项式 sys=tf(num,den); %高阶系统建模 step(sys); %绘制高阶系统的单位阶跃响应曲线 grid; %添加栅格 title('单位阶跃响应'); %标注标题 xlabel('t'); ylabel('c(t)'); %标注横、纵坐标 轴

(整理)matlab16常用计算方法.

常用计算方法 1．超越方程的求解 一超越方程为 x (2ln x – 3) -100 = 0 求超越方程的解。 [算法]方法一：用迭代算法。将方程改为 01002ln()3 x x =- 其中x 0是一个初始值，由此计算终值x 。取最大误差为e = 10-4，当| x - x 0| > e 时，就用x 的值换成x 0的值，重新进行计算；否则| x - x 0| < e 为止。 [程序]P1_1abs.m 如下。 %超越方程的迭代算法 clear %清除变量 x0=30; %初始值 xx=[]; %空向量 while 1 %无限循环 x=100/(2*log(x0)-3); %迭代运算 xx=[xx,x]; %连接结果 if length(xx)>1000,break ,end %如果项数太多则退出循环（暗示发散） if abs(x0-x)<1e-4,break ,end %当精度足够高时退出循环 x0=x; %替换初值 end %结束循环 figure %创建图形窗口 plot(xx,'.-','LineWidth',2,'MarkerSize',12)%画迭代线'.-'表示每个点用.来表示，再用线连接 grid on %加网格 fs=16; %字体大小 title('超越方程的迭代折线','fontsize',fs)%标题 xlabel('\itn','fontsize',fs) %x 标签 ylabel('\itx','fontsize',fs) %y 标签 text(length(xx),xx(end),num2str(xx(end)),'fontsize',fs)%显示结果 [图示]用下标作为自变量画迭代的折线。如P0_20_1图所示，当最大误差为10-4时，需要迭代19次才能达到精度，超越方程的解为27.539。 [算法]方法二：用求零函数和求解函数。将方程改为函数 100()2ln()3f x x x =-- MATLAB 求零函数为fzero ，fzero 函数的格式之一是 x = fzero(f,x0) 其中，f 表示求解的函数文件，x0是估计值。fzero 函数的格式之二是 x = fzero(f,[x1,x2])

matlab有限域上的运算

1 有限域基础知识 1.1 有限域（Galois域）的构造 令p为一个素数. 则对任意的一个正整数n，存在一个特征为p，元素个数为p n的有限域GF(p n). 注：任意一个有限域，其元素的个数一定为p n，其中p为一个素数（有限域的特征），n为一个正整数. 例1（有限域GF(p)）令p为一个素数，集合 GF(p)=Z p={0,1,2,…,p?1}. 在GF(p)上定义加法⊕和乘法⊙分别为模p加法和模p乘法，即任意的a,b∈GF(p)， a⊕b=(a+b)mod p, a⊙b=(a?b)mod p 则为一个有p个元素的有限域，其中零元素为0，单位元为1. 令a为GF(p)中的一个非零元素. 由于gcd(a,p)=1，因此，存在整数b,c，使得ab+pc=1. 由此得到a的逆元为a?1=b mod p. 域GF(p)称为一个素域(prime field).

例注1：给定a和p，例1中的等式ab+pc=1可以通过扩展的欧几里得除法得到，从而求得GF(p)中任意非零元素的逆元. 例2（有限域GF(p n)）从GF(p)出发，对任意正整数n，n≥2，我们可以构造元素元素个数为p n的有限域GF(p n)如下： 令g(x)为一个GF(p)上次数为n的不可约多项式，集合 GF(p n)=GF(p)[x]/?g(x)?={a0+a1x+a2x2+?+a n?1x n?1 | a i∈ GF(p),0≤i≤n?1} 在GF(p n)上定义加法⊕和乘法⊙分别为模g(x)加法和模g(x)乘法，即任意的a(x),b(x)∈GF(p n)， a(x)⊕b(x)=a(x)+b(x), a(x)⊙b(x)=(a(x)?b(x))mod g(x) 则为一个有p n个元素，特征为p的有限域，其中零元素为GF(p)中的0，单位元为GF(p)中的1. 令a(x)为GF(p n)中的一个非零元素. 由于gcd(a(x),g(x))=1，因此，存在GF(p)上的多项式b(x),c(x)，使得a(x)b(x)+g(x)c(x)=1. 由此得到 a(x)的逆元为a?1(x)=b(x)mod g(x). 域GF(p n)称为GF(p)的（n次）扩域(extension field)，而GF(p)称为GF(p n)的子域(subfield).

第三章 系统频率特性

第三章 系统频率特性 系统的时域分析是分析系统的直接方法，比较直观，但离开计算机仿真，分析高阶系统是困难的。系统频域分析是工程广为应用的系统分析和综合的间接方法。频率分析不仅可以了解系统频率特性，如截止频率、谐振频率等，而且可以间接了解系统时域特性，如快速性，稳定性等，为分析和设计系统提供更简便更可靠的方法。 本章首先阐明频率响应的特点，给出计算频率响应的方法，接着介绍Nyquist 图和Bode 图的绘制方法、系统的稳定裕度及系统时域性能指标计算。 3.1 频率响应和频率特性 3.1.1 一般概念 频率响应是指系统对正弦输入的稳态响应。考虑传递函数为G(s)的线性系统，若输入正弦信号 t X t x i i ωsin )(= (3.1-1) 根据微分方程解的理论，系统的稳态输出仍然为与输入信号同频率的正弦信号，只是其幅值和相位发生了变化。输出幅值正比于输入的幅值i X ，而且是输入正弦频率ω的函数。输出的相位与i X 无关，只与输入信号产生一个相位差?，且也是输入信号频率ω的函数。即线性系统的稳态输出为 )](sin[)()(00ω?ωω+=t X t x (3.1-2)

由此可知，输出信号与输入信号的幅值比是ω的函数，称为系统的幅频特性，记为)(ωA 。输出信号与输入信号相位差也是ω的函数，称为系统的相频特性，记为)(ω?。 幅频特性： )()()(0ωωωi X X A = (3.1-3) 相频特性： )()()(0ω?ω?ω?i -= (3.1-4) 频率特性是指系统在正弦信号作用下，稳态输出与输入之比对频率的关系特性，可表示为： )()()(0ωωωj X j X j G i = (3.1-5) 频率特性)(ωj G 是传递函数)(s G 的一种特殊形式。任何线性连续时间系统的频率特性都可由系统传递函数中的s 以ωj 代替而求得。 )(ωj G 有三种表示方法： )()()(ω?ωωj e A j G = (3.1-6) )()()(ωωωjV U j G += (3.1-7) )(sin )()cos()()(ω?ωωωωjA A j G += (3.1-8) 式中，实频特性： )(cos )()(ω?ωωA U = 虚频特性：

典型环节的单位阶跃响应

实验二 典型环节的单位阶跃响应 一、实验目的 1、根据对象的单位阶跃响应特性，掌握和深刻理解几种典型环节的特性以及它们特性参数的含义。 2、研究对象传递函数的零极点对系统动态特性的影响。 3、学习Matlab 的基本用法 ――求取阶跃响应、脉冲响应(step, impulse) ――基本做图方法(hold, plot) 二、实验内容 1、比例环节 求取K s G )(在不同比例系数K 下的单位阶跃响应，说明比例系数对系统动态过程的影响。 0.10.20.30.40.50.60.70.80.91 G(s)=K,在不同比例系数K 下的单位阶跃响应 Time (sec) A m p l i t u d e 由上图可以看出： 因为G （s ）=K ，所以被控对象是一个单纯的比例系统。随着K 的增加，系统的终值是输入信号的K 倍。 2、一阶惯性环节

（1） 求取1 )(+= Ts K s G 的单位阶跃响应，其中放大倍数K =２，时间常数T =２。 1)(+= Ts K s G 的单位阶跃响应如下图： 024681012 0.20.40.60.811.2 1.41.61.8 2G(s)=2/(2s+1)的单位阶跃响应 Time (sec) A m p l i t u d e

（2） 求取1 22 )(+= s s G 的单位脉冲响应，可否用step 命令求取它的脉冲响应？ 122 )(+= s s G 的单位脉冲响应如下图： 024681012 0.10.20.30.40.50.6 0.70.80.9 1G(s)=2/(2s+1)的单位m 脉冲响应 Time (sec) A m p l i t u d e 把传递函数乘以s 再求其单位阶跃响应，就可获得乘s 前的传递函数的脉冲响应。如下图： 024681012 0.10.20.30.40.50.6 0.70.80.9 1G(s)=2*s/(2s+1)的单位m 阶跃响应 Tim e (sec) A m p l i t u d e

(整理)二阶系统的阶跃响应.

实验一 一、二阶系统的阶跃响应 实验报告 ＿＿＿系＿＿专业＿＿＿班级 学号＿＿＿姓名＿＿＿成绩＿＿＿指导教师＿＿一、实验目的 1、学习实验系统的使用方法。 2、学习构成一阶系统（惯性环节）、二阶系统的模拟电路，分别推导其传递函数。了解电路参数对环节特性的影响。 3、研究一阶系统的时间常数T 对系统动态性能的影响。 4、研究二阶系统的特征参数，阻尼比ξ和无阻尼自然频率n ω对系统动态性能的影响。 二、实验仪器 1、EL-AT-II 型自动控制系统实验箱一台 2、计算机一台 三、实验内容 （一） 构成下述一阶系统（惯性环节）的模拟电路，并测量其阶跃响应。 惯性环节的模拟电路及其传递函数如图1-1。 （二）构成下述二阶系统的模拟电路，并测量其阶跃响应。 典型二阶系统的闭环传递函数为 ()2222n n n s s s ωζωω?++= （1） 其中ζ和n ω对系统的动态品质有决定的影响。 图1-1 一阶系统模拟电路图 R1 R2

构成图1-2典型二阶系统的模拟电路，并测量其阶跃响应： 电路的结构图如图 1-3 系统闭环传递函数为 ()()()()2 2 2/1//11/2T S T K s T s U S U s ++==? 式中 T=RC ，K=R2/R1。 比较（1）、（2）二式，可得 n ω=1/T=1/RC ξ=K/2=R2/2R1 （3） 由（3）式可知，改变比值R2/R1，可以改变二阶系统的阻尼比。改变RC 值可以改变无阻尼自然频率n ω。 今取R1=200K ，R2=0K Ω，50K Ω，100K Ω和200K Ω，可得实验所需的阻尼比。图1-2 二阶系统模拟电路图 图1-3 二阶系统结构图 R2

闭环零点对二阶系统单位阶跃响应的影响

闭环零点对二阶系统单位阶跃响应的影响作者： 单位： 邮编： 摘要 在工程上电路中出现两个储能元件时便构成了二阶系统。由于欠阻尼二阶系统最具有实际意义，并且二阶系统往往需要满足工程最佳参数的要求，然而仅仅通过改变开环放大系数从而满足工程要求则可能会出现系统稳态误差增大的现象，设置具有闭环零点的二阶系统既可以达到满足工程所需的阻尼比，又可保证系统稳态精度。 在全面的分析了二阶系统之后，得出二阶系统的动态变化，由此引入带有零点的二阶系统，并分析了在欠阻尼状态下二阶系统的单位阶跃响应，并分析了其上升时间、峰值时间、调节时间、最大超调量，与没有零点的二阶系统进行了动态特性的对比。在此基础上分析了零点位置变化对二阶系统的影响。得到了重要结论。 关键字：二阶系统上升时间峰值时间调节时间最大超调量

0 引言 在已经知道了二阶系统的动态特性的基础之上，进一步研究具有闭环零点的二阶系统。并通过对比二阶系统和具有闭环零点的二阶系统，得出一定的结论。讨论当零点移动时对动态特性的影响。对满足工程所需的阻尼比，保证系统稳态精度具有重要作用。 1 二阶系统 用二阶微分方程描述的系统成为二阶系统。 等效开环传递函数方框图： 其闭环传递函数方框图： 其中n ω无阻尼自然振荡角频率，ξ为阻尼比。 W B （s ）=2n 22n 2n s s ＋ωξω＋ω （1-1） 二阶系统的特征方程为： 2n 22n s s ＋ωξω＋=0 两根为S 1,2=12n n －ξω－ξω 二阶系统极点分布图：

1、当ξ>1时，（过阻尼） 2、当0<ξ<1时，（欠阻尼） 3、当ξ=1时，（临界阻尼） 4、当ξ=0时，（无阻尼） 5、当ξ<0时，（发散振荡） 在不同的阻尼比时，二阶系统的动态响应有很大的差别，因此阻尼比ξ是二阶系统的重要参数，当ξ<0时系统不可以正常工作，而在ξ>1时，系统动态响应进行得太慢，所以对二阶系统来说欠阻尼是最有实际意义的。