自动控制大作业

自动控制原理课程大作业

班级：021111

成员：02111005 付裕深 02111006 高慧聪 02111004 崔贵福

2014 年《自动控制技术》课程大作业

一、课程习题

1．带飞球式调节器的瓦特蒸汽机是近代工业革命兴起的标志，同样也是一个典型的反馈控制系统，请画出该系统的组成框图，并注明下列器件的位置，并说明与每个信号相关的装置。 受控过程

过程要求的输出信号 传感器 执行机构

执行机构的输出信号 调节器

调节器输出信号 参考信号 误差信号

图1.1瓦特离心式调速器示意图

解：

受控过程 蒸汽机运行 过程要求的输出信号 转速 传感器 履带 执行机构 蒸汽阀

执行机构的输出信号 蒸汽推力 调节器 调速器

调节器输出信号 调速器转速 参考信号 规定转速 误差信号 转速偏差

2.倒立摆控制系统是一个典型的非线性、强耦合、多变量和不稳定系统。倒立摆的控制问题就是给连接摆杆的小车施加控制力，使摆杆尽快的达到一个平衡位置，并且保证不出现过大的振荡。当摆杆到达期望的平衡位置后，系统能克服随机扰动保持在平衡点。

如图2 所示是一个简单的一阶倒立摆系统，这里忽略空气阻力和各种次要的摩擦力，将倒立摆系统看做是一个由小车和均匀刚性杆组成的系统。假设系统初始状态时，摆杆垂直于小车处于平衡状态，此时摆杆受到冲激信号作用产生一个微小的偏移。建立此时系统的控制系统数学模型。

图1.2 小车-单摆系统示意图

解：

1. 一阶倒立摆的微分方程模型

对一阶倒立摆系统中的小车和摆杆进行受力分析，其中，N和 P为小车与摆杆相互作用力的水平和垂直方向的分量。

图 1-2 小车及摆杆受力图

分析小车水平方向所受的合力，可以得到以下方程：

（1-1）

由摆杆水平方向的受力进行分析可以得到下面等式：

（1-2）

即：

把这个等式代入式(1-1)中，就得到系统的第一个运动方程：

（1-4）

为了推出系统的第二个运动方程，我们对摆杆垂直方向上的合力进行分析，可以得到下面方程：

（1-3）

（1-5）

即：

（1-6）

力矩平衡方程如下：

（1-7）

由于

面有负号。

合并这两个方程，约去 P和 N，得到第二个运动方程：

（1-8）

设

，（φ是摆杆与垂直向上方向之间的夹角），假设φ <<1弧

。用u代表所以等式前

度， 则可以进行近似处理：

被控对象的输入力F，利用上述近似进行线性化得直线一阶倒立摆的微分方程为：

（1-9）

2. 一阶倒立摆的传递函数模型

对式(1-9)进行拉普拉斯变换，得：

（2-1）

注意：推导传递函数时假设初始条件为 0。

由于输出为角度φ，求解方程组的第一个方程，可得：

（2-2）

或

（2-3）

如果令

，则有：

（2-4）

把上式代入方程组（2-1）的第二个方程，得：

（2-5）

整理后得到传递函数：

（2-6）

其中

。

3. 一阶倒立摆的状态空间模型

设系统状态空间方程为：

（3-1）

方程组（2-9）对

解代数方程，得到解如下：

（3-1）

整理后得到系统状态空间方程：

（3-2）

（3-3）

摆杆的惯量为

，代入（1-9）的第一个方程为：

得：

化简得：

设

，

则有：

（3-4）

（3-5）

3.设计用来保持飞机俯仰姿态角的自动驾驶仪的方框图如图3所示

图1.3飞机俯仰控制系统示意图

升降角和俯仰姿态角之间的传递函数为：

（s）50( +1)( +2)

= ( )=

自动驾驶仪控制者根据以下传递函数用俯仰姿态角误差来调整飞机升降：

（s） ( +3)

=D(s)= 计算输入参考为单位阶跃变化时K的值，使得系统超调量小于10%，上升时间

小于0.5s. 解：

系统根轨迹为：

编程计算： clc; clear; close all; mn=0; mx=1; ka=0.01; kf=10; dk=0.01; K=ka:dk:kf;

b=50*conv([1,3],conv([1,1],[1,2]));

a=conv([1,10],conv([1,5,40],[1,0.03,0.06])); sys=tf(b,a); sys1=K*sys; j=1;k=1;

for i=1:length(K)

sys2(i)=feedback(sys1(i),1,-1); [Y,T]=step(sys2(i));

post1=(find(Y>mx*Y(end))); tr(i)=T(2)*(post1(1));

or(i)=(max(Y)-Y(end))/Y(end);

if tr(i)<0.5&or(i)<0.1 k1(j)=i;j=j+1; end end if j>1

mid=ceil(length(k1)/2); subplot(2,2,1);

pzmap(sys1(k1(mid))); title('开环零极点分布图'); subplot(2,2,2); rlocus(sys);

title('根轨迹曲线') grid on; subplot(2,2,3);

pzmap(sys2(k1(mid))); title('闭环零极点分布图'); subplot(2,2,4);

step(sys2(k1(mid)));

title(['阶跃响应\delta=',num2str(100*or(k1(mid))),'% t_r=',num2str(tr(k1(mid)))]); else disp('不存在满足指标的K值'); end

K=2.6；

4.我们希望设计一个汽车的速度自动控制系统,假定(a)汽车的质量为

1500kg;(b)加速度计提供控制信号U,即当其角速度变化1度时,便在汽车上施 加10N的力;(c)空气的摩擦阻力正比于速度,其比例系数为10N s/m。 1） 试求出从输入U到汽车转速间的传递函数 2） 假设转速的变化由下式给出：

10.05

( )+ ( )

其中V的单位为米/秒，U的单位为度，W为公路等级。设计一个比例控制器U =-kpV，使转速误差小于1m/s，此时路面的等级为常值的2%。 3） 试讨论当对系统施加积分控制时，会带来什么好处。 4）若该系统在纯积分控制作用下，适当选择反馈增益使系统处于临界阻尼状态。 解： （1）

( )

V(s)=

（2）

程序代码：

V(s)=

clc; clear; close all; k=0.1:5; j=1;

for i=1:length(k)

b=0.05*0.02*[1,0.002];

a=conv([1,k(i)+0.002],[1,0.02]); sys(i)=tf(b,a); [y,t] = impulse(sys); if max(y)<1

post(j)=i;j=j+1; end end

if length(post)>0

impulse(sys(post(1))); else

disp('无满足要求的K'); end

运行结果：

K=0.3 （3）

对系统施加积分控制会提升系统类型号减小稳态误差。 （4）

特征方程：573 2+10 + =0 临界阻尼方程有重根 =0.0436

5.假设一种位置伺服控制系统中受控对象的传递函数为：

10

G(s)=在单位反馈结构中设计传递函数为 D(s)的串联补偿，使以下闭环性能指标 得到满足：

参考阶跃输入响应的超调量不大于16%； 参考阶跃输入响应的上升时间不超过0.4秒。 单位斜坡输入的稳态误差小于0.02

1）设计超前补偿环节使系统满足动态响应指标。

2）如果 D(s)为比例控制器，其值为kp，速度常数Kv为多少？

3）设计一个滞后补 …… 此处隐藏：3215字，全部文档内容请下载后查看。喜欢就下载吧 ……