H-inf H无穷控制.ppt

H 设计的m函数 目录 广义对象的求取连续系统H 设计的m函数离散系统H 设计的m函数 标准的H 问题的框图如图 所示 图中G为系统的广义对象 K为控制器 图1 求解图 中对象G的方法有两种 1 m函数调用法2 直接求取法 一 广义对象的求取 图2加权灵敏度问题 下面通过图2所示的加权灵敏度问题的例子来看一下如何通过m函数调用来求取系统的广义对象G 1 m函数调用法 系统除去控制器K以外的部分就是广义对象G 它是两入两出的 输入信号是w和u 输出信号是z和y 可用传递函数表示为 设图2中的对象P和灵敏度权函数W分别为 将参数代入 可以得到广义对象G为 G送进去以后 调用下面的三个m函数 就可以得到广义对象G的状态空间实现 A B C D ssdata sys sys minreal ss G 图2加权灵敏度问题 G通过下面的函数送进去 G tf 0100 11 tf 1002000 1224120 1 tf 120 12120 G ltisys A B C D 这个G就是我们求解问题时所用的G 它是这样送进去的 用上面的函数调用法来求取G的状态实现 是非常简单的 但是从上面的结果可以看出 用这种方法得到的状态变量纯粹是数值上的运算 脱离了物理概念 本例中得到的广义对象G 图2加权灵敏度问题 根据结果只能知道这个广义对象的输入输出之间的关系 这几个状态变量之间的关系与实际的物理系统之间的状态没有直接联系 没有物理意义 下面我们仍用上面的例子 用直接建立状态变量的方法来求取广义对象G的状态空间实现 A B C D 首先来求对象P的状态空间实现 设被控对象P的状态变量为x1和x2 根据P的传递函数 可以得到如下的状态方程 2 直接求取法 设权函数W的状态变量为x3 根据W的传递函数 可以得到权函数W的状态空间实现 根据图2中各信号的关系 进一步可以得到广义对象G的状态空间实现为 图2加权灵敏度问题 前面讲的这部分内容是关于广义对象G如何送进去 这里我们讲了两种方法 1 m函数调用法 2 直接求取法 接下来要讲的是第二部分的内容 连续系统H 设计的m函数 1 函数hinfsyn 该函数用来计算系统的H 控制器k 函数的调用形式为 k g gfin hinfsyn G nmeas ncon gmin gmax tol 该函数用的是 DGKF文献 中的算法 1 Doyle J C K Glover P Khargonekar andB Francis State spacesolutionstostandardH2andH controlproblems IEEETransactionsonAutomaticControl vol 34 no 8 pp 831 847 August1989 2 Glover K andJ C Doyle State spaceformulaeforallstabilizingcontrollersthatsatisfyanH normboundandrelationstorisksensitivity SystemsandControlLetters vol 11 pp 167 172 1988 二 连续系统H 设计的m函数 该函数用来计算系统的H 控制器k 函数的调用形式为 k g gfin hinfsyn G nmeas ncon gmin gmax tol 其中输入变量中的G为如下定义的两入两出的广义对象 也是我们第一部分内容里所讲的用G ltisys A B C D 送进去的G 1 函数hinfsyn G 系统的广义对象 nmeas 连接到控制器的测量输出的个数 ncon 控制输入的个数 gmin 的下界 gmax 的上界 tol 的迭代精度 k H 最优控制器 g 闭环控制系统 gfin 最终的值 k g gfin hinfsyn G nmeas ncon gmin gmax tol 算例 PS T混合灵敏度问题 本例的H 问题是要求解如下的有约束的优化设计问题 K 图3PS T问题 图3中参数如下 图中除去K以外的部分就是广义对象G 按照我们第一部分内容所讲的方法把参数送进去以后 得到系统广义对象G的状态空间实现矩阵如下 由于调用函数hinfsyn时对象要满足假设中秩的要求 设计中取Dp 10 6 以后第4章讲DGKF法时还要提到 广义对象G由下面的函数送进去 G ltisys A B C D 本例中函数的调用形式如下 hinfsyn G nmeas ncon gmin gmax tol k g gfin hinfsyn G 1 1 0 1 2 0 0001 函数调用中的迭代过程如下 gammahamx eigxinf eighamy eigyinf eignrho xyp f2 0006 8e 0001 2e 0161 0e 003 1 9e 0170 0000p1 0506 7e 0001 2e 0161 0e 003 1 9e 0170 0000p0 5756 6e 0001 2e 0161 0e 003 1 9e 0170 0000p 设计中权函数W1中的是可变的 要取尽可能的最大值 这里给出的是当取1000时的迭代过程 0 3386 2e 0001 2e 0161 0e 003 1 9e 0170 0000p0 2195 3e 000 3 3e 003 1 0e 003 1 9e 0170 0000f0 2785 9e 000 2 9e 002 1 0e 003 1 9e 0170 0000f0 3086 1e 000 4 4e 001 1 0e 003 1 9e 0170 0000f0 3236 1e 0001 2e 0161 0e 003 1 9e 0170 0000p0 3156 1e 0001 2e 0161 0e 003 1 9e 0170 0000p0 3126 1e 0001 2e 0161 0e 003 1 9e 0170 0000p0 3106 1e 000 1 3e 000 1 0e 003 1 9e 0170 0000f0 3116 1e 000 1 6e 002 1 0e 003 1 9e 0170 0000f0 3116 1e 0001 2e 0161 0e 003 1 9e 0170 0000p0 3116 1e 0001 2e 0161 0e 003 1 9e 0170 0000p0 3116 1e 0001 2e 0161 0e 003 1 9e 0170 0000p0 3116 1e 0001 3e 0161 0e 003 1 9e 0170 0000pGammavalueachieved 0 3107 1 逐渐增大 当增大到100000时 这就是最终的设计结果 函数调用中的迭代过程如下 gammahamx eigxinf eighamy eigyinf eignrho xyp f2 0002 1e 0014 6e 0131 0e 0030 0e 0000 0000p1 0501 9e 0014 7e 0131 0e 0030 0e 0000 0000p0 5751 4e 001 9 1e 004 1 0e 0030 0e 0000 0000f0 8131 8e 001 3 5e 000 1 0e 0030 0e 0000 0000f0 9311 9e 001 1 4e 001 1 0e 0030 0e 0000 0000f0 9911 9e 001 7 9e 001 1 0e 0030 0e 0000 0000f1 0201 9e 0014 7e 0131 0e 0030 0e 0000 0000p1 0051 9e 001 1 1e 003 1 0e 0030 0e 0000 0000f1 0131 9e 0014 7e 0131 0e 0030 0e 0000 0000p 1 0091 9e 0014 7e 0131 0e 0030 0e 0000 0000p1 0071 9e 0014 7e 0131 0e 0030 0e 0000 0000p1 0061 9e 001 4 2e 003 1 0e 0030 0e 0000 0000f1 0071 9e 0014 7e 0131 0e 0030 0e 0000 0000p1 0071 9e 001 1 6e 004 1 0e 0030 0e 0000 0000f1 0071 9e 0014 7e 0131 0e 0030 0e 0000 0000p1 0071 9e 001 5 1e 004 1 0e 0030 0e 0000 0000fGammavalueachieved 1 0067 设计所得的H 控制器 设计所得的闭环系统的奇异值Bode图如图4所示 图4闭环系统奇异值Bode图 2 函数hinf 函数的调用形式为 sscp sscl hinf G ssu 该函数用的是下面文献中的算法 对于D11不为0的情形 可以用该函数求解 M G Safonov D J N LimebeerandR Y Chiang SimplifyingtheH TheoryviaLoopShifting MatrixPencilandDescriptorConcepts Int J Contr vol 50 no 6 pp 2467 2488 1989 函数的输入变量G为如下定义的广义对象 图5 函数的输入变量中的ssu对应的就是图5中的U s 是个可调参数 一般都取默认值0 此时所求得的H 控制器是中心控制器 sscp sscl hinf G ssu 输出变量中的sscp表示控制器F s sscl表示闭环传递函数 算例 S KS T问题 图中参数如下 图6 本例的H 问题是要求解如下的有约束的优化问题 调用函数hinf时 其输入变量G有自己的调用形式 要用如下的几个函数调用来送进去 W1 900 15 0 010 020 01 W2 0 0001 W3 0 11 3 16 3003 16 P ss ap bp cp dp G augtf P W1 W2 W3 其中ap bp cp dp为对象P的状态空间实现 sscp sscl hinf G ssu 图6 本例中函数hinf的调用形式为 图7闭环系统的奇异值Bode图 sscp sscl hinf G 当时 闭环系统的H 范数为0 8916 当时 闭环系统的H 范数为0 9998 设计所得的H 控制器为 hinf G ssU 该函数用的是下列文献中的算法 1 Doyle J C K Glover P Khargonekar andB Francis State spacesolutionstostandardH2andH controlproblems IEEETransactionsonAutomaticControl vol 34 no 8 pp 831 847 August1989 2 P Gahinet A Nemirovskii A J Laub M Chilali TheLMIControlToolbox Proc oftheIEEEConf onDec andControl 1994 2038 2041 3 函数hinfric 函数的调用形式为 其中输入变量中的G为如下定义的广义对象 输入变量中的r是2维列向量 r 1 表示量测输出的个数 r 2 表示控制输入的个数 输出变量gopt表示最优的H 性能 输出变量k表示H 中心控制器 gopt k hinfric G r 该函数用的是基于Riccati方程的算法 是解析法 其最优控制器是在任意阶控制器中寻优的 所以优化效果比较好 函数对对象的秩的要求比较低 算例 S KS T问题 图中参数如下 图8 本例的H 问题是要求解如下的有约束的优化问题 函数调用时的广义对象G由下面的函数送进去 hinfric G r G ltisys A B C D gopt k hinfric G 1 1 H 优化设计所得的闭环系统的最终值为0 1011 所得的H 中心控制器为 本例中函数的调用形式如下 图9加权闭环系统的奇异值Bode图 该函数用的是下面文献中的算法 是基于LMI的算法 对对象的秩的要求也比较低 1 P Gahinet P Apkarian ALinearMatrixInequalityApproachtoH Control Int J ofRobustandNonlinearcontrol 1994 4 421 448 函数的调用形式为 其中输入变量中的G为如下定义的广义对象 输入变量中的r是2维列向量 r 1 表示量测输出的个数 r 2 表示控制输入的个数 输出变量gopt表示最优的H 性能 输出变量k表示H 中心控制器 gopt k hinflmi G r 4 函数hinflmi 算例 S KS T问题 图10 图中参数如下 本例的H 问题是要求解如下的有约束的优化问题 函数调用时的广义对象G由下面的函数送进去 hinflmi G r G ltisys A B C D gopt k hinflmi G 1 1 H 优化设计所得的闭环系统的最终值为0 1064 所得的H 中心控制器为 本例中函数的调用形式如下 图11加权闭环系统的奇异值Bode图 5 函数hinfnorm 该函数用来计算系统的H 范数 函数的调用形式为 out hinfnorm sys tol iiloc sys 为系统矩阵 当已知系统的传递函数T时 调用如下的函数可得到sys A B C D ssdata T sys ltisys A B C D out hinfnorm sys tol iiloc tol 为H 范数的上下界之间的相对精度 iiloc 为假定的范数值所对应的初始频率点 out是一个的行向量 分别表示的下界 上界以及下界所对应的频率 算例 首先调用函数ltisys将系统矩阵sys送进去 sys ltisys A B C D 如下调用函数hinfnormhinfnorm sys tol iiloc out hinfnorm sys 0 000001 0 求得的out 110 设系统的状态空间实现为 下界 上界 下界对应频率 6 函数norminf 该函数的调用形式为 gain peakf norminf sys tol 输入变量中的sys为系统矩阵 tol为相对精度 输出变量gain为求得的峰值增益 peakf为该范数值对应的频率 1 函数dhfsyn 对应的连续函数为hinfsyn 的调用形式 k g gfin dhfsyn G nmeas ncon gmin gmax tol h 输入变量中的h为采样时间 G为离散