最少拍无纹波控制器的设计

目录

前言： (1)

1 课题简介 (2)

课程设计目的 (2)

课程设计内容 (2)

2 最小拍无纹波系统控制算法设计 (3)

设计原理 (3)

算法实现 (4)

3 最小拍无纹波控制软件编程设计 (5)

运用simulink进行仿真 (5)

程序仿真 (8)

4 结果分析 (9)

5最少拍无纹波控制系统对典型输入的适应性问题 (9)

6设计总结 (10)

参考文献 (11)

最少拍无纹波控制器的设计3

摘要：本实验介绍了对一阶惯性环节控制对象,采用最少拍无纹波控制算法设计的一种数字控制器，《计算机控制技术》是一门理论性、实用性和实践性都很强的课程，课程设计环节应占有更加重要的地位。计算机控制技术的课程设计是一个综合运用知识的过程，它需要控制理论、程序设计、硬件电路设计等方面的知识融合。通过课程设计，加深对学生控制算法设计的认识，学会控制算法的实际应用，使学生从整体上了解计算机控制系统的实际组成，掌握计算机控制系统的整体设计方法和设计步骤，编程调试，为从事计算机控制系统的理论设计和系统的调试工作打下基础。

关键词：最少拍无纹波控制；控制；计算机控制；

；

前言：本实验通过对最少拍无纹波控制器的设计及仿真了解和掌握了最少拍无纹波设计及有纹波设计。

首先，通过学习和搜集相关书籍资料了解和掌握了最少拍控制器的设计原理，从而分别根据单位阶跃信号输入和单位速度信号输入情况，设计了不同的最少拍无纹波控制器，并采用Simulink进行了仿真，同时又通过matlab程序验证了仿真结果的正确性。

其次，我们以单位速度信号输入为例，比较了有纹波和无纹波控制器的区别，最终得出结论：最少拍无纹波调整时间较长，但精度较高；

最后，我们通过选择不同的输入信号对同一个最少拍无纹波控制器进行仿真，研究了最少拍无纹波控制系统对典型输入的适应性问题，最终发现根据某一种输入信号情况设计的无纹波控制器可适用于较低阶的输入信号情况，但不适用于更高阶的输入信号情况。

1 课题简介

课程设计目的

1. 采用零阶保持器离散化，采样周期取秒；

2. 利用matlab 软件中的c2d()函数、solve()等函数求解)(z D 的参数；

3. 绘制出单位加速度输入信号闭环系统的仿真曲线；

4. 分析单位阶跃输入信号和单位速度输入信号时闭环系统动态性能和稳态性能，并绘制出系统的响应曲线；

5. 课程设计报告中要给出系统控制算法的程序流程图； 完成课程设计文档、按照课程设计要求撰写课程设计报告。 课程设计内容

下面以一个具体实例介绍最少拍系统的设计和仿真。 如图1所示的采样-数字控制系统，

图1 离散控制系统结构图

其中对象8)

+(s 4)+(s 1)+(s s 2)+(s 68)(2

s G

选择采样周期T=，试设计无纹波最少拍控制器，并分析仿真结果 1. 分别在单位阶跃/单位速度输入下设计无纹波有限拍控制器

2. 在Simulink 仿真环境画出仿真框图及得出仿真结果，画出数字控制器和系统输出波形。

3. 与有纹波系统进行对比分析（选用单位速度输入进行对比分析即可）

4. 探讨最少拍无纹波控制系统对典型输入的适应性问题

2 最小拍无纹波系统控制算法设计

设计原理

最小拍控制的闭环脉冲传递函数要求有如下的形式：

N N z z z z ---+++=Φφφφ............)(2211

这一形式表明经历有限个采样周期后输出能跟上输入的变化，系统在采样点没有静差。

根据z 变换的终值定理和系统的稳态误差的情况，要求系统的

)()1()(1)(1z F z z z q e --=Φ-=Φ

即有)()1(1)(1)(1z F z z z q

e ---=Φ-=Φ

这里F(z)是关于1

-z 的待定系数多项式。显然，为了使)(z Φ能够实现，F(z)首项应为1，即

p p z z f z f z F ---++++=φ............1)(2211

因此最少拍控制器D(z)为q

q

z z z G z z z z G D )

1)(()1(1)(1)()(111-----=Φ-Φ=

图2 控制原理图

最小拍无纹波控制系统要求在非采样时间的时候也没有偏差，因此必须满足：

①对阶跃输入，当t ≥NT 时，有y(t)=常数。 ②对速度输入，当t ≥NT 时，有y'(t)=常数。 ③对加速度输入，当t ≥NT 时，有y''(t)=常数。

因此，设计最小拍无纹波控制器时，对速度输入函数进行设计，则Gc(s)必须至少有一个积分环节，使得控制信号u(k)为常值时，Gc(s)的稳态输出是所要求的速度函数。同理，若针对加速度输入函数设计的无纹波控制器，则Gc(s)中必须至少有两个积分环节。

最小拍控制的广义对象含有D 个采样周期的纯滞后d

z z A z B z G -=

)

()()(所以)()()()()

()

()(z R z z R z A z B z z z U d μ

Φ=Φ=

-其中)()()()(z A z B z z z d -Φ=Φμ。 要使控制信号u(k)在稳态过程中为常数或0，那么)

(z μΦ只能是关于1-z 的有限多项式。

因此

)

()]1([)()()(211

2z F z b z z F z B z z i i d

d

-=--∏-==Φω

w 为G(z)所有零点数(包括单位圆内、单位圆上以及单位圆外的零点)。w

b b b ,...,21为其

所有零点。 算法实现

2.2.1单位阶跃输入

（1）带零阶保持器的广义被控对象为G(s) 通过matlab ，z 变换程序为

num=[0 0 0 0 68 136] den=[1 13 44 32 0 0] Gs=tf(num,den); Gz=c2d(Gs,,'zoh'); zsys=zpk(Gz)

num = 0 0 0 0 68 136 den = 1 13 44 32 0 0 Zero/pole/gain: (z+ (z+ (z+

------------------------------------------------------ (z-1)^2

Sampling time: 即

(z+ (z+ (z+

G(z)=------------------------------------------------------

(z-1)^2