直流电机PI控制器设计与稳态性能分析(1)

  • 格式:docx
  • 大小:1.11 MB
  • 文档页数:16

目 录 1 计算直流电机的系统微分方程…………………………………………………… 1 2 计算kP和kI参数值 ………………………………………………………………… 2 3 计算在不同输入情况下,系统的稳态误差 …………………………………… 3 3.1单位阶跃参考输入时系统的稳态误差 …………………………………………… 3 3.2单位斜坡参考输入时系统的稳态误差 …………………………………………… 4 3.3单位阶跃扰动输入时系统的稳态误差 …………………………………………… 4 3.4单位斜坡扰动输入时系统的稳态误差 …………………………………………… 5 4 用MATLAB验证各输入条件下的稳态误差 …………………………………… 6

4.1验证单位阶跃输入时系统的稳态误差 …………………………………………… 6 4.2验证单位斜坡输入时系统的稳态误差 …………………………………………… 7 4.3验证单位阶跃扰动输入时系统的稳态误差 ……………………………………… 9 4.4验证单位阶跃扰动输入时系统的稳态误差 ……………………………………… 10

5 课程设计小结………………………………………………………………………… 12 6 心得体会……………………………………………………………………………… 12 参考文献…………………………………………………………………………………… 13 1

1、 写出以aV和W为输入的直流电机系统微分方程 图1 如图1所示,R为系统给定输入,W为系统扰动输入,由题意可知:

化简得:

所以所求的系统微分方程式为:

Gc(s) 300

031s

R Y av e +

- +

1200

W -

1(300*1200)30aVWYs

(3001200)(30)aVWYs3030012000adyyvwdt 2 230300(30300)4*3002ppIkkks

2、 已知特征方程的根包括3030j,计算kP和kI的值。 由题目已知特征方程的部分根,可以先求出系统的闭环传递函数,写出特征方程,再将特征方程根带入方程求得方程系数。具体过程如下:

① 先求出CGS的表达式: 由)(0tIpaedtkekv得

② 由()CGs算得系统的开环传递函数为:

③ 再由开环传递函数写出闭环特征方程:

④ 求得特征方程根为:

⑤ 将此根与3030j比较得:

**/pIaIcpkEskEssVskGskEsEss



3001()()30030(30)pIIpkskkGskssss



2()(30300)3000pIDssksk

0.16pIkk 3

3、 计算在不同输入情况下,系统的稳态误差 系统稳态误差的定义:当系统的过渡过程结束以后,就进入了稳态,而系统的实际输出与期望输出的偏差量称为稳态误差。稳态误差描述了控制系统的控制精度。 稳态误差产生的原因:1.组成系统的元件不完善,例如静摩擦、间隙、不灵敏区以及放大器的零点、老化或变质等。这方面引起的误差通常称为静差,消除静差可以通过优化元件来解决;2.系统结构造成的。消除这个误差的方法只能是改变系统结构。 控制系统还经常处于各种扰动作用之下,给定输入作用产生的误差称为系统给定误差,而扰动作用产生的误差称为系统扰动误差。 系统在参考输入和扰动输入作用下的误差信号的拉氏变换为: 21212

()()1()()()1()()()1()()()GsHsEsRsNsGsGsHsGsGsHs



定义 (3.1)

为扰动误差传递函数。

3.1 单位阶跃参考输入时系统的稳态误差 当输入信号1(s)=,()0sRNs时

则:

其中: 故单位阶跃输入时系统的稳态误差为: 101sspeK

111()()1()1()EsRsGsGss



0001()lim()lim1()111lim()1ssssPsesEsGsGsK





00300lim()lim(30)pIpsskskKGsss



212

()()()1()()()eNGsHssGsGsHs

 4

3.2 单位斜坡参考输入时系统的稳态误差 当输入信号21(s)=,()0sRNs时

其中: 故单位斜坡输入时系统的稳态误差为:

3.3单位阶跃扰动输入时系统的稳态误差 由题目可知,在式3.1中,()()NsWs,当()0Rs时,系统误差信号:

222()()()()()1300()()130()3001(30)()601800eWC

pI

EssWsGsWsGsGssWsksksssWsss











因为1200()Wss,故误差信号:

21200()601800Esss



2200000011()1()11()lim()lim1()11lim()limlim()11lim()sssssssvsEsGssesEssGssssGsssGssGsK











1()1600.017ssveK



00300()lim()lim(30)60pIvsskskKsEssss

 5

由终值定理求单位阶跃扰动输入时系统的稳态误差为: 2001200()lim()lim6018000sssssesEsss



3.4单位斜坡扰动输入时系统的稳态误差 同上一节,系统的误差信号为:

由于21200()Wss,故误差信号: 21200()(601800)Essss



由终值定理得: 020()lim()1200lim60180023ssssessEsss



2()()601800sEsWsss 6

4、 用MATLAB验证各输入条件下的稳态误差 4.1 验证单位阶跃输入时系统的稳态误差 由第三节的计算可以得到系统的闭环传递函数如下所示:

2()301800()1()601800GsssGsss



在MATLAB中输入如下程序,绘出系统的单位阶跃响应: 程序代码: num=[30,1800]; %系统传递函数分子多项式 den=[1,60,1800]; %系统传递函数分母多项式 step(num,den) %系统单位阶跃响应 grid on %绘制网格 绘出的单位阶跃输入响应图像如图2

图2 7

程序代码原理:step为MATLAB自带的单位阶跃响应函数,num与den代码为传递函数表达式。grid代码则在图像中绘制网格。 从图中稳态输出取参考点,由图中该点数据可以看出,系统的最终稳态输出为0.997,即C(s)=0,997,故可得到系统的误差信号E(s)()C(s)0.003Rs,因此系统稳态误差

与理论计算值一致。 4.2 验证单位斜坡输入时系统的稳态误差 系统的闭环传递函数为:

单位斜坡输入信号为21()Rss,编写MATLAB程序代码如下: num=[30,1800]; %系统传递函数分子多项式

den=[1,60,1800]; %系统传递函数分母多项式 sys=tf(num,den); %定义sys函数 R=tf(1,[1,0,0]); %定义R C=sys*R; impulse(C) %绘制C的波形 hold on % impulse(R) %绘制R的波形 grid on

2()301800()1()601800GsssGsss



0()lim()0sssesEs 8

程序代码原理:R=tf(1,[1,0,0])意思是21Rs,impulse代码的含义是单位冲激响应。num与den代码为传递函数表达式。grid代码则在图像中绘制网格。代码中绘制的C就是输出信号。 绘出的单位斜坡输入响应与输入信号的图像如图3所示:

图3 稳态输出局部放大图如图4所示:

图4