哈工大现代控制理论实验报告一

Harbin Institute of Technology

现代控制理论基础上机实验一

亚微米超精密车床振动控制系统

的状态空间法设计

院系：航天学院控制科学与工程系

专业：探测制导与控制技术

姓名：

班号：

学号：

指导教师：史小平

哈尔滨工业大学

2015年5月26日

目录

一、工程背景介绍及物理描述 (3)

1.1 工程背景介绍 (3)

1.2实验目的 (3)

1.3工程背景的物理描述 (3)

二．闭环系统的性能指标 (5)

三．实际给定参数 (6)

四．车床振动系统的开环状态空间模型 (6)

五．状态反馈控制律的设计过程 (7)

六. 闭环系统数字仿真的MATLAB编程 (8)

6.1源程序 (8)

6.2 运行截图 (9)

七. 实验结论及心得 (10)

7.1实验结论 (10)

7.2 心得体会 (11)

一、工程背景介绍及物理描述

1.1 工程背景介绍

超精密机床是实现超精密加工的关键设备，而环境振动又是影响超精密加工精度的重要因素。为了充分隔离基础振动对超精密机床的影响，目前国内外均采用空气弹簧作为隔振元件，并取得了一定的效果，但是这属于被动隔振，这类隔振系统的固有频率一般在2Hz左右。

这种被动隔振方法难以满足超精密加工对隔振系统的要求。为了解决这个问题，有必要研究被动隔振和主动隔振控制相结合的混合控制技术。其中，主动隔振控制系统采用状态空间法设计，这就是本次上机实验的工程背景。

1.2实验目的

通过本次上机实验，熟练掌握：

1. 控制系统机理建模；

2. 时域性能指标与极点配置的关系；

3. 状态反馈控制律设计；

4. MATLAB语言的应用。

四个知识点。

1.3工程背景的物理描述

上图表示了亚微米超精密车床隔振控制系统的结构原理，其中被动隔振元件为空气弹簧，主动隔振元件为采用状态反馈控制策略的电磁作动器。

上图表示一个单自由度振动系统，空气弹簧具有一般弹性支承的低通滤波特性，其主要作用是隔离较高频率的基础振动，并支承机床系统；主动隔振系统具有高通滤波特性，其主要作用是有效地隔离较低频率的基础振动。主、被动隔振系统相结合可有效地隔离整个频率范围内的振动。

床身质量的运动方程为：

p a 0

ms F F ++=

p

F ——空气弹簧所产生的被动控制力； a

F ——作动器所产生的主动控制力。

假设空气弹簧内为绝热过程，则被动控制力可以表示为：

p 0r r r e e

{1[/()]}n F cy k y p V V A y A =++-+

r

V ——标准压力下的空气弹簧体积；

y s s =-——相对位移（被控制量）；

r

p ——空气弹簧的参考压力； r

A ——参考压力下单一弹簧的面积；

e r

4A A =——参考压力下空气弹簧的总面积；

n ——绝热系数。

电磁作动器的主动控制力与电枢电流、磁场的磁通量密度及永久磁铁和电磁铁之间的间隙面积有关，这一关系具有强非线性。

由于系统工作在微振动状况，且在低于作动器截止频率的低频范围内，因此主动控制力可近似线性化地表示为：

a e a F k I =

e k ——力-电流转换系数； a

I ——电枢电流。

其中，电枢电流

a

I 满足微分方程：

a a a (,)()

LI RI E I y u t ++=

L ——控制回路电枢电感系数； R ——控制回路电枢电阻； E ——控制回路反电动势；

u ——控制电压。

二．闭环系统的性能指标

闭环系统单位阶跃响应的超调量不大于5%，过渡过程时间不大于0.5秒（0.02?=）。

三．实际给定参数

=1200N/m m=120kg =980N/A c=0.2 R=300 L=0.95H

四．车床振动系统的开环状态空间模型

首先假定0

s 为常数，将式0

y s s =-两边求关于时间的二阶导数可得：

()p a 1

y s F F m ==-+

()0r r r e e e a 1

{1[/()]}n cy k y p V V A y A k I m =-

++-++ （5） 记为

()0e a 1

y cy k y k I m ω=-

+++ （6）

其中

r r r e e

{1[/()]}n p V V A y A ω=-+。

式（6）两边求导数得：

()0e a 1

y cy k y k I m ω=-

+++ （7）

由式（6）可得：

0a e

my cy k y I k ω

+++=-

（8）

由式（7）可得：

0a e

my cy k y I k ω

+++=-

（9）

将式（8）和（9）代入式（4）可得：

00a e e

(,)()

my cy k y my cy k y L

R E I y u t k k ωω

++++++--+=

()()00e a e (,)()

Lmy Lc Rm y Lk Rc y Rk y L R k E I y k u t ωω+++++++-=-将非线性项

e a (,)

L R k E I y ωω+-视为干扰信号，略去不计，可得线性化模型：

()()00e ()

Lmy Lc Rm y Lk Rc y Rk y k u t +++++=- （10）

五．状态反馈控制律的设计过程

（1）根据带入实际参数后的线性化模型列写系统状态空间表达式： 令 则

u y = x

（2）设计一状态反馈控制律，使对y 的控制满足性能指标

由性能指标要求的闭环系统单位阶跃响应的超调量不大于5%，过渡过程时间不大于0.5秒（0.02?=），取阻尼比为0.8， ，进一步的，将闭环极点配置如下：

设状态反馈控制律

，则

解得

故加入控制律后的状态方程为

u y = x

六. 闭环系统数字仿真的MATLAB 编程

给定初始条件：

51(0)610m

x -=?，

52(0)210m/s

x -=?，

52

3(0)0.810m/s x -=-?

下面进行仿真分析。

6.1源程序

第一个文件simu046.m function dx=simu046(t,x)

A=[0,1,0; 0,0,1; -5850,-940,-63]; dx=A*x;

第二个文件do_simu046.m

[t,x]=ode45('simu046',[0,1],[6*10^-5,2*10^-5,-0.8*10^-5]); subplot(3,1,1); plot(t,x(:,1),'r-'); legend('x_1'); grid;

subplot(3,1,2); plot(t,x(:,2),'b-'); legend('x_2'); grid;

subplot(3,1,3);

plot(t,x(:,3),'k-');

legend('x_3');

grid;

阶跃响应曲线绘制：

A=[0,1,0; 0,0,1; -5850,-940,-63];

B=[0;0;-8.6];

C=[1 0 0];

G1=ss(A,B,C,0);

G2=-tf(G1)

t=0:0.02:5;

p=step(G2,t);

plot(t,p);grid;

xlabel('t-sec');

6.2 运行截图

阶跃响应曲线：

七. 实验结论及心得

7.1实验结论

闭环系统单位阶跃响应的超调量不大于5%，过渡过程时间不大于0.5秒

?=）。

（0.02

下面进行验证，代码如下：

A=[0,1,0; 0,0,1; -5850,-940,-63];

B=[0;0;-8.6];

C=[1 0 0];

G1=ss(A,B,C,0);

G2=-tf(G1)

[y, t, x]=step(G2);

mp=max(y);

yss=y(cs)

ct=(mp - yss)/yss

运行结果如下：

ct =

0.0159

即超调量为1.59%，满足指标要求。再利用plot界面的data cursor，得到调整时间约为0.32s，同样满足指标要求，状态反馈设计起到了预期的效果。

7.2 心得体会

经过本次上机实验，我认识到现代控制理论这门课并不只是理论分析，它同样具有工程实践的理论指导意义，经过状态反馈控制率的设计，实现了对现有系统性能指标的改良。通过实验进一步理解掌握了状态空间法的设计步骤，充分将所学的知识应用到了实际的系统设计中，同时也第一次用MATLAB 建立系统的状态空间模型，以及性能指标分析，收获很多。

现代控制理论实验报告

实验报告 ( 2016-2017年度第二学期) 名称：《现代控制理论基础》 题目：状态空间模型分析 院系：控制科学与工程学院 班级： ___ 学号： __ 学生姓名： ______ 指导教师： _______ 成绩： 日期： 2017年 4月 15日

线控实验报告 一、实验目的： l.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验内容 1 第一题：已知某系统的传递函数为G (s) S23S2 求解下列问题： （1）用 matlab 表示系统传递函数 num=[1]; den=[1 3 2]; sys=tf(num,den); sys1=zpk([],[-1 -2],1); 结果： sys = 1 ------------- s^2 + 3 s + 2 sys1 = 1 ----------- (s+1) (s+2) （2）求该系统状态空间表达式： [A1,B1,C1,D1]=tf2ss(num,den); A = -3-2 10 B = 1 C = 0 1

第二题：已知某系统的状态空间表达式为： 321 A ,B,C 01：10 求解下列问题： （1）求该系统的传递函数矩阵： （2）该系统的能观性和能空性： （3）求该系统的对角标准型： （4）求该系统能控标准型： （5）求该系统能观标准型： （6）求该系统的单位阶跃状态响应以及零输入响应：解题过程： 程序： A=[-3 -2;1 0];B=[1 0]';C=[0 1];D=0; [num,den]=ss2tf(A,B,C,D); co=ctrb(A,B); t1=rank(co); ob=obsv(A,C); t2=rank(ob); [At,Bt,Ct,Dt,T]=canon(A,B,C,D, 'modal' ); [Ac,Bc,Cc,Dc,Tc]=canon(A,B,C,D, 'companion' ); Ao=Ac'; Bo=Cc'; Co=Bc'; 结果： （1） num = 0 01 den = 1 32 （2）能控判别矩阵为： co = 1-3 0 1 能控判别矩阵的秩为： t1 = 2 故系统能控。 （3）能观判别矩阵为： ob = 0 1

哈工大单片机实验报告(上传)

H a r b i n I n s t i t u t e o f T e c h n o l o g y 单片机原理与应用 实验报告 学生姓名： 学号： 班级： 专业： 任课教师： 所在单位： 2013年5月

软件实验 在软件实验部分，通过实验程序的调试，使学生熟悉MCS-51的指令系统，了解程序设计过程，掌握汇编语言设计方法以及如何使用实验系统提供的调试手段来排除程序错误。 实验一清零程序 一、实验目的 掌握汇编语言设计和调试方法，熟悉键盘操作。 二、实验内容 把2000~20FFh的内容清零。 三、程序框图 四、实验过程 实验中利用MOVX语句，将外部存储器指定内容清零。利用数据指针DPTR完成数据传送工作。程序采用用循环结构完成，R0移动单元的个数，可用CJNE比较语句判断循环是否结束。 五、实验结果及分析 清零前清零后

【问题回答】清零前2000H~20FFH中为内存里的随机数，清零后全变为0。 六、实验源程序 AJMP MAIN ORG 0640H MAIN: MOV R0, #00H MOV DPL, #00H MOV DPH, #20H LOOP: MOV A, #00H MOVX @DPTR, A INC DPTR INC R0 CJNE R0, #0FFH, LOOP MOVX @DPTR, A END 实验二拆字程序 一、实验目的 掌握汇编语言设计和调试方法。 二、实验内容 把2000h的内容拆开，高位送2001h低位，低位送2002h低位，2001h、2002h高位清零，一般本程序用于把数据送显示缓冲区时用。 三、程序框图 四、实验过程 将寄存器中内容送入2000H，分别将高低四位移到低位，将高四位置零然后移入2001H 和2002H中。利用MOVX语句、DPTR指针可实现数据的传送，利用高低四位交换语句SWAP和与语句ANL可进行对高低位的清零。

现代控制理论课程报告

- 现代控制理论课程总结 学习心得 从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程，在刚拿到课本的时候，没上张老师的课之前，咋一看，会认为开课的内容会是上学期学的控制理论基础的累赘或者简单的重复，更甚至我还以为是线性代数的复现呢！根本没有和现代控制论联系到一起。但后面随着老师讲课的风格的深入浅出，循循善诱，发现和自己想象的恰恰相反，张老师以她特有的讲课风格，精心准备的 ppt 课件，向我们展示了现代控制理论发展过程，以及该掌握内容的方方面面，个人觉得，我们不仅掌握了现代控制理论的理论知识，更重要的是学会了掌握这门知识的严谨的逻辑思维和科学的学习方法，对以后学习其他知识及在工作上的需要大有裨益，总之学习了这门课让我受益匪浅。 由于我们学习这门课的课时不是很多，并结合我们学生学习的需求及所要掌握的课程深入程度，张老师根据我们教学安排需要，我们这学期学习的内容主要有： 1.绪论；2.控制系统的状态表达式；3.控制系统状态表达式的解；4.线性系统的能空性和能观性；5.线性定常系统的综合。而状态变量和状态空间表达式、状态转移矩阵、系统的能控性与能观性以及线性定常系统的综合是本门课程的主要学习内容。当然学习的内容还包括老师根据多年教学经验及对该学科的研究的一些深入见解。 在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的必修课。 经典控制理论的特点 经典控制理论以拉氏变换为数学工具，以单输入－单输出的线性定常系统为主要的研究对象。将描述系统的微分方程或差分方程变换到复数域中，得到系统的传递函数，并以此作为基础在频率域中对系统进行分析和设计，确定控制器的结构和参数。通常是采用反馈控制，构成所谓闭环控制系统。经典控制理论具有明显的局限性，突出的是难以有效地应用于时变系统、多变量系统，也难以揭示系统更为深刻的特性。当把这种理论推广到更为复杂的系统时，经典控制理论就显得无能为力了，这是因为它的以下几个特点所决定。 [ 1．经典控制理论只限于研究线性定常系统，即使对最简单的非线性系统也是无法处理的；这就从本质上忽略了系统结构的内在特性，也不能处理输入和输出皆大于1的系统。实际上，大多数工程对象都是多输入－多输出系统，尽管人们做了很多尝试，但是，用经典控制理论设计这类系统都没有得到满意的结果；2．经典控制理论采用试探法设计系统。即根据经验选用合适的、简单的、工程上易于实现的控制器，然后对系统进行分析，直至找到满意的结果为止。虽然这

现代控制理论实验

华北电力大学 实验报告| | 实验名称状态空间模型分析 课程名称现代控制理论 | | 专业班级：自动化1201 学生姓名：马铭远 学号：2 成绩： 指导教师：刘鑫屏实验日期：4月25日

状态空间模型分析 一、实验目的 1.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验仪器与软件 1. MATLAB7.6 环境 三、实验内容 1 、模型转换 图 1、模型转换示意图及所用命令 传递函数一般形式： MATLAB 表示为： G=tf(num,den)，，其中 num,den 分别是上式中分子，分母系数矩阵。 零极点形式： MATLAB 表示为：G=zpk(Z,P,K) ，其中 Z，P ，K 分别表示上式中的零点矩阵，极点矩阵和增益。 传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)； 状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu 表示对系统的第 iu 个输入量求传递函数；对单输入 iu 为 1。

例1：已知系统的传递函数为G（S）= 2 2 3 24 11611 s s s s s ++ +++ ,利用matlab将传递函数 和状态空间相互转换。 解：1.传递函数转换为状态空间模型： NUM=[1 2 4];DEN=[1 11 6 11]; [A,B,C,D] = tf2ss(NUM,DEN) 2.状态空间模型转换为传递函数： A=[-11 -6 -11;1 0 0;0 1 0];B=[1;0;0];C=[1 2 4];D=[0];iu=1; [NUM,DEN] = ss2tf(A,B,C,D,iu); G=tf(NUM,DEN) 2 、状态方程状态解和输出解 单位阶跃输入作用下的状态响应： G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x). 零输入响应 [y,t,x]=initial(G,x0)其中，x0 为状态初值。

哈工大_控制系统实践_磁悬浮实验报告

研究生自动控制专业实验 地点：A区主楼518房间 姓名：实验日期：年月日斑号：学号：机组编号： 同组人：成绩：教师签字：磁悬浮小球系统 实验报告 主编：钱玉恒，杨亚非 哈工大航天学院控制科学实验室

磁悬浮小球控制系统实验报告 一、实验内容 1、熟悉磁悬浮球控制系统的结构和原理； 2、了解磁悬浮物理模型建模与控制器设计； 3、掌握根轨迹控制实验设计与仿真； 4、掌握频率响应控制实验与仿真； 5、掌握PID控制器设计实验与仿真； 6、实验PID控制器的实物系统调试； 二、实验设备 1、磁悬浮球控制系统一套 磁悬浮球控制系统包括磁悬浮小球控制器、磁悬浮小球实验装置等组成。在控制器的前部设有操作面板，操作面板上有起动/停止开关，控制器的后部有电源开关。 磁悬浮球控制系统计算机部分 磁悬浮球控制系统计算机部分主要有计算机、1711控制卡等； 三、实验步骤 1、系统实验的线路连接 磁悬浮小球控制器与计算机、磁悬浮小球实验装置全部采用标准线连接，电源部分有标准电源线，考虑实验设备的使用便利，在试验前，实验装置的线路已经连接完毕。 2、启动实验装置 通电之前，请详细检察电源等连线是否正确，确认无误后，可接通控制器电源，随后起动计算机和控制器，在编程和仿真情况下，不要启动控制器。 系统实验的参数调试

根据仿真的数据及控制规则进行参数调试（根轨迹、频率、PID 等），直到获得较理想参数为止。 四、实验要求 1、学生上机前要求 学生在实际上机调试之前，必须用自己的计算机，对系统的仿真全部做完，并且经过老师的检查许可后，才能申请上机调试。 学生必须交实验报告后才能上机调试。 2、学生上机要求 上机的同学要按照要求进行实验，不得有违反操作规程的现象，严格遵守实验室的有关规定。 五、系统建模思考题 1、系统模型线性化处理是否合理，写出推理过程？ 合理，推理过程： 由级数理论，将非线性函数展开为泰勒级数。由此证明，在平衡点)x ,(i 00对 系统进行线性化处理是可行的。 对式2x i K x i F )(),(=作泰勒级数展开，省略高阶项可得： )x -)(x x ,(i F )i -)(i x ,(i F )x ,F(i x)F(i,000x 000i 00++= )x -(x K )i -(i K )x ,F(i x)F(i,0x 0i 00++= 平衡点小球电磁力和重力平衡，有 (,)+=F i x mg 0 |,δδ===00 i 00 i i x x F(i,x) F(i ,x )i ；|,δδ===00x 00i i x x F(i,x)F (i ,x )x 对2 i F(i,x )K()x =求偏导数得：

现代控制理论实验报告

现代控制理论实验报告

实验一系统能控性与能观性分析 一、实验目的 1.理解系统的能控和可观性。 二、实验设备 1.THBCC-1型信号与系统·控制理论及计算机控制技术实验平台； 三、实验容 二阶系统能控性和能观性的分析 四、实验原理 系统的能控性是指输入信号u对各状态变量x的控制能力，如果对于系统任意的初始状态，可以找到一个容许的输入量，在有限的时间把系统所有的状态引向状态空间的坐标原点，则称系统是能控的。 对于图21-1所示的电路系统，设iL和uc分别为系统的两个状态变量，如果电桥中 则输入电压ur能控制iL和uc状态变量的变化，此时，状态是能控的。反之，当 时，电桥中的A点和B点的电位始终相等，因而uc不受输入ur的控制，ur只能改变iL的大小，故系统不能控。 系统的能观性是指由系统的输出量确定所有初始状态的能力，如果在有限的时间根据系统的输出能唯一地确定系统的初始状态，则称系统能观。为了说明图21-1所示电路的能观性，分别列出电桥不平衡和平衡时的状态空间表达式： 平衡时：

由式（2）可知，状态变量iL和uc没有耦合关系，外施信号u只能控制iL的变化，不会改变uc的大小，所以uc不能控。基于输出是uc，而uc与iL无关连，即输出uc中不含有iL的信息，因此对uc的检测不能确定iL。反之式（1）中iL与uc有耦合关系，即ur的改变将同时控制iL和uc的大小。由于iL与uc的耦合关系，因而输出uc的检测，能得到iL 的信息，即根据uc的观测能确定iL（ω） 五、实验步骤 1．用2号导线将该单元中的一端接到阶跃信号发生器中输出2上，另一端接到地上。将阶跃信号发生器选择负输出。 2．将短路帽接到2K处，调节RP2，将Uab和Ucd的数据填在下面的表格中。然后将阶跃信号发生器选择正输出使调节RP1，记录Uab和Ucd。此时为非能控系统，Uab和Ucd没有关系（Ucd始终为0）。 3．将短路帽分别接到1K、3K处，重复上面的实验。 六、实验结果 表20－1Uab与Ucd的关系 Uab Ucd

哈工大计算机网络实验报告之五

计算机网络课程实验报告 实验5：利用Ethereal分析TCP、UDP、ICMP协议 继续学习Ethereal的使用； 利用Ethereal分析TCP、UDP和ICMP协议。 TCP协议采用了哪些机制保证可靠数据传输。（3分） 数据重传和数据确认应答机制 Traceroute的工作过程，用自己的话来描述，200字以内，超过酌情扣分。 （4分） 构造数据包，来检查到达一个主机时经过了哪些路由。主机发送给目的地址的数据包的TTL是从1逐个递增的，而数据包每到达一个路由器，它的TTL值就会减1，当TTL减到0时，该数据包被取消，传回一个数据包给主机，我们就能捕获这个路由器的IP地址了。如果收到"超时错"，表示刚刚到达的是路由器，而如果收到的是"端口不可达" 错误，表示刚刚到达的就是目的主机，路由跟踪完成，程序结束。 阐述一下为什么应用程序开发者会选择将应用程序运行在UDP而不是TCP 之上？（3分） UDP没有拥塞控制机制，发送方可以以任何速率向下层注入数据。很多实时应用是

可以容忍一定的数据丢失的，同时又对速率有很高要求（比如在线视频播放），这时开发者会倾向选择UDP协议，避免使用TCP协议的拥塞控制机制产生的分组开销。 实验过程： 使用Ethereal分析TCP协议: (15分)得分：抓取本机与https://www.sodocs.net/doc/605776918.html,/ethereal-labs/alice.txt通信过程中的网络数据包。根据操作思考以下问题： 客户服务器之间用于初始化TCP连接的TCP SYN报文段的序号（sequence number）是多少？在该报文段中，是用什么来标示该报文段是SYN报文段的？ Seq=0 Flags中的syn位为1，ack位为0，说明是syn报文段 服务器向客户端发送的SYNACK报文段序号是多少？该报文段中，Acknowledgement字段的值是多少？https://www.sodocs.net/doc/605776918.html,服务器是如何决定此值 的？在该报文段中，是用什么来标示该报文段是SYNACK报文段的？ Seq=0 Ack=1,服务器根据客户端发送的SYN报文的Seq值加一后得到此值 Flags中的Ack和Syn位都为1，所以是SYNACK报文

现代控制理论综合设计报告—你懂得

《现代控制理论综合设计报告》 问题重述： 图示为单倒立摆系统的原理图，其中摆的长度l=1m，质量m=0.1kg，通过铰链安装小车上，小车质量M=1kg，重力加速度g=9.8m/s2。控制的目的是当小车在水平方向上运动时，将倒立摆保持在垂直位置上。 分别列写小车水平方向的力平衡方程和摆的转矩平衡方程，通过近似线性化处理建立系统的状态空间表达式； 绘制带状态观测器状态反馈系统的模拟仿真图，要求系统期望的特征值为：-1，-2，-1+j，-1-j；状态观测器的特征值为：-2，-3，-2+j，-2-j； 根据模拟仿真图，分别绘制系统综合前后的零输入响应曲线 本文的仿真实验亮点如下： ●对单倒立摆进行传统的传递函数、状态空间建模，全面分析了单倒立摆的物理性质。 ●在物理模型建立时，强调了角速度θ不能近似为0。 ●建立状态空间表达时，选择位移x和角度θ作为输出，是一个多输出系统。但增加了状 态观测器设计的复杂度。 ●在摆运动过程中，初始扰动角θ可达60度左右；而且调节过程中，倒立摆θ在（-90,90） 范围内变化，符合实际情况。 ●在仿真波形图中，展示了状态观测器的跟踪过程，体现了其在反馈控制中起到的作用。 ●在初始扰动60度下，分别在原始系统、状态反馈系统、带状态观测器反馈系统，进行 了零输入响应、阶跃输入响应的仿真实验。 ●解释了带状态观测器反馈时，阶跃输入，但系统前1秒处于稳态的现象的原因。

1单级倒立摆数学模型的建立 倒立摆系统是一个典型的非线性、强耦合、多变量和不稳定系统，作为控制系统的被控对象，许多抽象的控制概念都可以通过倒立摆直观地表现出来。本设计是以一阶倒立摆为被控对象来进行设计的。 传递函数法：对SISO 系统进行分析设计，在这个系统中θ作为输出，因为它比较直观，作用力u 作为输入。 状态空间法：状态空间法可以进行单输入多输出系统设计，因此在这个实验中，我们将尝试同时对摆杆角度和小车位置进行控制，并给小车加一个阶跃输入信号。 本文利用Matlab ，对系统的传递函数和状态空间进行分析，并用指令计算状态空间的各种矩阵，仿真系统的开环阶跃响应。Matlab 将会给出系统状态空间方程的A,B,C 和D 矩阵，并绘出在给定输入为阶跃信号时系统的响应曲线。 在忽略了空气阻力、各种摩擦之后，可将直线一级倒立摆系统抽象成小车和匀质杆组成的系统。 假设系统内部各相关参数为： φ和θ都表示摆杆与垂直向上方向的夹角 l L 、都表示 摆杆长度 1m M 小车质量 1kg m 摆杆质量 0.1kg x 小车位置 单倒立摆系统力的平衡方程分析 小车、摆杆力的分析图如下所示： 小车的平衡方程：u H Mx -= 摆杆的X 轴方向力的平衡方程：2 2(sin )d H m x l dt θ=+ 摆杆Y 轴方向，力的平衡方程：2 2(lcos )d V mg m dt θ-= 摆杆的转矩平衡方程：sin cos VL HL I θθθ-= 选择摆杆的质心在端点处，则惯性惯量2 12ml I = 方程的线性化处理 当θ很小时，可对方程进行线性化。由于控制的目的当小车在水平方向上运动时，将倒立摆保持在垂直位置上。在施加合适的外力下，θ比较小，接近于0，sin ,cos 1θθθ→→,对以 上方程进行线性化。但要注意的是，θ不能约等于0，因为摆杆的角速度在实际情况中是比较快的。但对以上方程先求导会产生θ及其平方项，但这些项都和sin θ相乘，于是这些项还是约等于0。另外，如果先线性化，再求导，则不会产生以上需要考虑的问题。线性化后方程如下：

现代控制理论实验报告河南工业大学

河南工业大学 现代控制理论实验报告姓名：朱建勇 班级：自动1306 学号：201323020601

现代控制理论 实验报告 专业: 自动化 班级: 自动1306 姓名: 朱建勇 学号: 201323020601 成绩评定: 一、实验题目： 线性系统状态空间表达式的建立以及线性变换 二、实验目的 1. 掌握线性定常系统的状态空间表达式。学会在MATLAB 中建立状态空间模型的方法。 2. 掌握传递函数与状态空间表达式之间相互转换的方法。学会用MATLAB 实现不同模型之 间的相互转换。 3. 熟悉系统的连接。学会用MATLAB 确定整个系统的状态空间表达式和传递函数。 4. 掌握状态空间表达式的相似变换。掌握将状态空间表达式转换为对角标准型、约当标准 型、能控标准型和能观测标准型的方法。学会用MATLAB 进行线性变换。 三、实验仪器 个人笔记本电脑 Matlab R2014a 软件 四、实验内容 1. 已知系统的传递函数 (a) ) 3()1(4)(2++=s s s s G

(b) 3486)(22++++=s s s s s G

(c) 6 1161)(232+++++=z z z z z z G （1）建立系统的TF 或ZPK 模型。 （2）将给定传递函数用函数ss( )转换为状态空间表达式。再将得到的状态空间表达式用函 数tf( )转换为传递函数，并与原传递函数进行比较。 （3）将给定传递函数用函数jordants( )转换为对角标准型或约当标准型。再将得到的对角 标准型或约当标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。 （4）将给定传递函数用函数ctrlts( )转换为能控标准型和能观测标准型。再将得到的能控标 准型和能观测标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。

哈工大天线实验报告

Harbin Institute of Technology 天线原理实验报告 课程名称：天线原理 班级： 姓名： 学号： 同组人： 指导教师： 实验时间： 实验成绩： 注：本报告仅供参考 哈尔滨工业大学

一、实验目的 1. 掌握喇叭天线的原理。 2. 掌握天线方向图等电参数的意义。 3. 掌握天线测试方法。 二、实验原理 1. 天线电参数 (1).发射天线电参数 a.方向图：天线的辐射电磁场在固定距离上随空间角坐标分布的图形。 b.方向性系数：在相同辐射功率，相同距离情况下，天线在该方向上的辐射功率密度Smax与无方向性天线在该方向上的辐射功率密度S0之比值。 c.有效长度：在保持该天线最大辐射场强不变的条件下，假设天线上的电流均匀分布时的等效长度。 d.天线效率：表征天线将高频电流或导波能量转换为无线电波能量的有效程度。 e.天线增益：在相同输入功率、相同距离条件下，天线在最大辐射方向上的功率密度Smax与无方向性天线在该方向上的功率密度S0之比值。 f.输入阻抗：天线输入端呈现的阻抗值。 g.极化：天线的极化是指该天线在给定空间方向上远区无线电波的极化。 h.频带宽度：天线电参数保持在规定的技术要求范围内的工作频率范围。 (2).接收天线电参数:除了上述参数以外，接收天线还有一些特有的电参数：等效面积和等效噪声温度。 a.等效面积：天线的极化与来波极化匹配，且负载与天线阻抗共轭匹配的最佳状态下，天线在该方向上所接收的功率与入射电波功率密度之比。 b.等效噪声温度：描述天线向接收机输送噪声功率的参数。 2. 喇叭天线 由逐渐张开的波导构成，是一种应用广泛的微波天线。按口径形状可分为矩形喇叭天线与圆形喇叭天线等。波导终端开口原则上可构成波导辐射器，由于口径尺寸小，产生的波束过宽；另外，波导终端尺寸的突变除产生高次模外，反射较大，与波导匹配不良。为改善这种情况，可使波导尺寸加大，以便减少反射，又可在较大口径上使波束变窄。 (1).H面扇形喇叭：若保持矩形波导窄边尺寸不变，逐渐张开宽边可得H面扇

现代控制理论实验报告

现代控制理论实验报告 组员： 院系：信息工程学院 专业： 指导老师： 年月日

实验1 系统的传递函数阵和状态空间表达式的转换 [实验要求] 应用MATLAB 对系统仿照[例]编程，求系统的A 、B 、C 、阵；然后再仿照[例]进行验证。并写出实验报告。 [实验目的] 1、学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法； 2、通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数相互转换方法。 [实验内容] 1 设系统的模型如式示。 p m n R y R u R x D Cx y Bu Ax x ∈∈∈?? ?+=+=& 其中A 为n ×n 维系数矩阵、B 为n ×m 维输入矩阵 C 为p ×n 维输出矩阵，D 为传递阵，一般情况下为0，只有n 和m 维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式示。 D B A SI C s den s num s G +-== -1)() () (()( 式中，)(s num 表示传递函数阵的分子阵，其维数是p ×m ；)(s den 表示传递函数阵的按s 降幂排列的分母。 2 实验步骤 ① 根据所给系统的传递函数或（A 、B 、C 阵），依据系统的传递函数阵和状态空间表达式之间的关系如式，采用MATLA 的编程。注意：ss2tf 和tf2ss 是互为逆转换的指令； ② 在MATLA 界面下调试程序，并检查是否运行正确。 ③ [] 已知SISO 系统的状态空间表达式为，求系统的传递函数。

, 2010050010000100001 0432143 21u x x x x x x x x ? ? ??? ? ??????-+????????????????????????-=????????????&&&&[]??? ? ? ???????=43210001x x x x y 程序: A=[0 1 0 0;0 0 -1 0;0 0 0 1;0 0 5 0]; B=[0;1;0;-2]; C=[1 0 0 0]; D=0; [num,den]=ss2tf(A,B,C,D,1) 程序运行结果: num = 0 den = 0 0 0 从程序运行结果得到:系统的传递函数为: 2 4253 )(s s s S G --= ④ [] 从系统的传递函数式求状态空间表达式。 程序: num =[0 0 1 0 -3]; den =[1 0 -5 0 0]; [A,B,C,D]=tf2ss(num,den) 程序运行结果: A = 0 5 0 0 1 0 0 0 0 1 0 0

现代控制理论课程报告

现代控制理论课程总结 学习心得 从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程，在刚拿到课本的时候，没上张老师的课之前，咋一看，会认为开课的内容会是上学期学的控制理论基础的累赘或者简单的重复，更甚至我还以为是线性代数的复现呢！根本没有和现代控制论联系到一起。但后面随着老师讲课的风格的深入浅出，循循善诱，发现和自己想象的恰恰相反，张老师以她特有的讲课风格，精心准备的ppt 课件，向我们展示了现代控制理论发展过程，以及该掌握内容的方方面面，个人觉得，我们不仅掌握了现代控制理论的理论知识，更重要的是学会了掌握这门知识的严谨的逻辑思维和科学的学习方法，对以后学习其他知识及在工作上的需要大有裨益，总之学习了这门课让我受益匪浅。 由于我们学习这门课的课时不是很多，并结合我们学生学习的需求及所要掌握的课程深入程度，张老师根据我们教学安排需要，我们这学期学习的内容主要有：1.绪论；2.控制系统的状态表达式；3.控制系统状态表达式的解；4.线性系统的能空性和能观性；5.线性定常系统的综合。而状态变量和状态空间表达式、状态转移矩阵、系统的能控性与能观性以及线性定常系统的综合是本门课程的主要学习内容。当然学习的内容还包括老师根据多年教学经验及对该学科的研究的一些深入见解。 在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的必修课。 经典控制理论的特点 经典控制理论以拉氏变换为数学工具，以单输入－单输出的线性定常系统为主要的研究对象。将描述系统的微分方程或差分方程变换到复数域中，得到系统的传递函数，并以此作为基础在频率域中对系统进行分析和设计，确定控制器的结构和参数。通常是采用反馈控制，构成所谓闭环控制系统。经典控制理论具有明显的局限性，突出的是难以有效地应用于时变系统、多变量系统，也难以揭示系统更为深刻的特性。当把这种理论推广到更为复杂的系统时，经典控制理论就显得无能为力了，这是因为它的以下几个特点所决定。 1．经典控制理论只限于研究线性定常系统，即使对最简单的非线性系统也是无法处理的；这就从本质上忽略了系统结构的内在特性，也不能处理输入和输出皆大于1的系统。实际上，大多数工程对象都是多输入－多输出系统，尽管人们做了很多尝试，但是，用经典控制理论设计这类系统都没有得到满意的结果；2．经典控制理论采用试探法设计系统。即根据经验选用合适的、简单的、工程上易于实现的控制器，然后对系统进行分析，直至找到满意的结果为止。虽然这种设计方法具有实用等很多完整，从而促使现代控制理论的发展：对经典理论的精确化、数学化及理论化。优点，但是，在推理上却是不能令人满意的，效果也

(完整版)哈工大-数值分析上机实验报告

实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。 前言：（目的和意义） 掌握二分法与Newton法的基本原理和应用。 数学原理： 对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。 对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。

Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式 产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。 程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x); 写成如上形式即可，下面给出主程序。 二分法源程序： clear %%%给定求解区间 b=1.5; a=0;

%%%误差 R=1; k=0;%迭代次数初值 while (R>5e-6) ; c=(a+b)/2; if f12(a)*f12(c)>0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);

现代控制理论实验报告3

实验三 利用MATLAB 导出连续状态空间模型的离散化模型 实验目的： 1、基于对象的一个连续时间状态空间模型，导出其相应的离散化状态空间模型； 2、通过编程、上机调试，掌握离散系统运动分析方法。 实验原理： 给定一个连续时间系统的状态空间模型： ()()()()()() x t Ax t Bu t y t Cx t Du t =+=+ （3.1） 状态空间模型（3.1）的输入信号()u t 具有以下特性： ()(),u t u kT kT t kT T =≤≤+ （3.2） 已知第k 个采样时刻的状态()x kT 和第k 个采样时刻到第1k +个采样时刻间的输入()()u t u kT =，可得第1k +个采样时刻(1)k T +处的状态 (1)((1))((1))()((1))()k T kT x k T k T kT x kT k T Bu d τττ++=Φ+-+Φ+-? （3.3） 其中： ((1))((1))A k T kT AT k T kT e e +-Φ+-== ((1))((1))A k T k T e ττ+-Φ+-= 由于输入信号在两个采样时刻之间都取常值，故对式（3.3）中的积分式进行一个时间变量替换(1)k T στ=+-后，可得 0((1))()()()AT A x k T e x kT e d Bu kT τ σσ+=+? （3.4） 另一方面，以周期T 对输出方程进行采样，得到 ()()()y kT Cx kT Du kT =+ 在周期采样的情况下，用k 来表示第k 个采样时刻kT 。因此，连续时间状态空间模型

（3.1）的离散化方程可以写成 (1)()()()()()()() x k G T x k H T u k y k Cx k Du k +=+=+ （3.5） 其中： 0()()()AT A G T e H T e d B τσσ==? （3.6） 已知系统的连续时间状态空间模型，MATLAB 提供了计算离散化状态空间模型中状态矩阵和输入矩阵的函数： [G ,H]=c2d(A,B,T) 其中的T 是离散化模型的采样周期。 实验步骤 1、导出连续状态空间模型的离散化模型，采用MA TLAB 的m-文件编程； 2、在MA TLAB 界面下调试程序，并检查是否运行正确。 例3.1 已知一个连续系统的状态方程是 010()()()2541x t x t u t ????=+????--???? 若取采样周期0.05T =秒，试求相应的离散化状态空间模型。 编写和执行以下的m-文件： A=[0 1;-25 –4]; B=[0;1]; [G ,H]=c2d(A,B,0.05) 得到 G= 0.9709 0.0448 -1.1212 0.7915 H= 0.0012 0.0448 因此，所求的离散化状态空间模型是 0.97090.04480.0012(1)()()1.12120.79150.0448x k x k u k ????+=+????-????

哈工大威海计算机网络实验报告1资料

计算机网络与通信实验报告（一）学号姓名班级报告日期 2015.04.15 实验内容网络常用命令的使用 实验目的1.熟悉网络命令的使用，例如ping，tracert，netstat，ipconfig等，对结果进行分析判断。 2.熟悉dns的层次查询，以及smtp协议。 实验预备知识结合实验报告相关知识以及老师课堂演示、笔记。 实验过程描述1.按照实验报告步骤所指，一步步熟悉ping tracert ipconfig 等网络命令，并对结果进行相应分析、截图。 2.Dns层次查询时，首先网上搜索全球13 个根域名服务器的ip，选择其中一个ip 对学校主页https://www.sodocs.net/doc/605776918.html, 进行层次分析，依次进行cn https://www.sodocs.net/doc/605776918.html, https://www.sodocs.net/doc/605776918.html, https://www.sodocs.net/doc/605776918.html, 的域名分析，最终得到主页ip，然后使用ping命令ping得主页ip 相比较，结果一致，查询成功。 3.熟悉掌握SMTP协议。Dos 命令下依次输入telnet相关命令，并使用事先转换成base64 的用户名、密码登陆邮箱。登陆成功后给自己的邮箱发送信息，最后退出。操作、邮箱截图如下。 实验结果见表格下方截图。 实验当中问题及解决方法1、telnet命令刚开始dos无法识别，属于不认识的命令。上网查询资料后，在控制面板中设置后成功解决。 2、熟悉SMTP协议时，telnet 登陆邮箱并发送信件，期间出现好多错误，比如单词拼写错误，指令错误。重复多次后最终成功实现。 成绩（教师打分）优秀良好及格不及格

实验相关截图 一、网络命令的使用 1.ping 命令

2.tracert 命令

利用MATLAB设计状态观测器—现代控制理论实验报告

实验六利用MATLAB设计状态观测器 ******* 学号 1121*****

实验目的： 1、学习观测器设计算法； 2、通过编程、上机调试，掌握基于观测器的输出反馈控制系统设计方法。 实验原理： 1、全阶观测器模型： () ()x Ax Bu L y Cx A LC x Bu Ly =++-=-++ 由极点配置和观测器设计问题的对偶关系，也可以应用MATLAB 中极点配置的函数来确定所需要的观测器增益矩阵。例如，对于单输入单输出系统，观测器的增益矩阵可以由函数 L=(acker(A ’,C ’,V))’ 得到。其中的V 是由期望的观测器极点所构成的向量。类似的，也可以用 L=(place(A ’,C ’,V))’ 来确定一般系统的观测器矩阵，但这里要求V 不包含相同的极点。 2、降阶观测器模型： ???w Aw By Fu =++ b x w Ly =+ 基于降阶观测器的输出反馈控制器是： ????()[()]()b a b b a b w A FK w B F K K L y u K w K K L y =-+-+=--+ 对于降阶观测器的设计，使用MATLAB 软件中的函数 L=(acker(Abb’,Aab’,V))’ 或 L=(place(Abb’,Aab’,V))’ 可以得到观测器的增益矩阵L 。其中的V 是由降阶观测器的期望极点所组成的向量。 实验要求 1．在运行以上例程序的基础上，考虑图6.3所示的调节器系统，试针对被控对象设计基于全阶观测器和降 阶观测器的输出反馈控制器。设极点配置部分希望的闭环极点是1,22j λ=-± （a ） 对于全阶观测器，1 8μ=-和 28μ=-； （b ） 对于降阶观测器，8μ=-。 比较系统对下列指定初始条件的响应： （a ） 对于全阶观测器： 1212(0)1,(0)0,(0)1,(0)0x x e e ==== （b ） 对于降阶观测器： 121(0)1,(0)0,(0)1x x e === 进一步比较两个系统的带宽。

现代控制理论结课论文

现代控制理论方法综述 研电1610秦晓 1162201332 摘要：本文将控制理论方法分为现代控制理论基础，线性最优控制，非线性最优控制三大部分，查阅文献，综述了每一部分中的经典控制方法，以及每种控制方法的优缺点和在工业中的应用，最后提出了目前在现代控制理论中依旧存在的问题。 1.引言 电力系统是一个复杂的非线性动态大系统，对于这个规模庞大的系统，研究其运行的动态特性进而构建先进的安全控制系统是极富挑战性的课题。同时，各种新技术的应用，一方面增强了系统的调控能力和经济效益，另一方面也极大的增加了电网控制的复杂性，对电力系统的安全稳定运行提出了更严格的要求。因此，改善与提高我国电力系统的动态品质、安全稳定和经济性成为了电力工作者的首要任务。提高电力系统稳定性的最经济和最有效的手段之一是采用先进的控制理论和方法。在过去的时间里，电力工作者们为改进与发展电力系统控制技术进行了大量研究。本文主要梳理总结电力系统在现代控制方面的研究成果，分析了电力系统控制技术的发展趋势，并总结了目前现代控制理论还需要解决的问题。 2.现代控制的基础 现代控制理论的基础是经典控制理论，在20世纪20年代到50年代间，为了满足第二次世界大战前后军事技术和工业发展的需求，经典控制理论有了飞速的发展。经典控制理论主要研究线性时不变、单输入单输出的控制问题。在分析和设计大型反馈控制系统时，经典控制论主要采用频域法，其中以 Nyquist 判据、Bode 图和根轨迹法最为广泛[1~2]。经典控制理论的设计目标是使闭环系统特征方程的特征根全部位于左半开平面上。上述设计目标可以描述为一类无目标函数的优化问题，即约束满足问题。由于使系统稳定的控制器解并不唯一，所以根据经典控制理论设计的 PID 控制器往往带有较大的冗余性[3]。也正是由于经典控制理论设计目标及方向简单明确，计算方便，特别适合需要依赖工程经验或现场测试进行控制器设计的系统，所以至今仍在工业中广泛应用。 在上世纪70年代以前，经典控制是电力系统控制的主流。如发电机励磁控制AVR主要采用单变量反馈方式，即采用发电机端电压偏差作为反馈量的 PID 控制方式。随着发电技术的进步和电力系统自身规模的增长，人们逐渐发现这种单输入控制方式难以满足电力系统对抑制振荡和提高稳定极限方面的要求。最早报道的互联电力系统低频振荡发生于20世纪60年代，北美MAPP的西北联合系统和西南联合系统进行互联试运行时发生了低频振荡，造成联络线过流跳闸[4]。之后，随着大容量机组的不断投运，以及快速、高放大倍数励磁系统越来越广泛的使用，使得低频振荡现象在世界各国大型互联电网中时有发生，这对电网安全产生了严重威胁。为解决这个问题，文献[5]采用转速偏差作为附加反馈与AVR并联，发展出PSS+AVR的励磁控制方式。进入21世纪以来，我国电网互联程度不断提高，系统中出现了

《现代控制理论》实验报告

. 现代控制理论实验报告 组员： 院系：信息工程学院 专业： 指导老师： 年月日

实验1 系统的传递函数阵和状态空间表达式的转换 [实验要求] 应用MATLAB 对系统仿照[例1.2]编程，求系统的A 、B 、C 、阵；然后再仿照[例1.3]进行验证。并写出实验报告。 [实验目的] 1、学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法； 2、通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数相互转换方法。 [实验内容] 1 设系统的模型如式(1.1)示。 p m n R y R u R x D Cx y Bu Ax x ∈∈∈?? ?+=+= (1.1) 其中A 为n ×n 维系数矩阵、B 为n ×m 维输入矩阵 C 为p ×n 维输出矩阵，D 为传递阵，一般情况下为0，只有n 和m 维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式(1.2)示。 D B A SI C s den s num s G +-== -1)() () (()( (1.2) 式(1.2)中，)(s num 表示传递函数阵的分子阵，其维数是p ×m ；)(s den 表示传递函数阵的按s 降幂排列的分母。 2 实验步骤 ① 根据所给系统的传递函数或（A 、B 、C 阵），依据系统的传递函数阵和状态空间表达式之间的关系如式(1.2)，采用MATLA 的file.m 编程。注意：ss2tf 和tf2ss 是互为逆转换的指令； ② 在MATLA 界面下调试程序，并检查是否运行正确。 ③ [1.1] 已知SISO 系统的状态空间表达式为(1.3)，求系统的传递函数。

中南大学现代控制理论实验报告

中南大学 现代控制理论实验报告 指导老师：年晓红、郭宇骞 姓名： 学号： 专业班级： 实验日期： 2015.6.11 学院：信息科学与工程学院

实验1 用MATLAB分析状态空间模型 1、实验设备 PC计算机1台，MATLAB软件1套。 2、实验目的 ①学习系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数 相互转换的方法； ②通过编程、上机调试，掌握系统状态空间表达式与传递函数相互转换方法学 习系统齐次、非齐次状态方程求解的方法，计算矩阵指数，求状态响应； ③通过编程、上机调试，掌握求解系统状态方程的方法，学会绘制状态响应曲 线； ④掌握利用MATLAB导出连续状态空间模型的离散化模型的方法。 3、实验原理说明 参考教材P56~59“2.7 用MATLAB分析状态空间模型” 参考教材P99~101“3.8 利用MATLAB求解系统的状态方程” 4、实验步骤 ①根据所给系统的传递函数或A、B、C矩阵，依据系统的传递函数阵和状 态空间表达式之间的关系式，采用MATLAB编程。 ②在MATLAB界面下调试程序，并检查是否运行正确。 ③根据所给系统的状态方程，依据系统状态方程的解的表达式，采用 MATLAB编程。

④ 在MATLAB 界面下调试程序，并检查是否运行正确。 5、实验习题 题1.1 已知SISO 系统的传递函数为 2432 58 ()2639 s s g s s s s s ++=++++ （1）将其输入到MATLAB 工作空间； （2）获得系统的状态空间模型。 解： (1) num=[1,5,8] ; den=[1,2,6,3,9] ; G=tf(num , den) Transfer function: s^2 + 5 s + 8 ----------------------------- s^4 + 2 s^3 + 6 s^2 + 3 s + 9 (2) G1=ss(G) a = x1 x2 x3 x4 x1 -2 -1.5 -0.75 -2.25 x2 4 0 0 0 x3 0 1 0 0