信息 与控 制
年 第 期
自适应 鲁棒控 制算 法
刘新平 刘宏 才
日匕京科技大 学 自动化系
叫
摘要 本文根据 内模 原理 提 出了 三种新 的 自适应鲁榨 控 制器 , 克服 了
等人 提 出的鲁棒 控 制
器 的缺点 本文提 出 的算法适用于 跟踪任何参考信号 , 当系统的 阶次发生 变化时均 能跟 踪与柯节 具有 较强
的鲁棒性
关锭词 自适应控 制 , 自适应鲁棒 控制 , 鲁棒 自适应控制 反诩控制
前言
应用 自适应控制 时 常做 一 些假设 , 例如 已 知被控 系 统的 阶次 、 时延及有关扰动特性等 。 但实 际系统
的复杂性 , 加上工 况和 环境 的变化 , 系统特性 的 非线性和 随 机干扰 , 就可 能使原系统发 散 , 控制算法 失
等人 认为其根本反 因 在于 现存 的 自适应算 法都没有从鲁 棒控制器 出发 , 因此要 获得各
效 。
棒性强 的 白适 应控 制算法 , 应先考虑 一个具有 各棒 性 的控 制器 , 然后 再设法使它 自适应 。 所 以 , 鲁棒性
自适 应控 制算法就可 由一个具有各棒性 的参数估计器加上一 个具有 鲁 棒性 的控制器组成 。
和
勺 提 出 了 内模原理 , 他 们认为各 棒 性控 制应是反 馈 控制 , 在反馈 回 路 中应含
有参 考 输 人 和 千扰 的动 态 模 态 。
‘ 〕 指 出 鲁 棒 控 制 器 可 分 为伺 服 机 构 和 稳 定 补 偿 机 构 两 部
分 当系 统存 在未 建模 误差 或外 界 千扰 时稳定补偿 机构 维 持整个 系统 的稳定性 而伺服机构保持参考信号
, 提 出 如 果 采 川 非 最 小 系 统 能 增 强 系 统 的 鲁 棒 性 , 基 于 上 述 思 想
的 跟 踪 和 训 节 。
等人 ’ 在 此 基 础 土 进 步 改
一 种 新 的 自适 应 鲁 棒 控 制 策 略
进 , 提 出 种 仓 棒 控 制 器 然 后 川 确 定性 等价 原理 再使它 自适应 化 , 并采川 死 以法进行参数估 计 , 他们
证明 厂算 法 的 收敛 性 , 并成 功地 应 用 几蒸馏塔 的控 制上
〔 〕提 出
。 和
翻
。 。
和
但
等人 提 出 的算 法 有 不 足 之 处 第 , 要求 被估计 后 的系统 总是 逆 稳定 的 , 只 适川 于 最小相
二 ,
位 系统 , 第 二 , 内模原理 需 要 】个多 顶 式 川 一 ’ , 它需 满足 两个条 件
,, 为 给定 信号 , 叫 为 系 统 不 可 测 的 确 定 性 扰 , 实 际上 要找 到 这样 的
对有 些给 定信号 不 适 用 , 第三 , 该算 法是从考 虑 系 统建 模误 差 出发 而未考虑 时 延 变化 的影 响 。
一 ’凡 一 。, ‘ 一 ,
一 ’是 十 分 困难 的 , 该算 法
本文 灯先 提 出
一 种 适 合 」跟踪 任何 参考 信号 的 新算法 , 且 厂 ’ 多项 式选择容 易 , 然后 再 将它推
等人 提 出的算法 改进 , 使适
。 〔砂 多 步 预 测
的思 想 把
广到 非 最小 相 位系 统 。 最 后 利 用
川 几 寸延 变化 的情 况
各棒控 制 器 算法 一
宜 系统 模型
假设 一
‘ 一个 单输人 一单输 出系统 如下
, 二 一 ’ , 、卜 。 。 一 ’ 一 ‘ · 、卜 望
,
,
一 ’ 一
。
一 ’ 一
, 一 。
一 ’
·
△
一 一
一 , 为后移算子
试 分别 为系统 的输出和输人 ,
袱 为随机干扰 , 。 , 。 分别为 刃 ,
一 ‘ , 凡 一 , 的阶次 , 为
式 中
只
纯时延 ‘
①收 到 本 文 的 时 间是 卫 年吕月 日 。
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
自适应鲁 棒 控制算 法
控制器的结构
根据 内模原理 , 控 制器 的结构 由两部分组成
二 丝生 」丝且
一 ’
叮
△
厂 ’
式中
一 , 为一稳定的多项 式
三笠生一 , 甲 , 尸气 扭冲
。 为 控 制误 差 , 。 二 , 卜夕 ,
以 为参考信号 , 只 为系统 的实 际输 出 , 叭 为一 辅助 信号 ,
它是 由最优性 能指 标求得 的
刀 一 ’
二 十
万 一 ’ 十 · 十
为
当控 制过程 的 确定 性 干 扰 为 未 知时 , 则选
过程的确定性 干扰 已知 时 , 则选 ‘ , 满 足 厂 ’
一 , 的 阶次
一 , 为一稳 定 多项 式 , 也可选
幻为 确定性干 扰 。
二 ,
一 , 二 卜 一 , 当控 制
由式 可 知 , 控 制器 是 由误差 驱 动 的 , 等式 右边第 一项 为右和及机构 , 第 二项 为稳定补偿机构 , 系
统的框 图如 图 所示
自自适 应稳稳
定定 补 偿 机 构构
辅助 系统
辅助 信号 爪 的获取是 关键 , 可认为控 制误差 以 是 山 。 产生 的 一卜面来 推 导 叹幻与 叮 的关系
图
自适应 各 棒 控 制 系 统框 图
一
,
£
二 一
夕
一 夕
二
夕
一
,
式
式 中
人 为 。 一 , 多项 式 , 省 略 了 厂 , 算子
刀 为 刀 一 , 多 项式 , 省 略 了 一 , 算子 由式 和 式 得
侧 刀 。 、卜 付 一 ‘ 、 旦丛且
·
用
凸
夕
冲
由式 和式 得
式中
注 。 七 一 刀 一 。 无
, 无 一 △注
曰
一 甘 △
一
刀
·
,
式 即为辅助系统 , 收敛性要 求式 为最 小相位系统 。
性能指标的极 小化
斌 是根据下列 性 能指 标极 小化来求得
一 于【尸 。 、 ‘ ’ 又乃
△ 无 ’
’
·
一 、
尸 一 ‘ 一 二了 , 仄 一 , 为一 稳定多项式 , 又为控制 因子 , 可 取其为常数或多项 式 。
为了得到 尸。 丸十的的预报值 , 先考 虑下 列恒等式
一
— 巴 仃 一二丁一 一 一
八
尸
习
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
城一中,
一勺协︸沪户习
信 弘 与控 制
年 第 期
, 一 ‘
一 石 万 十 。 二 十众大
矿 一 , 凡十 。矛
。 为 的 阶 次
肠 、、 一‘ 二 十
蜘 分别为 习 一 ,
性式 和 式 得
“
, 献 ’ 的 阶次 。
一 臀 , 十 “
十 刃 一 △“ 御 ,
司 十 命
因此可 得 凡 十的的预报值为
£
十
写成 向量形式
二 缪 。 共 。 一 △、
诬沪
夕
, 、 司
。 ’ 仕
、 二 。 , 仕
式中
油工、七卜令价︻‘少
小‘‘阴百‘£卜与·二又丫厂︼
价 、 丁 一 。 、 ,
“ 一 ‘。 … ,
。 · 。 , 。 、一 二 , 上
,
£
,
“ 一
。 、 , …‘ ,
, … ,
, …,
上 。 一 力 ,
才
, …
夕 ,
一 一
为 召 的 阶次 。 夕 为 的阶次 , 为 △ 用 的阶次 , 为 一 , 为 的阶次 。
一
一 无 一 , 一 , 一 。
一
一
预报误差 仪 十 刃 二 广 ’城 十的
由式 可翅卜 步 化简性 能 指标 , 得
‘
汀
。 又’ 尹 △。
。 无 一 。
王 止二 弋可 求 出 叮伏 , 然后 代人 式 且 可求 出 仕 。
白二返应控 制算 法 的 少骤
竹 根据 式 辨 识
再将 孟代 人 式 求 出 。 ’ 无 司
式 可 求
将 即 代入 式 求 。
,
手亏仁复 。
白适应 树 科棒控制 器
名 下 歹 条 件满 足 时 , 即可得 到 自适应控 制器 。
系统模 耳健为
月 ,
一 ’ 一
一 ’
一 , 凡 一 , 一 方。
一
一, 又 即最小方差控 制
外部 输人信号 为一 恒定值 , △凡 二 仅
选 取 仕一 ’ 一 卜 一 ’。
则 至式 得
。二 一
介 〔石斌 一 ,琪厂 △
天 十
一 ’ 厂
£
口 。 、 一
。 仅“为传统 的
形 式 。
扩展 的拜棒控 制器 算法二
上 一 节的算法 需 满足三 个条件
均 辅助 系统 是逆稳定 的
月 。仔 ” , 刀 、 一 ’ 和 办 一 , 没有公 囚子 ,
即互逆
一 , 是稳 足 的 。
参
一 二 一 ’
卜 ,
—几 一一 —一一
十
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
自适 应答 榨 控 制算 法
条 件 过于苛刻 , 本节利用 极 点配置 方法来求 叮
对辅助系统 进 一步改造 便之扩展为最小 位 系统
设广 义输 出
为
无 一 £ 一 ’ 七
一 甘 。
叮
式 中
仕一 , 为一稳定多项 式
适 当选 取
,
,
厂‘ 使
一 , 为所希望 的稳定多项 式 ,
为简单起见 , 多项 式 中算 了 扩 ’都右 略
£ 一 , 刀 一 ,
一 , 月 一 ,
一 ,
由式
和 得
月 丸
一 ‘
叮
一 △ ,
夕 ,
、
令万
声 十 月
则式 为
万 、 一 一 ‘万。 、 亡 、 一 £△月 刀 ,
上式 为 扩展 的 系 统 , 通 过适 当 的 选 取 , 使万稳定 , 即土式 为 最小相 位 系统
为 了求得 扩展 系统 的 预 报值 , 先考 虑 卜列 恒等 式
亡 吞万十 户 一
扣中幻式
‘ 一 诃。 十 万,
一 ’
十 …
一
二
十
二
‘
万
万 , 万‘ 分别 为多项式万,
由式 和 一 得
亡的阶次
十 …
一 , 广
几
元 一 ‘
一
‘ 十 的 一 臀叮
一 乌 △, 。 , 、 过
下
司 一
︸﹄
则
司广义输 出的预 报值为
二 臀, “ , 答‘“ ,
同样上式也可写成 向量形式
定义广义输出误差 产 和估计误差 引 如下
一︸
△
夕 ,
‘ 无 二 无 一五 一 , 夕 , 无
叔
一洲无 一司
共 为参考信号
叔 一 , 为可选择 的稳定多 项式
控制辅助信号 刀 由下式给 出
仿第二节方法 同样可求 出
本算法适用于 非最小相位系统 , 但需 满足
丑 扩 ,
,
司
石 一 , 是稳定多项 式
离线选择 一 ’和
一 ‘ , 使 一 ’刀 一 ,
一 , 注 一 ’为一稳定多项式
本算法的缺点是选择的多项 式多 ,
, 的选择较困难
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
信息与控 制
年 第 期
多步预 测的 自适应 鲁棒控制 算法三
等人 闭 提 出 了一种 多步预测
当系 统 的时 延 发 生 较大 变化时 , 上述两种 算 法是不适 用 的 ,
控 制思 想 , 其关键是采用 滚动 式性能指标 , 在时 刻 求得 的控制量 就考虑 了 十 时刻后 步 的
控制误差 由此增强 了鲁棒 性 ,
控制适用 于变时延 、 变结构和非最刁湘 位系统
算法 三与算法一不 同之处在 于性能指标 的不 同
考虑如 下性能指标
一
艺 那 、
’ 艺 以刀 △ 无一
刀 ’
其中
为预测长度
为控制长度 ,
为了求 那 十 多步 的预报值 , 将辅助系统 式重写如下
, 伏 一 △注 刀 , 灸
。
一 刀 一 ‘ , 无
考虑下列 恒等式
式中
以式 和 得
一吸一一,
一一
凡, 十 一 二一二
尹
止二 ,
乓二 司切 ,
一
一 ,
十 , 。 气卜 一 一 什
二 十奋 一‘
一 , , 碑 才 一 。 一
£
、 望 ” 一
‘
,
‘ 十 万乞“凡
一 △刁
, 夕 ,
由上式得 八
的估计值
£
十 ‘ 一方 刀 一
一七 ‘ 叮 一 △
气 口 ”
夕 ,
了占计误差 亡 二 拼
令
将 预报值写 成矩阵形式
产
,
, 二 心 十 , 一
一
”
一 一 一
一 ,
扮一 ‘
一一丫
一
叮
试 。
才 一,
,
,
斗一
十 无
一
万 一 洲 一
‘ 拟 一 ,
才 “ 一泌 一
“ 一 ,澎 一 万
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
自适 应各 棒 控 制 算法
一
丫
”‘
‘汀 ,
,
‘
‘ … … …十 … … …
一 ’ 」
一 “‘
一 ’ 川 瓜 ,
定义 当 户
爪 刀
,
,
二尸 ,
二七 ‘
嵘
。
一 ‘ △汉 ,
,
, 、
才一 干 一
艺 , ,
一 咨一 寸
一 ’。,
一 口 当
一
。
一 , △ , 夕 , 十
当 ,
一
,
· ,
一
瓜 , 均 为 时刻 的 已 知 值
令 。 ’
夕式
为
现 在 考察性 能 指 标
。 ‘
十
, 一
, 。
叮‘ ” 叮’
, … ,
叮‘
一 今 ,
, …,
心 ,
一
一
一
丁
材
‘ 刀 一 万 一 ’
’
己 万 一
气 十 一 爪 一 ’ 二 勺 十 一 , 一
、 。‘ 十 瓜
卜 份
一 万于 艺 。 、
’ 艺 △。 、
一
十 又。 、
一
为 了将性 能指标 化成关于 叮
‘ 的形 式 , 对土式 中 的第三项 和 第 三项 也写 成矩阵形式 , 并经定 义 可得
第二项 尸 △叹 十 司 为
式 中
£
卜二一
二 。, 瓜
。 。‘ 爪
。
。
, …, 尸。 尸 刀一
,
。犷 。 一 ,
…, 尸叹
一
‘几
︸
万
封
‘ 万 一 万 一 才一
“ ’
‘ 万 一
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
亨厂之、 控 制 , 年 第 期
几‘、︸
脚
。
‘ 万 一 万 一 卜 才
“ ‘
‘ 万 一 ,
‘ 二
瓜 , 一 叽 ,
, 哟 , 一 吞,
,
万 一
丁
,
,
,
,
,
吞,
£
少、
净一 十 一 。
艺 , ,
一 场,
一
‘一 一
艺 , , 一 ’丫
二 , 一 , 口 刃‘
, 为 时刻的 已知值
石,
一 刃
一, 廿
︸
,一
£ 一 哟 一 , △ ,
,
一
‘
一
一 内
一, 才
一
一
《 一 哟 一 , △ , 夕
一 十
£ 一
一
‘
£
…咬
叉 ‘,
才
十
一 ’。‘ 、
一
二叨二动
£ 一
一 , △ , 夕 ,
一
爪 ,
一 一
艺 , ,
一 ‘一 沙
一 , 。,
一
十
£ 一
一
△月 ,
,
一
续
一
将式
中的第三项 叭 一司也写成 矩 阵形式 , 然后 也 可得
式中
。犷一 、 。‘ 耘
二 【刀 ,
叮 刃二
。
当
付
们
︸
·· …””
侃 一 呱 ,
。 」
, …, 吞 ,
一
了
这样性能指标 可 写成
, 一 ’ 了 。 ’ 仕 鉴 。梦 。
£犷 。扩
产目 ,
一 ‘ ,
,
, 几
, 。,
「 ,
,
,
二 。, 瓜 爪 爪
,
‘ 。
。 。,
寿 瓜, 几 ,
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
自适应 鲁 榨 控制 算法
刀 “ 万
劝十 们
石 石 痴 十扁
一
, 。‘ 几 了
, 。· 几
, 丫 几 ‘ , 。‘ 吞
令
对式 〔, , 黔一
得辅助 信号 刀 产 为
, ‘ 一 二 ,
。二 、 一 ’ 。二几 。二,沁
由式 可求得 时刻 的 丫 , 由 叮二 叮 ‘ 可得 叮 , 代人 式 可 得 。
本算法适用于变结构 、 变时延 和最小相位系统 如仿照 算法二对辅助系统进 行改造 , 仍 可适 川 于非
最小相位系统
本算 法计算量 大 , 当 , 值较大时 , 计算 时间更长
仿真结果
对参考信号 的跟踪 性能
某加 热炉 为最小相位系 统 , 其模型如 下
二 夕 一 一
夕 一
,
。 一
一
当 给 定 信 号 为 阶 跃 、 斜 坡 、 抛 物 线 、 正 弦 、 方 波 及 指 数 等 ,
一 与二
一 一 一 ’
一
又二
一 ’ 一 一 一 , , 给定 信号 的幅度最大 为
用 递推最 小二乘法拈计参 数 , 噪 声方差 护 二
算法 一 的
仿真结果 如 图 所示 限于 篇幅仅 以方波为例 , 该算法对任何参考信号均 能跟 踪
对 系 统结构变化 时 的性能
在
步 以前 系统结构为实 际值 , 以后 系统结 构变为 。“
, 二 算法一仍能跟踪调 节 , 仿真 结
果如 图 所示
一
图 给定信号为方波时 的仿真 曲线 ’
图 结构变化 时 的仿真 曲线
反 非最小相位系统 算法二 性能
设一非最小相位系统如下
二
夕
,
一
一
夕 一
。 一
一
一 ,
一
取
方 波信号 幅值 为 , 噪 声方 差 尹
一 , , 天 一 , 二 一
一 ,
, 用递 推最小 二 乘法估 计系 统参数 算法 二仍 能跟 踪 调
节 , 仿真结果如 图 所示
© 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net