-
摘要: 随着指纹识别技术的广泛应用,大量指纹图像需要被收集和存储.在指纹识别系统中,对于大容量的指纹数据库,指纹图像必须经过压缩后存储以减少存储空间,本文提出了基于自适应稀疏变换的指纹图像压缩算法.该算法在离线状态下提取指纹图像特征训练超完备字典;在编码过程中,首先利用差分预测编码和稀疏变换将待压缩指纹图像转换到稀疏域,然后对直流系数和稀疏表达系数进行量化和熵编码,从而实现图像信息的压缩.实验表明,在中低码率段,本文算法相比于JPEG、JPEG2000和WSQ等主流压缩算法表现出更优越的率失真性能;在相同码率时,本文算法生成的压缩图像的主观视觉效果更好,指纹识别率更高.Abstract: With the wide application of fingerprint identification technology, a large number of fingerprint images need to be collected and stored. In fingerprint identification, as for the fingerprint database with large-capacity, the fingerprint images must be stored after compression to reduce the storage space. In this paper, a fingerprint image compression algorithm based on adaptive sparse transformation is proposed. The feature of the fingerprint image is extracted offline to train the over-complete dictionary. In the encoding process, the fingerprint image to be compressed is converted to sparse domain by utilizing the differential predictive coding and sparse transformation in the first place; after that the DC coefficients and the sparse coefficients are quantized and entropy coded to achieve the compression of the image information. Experimental results show that the proposed algorithm outperforms the mainstream compression methods, such as JPEG, JPEG2000 and WSQ, in terms of ratio-distortion performance of decoded fingerprint image, especially at low to medium bit rates. At the same bit rate, the compression image generated by the proposed algorithm exhibits better subjective visual effect and higher fingerprint recognition rate.
-
20 世纪 50 年代以来,基于模型的控制理论得到了快速发展与完善.其建立控制系统的过程要经历三个阶段,依次是建立模型、分析模型和依靠模型设计控制律[1].但随着受控对象越来越复杂,如何有效地对受控对象进行建模已成为一个棘手的问题.原因是对于一个系统,如果所建数学模型过于复杂,则难以设计控制律或所得控制律工程上难以实现[2-3].如果过于简单,则难以反映实际系统的动态特性,基于该模型设计的控制律在实际应用中难以取得理想的控制效果. 其次,建立系统模型的方法主要有机理建模和系统辨识两种方法,无论哪种方法,所建模型都是对真实系统的逼近,实际系统一定会存在影响控制系统鲁棒性的未建模动态和其他不确定性因素[4-5].
目前,主要存在两大类方法用以解决上述问题.第一大类方法是基于数据驱动的控制方法. 较为典型的有虚拟参考反馈整定(Virtual reference feedback tuning,VRFT)、同步扰动随机逼近(Simultaneous perturbation stochastic approximation,SPSA)和无模型自适应控制 (Model free adaptive control,MFAC)三种方法. VRFT 是 Guardabassi等 在 2000 年提出的,其特点是根据离线数据设计控制律. 该方法基于模型参考自适应,通过参数辨识直接得到控制器参数,使闭环系统动态特性逼近参考模型[6]. 由于实际系统的动态信息无法通过一次激励完全得到,所以运用该方法设计所得的控制器很难保证闭环系统的稳定性. SPSA 是Spall在 1993 年提出的,其特点是迭代辨识系统参数,控制效果易受系统结构变化或参数变化的影响.使用该方法同样难以保证闭环系统的稳定性[7]. MFAC 方法是侯忠生在 1994年提出的一种无模型自适应控制方法[8], 其特点是控制器的设计不需要系统的模型信息,而是基于输入输出数据直接计算得到控制量.文献[5, 9-10]证明了紧格式和偏格式动态线性化控制方法闭环系统的稳定性.第二大类方法是吴宏鑫等提出的特征建模方法.特征建模方法广泛应用于航天器和工业控制中,是一种综合考虑对象动力学特征和控制性能要求的建模方法[2].该方法与文献[5]提出的全格式动态线 性化模型有异曲同工之处.不同的是,特征建模方法强调在采样周期足够小的情况下,特征模型与原模型的输入输出具有等价性[2],而基于全格式的动态线性化方法则考虑存在时变线性系统与原系统等价[5, 9-10].
侯忠生提出的无模型自适应控制,因其理论严密和计算量较小,在工程上有广阔的应用前景. 经过多年的发展,侯忠生分别针对单输入单输出系统和多输入多输出系统设计了基于紧格式、偏格式、全格式的无模型自适应控制. 在文献[11] 中,提出了控制器的动态线性化方法.在文献[12]1中,提出了一种新的PPD(Pseudo-partial derivative)参数估计方法.但以上方法均未考虑实际物理系统中可能出现的执行器饱和问题.由于执行器的执行能力都是有限的,而执行器饱和这种非线性环节可能会导致控制效果严重恶化甚至丧失闭环稳定性[13].具体来讲,忽略可能出现的执行器饱和,将无法发挥其最大的控制能力.同时在执行器饱和后,作用在系统上的实际控制量并不取决于控制器计算的数值,故用于模型参数辨识的输入输出就可能出现明显偏差,严重时会造成系统的不稳定.文献[12]虽设置了控制量的速率饱和,但没有对可能出现的饱和做任何优化,同时未考虑执行器的幅值饱和.
借鉴文献[5] 中 PJM(Pseudo Jacobian matrix)辨识技术,本文针对执行器饱和问题,提出了基于紧格式改进无 模型自适应控制方法,并证明了该控制方法的闭环稳定性. 进一步,对蒸馏塔 Wood/Berry的仿真实验表明了该方法具有更强的跟踪能力 且对初始值不敏感.本文对执行器同时设置了速率饱和和位置饱和并进行了优化.
1本文与文献[12]的区别主要在于本文提出的方法考虑了位置饱和对系统造成的影响,并对可能出现的执行器饱和问题进行了优化.在仿真时,文献[12]并没有考虑不同初始参数对系统造成的影响.而本文用仿真结果说明改进后的无模型自适应方法具有初值不敏感性.
1. 问题描述
本节简述文献[5]提出的无模型自适应控制方法的相关概念以及计算思路. 在此基础上,分析已有方法的不足并提出本文要解决的问题.
考虑如下多输入多输出离散系统:
$\begin{array}{*{35}{l}} y(k+1)= & \ f(y(k),\cdots ,y(k-{{n}_{y}}), \\ {} & \ u(k),\cdots ,u(k-{{n}_{u}})) \\ \end{array}$
(1) 其中, ${{ u}(k)\in{\bf R}^m}$ , ${{ y}(k)\in{\bf R}^n}$ ,它们分别是系统 ${k}$ 时刻的输入和输出向量. ${n_y}$ 和 ${n_u}$ 是未知整数. ${{ f}(\cdot)}$ 是一个未知非线性函数. 假设 ${{f}(\cdot)}$ 关于 ${{ u}(k)}$ 的偏导数连续,且系统 (1) 满足广义Lipschitz 假设. 则可以得到定理1[5].
定理1. 对满足广义Lipschitz假设且 ${{ f}(\cdot)}$ 关于 ${{ u}(k)}$ 偏导数连续的非线性系统(1),当 ${\|\Delta{ u}(k)\|}\neq0$ 时,一定会存在一个被称为PJM (Pseudo Jacobian matrix) 的时变参数 ${\Phi_c(k)\in{\bf R}^{n× m}}$ ,使系统(1) 可转化为如下紧格式动态线性化数据模型(Compact form dynamic linearization,CFDL):
$align\Delta y(k+1)={{\Phi }_{c}}(k)\Delta u(k)$
(2) 其中
$\begin{align} & {{\Phi }_{c}}(k)= \\ & \qquad \left[ \begin{matrix} {{\phi }_{11}}(k) & {{\phi }_{12}}(k) & \cdots & {{\phi }_{1m}}(k) \\ {{\phi }_{21}}(k) & {{\phi }_{22}}(k) & \cdots & {{\phi }_{2m}}(k) \\ \vdots & \vdots & \ddots & \vdots \\ {{\phi }_{n1}}(k) & {{\phi }_{n2}}(k) & \cdots & {{\phi }_{mm}}(k) \\ \end{matrix} \right]\in {{\mathbf{R}}^{n\times m}} \\ \end{align}$
(3) 且对任意时刻 $k$, $\|{\Phi}_c(k)\|$ 有界.
基于紧格式动态线性化的无模型自适应控制是运用参数辨识的方法,动态计算 PJM 时变参数的数值,并在此基础上进行控制的算法.详细的计算步骤可参见文献 [9-14]. 在现有的方法中,仅文献 [12]考虑了不完全的执行器饱和问题. 对于实际的物理系统,执行器的执行能力都是有限的,具体体现在控制量的幅度与变化速率的有界性.执行器的执行能力可完整表示为
$\left\{ \begin{array}{*{35}{l}} \Delta {{u}_{\min }}(k)\le \Delta u(k)\le \Delta {{u}_{\max }}(k) \\ {{u}_{L}}\le u(k)\le {{u}_{U}} \\ \end{array} \right.$
(4) 其中, $\Delta{ u}_{\min}$ 和 $\Delta{ u}_{\max}$ 分别表示控制量变化速率的最小和最大值, ${ u}_L$ 和 ${ u}_U$ 分别表示控制幅度的最小和最大值.
2. 无模型自适应控制的执行器饱和优化
本节提出基于紧格式的优化无模型自适应控制,通过综合分析执行器的执行能力给出控制算法.
2.1 改进控制器设计
在辨识 PJM 参数时,采用如下计算方法[5]:
$\begin{align} & {{{\hat{\Phi }}}_{c}}(k)={{{\hat{\Phi }}}_{c}}(k-1)+ \\ & \frac{\eta (\Delta y(k)-{{{\hat{\Phi }}}_{c}}(k-1)\Delta u(k-1))\Delta {{u}^{\text{T}}}(k-1)}{\mu +\|\Delta u(k-1){{\|}^{2}}} \\ \end{align}$
(5) 其中, $\mu>0$ , $\eta\in(0,2]$ . 则可知 PJM时变参数 $\hat{\Phi}_c(k)$ 有界.
进一步对 $y(k+1) $ 进行滤波[12]
$\left\{ \begin{array}{*{35}{l}} \overline{y}(k+1)=\widehat{y}(k)+{{{\hat{\Phi }}}_{c}}(k)\Delta -u(k) \\ \widehat{y}(k+1)=\overline{y}(k+1)+K(-y(k+1)-\overline{y}(k+1)) \\ \end{array} \right.$
(6) 其中, $K\in[0,1]$ 为常数, $\bar{{ y}}(k+1) $ 为 ${ y}(k+1) $ 的预报值, $\hat{{ y}}(k+1) $ 为 ${ y}(k+1) $ 的滤波值.滤波值是控制量 ${ u}(k)$ 作用于控制对象后,测量得到 ${ y}(k+1) $ ,然后再进行计算得到的.
控制算法的设计目标是让预报值 $\bar{{ y}}(k+1) $ 尽可能地跟踪参考输入 ${ y}^*(k)$ ,同时控制量符合实际要求. 为达成此目标,考虑如下控制输入准则函数
$\begin{align} & J(\Delta u(k))= \\ & \qquad \|{{y}^{*}}(k+1)-\overline{y}(k+1){{\|}^{2}}+\lambda \|\Delta u(k){{\|}^{2}} \\ \end{align}$
(7) 其中,λ>0 是一个权重因子,用来限制控制量的变化幅度.
对于执行器执行能力(4),设
$\left\{ \begin{array}{*{35}{l}} {{u}_{LB}}(k)=\max \{\Delta {{u}_{\min }}(k),{{u}_{L}}-u(k-1)\} \\ {{u}_{UB}}(k)=\min \{\Delta {{u}_{\max }}(k),{{u}_{U}}-u(k-1)\} \\ \end{array} \right.$
(8) 则在 $k$ 时刻控制量增量范围可表示为
$\left[ \begin{matrix} -I \\ I \\ \end{matrix} \right]\Delta u(k)\le \left[ \begin{matrix} -{{u}_{LB}}(k) \\ {{u}_{UB}}(k) \\ \end{matrix} \right]$
(9) 联立式 (7) 和式 (9),并做如下变换:
$\left\{ \begin{array}{*{35}{l}} x=\Delta u(k) \\ E=2({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I) \\ F=-2\hat{\Phi }_{c}^{\text{T}}(k)({{y}^{*}}(k+1)-\overline{y}(k)) \\ G=\|{{y}^{*}}(k+1)-\widehat{y}(k){{\|}^{2}} \\ M=\left[ \begin{matrix} -I \\ I \\ \end{matrix} \right] \\ [4mm]\gamma =\left[ \begin{matrix} -{{u}_{LB}}(k) \\ {{u}_{UB}}(k) \\ \end{matrix} \right] \\ \end{array} \right.$
(10) 忽略常量 $G$ ,其等价于求解:
$\left\{ \begin{array}{*{35}{l}} J(x)=\frac{1}{2}{{x}^{\text{T}}}Ex+{{x}^{\text{T}}}F \\ Mx\le \gamma \\ \end{array} \right.$
(11) 当x取何值时,J的数值最小. 对于极值问题(11),采用比较成熟的 Hildreth 二次编程方法求解[15].
定义 ${H}$ 、 ${K}$ 矩阵和 ${ \chi}$ 向量如下:
$\left\{ \begin{array}{*{35}{l}} H=M{{E}^{-1}}{{M}^{\text{T}}} \\ K=\gamma +M{{E}^{-1}}F \\ \chi \le 0 \\ \end{array} \right.$
(12) 其中, $\chi \in {{\mathbf{R}}^{n\times 1}}$ 是拉格朗日乘子.采用如下步骤迭代求改进解:
$\begin{align} & \chi \leftarrow 0 \\ & iter\leftarrow 1 \\ & ite{{r}_{\max }}\leftarrow 20 \\ & \epsilon \leftarrow {{\epsilon }_{0}} \\ & \text{while}(\|{{\chi }^{iter}}-{{\chi }^{iter-1}}\|<\epsilon \ \text{or}\ iter<ite{{r}_{\max }})\ \\ & \{\text{for}\ i=1:n \\ & \{\omega \leftarrow -\frac{K(i)+\sum\limits_{j=1}^{i-1}{H}(i,j)\chi (j)+\sum\limits_{j=i+1}^{n}{H}(i,j)\chi (j)}{H(i,i)} \\ & \chi (i)\leftarrow \max (0,\omega ) \\ & \ iter\leftarrow iter+1 \\ & \ \} \\ & \ \} \\ & x=-{{E}^{-1}}(F+{{M}^{\text{T}}}\chi ) \\ \end{align}$
(13) 其中,迭代次数 $iter_{\max}$ 可变.从以上的计算过程中可以看出,该方法基于求解原问题的对偶问题,通过迭代逼近拉格朗日乘子来得到最终解.在整个运算过程中没有矩阵求逆运算,因而计算量小、易于编程实现.由文献[15]可知,当不等式可行解不为空时,该算法收敛.同时易知没有约束式(9) 时,式(7) 的最优解为 ${ x}=-E^{-1}F$ .结合式 $(\ref{xiuci})$ 可知,增加约束后的最优解相当于在原有的基础上添加了一个修正项,以满足不等式约束(9) .
2.2 改进控制算法闭环稳定性证明
为了严格分析改进算法的闭环稳定性,做如下假设:
假设1.存在充分大的 λ使得 $\Phi{E}^{-1}({F}+{M}^{\rm T}{ \chi})$ 正定.
如果假设 1 不能满足,则说明该改进算法无法保证此系统在遇到执行器饱和问题时的闭环稳定性,需要更复杂的控制算法才能进行有效控制.
基于假设 1可以证明定理2.
定理2. 对于非线性系统 $(1)$ ,当 $\Delta{ u}(k)$ 幅值有界时,在满足假设 1条件下,由辨识方案(1)和本文提出的迭代算法具有如下性质: 当 ${ y}^*(k+1) $ $=$ ${ y}^*$ $=$ ${\rm {const}}$ 时,存在一个正数 ${{\lambda }_{\min }}>0$ ,使得当 $\lambda \ge {{\lambda }_{\min }}$ 时,有:
1) 系统跟踪误差序列是有界的,即 $\|{ y}(k+1) -{ y}^*\|$ 有界.
2) 闭环系统是 BIBO (Bounded-input bounded-output)稳定的,即输出序列 $\{{ y}(k)\}$ 和输入序列 $\{{ u}(k)\}$ 是有界的.
证明.
1) 证明 $\|\tilde{{ y}}(k)\|=\|{ y}(k)-\hat{{ y}}(k)\|$ 有界. 由文献[15]中的定理, $\|\Phi_c(k)-\hat{\Phi}_c(k)\|$ 有界,且 $\Delta{ u}(k)$ 有界,不妨设
$\|{{\Phi }_{c}}(k)-{{\hat{\Phi }}_{c}}(k)\|\|\Delta u(k)\|\le b$
(14) 则有
$\begin{align} & \|\widetilde{y}(k+1)\|=\|y(k+1)-\widehat{y}(k+1)\|= \\ & (1-K)\|\widetilde{y}(k)+{{{\tilde{\Phi }}}_{c}}(k)\Delta u(k)\|\le \\ & (1-K)(\|\widetilde{y}(k)\|+\|{{{\tilde{\Phi }}}_{c}}(k)\Delta u(k)\|)\le \\ & (1-K)\|\widetilde{y}(k)\|+(1-K)b\le \cdots \le \\ & {{(1-K)}^{k}}\|\widetilde{y}(1)\|+\frac{(1-k)b(1-{{(1-K)}^{k}})}{K} \\ \end{align}$
(15) 故 $\|\tilde{{ y}}(k)\|$ 有界,设为 $c$ ,即 $\|\tilde{{ y}}(k)\|\leq{c}$ .
2) 证明存在 λ,使得:
$I-{{\Phi }_{c}}(k)S(k){{({{\hat{\Phi }}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)$
(16) 的特征根的绝对值都小于1.其中, ${S}(k)$ 为一个对角矩阵, $S_i(k)$ 表示矩阵 $S(k)$ 中第 $k$ 个对角元.另记 ${\delta}(k)$ 为 $m$ 维列向量, $\epsilon_0$ 为某一大于0的正数.且 $S(k)$ , ${ \delta}(k)$ 满足
$\begin{align} & \text{for}\ i=1:m \\ & \text{if}\ {{({{E}^{-1}}F)}_{i}}\le {{\epsilon }_{0}}\ \text{then}\ \\ & {{S}_{i}}(k)\leftarrow \frac{{{({{E}^{-1}}(F+{{M}^{\text{T}}}\chi ))}_{i}}}{{{\epsilon }_{0}}} \\ & \ {{\delta }_{i}}(k)\leftarrow \frac{{{(\mathbf{x})}_{i}}}{{{S}_{i}}}-{{({{E}^{-1}}F)}_{i}}\ \\ & \text{else}\ \ \\ & {{S}_{i}}(k)\leftarrow \frac{{{({{E}^{-1}}(F+{{M}^{\text{T}}}\chi ))}_{i}}}{{{({{E}^{-1}}F)}_{i}}} \\ & \ \ {{\delta }_{i}}(k)\leftarrow 0 \\ \end{align}$
则易知
$x=S(k){{E}^{-1}}F+S(k)\delta (k)$
(17) 用 $\|\cdot\|_2$ 表示矩阵的谱范数,用 $\rho(\cdot)$ 表示矩阵的谱.由于 $\Phi_c(k)$ , ${S}(k)$ , $\hat{\Phi}_c^{\rm T}(k)$ , ${\delta}(k)$ 均有界,则不妨设
$\begin{align} & \|{{\Phi }_{c}}(k){{\|}_{2}}\|S(k){{\|}_{2}}\|\hat{\Phi }_{c}^{\text{T}}(k){{\|}_{2}}\le {{e}_{1}} \\ & \|S(k){{\|}_{2}}\|\hat{\Phi }_{c}^{\text{T}}(k){{\|}_{2}}\le {{e}_{2}} \\ & \|{{\Phi }_{c}}(k){{\|}_{2}}\|S(k)\delta (k){{\|}_{2}}\le {{e}_{3}} \\ \end{align}$
(18) 取 $e_0\!=\!\max\{e_1,e_2\}$ , $λ_{\min}\!=\!e_0+\rho(\hat{\Phi}_c(k)\hat{\Phi}_c^{\rm T}(k))$ ,因为 $(\hat{\Phi}_c\hat{\Phi}_c^{\rm T}(k)+λ{I})^{-1}$ 为对称矩阵,且当 $λ>λ_{\min}$ 时有
$\begin{align} & \|{{({{{\hat{\Phi }}}_{c}}\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}{{\|}_{2}}= \\ & \qquad \rho ({{({{{\hat{\Phi }}}_{c}}\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}})<\frac{1}{{{e}_{0}}} \\ \end{align}$
(19) 则进一步有
$\begin{align} & \rho ({{\Phi }_{c}}(k)S(k){{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k))\le \\ & \qquad \|{{\Phi }_{c}}(k){{\|}_{2}}\|S(k){{\|}_{2}}\|{{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}{{\|}_{2}}\times \\ & \qquad \|\hat{\Phi }_{c}^{\text{T}}(k){{\|}_{2}}<1 \\ \end{align}$
(20) 由假设1可知,当 $\lambda >{{\lambda }_{\min }}$ 时,式 (16) 的特征根的绝对值均小于1.
3) 证明跟踪误差有界.由步骤2知
$\rho (I-{{\Phi }_{c}}(k)S(k){{({{\hat{\Phi }}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k))<1$
(21) 则存在足够小 $\epsilon$ 和范数 $\|\cdot\|$ ,使得
$\begin{align} & \|I-{{\Phi }_{c}}(k)S(k){{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)\|\le \\ & \qquad \epsilon +\rho (I-{{\Phi }_{c}}(k)S(k)\text{ }\lambda \\ & \qquad {{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k))<1 \\ \end{align}$
(22) 对任意 $k$ 取
$\begin{align} & \|I-{{\Phi }_{c}}(k)S(k){{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)\|< \\ & {{d}_{1}}<1 \\ & \|{{\Phi }_{c}}(k)S(k){{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)\|<{{d}_{2}} \\ \end{align}$
(23) 则
$\begin{align} & \|e(k+1)\|= \\ & \|{{y}^{*}}(k+1)-y(k+1)\|= \\ & \|{{y}^{*}}(k+1)-y(k)-{{\Phi }_{c}}(k)S(k)\times \\ & {{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)({{y}^{*}}(k)- \\ & \widehat{y}(k))-{{\Phi }_{c}}(k)S(k)\delta (k)\|\le \| \\ & I-{{\Phi }_{c}}(k)S(k){{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}} \\ & \hat{\Phi }_{c}^{\text{T}}(k)\|+\|{{y}^{*}}(k)-y(k)\|+\|{{\Phi }_{c}}(k)S(k)\times \\ & {{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)\|\|\widetilde{y}(k)\|+\| \\ & {{\Phi }_{c}}(k)S(k)\delta (k)\|\le \\ & {{d}_{1}}\|e(k)\|+{{d}_{2}}c+{{e}_{3}}\le \cdots \le \\ & d_{1}^{k}\|e(1)\|+\frac{({{d}_{2}}c+{{e}_{3}})(1-d_{1}^{k})}{1-{{d}_{1}}} \\ \end{align}$
则 $\|{ e}(k)\|$ 有界,设其界为 $f_0$ ,即 $\|{ e}(k)\|\leq{f_0}$ .
4) 由 ${ y}^*(k)$ 有界知 ${ y}(k)$ 有界.又因为在求解 ${ x}$ 的过程中,加入了控制量位置饱和限制,所以 $\|{ u}\|$ 满足约束条件,故有界.进而可知 ${ u}(k)$ 有界. $\Box$
推论1. 如果不使用滤波器,即 $\hat{{ y}}(k)={ y}(k)$ .且当 ${ \delta}(k)\equiv 0$ 成立时,定理2中的结论可加强为:
1) 系统跟踪误差序列收敛, ${{\lim }_{k \to \infty } }\|{ y}(k+1) -{ y}^*\|=0$ ;
2) 闭环系统是 BIBO 稳定的,即输出序列 $\{{ y}(k)\}$ 和输入序列 $\{{ u}(k)\}$ 是有界的.
证明.
$\begin{align} & \|e(k+1)\|\le \ \\ & \|{{y}^{*}}(k+1)-y(k)-{{\Phi }_{c}}(k)S(k) \\ & \times \ {{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)({{y}^{*}}(k)-\widehat{y}(k))\|\le \ \\ & \|I-{{\Phi }_{c}}(k)S(k){{({{{\hat{\Phi }}}_{c}}(k)\hat{\Phi }_{c}^{\text{T}}(k)+\lambda I)}^{-1}}\hat{\Phi }_{c}^{\text{T}}(k)\|\times \\ & \|{{y}^{*}}(k)-y(k)\|\le \\ & \ {{d}_{1}}\|e(k)\|\cdots \le {{d}_{1}}^{k}\left\| e(1) \right\| \\ \end{align}$
(24) 故知 ${{\lim }_{k \to \infty } }\|{ y}(k+1) -{ y}^*\|=0$ ,结论1) 得证. 结论2) 的证明类似于定理2中结论2) 的证明. $\square$
推论1中 ${ \delta}(k)\equiv0$ 的物理含义是在没有式(9) 约束的前提下,仅根据式(7) 计算所得的 ${\Delta{ u}}$ 不为零.在实际系统中,由于噪声和数值计算误差的存在,求解式(7) 所得的数值将很少为零,所以该假设在实际系统中有一定意义.
3. 仿真实验
蒸馏塔被广泛应用于化学工业.但在蒸馏塔内出现的延迟给控制律设计造成了很大困难.采用无模型自适应控制可进行有效的跟踪控制[12].在本仿真中,使用了 Wood/Berry蒸馏塔模型,如图 1所示[16].其中 $u_1$ 代表回流速率(IB/ $\min$ ), $u_2$ 代表蒸汽流量 (IB/ $\min$ ), $y_1$ 代表上部成份 (mol $%$ methanol), $y_2$ 代表底层成份 (mol $%$ methanol).选择如下离散系统作为 Wood/Berry 蒸馏塔模型.
$\left\{ \begin{array}{*{35}{l}} {{y}_{1}}(z)=\frac{0.7665}{z-0.9419}{{u}_{1}}(z)+\frac{0.9{{z}^{-2}}}{z-0.9535}{{u}_{2}}(z) \\ [4mm]{{y}_{2}}(z)=\frac{0.6055{{z}^{-6}}}{z-0.9124}{{u}_{1}}(z)+\frac{1.3472{{z}^{-2}}}{z-0.90311}{{u}_{2}}(z) \\ \end{array} \right.$
期望输出信号为
$\begin{align} & y_{1}^{*}(k)=\left\{ \begin{array}{*{35}{l}} 40, & k\le 1000 \\ 90, & k>1000 \\ \end{array} \right. \\ & y_{2}^{*}(k)=\left\{ \begin{array}{*{35}{l}} 40, & k\le 1000 \\ 85, & k>1000 \\ \end{array} \right. \\ \end{align}$
(25) 为了对比在出现执行器饱和时原方法与改进后方法的控制效果,现描述改进前 ${\hat\Phi_c (k)}$ 的计算方法如下[12]:
${{{\hat{y}}}_{i}}(k+1)={{{\hat{y}}}_{i}}(k)+\Delta {{u}^{\text{T}}}(k)\hat{\phi }_{i}^{\text{T}}(k)+{{k}_{i}}{{{\tilde{y}}}_{i}}(k)$
(26) $\begin{align} & \hat{\phi }_{i}^{\text{T}}(k+1)=\hat{\phi }_{i}^{\text{T}}(k)+2\Delta u(k){{({{\left\| \Delta u(k) \right\|}^{2}}+{{\mu }_{i}})}^{-1}}\times \\ & \ ({{{\tilde{y}}}_{i}}(k+1)-{{F}_{i}}{{{\tilde{y}}}_{i}}(k)) \\ \end{align}$
(27) 其中, ${{\hat y}_i}(k)$ 为第 $i$ 个输出分量的估计值, ${{\tilde y}_i}(k)=$ ${y_i}(k)-{{\hat y}_i}(k)$ 为相应的估计误差. ${F_i} = 1 -{k_i}$ , ${k_i}$ 是 $K$ 矩阵对角线上相应的元素. $ {{\hat {\phi}}_i}^{\rm T}(k)$ 为 ${\hat \Phi_c(k)}$ 矩阵的第 $i$ 个行向量.改进前 ${ u}(k)$ 的计算方法如下[12]:
$\begin{align} & u(k)=u(k-1)+{{{\hat{\Phi }}}_{c}}^{\text{T}}(k){{\left[ {{{\hat{\Phi }}}_{c}}(k){{{\hat{\Phi }}}_{c}}^{\text{T}}(k)+\alpha \right]}^{-1}}\times \\ & \ [{{y}^{*}}(k+1)-\hat{y}(k)-K\tilde{y}(k),\ \ \left\| \Delta u(k) \right\|\le \delta \\ \end{align}$
(28) $u(k)=u(k-1)+\delta \text{sgn}(\Delta u(k)),\ \Delta u(k)$
(29) 在本文仿真过程中, $\alpha={\rm{diag}}\{0.003,0.0015\}$ 且 $\delta =0.02$ .改进前与改进后算法共用以下参数并取相同数值:采样周期为 $T_s=1$ s, $K = {\rm{diag}}\{0.9$ , $0.9\}$ , $\mu_1=\mu_2=9$ ,PJM 参数初始值为 ${{\hat{\Phi }}_{c}}(0)=\left[ {matrix} 910 & 750 \\ 450 & 520 \\ {matrix} \right]$.
同时限定被控系统执行器执行能力为
$\left\{ \begin{array}{*{35}{l}} 0\le {{u}_{1}}\le 1, & -0.02\le \Delta {{u}_{1}}\le 0.02 \\ 0\le {{u}_{2}}\le 4, & -0.02\le \Delta {{u}_{2}}\le 0.02 \\ \end{array} \right.$
(30) 在此基础上,分别对改进前后的算法进行仿真可得图 2.
2原控制方法之所以会出现一个控制输入饱和,而另一个却接近于0,是因为进入稳态后,该控制算法不断地试图通过增加 ${u_2}$ 来消除误差.但受执行器饱和的影响,系统的状态并没有改变,最终导致了净差的产生.
为了对比优化前后控制算法对初始参数的敏感度,设初始参数 $\hat{\Phi}_c(0) ={bmatrix}1000& 1000\\1000& 1000{bmatrix}$ .再进行仿真后可得图 3.在没有限制执行能力时,无模型自适应控制可以较好地跟踪信号[12].但对比图 3可以发现,改进前无模型自适应控制在遇到执行器饱和后出现了明显的静差.这是由于控制器在计算控制输出时并没有考虑执行器的执行能力,导致两个执行器中有一个饱和,而另一个却没有发挥应有的控制能力.进一步分析图 3可知,在初始参数有摄动的情况下,改进前的控制算法出现了剧烈的抖动,并基本失去了跟踪参考输入的能力.这是由于在原算法的计算过程中,没有考虑执行器的实际情况,间接导致了系统无法校正PJM参数的误差,最终引起系统剧烈抖动并产生了较大的净差. 相比之下,改进后无模型自适应控制算法充分考虑了执行器的执行能力,可有效地跟踪参考输入,显示了算法的有效性.
4. 结论
本文针对传统无模型自适应控制算法无法应对执行器饱和的问题,提出了一种改进无模型自适应控制算法.并对该算法的闭环稳定性进行了严格证明.该算法具有实现简单、计算量小的优点. 结合蒸馏塔 Wood/Berry模型对比了算法改进前后的控制效果.仿真结果说明,改进算法相比传统方法具有跟踪能力强和对初始参数依赖弱的优点,能有效处理执行器饱和问题.
-
表 1 四种压缩算法的时间复杂度比较 (s)
Table 1 The comparison of time complexity about 4 compression algorithms (s)
码率 (bpp) 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 JPEG 0.16 0.16 0.16 0.16 0.17 0.17 0.17 0.17 0.17 JPEG2000 0.08 0.08 0.08 0.08 0.07 0.08 0.08 0.08 0.08 K-SVD-SR 1.26 1.74 2.20 2.64 3.19 3.61 4.13 4.49 4.89 本文算法 2.13 2.80 3.26 2.88 4.53 4.95 5.83 5.76 6.84 -
[1] Pennebaker W B, Mitchell J L. JPEG:Still Image Data Compression Standard. US:Springer, 1993. http://cn.bing.com/academic/profile?id=1540900967&encoded=0&v=paper_preview&mkt=zh-cn [2] Marcellin M W, Gormish M J, Bilgin A, Boliek M P. An overview of JPEG-2000. In:Proceedings of the 2000 Data Compression Conference. Snowbird, UT:IEEE, 2000.523-541 [3] Bradley J N, Brislawn C M, Hopper T. FBI wavelet/scalar quantization standard for gray-scale fingerprint image compression. In:Proceedings of the SPIE 1961, Visual Information Processing Ⅱ. Orlando, FL:SPIE, 1993.293-304 [4] Skodras A, Christopoulos C, Ebrahimi T. The JPEG 2000 still image compression standard. IEEE Signal Processing Magazine, 2001, 18(5):36-58 doi: 10.1109/79.952804 [5] Shao G, Wu Y, Yong A, Liu X, Guo T. Fingerprint compression based on sparse representation. IEEE Transactions on Image Processing, 2014, 23(2):489-501 doi: 10.1109/TIP.2013.2287996 [6] Olshausen B A, Field D J. Sparse coding with an overcomplete basis set:a strategy employed by V1? Vision Research, 1997, 37(23):3311-3325 doi: 10.1016/S0042-6989(97)00169-7 [7] Emmanuel B, Mu'Azu M, Sani S, Garba S. A review of wavelet-based image processing methods for fingerprint compression in biometric application. British Journal of Mathematics and Computer Science, 2014, 4(19):2781-2798 doi: 10.9734/BJMCS [8] Qian C, Xu Z. Robust visual tracking via sparse representation under subclass discriminant constraint. IEEE Transactions on Circuits and Systems for Video Technology, 2016, 26(7):1293-1307 doi: 10.1109/TCSVT.2015.2424091 [9] Sun B, Liu Z, Sun Y, Su F, Cao L, Zhang H. Multiple objects tracking and identification based on sparse representation in surveillance video. In:Proceedings of the 2015 IEEE International Conference on Multimedia Big Data (BigMM). Beijing, China:IEEE, 2015.268-271 [10] Cheng M, Wang C, Li J. Single-image super-resolution in RGB space via group sparse representation. Iet Image Processing, 2015, 9(6):461-467 doi: 10.1049/iet-ipr.2014.0313 [11] Tropp J A. Greed is good:algorithmic results for sparse approximation. IEEE Transactions on Information Theory, 2004, 50(10):2231-2242 doi: 10.1109/TIT.2004.834793 [12] Mallat S G, Zhang Z F. Matching pursuits with time-frequency dictionaries. IEEE Transactions on Signal Processing, 1993, 41(12):3397-3415 doi: 10.1109/78.258082 [13] Pati Y C, Rezaiifar R, Krishnaprasad P S. Orthogonal matching pursuit:recursive function approximation with applications to wavelet decomposition. In:Proceedings of the 1993 Conference Record of the 27th Asilomar Conference on Signals, Systems, and Computers. Pacific Grove, CA:IEEE, 1993.40-44 http://cn.bing.com/academic/profile?id=185716565&encoded=0&v=paper_preview&mkt=zh-cn [14] Gharavi-Alkhansari M, Huang T S. A fast orthogonal matching pursuit algorithm. In:Proceedings of the 1998 IEEE International Conference on Acoustics, Speech, and Signal Processing. Seattle, WA:IEEE, 1998.1389-1392 http://dl.acm.org/citation.cfm?id=1892575 [15] Chen S S, Donoho D L, Saunders M A. Atomic decomposition by basis pursuit. SIAM Journal on Scientific Computing, 1998, 20(1):33-61 doi: 10.1137/S1064827596304010 [16] Zhu J Y, Wang Z Y, Zhong R, Qu S M. Dictionary based surveillance image compression. Journal of Visual Communication and Image Representation, 2015, 31:225-230 doi: 10.1016/j.jvcir.2015.07.002 [17] Setiawan A D, Suksmono A B, Mengko T L R, Gunawan H. Low-bitrate medical image compression. In:Proceedings of the ACA2011 IAPR Conference on Machine Vision Applications. Nara, Japan, 2011.544-547 [18] Xu J, Pi Y, Ming R. SAR image compression based on sparse representation. In:Proceedings of the 11th International Radar Symposium (IRS). Vilnius, Lithuania:IEEE, 2010.1-4 [19] Zhan X, Zhang R, Yin D, Huo C. SAR image compression using multiscale dictionary learning and sparse representation. IEEE Geoscience and Remote Sensing Letters, 2013, 10(5):1090-1094 doi: 10.1109/LGRS.2012.2230394 [20] Aharon M, Elad M, Bruckstein A. K-SVD:an algorithm for designing overcomplete dictionaries for sparse representation. IEEE Transactions on Signal Processing, 2006, 54(11):4311-4322 doi: 10.1109/TSP.2006.881199 [21] Anurakphanawan N, Lamsrichan P. Fingerprint recognition performance with WSQ, CAWDR, and JPEG2000 compression. In:Proceedings of the 6th International Conference of Information and Communication Technology for Embedded Systems (IC-ICTES). Hua-Hin, Thailand:IEEE, 2015.1-6 [22] Fingerprint images[Online], available:http://pan.baidu.com/s/1i3KxQZV,November27,2015 [23] Abraham J, Kwan P, Gao J B. Fingerprint matching using a hybrid shape and orientation descriptor. State of the Art in Biometrics. New York:InTech, 2011. 期刊类型引用(2)
1. 张跃中,肖敏,王璐,徐丰羽. 大规模超环神经网络分岔动力学. 自动化学报. 2022(04): 1129-1136 . 本站查看
2. 郭磊,吕欢,黄凤荣,石洪溢. 基于突触可塑性的无标度脉冲神经网络的动态特性研究. 生物医学工程学杂志. 2019(06): 902-910 . 百度学术
其他类型引用(2)
-