G11G12G GG2221K是双速率离散控制器应该被适当调整去满足相应的因果性、周期性和有限维性。 对于任意的T>0,Dr为连续空间上的延迟算子,Druc (t) = uc (tT);U为离散空间上的一步滞后算子;U2为离散空间上的一步超前算子。
定义1 如果(U2)q1KUp1=K成立,则称K为(p,q)-周期离散控制器。 定义2 如果连续系统G满足DrG=GDr,则称G为T-周期连续时变系统。 3.仿真算法 3.1仿真表达式
K是一个已知的定义(p,q)-周期的离散控制器,采样算子和保持算子如上所述,则HKS以T为周期的时变系统。如图1即可证明信号之间的关系,在已知既定的条件下下式成立:
第 2 页 共 9 页
DTHKSyc(t)DTuc(t)uc(tT)uc(tq1qh)uc[(kq1)qhr]Hu(kq1)HUq1u(k)HUq1Ky(k)HKUp1y(k)HKy(kp1)HKSyc(tp1ph)HKSyc(tT)HKSDTyc(t)DTHKSyc(t)HKSDTyc(t)DTHKSHKSDT0rqh 我们可以由定义2看到,HKS为T周期的时变系统。由于HKS的周期是T,因此同单速率系统类似,图1中输出与输入的关系可以表示为:
z[G11G12(IHKSG22)1HKSG21]w
或者是
z[G11G12H(IKSG22H)1HKSG21]w
在双速率采样控制系统输出与输入通道中,通过增加一个可见的采样器且保持快速,像在图1中显示的一样,这个可见快速采样器及保持器的采样周期均为T/n。 Wd是w以T/n为采样周期的采样信号,当输入信号的仿真时间为mph时,有: Wd=w(kT/n),k=0,1,…,mn/p1
zd与z的关系同上。显然,当n→∞时,wd=w,zd=z。为使离散时间序列的个数为正整数,n选为l的整数倍。研究图1所示系统的仿真,便可得到虚框中双速率采样控制系统连续输入输出信号的仿真结果。
图1中的zd=Snz,w=Hnwd,故由式(2)得
zdSnzSn[G11G12H(IKSG22H)1KSG21]Hnwd[G11nG12nSnH(IKSG22H)1KSHnG21n]wd其中G11n,G12n,G21n为对应于周期T/n的离散化。式(4)即为双速率采样控制系统的仿真表达式。
3.2 仿真表达式的计算
第 3 页 共 9 页
欲求表达式(4),首先要得到G11n, G12n,,G21n,,SnH,,SHn,以及(I-KSG22H)-1K等等变量 ,G11n, G12n, G21n 分别是连续传递函数G11, G12, G21以T为采样周期采样后的离散传递函数,均以计算。下面讨论SnH,SHn,(I-KSG22H)-1K的计算。
(1)计算SnH
图 2 SHn的输入与输出框图
图2中Hn的周期为T/n=lh/n,S的周期为ph,当x2 (0) = x1 (0), x2 (1) = x1 (n/p1), ..., x2 (m -1) = x1 ((m-1) n/p1), SHn =
m/pm/pm/p100000001000 000010mn/p(2)计算SnH
同理可得到SnH的表达式:
m10001000SnH2n/q 00010001
(3)计算(I-KSG22H)-1K
由离散采样以及离散算子的定义,有:
Ψ (k) = Φ (kp)Sp2l → l, Ψ = SpΦ
第 4 页 共 9 页
υ (kq + r) = Φ (k)Hq2l → l υ = HqΦ
R = 0,1, ..., q-1
可得:
SG22H = SrSAG22HhHq = SpG22dHq (5)
其中G22d可通过单一速率h离散化过程得到。对于 (I-KSpG22dHq)-1K而言,仍然是变量SpG22dHq和K的周期,本文应用离散提升算子将其变成为是不变系统,具体的过程如图3所示。此时仿真表达式(4)可表示为图4。
经离散提升后,周期时变环节 SpG22dHq 和K变成了时不变的Lp1SpG22dHqL-1 q1 和 Lq1K L-1 q1,具体的计算过程如下: (1)Lq1K L-1 q1 的计算
如果双速率控制器K的状态方程是
xk(k1)Axk(k)Bjy(p1kj)j0p11u(q1ki)Cixk(k)Di,jy(p1kj)
j0p11i0,1,,q11 与此同时,Lq1K L-1 q1的状态方程可以被标识为:
~~xk(k1)Axk(k)B~y(k) ~~~~u(k)CxkDy(k)k其中
图 3 用于提升离散信号的(I-KSG22H)-1K
第 5 页 共 9 页
图4 仿真示意
(2)Lp1SpG22dHqL-1 q1的计算
引理1:设P的状态变量为x,状态模型参数矩阵为[A、B、C、D],m、n和s 是满足如下关系式的正实数。对离散采样算子的系统状态变量可被表示成为一个状态模型。即:
~AA
其中
(j1)m1Di,jD[jm,(j1)m](im)特征函数X为:
rjmCAim1rB[0,im](r)
第 6 页 共 9 页
有上述结论,可由G22d求得Lp1SpG22dHqL-1 q1 状态空间模型矩阵。
综上所述,在图4中我们可以看到整个的仿真过程为:mph输入信号周期,然后:
4.仿真举例
图1中广义被控对象G为:
第 7 页 共 9 页
控制器K为:
采样周期: T1 = 2s, T2 = 3s, p = 2, q = 3, h = 1, p1 = 3, q1 = 2, l = 6, T = 6。令m=6,n依次令其等于4800,7200,9600,wd是单位阶跃输入信号。使用MATLAB 编程语言,并且进行系统仿真, 结果如图五所示:
5.结论
本文针对双速率采样控制系统的特点,应用离散提升技术、快速采样算子和快速
保持算子,研究双速率采样控制系统的仿真方法,并给出了具体的仿真步骤和方针实例。由于双速率控制器为时变控制器,所以有关双速率采样控制系统仿真精度的验证问题还
第 8 页 共 9 页
有待于进一步研究。
采样控制系统技术已经历十多年的发展,却存在着根本性的问题。尤其是自从采
用了提升技术,采样控制理论进入了一个新的发展阶段。由于能够计及采样时刻之间的性能,所以提升变换似乎已经成了采样控制系统分析和设计的唯一正确的方法,其应用也在逐步扩大,但是在现实设计中的应用却对其提出了更高的要求。对其提升技术本来是为了相关设计的需要而提出的,但很多现实情况不仅仅局限于个别领域。这就是采样控制系统的特殊性,尤其是在于其信号通道的结构上。采样控制系统的信号通道由两部分所构成,一个是连续通道,另一个是采样通道。采样控制系统提升后,其范数也不是完全等价的。考虑到这两个通道特点而提出的频率响应法也可以给出系统真实的频率响应诱导范数,将是采样控制系统分析和设计的正确方法。
第 9 页 共 9 页