文档介绍：该【自动控制设计(自动控制原理课程设计) 】是由【简爱】上传分享，文档一共【12】页，该文档可以免费在线阅读，需要了解更多关于【自动控制设计(自动控制原理课程设计) 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。工程施工文档
工程施工文档
工程施工文档
自动控制原理课程设计
本课程设计的目的着重于自动控制基本原理与设计方法的综合实际应用。主要内容包括:古典自动控制理论(PID)设计、现代控制理论状态观测器的设计、,掌握自动控制理论工程设计的基本方法和工具。
内容
某生产过程设备如图1所示,由液容为C1和C2的两个液箱组成,图中Q为稳态液体流量,为液箱A输入水流量对稳态值的微小变化,为液箱A到液箱B流量对稳态值的微小变化,为液箱B输出水流量对稳态值的微小变化,为液箱A的液位稳态值,为液箱A液面高度对其稳态值的微小变化,为液箱B的液位稳态值,为液箱B液面高度对其稳态值的微小变化,分别为A,B两液槽的出水管液阻。设为调节阀开度.
已知液箱A液位不可直接测量但可观,液箱B液位可直接测量。
图1某生产过程示意图
要求
工程施工文档
工程施工文档
工程施工文档
建立上述系统的数学模型;
对模型特性进行分析,时域指标计算,绘出bode,乃示图,阶跃反应曲线
对B容器的液位分别设计:P,PI,PD,PID控制器进行控制;
对原系统进行极点配置,将极点配置在-1+j和-1-j;(极点可以不一样)
设计一观测器,对液箱A的液位进行观测(此处可以不带极点配置);
如果要实现液位h2的控制,可采用什么方法,怎么更加有效?试之。
用MATLAB对上述设计分别进行仿真。
(提示:流量Q=液位h/液阻R,液箱的液容为液箱的横断面积,液阻R=液面差变化/流量变化.)
2 双容液位对象的数学模型的建立及MATLAB仿真过程一、对系统数学建模
如图一所示,被控参数的动态方程可由下面几个关系式导出:
液箱A:
液箱B:
消去中间变量,可得:
式中,——两液槽的容量系数
-—两液槽的出水端阻力
——第一个容积的时间常数
-—第二个容积的时间常数
_双容对象的放大系数
其传递函数为:
工程施工文档
工程施工文档
工程施工文档
,绘出bode,奈氏图,阶跃反应曲线
当输入为阶跃响应时的Matlab仿真:
令T1=T2=6;K=1
单位阶跃响应的MATLAB程序:
num1=[1];den1=[36121];
G1=tf(num1,den1);
figure(1);
step(G1);
xlabel('时间(sec)’);ylabel(’输出响应');title(’二阶系统单位阶跃响应');
step(G1,100);
运行结果如下:
阶跃反应曲线:
图1
c(∞)=1; c(tp)=1; tp=; td=10s; ts=45。5s;
最大超调量:δ(tp)=[c(tp)-c(∞)]/c(∞)*100%=0%
稳态误差分析:
开环传递函数,稳态误差;
用MATLAB绘制的奈氏图如下图2所示,其程序如下:
nyquist([1],conv([61],[61]))
工程施工文档
工程施工文档
工程施工文档
图2
在工程实践中,一般希望正相角裕度r为45~60,增益裕度KdB,即K。当系统为单位负反馈时的Bode图:用MATLAB绘制的奈氏图如下图3所示,其程序如下:
sys=tf([1],conv([61],[61]));margin(sys);figure
图3
三:对B容器的液位分别设计:P,PI,PD,PID控制器进行控制
PID控制的原理和特点
工程施工文档
工程施工文档
工程施工文档
(1)P控制:取P=9; I=0; D=0;
(2)PI控制:P=6,I=,D=0;
(3)PD控制:P=9,I=0,D=5;
工程施工文档
工程施工文档
工程施工文档
(4)PID控制:P=5,I=,D=4;
—1+j;-1-j
工程施工文档
工程施工文档
工程施工文档
根据传递函数
得微分方程
令
得状态方程
输出:
极点配置:
令K=1;T1=T2=2;
用MATLAB确定状态反馈矩阵K,使得系统闭环极点配置在(-1+j,-1-j),程序如下:A=[01;--1];
B=[0;1];
P=[—1+j;—1—j];
K=place(A,B,P)
运行结果为
K=
仿真:
工程施工文档
工程施工文档
工程施工文档
仿真图
,对液箱A的液位进行观测
建立状态观测器:
根据传递函数
得微分方程
令
得状态方程
输出:
全维观测器的建立:
令,得
工程施工文档
工程施工文档
工程施工文档
期望特征式:
对比1式和2式,得
得
所以全维状态观测器得方程是
工程施工文档
工程施工文档
工程施工文档
本实验中,需观测的状态为水箱A溶液的液位,
建立数学模型
R1=R2=1;c1=c2=1;
令状态观测器的极点为(-6—j,-6+j)
设计此给定系统状态观测器的MATLAB程序如下