1. 控制系统任务的物理描述
为了满足飞机品质的要求，飞机的纵向运动和横侧向运动都需要有能够连续工作的阻尼器，以用来调整飞机的飞行姿态，避免其出现不必要的俯仰和倾斜。维持飞机纵向运动的阻尼器称为俯仰阻尼器，维持飞机横侧向运动的阻尼器称为偏航阻尼器。本次课程大作业旨在通过运用Matlab 的经典控制系统设计工具对某型飞机偏航阻尼器进行控制系统的设计。
2. 控制系统对象的数学模型
巡航状态下，某型飞机侧向运动的状态空间模型为：
111
12131411122212223242122131
3233343132234142434441424()1()()()()2()()()3()()4t x t a a a a b b t x t a a a a b b u t a a a a b b u t x t t a a a a b b x t t x x x x ??????????????????????????????????=+????????????????????????
?????????? 111121314122122
2324234()()()()()()x t c c c c y t x t c c c c y t x t x t ??????????=??????????????
式中： 1()x t ：侧滑角（单位为rad ）
2()x t ：偏航角速度（单位为/rad s ）
3()x t ：滚转角速度（单位为/rad s ）
4()x t ：倾斜角（单位为rad ）
输入向量及输出向量分别为：
1()u t ：方向舵偏角（单位为rad ）
2()u t ：副翼偏角（单位为rad ）
1()y t ：偏航角速度（单位为/rad s ）
2()y t ：倾斜角（单位为rad ）
设飞机巡航飞行时的速度为0.8马赫，高度为40000英尺，此时模型的参数为：
0.055800.99680.08020.04150.598000.11500.031803.050000.38800.4650000.0805 1.00000A --????--??=??--????
0.007290.00000.475000.007750.15300.1430000B ????-??=?????
? 0 1.000000 1.0C ??=????
0000D ??=????
首先输入飞机状态空间模型参数。以及定义系统的状态变量、输入变量及输出变量，并建立状态空间模型。
在Matlab 命令窗口中输入如下命令：
>> A=[-0.0558 -0.9968 0.0802 0.0415; 0.5980 -0.1150 -0.0318 0; -3.0500
0.3880 -0.4650 0; 0 0.0805 1.000 0];
>> B=[0.00729 0.0000; -0.47500 0.00775; 0.15300 0.1430;0 0];
>> C=[0 1 0 0; 0 0 0 1];
>> D=[0 0; 0 0];
>> states={'beta','yaw','roll','phi'};
>> inputs={'rudder','aileron'};
>> outputs={'yaw rate','bank angle'};
>> sys=ss(A,B,C,D,'statename',states,'inputname',inputs,'outputname',outputs) 运行结果如图2-1所示：
图2- 1 状态空间模型
3.系统特性分析
根据前述系统的状态空间模型，首先分析系统的性能。
3.1.计算开环特征值
在Matlab中计算系统开环特征值，输入：
damp(sys)
所输入系统命令及运行结果如图3-1所示：
图3- 1 系统开环特征值
绘制零极点图，在Matlab的命令窗口中输入：
pzmap（sys）
运行结果如图3-2所示：
图3- 2 零极点图
由图可以看出，此模型含有接近虚轴的一对共轭极点，它们对应飞机的荷兰滚模态，此时，系统具有较小的阻尼，控制系统设计的目的是提高系统的阻尼比，
改善荷兰滚模态的阻尼特性。
3.2.计算系统的单位脉冲响应
在Matlab命令窗口输入：
impulse（sys）
运行后得到如图3-3所示的单位脉冲相应曲线：
图3- 3 单位脉冲对应曲线
由图可以看出，系统过渡过程振荡剧烈，飞机确实存在很小的阻尼，图中相应时间较长，而乘客及飞行员关心的是飞机在最初的几秒钟的行为，所以绘制飞机在最初的20s以内的单位脉冲响应曲线。
在Matlab命令窗口中输入：
impulse(sys,20)
所得响应时间为20秒的单位脉冲相应图形如下：
图3- 4 20s单位脉冲对应图形
由上图可以看出，飞机围绕非零倾角产生了震荡，因此在副翼脉冲信号作用下，飞机会发生改变。
图3- 5 bode图形
使用方向舵偏角作为控制输入，使用偏航角速度作为传感输入，为得到相应
的频率响应，在Matlab 命令窗口中输入如下命令：
s ys11=sys(‘yaw’,’rudder’);
bode(sys11)
运行后的bode 图如图3-5所示。
由图可以看出，方向舵的变化对小阻尼的荷兰滚模态具有明显的影响。
4. 控制系统的指标
一种比较合理的设计目标是确保自然频率s rad n /0.1<ω时，阻尼比30.0≥ξ。
5. 控制系统的设计
通过以上分析可知，只要通过改变系统的增益，就可确保系统性能得到改善。首先，应用根轨迹法确定合适的增益值。
在Matlab 命令窗口中输入：
rlocus(sys11)
运行后得到的曲线即为负反馈的根轨迹图。所得图形图5-1所示：
图5- 1 负反馈根轨迹图
由图可见，采用负反馈连接会使得系统立刻变得不稳定，为确保系统稳定，应当采用正反馈连接。
在Matlab窗口中进一步输入：
rlocus(-sys11)
sgrid
运行后得到正反馈的根轨迹图如图5-2所示：
图5- 2 正反馈根轨迹图
然后继续构成单输入单输出闭环反馈回路，在Matlab命令窗口中输入如下命令：
k=2.85；
cl11=feedback(sys11,-k);
运行后得到负反馈系统cl11，如图5-3所示：
图5- 3 负反馈系统
由下述的Matlab命令求取系统响应时间为20s的单位脉冲相应，并将其与前述的开环系统单位脉冲响应作比较。
图5- 4 响应曲线
在Matlab命令窗口中输入：
impulse (sys11,cl11,’o-‘,20)
运行后得到如图5-4所示的闭环系统的单位脉冲响应曲线。
由上图可以看出，与开环系统单位脉冲相应相比，闭环系统响应速度快，并且没有产生很大的震荡。
将全部多输入多输出模型构成闭合回路，分析在副翼输入信号作用下的响应。将系统由输入1连至输出1，构成反馈回路，在Matlab中输入如下指令：
cloop=feedback(sys,-k,1,1);
damp(cloop)
运行结果如图5-5所示：
图5- 5 反馈回路
图5- 6 脉冲响应曲线
绘制多输入多输出模型的脉冲响应曲线，在Matlab命令窗口输入：
北航最优化方法大作业参考
1 流量工程问题 1.1 问题重述 定义一个有向网络G=(N,E)，其中N是节点集，E是弧集。令A是网络G的点弧关联矩阵，即N×E阶矩阵，且第l列与弧里(I,j)对应，仅第i行元素为1，第j行元素为-1，其余元素为0。再令b m=(b m1,…,b mN)T，f m=(f m1,…,f mE)T，则可将等式约束表示成： Af m=b m 本算例为一经典TE算例。算例网络有7个节点和13条弧，每条弧的容量是5个单位。此外有四个需求量均为4个单位的源一目的对，具体的源节点、目的节点信息如图所示。这里为了简单，省区了未用到的弧。此外，弧上的数字表示弧的编号。此时，c=((5,5…,5)1 )T， ×13 根据上述四个约束条件，分别求得四个情况下的最优决策变量x=((x12,x13,…,x75)1× )。 13 图 1 网络拓扑和流量需求
1.2 7节点算例求解 1.2.1 算例1（b1=[4;-4;0;0;0;0;0]T） 转化为线性规划问题： Minimize c T x1 Subject to Ax1=b1 x1>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x1*=[4 0 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x1=20 1.2.2 算例2（b2=[4;0;-4;0;0;0;0]T） Minimize c T x2 Subject to Ax2=b2 X2>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x2*=[0 4 0 0 0 0 0 0 0 0 0 0 0]T 对应的最优值c T x2=20 1.2.3 算例3（b3=[0;-4;4;0;0;0;0]T） Minimize c T x3 Subject to Ax3=b3 X3>=0 利用Matlab编写对偶单纯形法程序，可求得: 最优解为x3*=[4 0 0 0 4 0 0 0 0 0 0 0 0]T 对应的最优值c T x3=40
欲索取更多考研资料，请上北京天问教育网站官网！ 自动控制原理综合 自动化科学与电气工程学院 2007年11月
931自动控制原理综合考试大纲（2008版） 一、考试组成 自动控制原理占90分; 理论力学占60分； 二、自动控制原理部分的考试大纲 （一）复习内容及基本要求 1．自动控制的一般概念 主要内容：自动控制的任务；基本控制方式：开环、闭环（反馈）控制；自动控制的性能要求：稳、快、准。 基本要求：反馈控制原理与动态过程的概念；由给定物理系统建原理方块图。 2．数学模型 主要内容：传递函数及动态结构图；典型环节的传递函数；结构图的等效变换、梅逊公式。 基本要求：典型环节的传递函数；闭环系统动态结构图的绘制；结构图的等效变换。 3．时域分析法 主要内容：典型响应及性能指标、一、二阶系统的分析与计算。系统稳定性的分析与计算：劳斯、古尔维茨判据。稳态误差的计算及一般规律。 基本要求：典型响应（以一、二系统的阶跃响应为主）及性能指标计算；系统参数对响应的影响；劳斯、古尔维茨判据的应用；系统稳态误差、终值定理的使用条件。 4．根轨迹法 主要内容：根轨迹的概念与根轨迹方程；根轨迹的绘制法则；广义根轨迹；零、极点分布与阶跃响应性能的关系；主导极点与偶极子。 基本要求：根轨迹法则（法则证明只需一般了解）及根轨迹的绘制；主导极点、偶极子等的概念；利用根轨迹估算阶跃响应的性能指标。 5．频率响应法 主要内容：线性系统的频率响应；典型环节的频率响应及开环频率响应；Nyquist稳定判据和对数频率稳定判据；稳定裕度及计算；闭环幅频与阶跃响应的关系，峰值及频宽的概念；开环频率响应与阶跃响应的关系，三频段（低频段，中频段和高频段）的分析方法。 基本要求：典型环节和开环系统频率响应曲线(Nyquist曲线和对数幅频、相频曲线)的绘制；系统稳定性判据（Nyquist判据和对数判据）；等M、等N圆图，尼柯尔斯图仅作一般了解；相稳定裕度和模稳定裕度的计算；明确最小相位和非最小相位系统的差别，明确截止频率和带宽的概念。 6．线性系统的校正方法 主要内容：系统设计问题概述；串联校正特性及作用：超前、滞后及PID；校正设计的频率法及根轨迹法；反馈校正的作用及计算要点；复合校正原理及其实现。 基本要求：校正装置的作用及频率法的应用；以串联校正为主，反馈校正为辅；以频率法为主，根轨迹法为辅；复合校正的应用。 7．线性连续系统的状态空间分析方法
现代控制理论 (主汽温对象模型) 班级: 学号: 姓名:
目录 一. 背景及模型建立 1.火电厂主汽温研究背景及意义 2.主汽温对象的特性 3.主汽温对象的数学模型 二.分析 1.状态空间表达 2.化为约当标准型状态空间表达式并进行分析 3.系统状态空间表达式的求解 4.系统的能控性和能观性 5.系统的输入输出传递函数 6.分析系统的开环稳定性 7.闭环系统的极点配置 8.全维状态观测器的设计 9.带状态观测器的状态反馈控制系统的状态变量图 10.带状态观测器的闭环状态反馈控制系统的分析 三.结束语 1.主要内容 2.问题及分析 3.评价
一.背景及模型建立 1.火电厂主汽温研究背景及意义 火电厂锅炉主汽温控制决定着机组生产的经济性和安全性。由于锅炉的蒸汽容量非常大、过热汽管道很长，主汽温调节对象往往具有大惯性和大延迟，导致锅炉主汽温控制存在很多方面的问题，影响机组的整个工作效率。主汽温系统是表征锅炉特性的重要指标之一，主汽温的稳定对于机组的安全运行至关重要。其重要性主要表现在以下几个方面： (1) 汽温过高会加速锅炉受热面以及蒸汽管道金属的蠕变，缩短其使用寿命。例如，12CrMoV 钢在585℃环境下可保证其应用强度的时间约为10万小时，而在 595℃时，其保证应用强度的时间可能仅仅是 3 万小时。而且一旦受热面严重超温，管道材料的强度将会急剧下降，最终可能会导致爆管。再者，汽温过高也会严重影响汽轮机的汽缸、汽门、前几级喷嘴和叶片、高压缸前轴承等部件的机械强度，从而导致设备损坏或者使用年限缩短。 (2) 汽温过低，会使得机组循环热效率降低，增大煤耗。根据理论估计可知：过热汽温每降低10℃，会使得煤耗平均增加0.2%。同时，汽温降低还会造成汽轮机尾部的蒸汽湿度增大，其后果是，不仅汽轮机内部热效率降低，而且会加速汽轮机末几级叶片的侵蚀。此外，汽温过低会增大汽轮机所受的轴向推力，不利于汽轮机的安全运行。 (3) 汽温变化过大会使得管材及有关部件产生疲劳，此外还将引起汽轮机汽缸的转子与汽缸的胀差变化，甚至产生剧烈振动，危及机组安全运行。 据以上所述，工艺上对汽温控制系统的质量要求非常严格，一般控制误差范围在±5℃。主汽温太高会缩短管道的使用寿命，太低又会降低机组效率。所以必须实现汽温系统的良好控制。而汽温被控对象往往具有大惯性、大延时、非线性，时变一系列的特性，造成对象的复杂性，增加了控制的难度。现代控制系统中有很多关于主汽温的控制方案，本文我们着重研究带状态观测器的状态反馈控制对主汽温的控制[1] 。 2.主汽温对象的特性 2.1主汽温对象的静态特性 主汽温被控对象的静态特性是指汽温随锅炉负荷变化的静态关系。过热器的传热形式、结构和布置将直接影响过热器的静态特性。现代大容量锅炉多采用对流过热器、辐射过热器和屏式过热器。对流过热器布置在450℃～1000℃烟气温度的烟道中，受烟气的横向和纵向冲刷，烟气以对流方式将热量传给管道。而辐射过热器则是直接吸收火焰和高温烟气的辐射能。屏式过热器布置在炉膛内上部
1.用薄钢板制造一体积5m 3，长度不小于4m ，无上盖的货箱，要求钢板耗量最小。确定货箱的长x 1、宽x 2和高x 3。试列出问题的数学模型。 解：min 32312122x x x x x x z ++= s.t 5321=x x x 41≥x 0,,321≥x x x 2.将下面的线性规划问题表示为标准型并用单纯形法求解 max f=x 1+2x 2+x 3 s ．t ．2x 1+x 2-x 3≤2 -2x 1+x 2-5x 3≥-6 4x 1+x 2+x 3≤6 x i ≥0 i=1，2，3 解：先化标准形： Min 321x x x z -+= 224321=+-+x x x x 6525321=++-x x x x 646321=+++x x x x 列成表格：
1 2 1 610011460105122001112----- 可见此表已具备1°，2°，3°三个特点，可采用单纯形法。首先从底行中选元素-1，由2/2,6/2,6/4最小者决定选第一行第一列的元素2，标以记号，迭代一次得 1 2 1 2102310401162010021212 11-------- 再从底行中选元素-2/3，和第二列正元素1/2，迭代一次得 1 2 12 32 30 210231040116201002121211- ------ 再从底行中选元素-3，和第二列正元素2，迭代一次得 4 2 3 3 410120280114042001112--- 再迭代一次得 10 2 30 2 10 6 221023 1010213000421021013-- 选取最优解：
一、选择题 1.下面关于建模和模型说法错误的是( C )。 A．无论是何种系统，其模型均可用来提示规律或因果关系。 B．建模实际上是通过数据、图表、数学表达式、程序、逻辑关系或各种方式的组合表示状态变量、输入变量、输出变量、参数之间的关系。 C．为设计控制器为目的建立模型只需要简练就可以了。 D．工程系统模型建模有两种途径，一是机理建模，二是系统辨识。 &&&&的类型是( B ) 。 2.系统()3()10() y t y t u t ++= A．集中参数、线性、动态系统。B．集中参数、非线性、动态系统。 C．非集中参数、线性、动态系统。D．集中参数、非线性、静态系统。 3.下面关于控制与控制系统说法错误的是( B )。 A．反馈闭环控制可以在一定程度上克服不确定性。 B．反馈闭环控制不可能克服系统参数摄动。 C．反馈闭环控制可在一定程度上克服外界扰动的影响。 D．控制系统在达到控制目的的同时，强调稳、快、准、鲁棒、资源少省。 x Pz说法错误的是( D )。 4.下面关于线性非奇异变换= A．非奇异变换阵P是同一个线性空间两组不同基之间的过渡矩阵。 B．对于线性定常系统，线性非奇异变换不改变系统的特征值。 C．对于线性定常系统，线性非奇异变换不改变系统的传递函数。 D．对于线性定常系统，线性非奇异变换不改变系统的状态空间描述。 5.下面关于稳定线性系统的响应说法正确的是( A )。 A．线性系统的响应包含两部分，一部是零状态响应，一部分是零输入响应。 B．线性系统的零状态响应是稳态响应的一部分。 C．线性系统暂态响应是零输入响应的一部分。 D．离零点最近的极点在输出响应中所表征的运动模态权值越大。 6.下面关于连续线性时不变系统的能控性与能观性说法正确的是( A ) 。 A．能控且能观的状态空间描述一定对应着某些传递函数阵的最小实现。 B．能控性是指存在受限控制使系统由任意初态转移到零状态的能力。 C．能观性表征的是状态反映输出的能力。 D．对控制输入的确定性扰动影响线性系统的能控性，不影响能观性。 7.下面关于系统Lyapunov稳定性说法正确的是( C ) 。
自动控制原理实验报告 一、实验名称：一、二阶系统的电子模拟及时域响应的动态测试 二、实验目的 1、了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系 2、学习在电子模拟机上建立典型环节系统模型的方法 3、学习阶跃响应的测试方法 三、实验内容 1、建立一阶系统的电子模型，观测并记录在不同时间常数T时的响应曲线，测定过渡过程时间T s 2、建立二阶系统电子模型，观测并记录不同阻尼比的响应曲线，并测定超调量及过渡过程时间T s 四、实验原理及实验数据 一阶系统 系统传递函数： 由电路图可得，取则K=1，T分别取：0.25, 0.5, 1 T 0.25 0.50 1.00 R2 0.25MΩ0.5M Ω1MΩ C 1μ1μ1μ T S 实测0.7930 1.5160 3.1050 T S 理论0.7473 1.4962 2.9927 阶跃响应曲线图1.1 图1.2 图1.3 误差计算与分析 （1）当T=0.25时，误差==6.12%； （2）当T=0.5时，误差==1.32%； （3）当T=1时，误差==3.58% 误差分析：由于T决定响应参数，而,在实验中R、C的取值上可能存在一定误差，另外，导线的连接上也存在一些误差以及干扰，使实验结果与理论值之间存在一定误差。但是本实验误差在较小范围内，响应曲线也反映了预期要求，所以本实验基本得到了预期结果。 实验结果说明 由本实验结果可看出，一阶系统阶跃响应是单调上升的指数曲线，特征有T确定，T越小，过度过程进行得越快，系统的快速性越好。 二阶系统 图1.1 图1.2 图1.3
系统传递函数： 令 二阶系统模拟线路 0.25 0.50 1.00 R4 210.5 C2 111 实测45.8% 16.9% 0.6% 理论44.5% 16.3% 0% T S实测13.9860 5.4895 4.8480 T S理论14.0065 5.3066 4.8243 阶跃响应曲线图2.1 图2.2 图2.3 注：T s理论根据matlab命令[os,ts,tr]=stepspecs(time,output,output(end),5)得出，否则误差较大。 误差计算及分析 1）当ξ=0.25时，超调量的相对误差= 调节时间的相对误差= 2）当ξ=0.5时，超调量的相对误差==3.7% 调节时间的相对误差==3.4% 4）当ξ=1时，超调量的绝对误差= 调节时间的相对误差==3.46% 误差分析：由于本试验中，用的参量比较多，有R1,R2,R3,R4;C1,C2;在它们的取值的实际调节中不免出现一些误差，误差再累加，导致最终结果出现了比较大的误差，另外，此实验用的导线要多一点，干扰和导线的传到误差也给实验结果造成了一定误差。但是在观察响应曲线方面，这些误差并不影响，这些曲线仍旧体现了它们本身应具有的特点，通过比较它们完全能够了解阶跃响应及其性能指标与系统参数之间的关系，不影响预期的效果。 实验结果说明 由本实验可以看出，当ωn一定时，超调量随着ξ的增加而减小，直到ξ达到某个值时没有了超调；而调节时间随ξ的增大，先减小，直到ξ达到某个值后又增大了。 经理论计算可知，当ξ=0.707时，调节时间最短，而此时的超调量也小于5%，此时的ξ为最佳阻尼比。此实验的ξ分布在0.707两侧，体现了超调量和调节时间随ξ的变化而变化的过程，达到了预期的效果。 图2.2 图2.1 图2.3
发动机空燃比控制器 引言：我主要从事自动化相关研究。这里介绍我曾经接触过的发动机空燃比控制器设计中的优化问题。 发动机空燃比控制器设计中的最优化问题 AFR =a f m m && (1) 空燃比由方程（1）定义，在发动机运行过程中如果控制AFR 稳定在14.7可以获 得最好的动力性能和排放性能。如果假设进入气缸的空气流量a m &可以由相关单元检测得到，则可以通过控制进入气缸的燃油流量f m &来实现空燃比的精确控制。由于实际发动机的燃油喷嘴并不是直接对气缸喷燃油，而是通过进气歧管喷燃油，这么做会在进 气歧管壁上液化形成油膜，因此不仅是喷嘴喷出的未液化部分燃油会进入气缸，油膜 蒸发部分燃油也会进入气缸，如方程（2）。这样如何更好的喷射燃油成为了一个问题。 1110101122211ττττ?? ?? -?? ??????????=+????????-????????????-???? ? ??? ?? ????????? ?f f f v X x x u x x X x y =x && (2) 其中12、,==ff fv x m x m &&=f y m &,=fi u m &这里面，表示油膜蒸发量ff m &、fv m &表示为液化部分燃油、fi m &表示喷嘴喷射的燃油，在τf 、τv 、X 都已知的情况下，由现代控制理论知识，根据系统的增广状态空间模型方程（3） 0000001 1 011011114.70ττττ????-?? ??????????=-+-??????????????? ??????????????? ?? ??=?????? f f v v a X X u +q q m y q x x x &&& (3) 其中()0 14.7?t a q = y -m &。由极点配置方法，只要设计控制器方程（4），就可以 使得y 无差的跟踪阶跃输入，那么y 也能较好的跟踪AFR *a m /&。 12-- u =K q K x (4) 这里面的12、K K 确定，可由主导极点概念降维成两个参数12C ,C ，虽然都是最终稳态无差，但是目标是使得瞬态过程中y 和阶跃输入y r 的差异尽可能的小。所以原问
2012年北航三系305（即导航方向）复试试题（回忆版） 先补充说明一下，“北航三系305（即导航方向）”的完整信息如下： 学校：北京航空航天大学 院系：自动化科学与电气工程学院 专业：控制科学与工程 方向：导航、制导与控制 北航三系305（即导航方向）复试总分500分，包括两部分：200分的笔试和300分的面试。 I、笔试 笔试总分200分，其中英语部分50分，专业课部分150分，考试时间2小时。 下面的试题是根据我自己的回忆，按题型分类总结的。 一、英语部分50分，题型为英汉互译。 1、英译汉，共3题，每题10分。 （1）一篇论文的摘要，介绍一种控制方法，约300个单词。 （2）关于通信方面的科普知识，约200个单词。 （3）概念介绍，关于经典控制理论与现代控制理论的对比，约100个单词。 2、汉译英，共2题，每题10分。 （1）“时变系统”与“时不变系统”的概念介绍，约100个汉字。 （2）“计算机仿真”的概念介绍，约100个汉字。 二、专业课部分150分，题型有判断题、选择题、填空题、问答题、计算题等。 1、判断题：关于电路原理，共10道小题。 2、选择题：关于自控原理，有6-7道小题。 3、填空题：在微分方程、传递函数、状态空间方程三种形式下，二阶系统的固有频率分别由什么量给出。 从第4题开始，以下题目均为大题形式： 4、给了一个弹簧振子的图，根据条件写出微分方程。 5、给了一个系统的电路图，让你补全负反馈部分的电路，并且画出系统的框图。
6、PID控制的特点是什么？比例、积分、微分对系统分别有什么影响？ 7、保持器的作用是什么？ 8、离散系统的稳定条件是什么？ 9、连续系统的稳定条件是什么？ 10、在阶跃输入下，为了研究被控对象，应记录响应的“过渡过程曲线”还是“稳定之后的曲线”？如果是正弦输入呢？ 11、画出一阶系统的阶跃响应曲线和脉冲响应曲线，说明一阶系统的特点。 12、机器语言、汇编语言、高级语言分别有什么特点？哪种语言执行速度最快？ 13、如果系统的时间常数为T，那么计算机控制系统的采样周期应该为多少？ 14、采样频率如何确定？它与计算机硬件有怎样的关系？ II、面试 面试总分300分，其中英语口试50分，专业综合面试250分。 面试的题目因人而异，随机性很强。 下面的面试题目是根据我自己的面试经历，按顺序记录下来的。 由于面试时头脑高度紧张，所以回忆得没有笔试题目那么完整，请谅解。 一、英语口试50分，所有的问答全部是英语。 1、进行1分钟左右的自我介绍 2、从若干张纸条中抽取一张，朗 读上面的英文，并翻译，约200个单词。 3、问答： （1）你为什么要考研？ （2）你为什么选择北航读研？ （3）你认为，线性系统与非线性系统的主要区别是什么？ 二、专业综合面试250分，又包括两部分：专业知识问答和其他问答，这两部分没有具体的分值比例和数量比例，我想主要根据面试老师的喜好以及你现场的表现。 1、专业知识问答： （1）什么是最小相位系统？ （2）什么是放大器的频率响应？ （3）在饱和状态下，三极管的极电结是什么状态？ （4）什么是组合逻辑电路？什么是时序逻辑电路？ （5）什么是竞争冒险现象？ 2、其他问答： （1）你在大学期间参加过科技、竞赛等活动吗？
现代控制理论大作业分析对象：汽车悬架系统 指导老师：周晓敏 专业：机械工程 姓名：白国星 学号：S2*******
1.建模 悬架是车轮或车桥与汽车承载部分之间具有弹性的连接装置的总称，具有传递载荷、缓和冲击、衰减振动以及调节汽车行驶中的车身位置等作用。传统汽车悬驾系统是被动悬驾，其参数不能改变，无法控制其对不同路面激励的响应，因此对不同路面的适应性较差。为提高汽车的行驶平顺性、操纵稳定性和制动性等性能，人们开始用主动悬架系统来代替传统的被动悬架系统。主动悬架系统能根据路面的情况通过一个动力装置改变悬挂架的参数，改善汽车的各方面性能。 对悬驾系统进行仿真计算首先要建立悬驾系统动力学模型，随后对所建立的模型进行仿真分析。为了简化模型，取汽车的一个车轮的悬驾系统进行研究，该模型可简化为一维二自由度的弹簧阻尼质量系统，图1所示为该模型的模拟图。 图1 悬架系统模型的模拟图 其中u为动力装置的作用力，w为路面位移，x1为车身位移，x2为悬驾位移，用车身位移来度量车身的振动情况，并视为系统的输出。路面状况以w为尺度，并视为系统的一个干扰输入。当汽车从平面落入坑时，w可用一个阶跃信
号来模拟。u 为主动悬架的作用力，它是系统的控制量。 进行受力分析，由牛顿第二规律可得车身悬架系统的动力学方程为: ()()()()() 1121212212122s s t m x K x x b x x u m x K x x b x x u K w x ?=-+-+?? =-+--+-??& &&&&&&& 设系统状态变量为: []1 2 12x x x x x =&& 则上面系统动力学方程可改写为状态空间表达式: x Ax Bu y Cx Du =+?? =+?& 其中: ()1 1 1 1222 200 100001s s s t s K K b b A m m m m K K K b b m m m m ????????--=????-+??-??? ? 12 200 001 01t B m K m m ?? ??????=????-???? []1000C = []00D = u u w ??=???? Matlab 系统模型程序代码： m1=800;m2=320;ks=10000;b=30000; kt=10*ks;
惯性导航基础课程大作业报告（一）光纤陀螺误差建模与分析 班级：111514 姓名： 学号 2014年5月26日
一.系统误差原理图 二.系统误差的分析 （一）漂移引起的系统误差 1. εx ，εy ，εz 对东向速度误差δVx 的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVx1=e*g*sin(L)/(Ws^2-Wie^2)*(sin(Wie*t)-Wie*sin(Ws*t)/Ws); mcVx2=e*((Ws^2-(Wie^2)*((cos(L))^2))/(Ws^2-Wie^2)*cos(Ws*t)-(Ws^2)*((sin(L))^2)*cos(Wi e*t)/(Ws^2-Wie^2)-(cos(L))^2); mcVx3=(sin(L))*(cos(L))*R*e*((Ws^2)*cos(Wie*t)/(Ws^2-Wie^2)-(Wie^2)*cos(Ws*t)/(Ws^2-Wi e^2)-1); plot(t,[mcVx1',mcVx2',mcVx3']); title('Ex,Ey,Ez 对Vx 的影响'); xlabel('时间t'); ylabel('Vx(t)'); 0,δλδL ,v v δδ
legend('Ex-mcVx1','Ey-mcVx2','Ez-mcVx3'); grid; axis square; 分析：εx，εy，εz对东向速度误差δVx均有地球自转周期的影响，εx，εy还会有舒勒周期分量的影响，其中，εy对δVx的影响较大。 2.εx，εy，εz对东向速度误差δVy的影响 clc;clear all; t=1:0.01:25; g=9.8; L=pi/180*39; Ws=2*pi/84.4*60; Wie=2*pi/24; R=g/(Ws)^2; e=0.1*180/pi; mcVy1=e*g*(cos(Wie*t)-cos(Ws*t))/(Ws^2-Wie^2); mcVy2=g*sin(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); mcVy3=g*cos(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t)); plot(t,[mcVy1',mcVy2',mcVy3']); title('Ex,Ey,Ez对Vy的影响'); xlabel('时间t'); ylabel('Vy(t)'); legend('Ex-mcVy1','Ey-mcVy2','Ez-mcVy3'); grid; axis square;
信息工程学院现代控制理论课程习题清单
3.有电路如图1-28所示。以电压U(t)为输入量，求以电感中的电流和电 容上的电压作为状态变量的状态方程，和以电阻 R 2上的电压作为输出 量的输出方程。 4.建立图P12所示系统的状态空间表达式。 M 2 1 f(t) 5.两输入u i ，U 2，两输出y i ，y 的系统，其模拟结构图如图 1-30所示, 练习题 ，输出为，试自选状态变量并列写出其状 2. 有电路如图所示，设输入为 态空间表达式。 C ri _ l- ------- s R 2 U i U ci L u A ------ — 2 R i
试求其状态空间表达式和传递函数阵。 6.系统的结构如图所示。以图中所标记的 x 1、x 2、x 3作为状态变量，推 导其状态空间表达式。 其中，u 、y 分别为系统的输入、 输出，1、 2 试求图中所示的电网络中，以电感 L i 、L 2上的支电流x i 、X 2作为状态 变量的状态空间表达式。这里 u 是恒流源的电流值，输出 y 是R 3上的 支路电压。 8. 已知系统的微分方程 y y 4y 5y 3u ，试列写出状态空间表达式。 9. 已知系统的微分方程 2y 3y u u ， 试列写出状态空间表达式。 10. 已知系统的微分方程 y 2y 3y 5y 5u 7u ，试列写出状态空间 表达式。 7. 3均为标量。
11. 系统的动态特性由下列微分方程描述 y 5 y 7 y 3y u 3u 2u 列写其相应的状态空间表达式，并画出相应的模拟结构图。 12. 已知系统传递函数 W(s) 坐 卫 2 ,试求出系统的约旦标准型 s(s 2)(s 3) 的实现，并画出相应的模拟结构图 13. 给定下列状态空间表达式 X 1 0 1 0 X 1 0 X 2 2 3 0 X 2 1 u X 3 1 1 3 X 3 2 X 1 y 0 0 1 x 2 X 3 (1)画出其模拟结构图；(2)求系统的传递函数 14. 已知下列传递函数，试用直接分解法建立其状态空间表达式，并画出状 态变量图。 15. 列写图所示系统的状态空间表达式。 16. 求下列矩阵的特征矢量 0 1 0 A 3 0 2 12 7 6 17. 将下列状态空间表达式化成约旦标准型(并联分解) (1)g(s ) s 3 s 1 3 2 s 6s 11s 6 ⑵ g(s ) s 2 2s 3 3 c 2 s 2s 3s 1
现代控制理论 直流电动机模型的分析 姓名：李志鑫 班级：测控1003 学号：201002030309
2 1直流电动机的介绍 1.1研究的意义 直流电机是现今工业上应用最广的电机之一，直流电机具有良好的调速特性、较大的启动转矩、功率大及响应快等优点。在伺服系统中应用的直流电机称为直流伺服电机，小功率的直流伺服电机往往应用在磁盘驱动器的驱动及打印机等计算机相关的设备中，大功率的伺服电机则往往应用在工业机器人系统和CNC铣床等大型工具上。[1] 1.2直流电动机的基本结构 直流电动机具有良好的启动、制动和调速特性，可以方便地在宽范围内实现无级调速，故多采用在对电动机的调速性能要求较高的生产设备中。 直流伺服电机的电枢控制：直流伺服电机一般包含3个组成部分： - 图1.1 ①磁极： 电机的定子部分，由磁极N—S级组成，可以是永久磁铁（此类称为永磁式直流伺服电机），也可以是绕在磁极上的激励线圈构成。 ②电枢： 电机的转子部分，为表面上绕有线圈的圆形铁芯，线圈与换向片焊接在一起。 ③电刷： 电机定子的一部分，当电枢转动时，电刷交替地与换向片接触在一起。 直流电动机的启动
电动机从静止状态过渡到稳速的过程叫启动过程。电机的启动性能有以下几点要求： 1）启动时电磁转矩要大，以利于克服启动时的阻转矩。 2）启动时电枢电流要尽可能的小。 3）电动机有较小的转动惯量和在加速过程中保持足够大的电磁转矩，以利于缩短启动时间。 直流电动机调速可以有： （1）改变电枢电源电压； （2）在电枢回路中串调节电阻； （3）改变磁通，即改变励磁回路的调节电阻Rf以改变励磁电流。 本文章所介绍的直流伺服电机，其中励磁电流保持常数，而有电枢电流进行控制。这种利用电枢电流对直流伺服电机的输出速度的控制称为直流伺服电机的电枢控制。如图1.2 Bm 电枢线路图1.2 ——定义为电枢电压（伏特）。 ——定义为电枢电流（安培）。 ——定义为电枢电阻（欧姆）。 ——定义为电枢电感（亨利）。 ——定义为反电动势（伏特）。 ——定义为励磁电流（安培）。 ——定义为电机产生的转矩（牛顿?米） ——定义为电机和反射到电机轴上的负载的等效粘带摩擦系数（牛顿?米∕度?秒） —定义为电机和反射到电机轴上的负载的等效转动惯量（千克?米）。 1.3建立数学模型 电机所产生的转矩，正比于电枢电流I与气隙磁通Φ的乘积，即： Φ (1-1) 而气隙磁通Φ又正比于激励电流，故式（1-1）改写为 (1-2)
1. 控制系统任务的物理描述 为了满足飞机品质的要求，飞机的纵向运动和横侧向运动都需要有能够连续工作的阻尼器，以用来调整飞机的飞行姿态，避免其出现不必要的俯仰和倾斜。维持飞机纵向运动的阻尼器称为俯仰阻尼器，维持飞机横侧向运动的阻尼器称为偏航阻尼器。本次课程大作业旨在通过运用Matlab 的经典控制系统设计工具对某型飞机偏航阻尼器进行控制系统的设计。 2. 控制系统对象的数学模型 巡航状态下，某型飞机侧向运动的状态空间模型为： 111 12131411122212223242122131 3233343132234142434441424()1()()()()2()()()3()()4t x t a a a a b b t x t a a a a b b u t a a a a b b u t x t t a a a a b b x t t x x x x ??????????????????????????????????=+???????????????????????? ?????????? 111121314122122 2324234()()()()()()x t c c c c y t x t c c c c y t x t x t ??????????=?????????????? 式中： 1()x t ：侧滑角（单位为rad ） 2()x t ：偏航角速度（单位为/rad s ） 3()x t ：滚转角速度（单位为/rad s ） 4()x t ：倾斜角（单位为rad ） 输入向量及输出向量分别为： 1()u t ：方向舵偏角（单位为rad ） 2()u t ：副翼偏角（单位为rad ）
现代控制理论大作业-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII
分类号：TH89 单位代码：10110 学号： 中北大学 综合调研报告题目: 磁盘驱动器读写磁头的定位控制 系别: 计算机科学与控制工程学院 专业年级: 电气工程与智能控制2014级 姓名: 何雨贾晨凌朱雨薇贾凯张钊中袁航 学号: 14070541 39/03/04/16/33/47 指导教师: 靳鸿教授崔建峰讲师 2017年5月7日
摘要 硬盘驱动器作为当今信息时代不可缺少的存储设备，在人们日常生活中正扮演着越来越重要的角色，同时它也成为信息时代科学技术飞速发展的助推器。然而，随着信息量的日益增长，人们对硬盘驱动器存储容量的要求越来越高。但另一方面由于传统硬盘驱动器的低带宽、低定位精度，导致磁头很难准确地定位在目标磁道中心位置，从而限制了存储容量的持续增加。 自IBM公司于1956年向全球展示第一台磁盘存储系统R.AMAC以来，随着存储介质、磁头、电机及半导体芯片等相关技术的不断发展，硬盘的存储容量成倍增长、读写速度不断提高。要保证可靠的读写性能，盘片的转速控制和磁头的定位控制问题具有重要意义。其中磁头的定位控制主要包括寻道控制与定位跟踪控制两个问题，如PID控制、自适应控制、模态切换控制等，这些控制方法大大提高了硬盘磁头伺服系统的性能。为达到更高的精度，磁头双级驱动模型成近年的研究热点，多种控制策略已有相关报道，但目前仍处于实验水平。 关键词： 磁盘驱动器；磁头；定位；控制 Abstract Hard disk drive (HDD), acted as requisite storage equipment in current information age,plays a more and more vital role in people’s daily life, and it becomes a roll booster in rapid development of science and technology. However, with the increase of information capacity, we put forward a severe request for HDD data storage capacity. Unfortunately, due to the low bandwidth, low positioning accuracy in conventional HDD, magnetic head is hard to be positioned onto the destination track center, thus it limits the continuing increase in storage capacity. Since IBM brought the first disk-the random access memory accounting machine(RAMAC) to market in 1956, the storage capacity and read/write speed have continuously increased along with the development of the techniques of media,read/write head, actuators and semiconducting chips. The problems of R/W head's settling control is definitely important in order to ensure the reliability of read and write performance. Track seeking and track following are two main stages of the hard disk servo system. Researchers have developed kinds of control strategies to implement the servo control from PID control to advanced control methods.Dual-stage actuator has attracted many researchers and engineers for its broaderbandwidth compared with single-stage actuator. Key Words：Hard Disk Drive;Heads; Location; Control
目标：使用带双步位移的QR 分解法求矩阵10*10[]ij A a =的全部特征值，并对其中的每一个实特征值求相应的特征向量。已知：sin(0.50.2)() 1.5cos( 1.2)(){i j i j ij i j i j a +≠+== (i,j=1,2, (10) 算法： 以上是程序运作的逻辑，其中具体的函数的算法，大部分都是数值分析课本上的逻辑，在这里特别写出矩阵A 的实特征值对应的一个特征向量的求法： ()[]()() []()[]()111111I 00000 i n n n B A I gause i n Q A I u Bu u λλ-?-?-=-?-?? ?-=????→=??????→= ?? ? 选主元的消元 检查知无重特征值 由于=0i A I λ- ，因此在经过选主元的高斯消元以后，i A I λ- 即B 的最后一行必然为零，左上方变 为n-1阶单位矩阵[]()()11I n n -?-，右上方变为n-1阶向量[]()11n Q ?-，然后令n u 1=-，则 ()1,2,,1j j u Q j n ==???-。
这样即求出所有A所有实特征值对应的一个特征向量。 #include 1、什么是对偶系统，从传递函数矩阵，特征多项式和能控、能观性说明互为对偶的两个系统之间的关系。 答：定义：如果两个系统满足A2=A1T，B2=C1T,C2=B1T，则称这两个系统互为对偶函数。互为对偶系统传递函数矩阵互为转置特征多项式相同，一个函数的能控性等价于另一个函数的能观性。 2、什么是状态观测器？简述构造状态观测器的原则。 答：系统的状态不易检测，以原系统的输入和输出为输入量构造，一动态系统，使其输出渐近于原系统状态，此动态系统为原系统的状态观测器。原则：（1）观测器应以原系统的输入和输出为输入量；（2）原系统完全能观或不能观于系统是渐近稳定的；（3）观测器的输出状态应以足够快速度超近于原系统状态；（4）有尽可能低的维数，以便于物理实现。 3、说明应用李氏第二法判断非线性系统稳定性基本思想和方法步骤和局限性。 答：基本思想：从能量观点分析平衡状态的稳定性。（1）如果系统受扰后，其运动总是伴随能量的减少，当达到平衡状态时，能量达到最小值，则此平衡状态渐近稳定：（2）如果系统不断从外界吸收能量，储能越来越大，那么这个平衡状态就是不稳定的：（3）如果系统的储能既不增加也不消耗，那么这个平衡状态时李亚普诺夫意义下的稳定。方法步骤：定义一个正定的标量函数V（x）作为虚构的广义能量函数，然后根据V（x）=dV（x）/dt的符号特征来判别系统的稳定性。局限性：李雅普诺夫函数V(x)的选取需要一定的经验和技巧。 4、举例说明系统状态稳定和输出稳定的关系。 答：关系：（1）状态稳定一定输出稳定，但输出稳定不一定状态稳定；（2）系统状态完全能观且能控=状态稳定与输出稳定等价。 举例： A的特征值 =-1 =1 所以状态不是渐进稳点的，W(s)的极点S=-1，所以输出稳点。 5、什么是实现问题？什么是最小实现？说明实现存在的条件。 答：（1）由系统的运动方程或传递函数建立SS表达式的问题叫做实现问题；（2）维数最小的实现方式时最小实现；（3）存在条件是m小于等于n。 6、从反馈属性、功能和工程实现说明状态反馈和输出反馈的优缺点。 答：（1）状态反馈为全属性反馈，输出反馈为部分信息反馈；（2）状态反馈在功能上优于输出反馈；（3）从工程上讲输出反馈优于状态反馈。 7、说明李氏第一法判断稳定性的基本思想和局限性。 答：（1）基本思想：将状态方程在平衡状态附近进行小偏差线性化，由系统矩阵的特征值判断系统稳定性。（2）局限性：对非线性系统，只能得出局部稳定性；系统虚轴上有特征值时不能判断稳定性。 8、简述线性时不变系统能控性定义，并说出两种判断能控性的方法。 答：（1）定义：如果存在一个分段连续的输入U(t)，能在有限时间区间{t0,tf}内，使系统由某一初始化状态x(t0)，转移到指定的任一终端状态x(tf)，则此状态是能控的。若系统所有状态都是能控的，则完全能控，否则不完全能控。（2）方法：约旦标准型判据，秩判据。 9、说明系统传递函数零、极点对消与系统能控能观性关系。现代控制理论
相关主题
文本预览