«上一篇
文章快速检索     高级检索
下一篇»
  智能系统学报  2018, Vol. 13 Issue (5): 707-715  DOI: 10.11992/tis.201704023
0

引用本文  

吕红丽, 宋玉晶, 段培永. 非线性布尔网络系统模糊建模与动态性能分析[J]. 智能系统学报, 2018, 13(5), 707-715. DOI: 10.11992/tis.201704023.
LYU Hongli, SONG Yujing, DUAN Peiyong. Fuzzy modeling and dynamic analysis of nonlinear Boolean networks systems[J]. CAAI Transactions on Intelligent Systems, 2018, 13(5), 707-715. DOI: 10.11992/tis.201704023.

基金项目

国家自然科学基金项目(61374187, 61403237).

通信作者

宋玉晶. E-mail:xiaosongyj@qq.com

作者简介

吕红丽,女,1978年生,副教授,主要研究方向为复杂系统建模控制与仿真、智能环境与网络化控制。发表学术论文20余篇;
宋玉晶,女,1990年生,硕士研究生,主要研究方向为智能环境与网络化控制;
段培永,男,1968年生,教授,博士生导师,主要研究方向为模糊控制、神经控制、数据挖掘、建筑智能环境与网络化控制研究。主持科研项目10余项,发表学术论文70余篇,被SCI、EI收录30余篇

文章历史

收稿日期:2017-04-19
网络出版日期:2018-04-25
非线性布尔网络系统模糊建模与动态性能分析
吕红丽1, 宋玉晶1, 段培永2    
1. 山东建筑大学 信息与电气工程学院,山东 济南 250101;
2. 山东师范大学 信息科学与工程学院,山东 济南 250014
摘要:针对非线性系统难以精确建模与动态性能分析的基本控制问题,基于模糊动态模型把布尔网络系统理论推广到非线性布尔网络系统,建立了模糊动态布尔网络控制系统的模型。引入模糊动态模型,对非线性布尔网络进行模糊建模,分别建立了非线性布尔网络系统的局部模型和全局模型。从系统的局部意义和全局意义上,对系统进行了能控性、能观性、稳定性等动态性能分析。最后,以多输入多输出的非线性布尔网络系统实例为具体研究对象,建立了系统的局部模型和全局模型,并对动态性能进行了仿真分析,得到了实验结果。实验结果表明,模糊动态布尔网络控制系统对非线性布尔网络系统的建模是有效的,动态性能分析是合理的,对模糊动态布尔网络控制系统的进一步分析有重要意义。
关键词模糊动态模型    布尔网络    半张量积    局部模型    全局模型    能控性    能观性    稳定性    
Fuzzy modeling and dynamic analysis of nonlinear Boolean networks systems
LYU Hongli1, SONG Yujing1, DUAN Peiyong2    
1. School of Information and Electrical Engineering, Shandong Jianzhu University, Ji’nan 250101, China;
2. School of Information Science and Engineering, Shandong Normal University, Ji’nan 250014, China
Abstract: Considering the difficulty in accurately modeling nonlinear systems and analyzing their dynamic properties, the Boolean network system theories are extended to the nonlinear Boolean network system based on a fuzzy dynamic model, establishing a model of fuzzy dynamic Boolean network control systems. The fuzzy dynamic model is introduced to build a fuzzy model of nonlinear Boolean network, establishing the local and global models of the nonlinear Boolean network systems. The dynamic properties, which include controllability, observability, and stability of the system, are analyzed from the local and global meanings of the system. Finally, a multi-input multi-output nonlinear Boolean network system is taken as a numerical example, and the local and global models of the system are established. The dynamic properties are simulated and analyzed, and the experimental results are obtained. The results show that the fuzzy dynamic Boolean network control system is effective in modeling nonlinear Boolean network systems and reasonably analyzes the dynamic properties, which is of great significance for further analysis of the fuzzy dynamic Boolean network control systems.
Key words: fuzzy dynamical model    Boolean network    semi-tensor product    local model    global model    controllability    observability    stability    

1965年,Zadeh[1]首次提出基于模糊集合的模糊数学概念。1974年,Mamdani[2]首先将模糊数学理论应用于工业控制中,实现了对锅炉和蒸汽机的模糊控制,使模糊控制从理论走向实际应用。1985年,日本学者Takagi等[3]又提出了以线性精确数学表达式为模糊规则后件的T-S模糊模型,将模糊控制系统与线性控制系统有效地结合起来。随后国内外学者基于T-S模糊模型,研究了非线性系统的大量控制问题,得到了丰富的模糊控制理论研究成果[4-9]。在T-S模糊模型的基础上,Feng等[5-6]提出了模糊动态模型(fuzzy dynamical model)。其主要思想是,构造一组线性模型,分别描述系统的动态特性,通过局部线性模型的加权组合得到系统的全局非线性动态模型。文献[10]中提出了一种离散时间模糊控制系统,将非线性离散时间系统作为模糊控制系统的后件部分,分析了系统的动态性能。

1969年,Kauffman[11]首先提出布尔网络模型,布尔网络是关于布尔状态变量的一种简单的逻辑动力系统,是当前学者专家们共同关心的热点问题。针对布尔网络研究缺少有效的数学工具问题,程代展教授在文献[12]中首次提出矩阵半张量积方法。这种方法将逻辑运算转换成代数运算,使得许多经典的处理量变过程的数学工具可直接用来分析逻辑动态系统。在文献[13]中,程代展教授将这种方法应用于布尔网络,将逻辑动态控制系统转化为普通离散时间系统,提出了一系列关于布尔网络的新理论。随后在文献[14-17]中研究了布尔网络系统的能控能观性等性质,形成了布尔控制网络分析设计的完整理论框架。之后,学者们在控制理论方面对线性布尔网络系统做了大量的深入研究[18-23],但是没有针对非线性布尔网络系统进行分析和研究。

为了解决非线性布尔网络的数学建模与分析问题,利用模糊动态模型的非线性特点,将模糊动态模型和布尔控制网络相结合,建立了模糊动态布尔网络系统的局部模型和全局模型,并且分别分析了系统局部模型和全局模型的能控性、能观性和稳定性。

1 预备知识 1.1 数学符号说明

为了叙述方便,文中用到的记号列表如下: ${{{M}}_{m \times n}}$ 表示所有 $m \times n$ 矩阵的集合, $ \otimes $ 表示矩阵的张量积; ${\rm Col}({{A}})({\rm Row}({{A}}))$ 为矩阵 ${{A}}$ 的列(行)集合; ${\rm Col}({{{A}}_i})({\rm Row}({{{A}}_i}))$ 为矩阵 ${{A}}$ 的第 $i$ 列(行);记 ${ \delta}_n^i$ 是单位矩阵 ${{{I}}_{{n}}}$ 的第 $i$ 列; ${{ {\varDelta}} _n} = \left\{ {{ \delta} _n^i|i = 1,2, \cdots ,n} \right\}$ ${{\cal{D}}_k} = \left\{ {0,1, \cdots } \right.,$ $k - 1\},$ $k \geqslant 2$ 。记逻辑变量:真 $ \sim T - 1$ ,假 $ \sim F - 0$ ,则 ${\cal{D}} = \left\{ {0,1} \right\}$ $f:{{\cal{D}}^n} \to \cal{D}$ 称为逻辑函数; ${{L}} \in {{\cal{M}}_{n \times r}}$ ,称 ${{L}}$ 为逻辑矩阵,如果 ${\rm Col}({{L}}) \subset {{ {\varDelta}} _m},$ $m \times r$ 维逻辑矩阵全体记为 ${{{L}}_{m \times r}}$ ;设矩阵 ${{L}} \in {{\cal{M}}_{n \times r}}$ 。其中, ${\rm Col}({{L}}) \subset {{ {\varDelta}} _n},$ ${{L}}$ 为逻辑矩阵,可简记为 ${{L}} = {\delta _n}\left[\!\!\!\!\! {\begin{array}{*{20}{c}} {{i_1}}\!\!&\!\!{{i_2}}\!\!&\!\! \cdots \!\!&\!\!{{i_r}} \end{array}} \!\!\!\!\! \right]$ ;矩阵 ${{{ B}}} \in {{\cal{M}}_{m \times n}}$ 为布尔矩阵,如果B $({b_{i,j}}) \in {\cal D},$ $m \times n$ 维布尔矩阵全体记为 ${{\cal{B}}_{m \times n}}$

1.2 模糊动态模型

模糊动态模型是非线性复杂系统模糊建模中一种典型的模糊建模方法。模糊动态模型每条规则的后件部分是一种状态空间形式的局部线性系统,因此对于局部的线性模型可以采用线性系统的理论体系去研究,然后通过模糊推理得到全局意义下的模糊控制系统的分析和设计。

一个m个输入n个输出的非线性系统,其模糊s控制器可表示为 $\varSigma \in f({y_1} \times {y_2}\cdots \times {y_n} \times {x_1} \times {x_2}\cdots \times {x_m})$ ,其中 $\{ {x_i}\} $ 是模糊控制器的输入,论域为 ${E_{{x_i}}}$ , $\{ {y_i}\} $ 是模糊控制器的输出,论域为 ${E_{{y_j}}}$ 。将模糊变量 ${x_i},{\kern 1pt} {\kern 1pt} {y_j}$ 按照隶属度函数进行模糊化, ${E_{{x_i}}} = \{ x_1^i,x_2^i, \cdots ,x_{{\alpha _i}}^i\} ,i = 1,2, \cdots ,m$ ,分别对应于“负大”“负中” $ \cdots $ 是基于隶属度的模糊集合; ${E_{{y_j}}} = \{ y_1^j,y_2^j, \cdots ,y_{{\beta _j}}^j\} ,$ $j = 1,2, \cdots ,n,$ 分别对应于“负大”“负中” $ \cdots $ 是基于隶属度的模糊集合。

假设总共有 $N$ 条模糊规则,则第 $k$ 条模糊规则为 ${R_k}$ , $k = 1,2, \cdots ,N$ ,其模糊规则为

$\begin{array}{c}{R_k}:\;{\rm{if}}\;{x_1}\;{\rm{is}}\;{A_{1k}}, {x_2}\;{\rm{is}}\;{A_{2k}} ,\cdots ,{x_m}\;{\rm{is}}\;{A_{mk}}\\{\rm{Then}}\left\{ \begin{array}{c}{y_{1k}} = p_{0k}^1 + p_{1k}^1{x_1} + \cdots + p_{ik}^1{x_i} + \cdots p_{mk}^1{x_m}\\ \vdots \\{y_{jk}} = p_{0k}^j + p_{1k}^j{x_1} + \cdots + p_{ik}^j{x_i} + \cdots p_{mk}^j{x_m}\\ \vdots \\{y_{nk}} = p_{0k}^n + p_{1k}^n{x_1} + \cdots + p_{ik}^n{x_i} + \cdots p_{mk}^n{x_m}\end{array} \right.\end{array}$ (1)

式中: ${A_{ik}}$ 是一个模糊集合; ${y_{jk}}$ 是第 $k$ 条模糊规则的第 $j$ 个输出, $j = 1,{\kern 1pt} 2, \cdots ,{\kern 1pt} n$ $p_{ik}^j$ 是第 $k$ 条模糊规则结论中第 $j$ 个输出的线性多项式函数中变量 ${x_i}$ 项的系数,一般为常数项,特别的 $p_{0k}^j$ 通常可以归一化为1。

对于第 $k$ 条规则 ${R_k}$ ,如果已知输入 ${x_1} = x_1^*,$ ${x_2} = x_2^*,{\kern 1pt} \cdots ,{\kern 1pt} {x_m} = x_m^*,{\kern 1pt} $ 则在结论部分的输出 ${y_{jk}}$ 可以由线性多项式函数计算得到:

$y_{jk}^* = p_{1k}^jx_1^* + p_{2k}^jx_2^* + \cdots + p_{mk}^jx_m^*$ (2)

每条规则的激活度 ${\mu _i}$

$\left\{ \begin{gathered} {\mu _1} = {\mu _{{A_{11}}}}(x_{_1}^*) \wedge {\mu _{{A_{21}}}}(x_{_2}^*) \wedge \cdots \wedge {\mu _{{A_{m1}}}}(x_{_m}^*) \\ \vdots \\ {\mu _k} = {\mu _{{A_{1k}}}}(x_{_1}^*) \wedge {\mu _{{A_{2k}}}}(x_{_2}^*) \wedge \cdots \wedge {\mu _{{A_{mk}}}}(x_{_m}^*) \\ \vdots \\ {\mu _N} = {\mu _{{A_{1N}}}}(x_{_1}^*) \wedge {\mu _{{A_{2N}}}}(x_{_2}^*) \wedge \cdots \wedge {\mu _{{A_{mN}}}}(x_{_m}^*) \\ \end{gathered} \right.$ (3)

式中: ${\mu _{{A_{ik}}}}(x_i^*)$ 表示论域中第 $i$ 个元素对 ${A_{ik}}$ 的隶属度, $ \wedge $ 是取小运算。

模糊动态模型的输出 ${y_j}$ 是由所有 $k$ 条规则( $k = 1,{\kern 1pt} 2, \cdots ,{\kern 1pt} N$ )的输出 $y_{jk}^*$ 加权平均得到的。模型的输出为

${Y_j} = \frac{{\sum\limits_{k = 1}^N {{\mu _k}} \cdot y_{_{jk}}^*}}{{\sum\limits_{k = 1}^N {{\mu _k}} }}$ (4)

一般取 $\sum\limits_{k = 1}^N {{\mu _k} = 1} $ ,则式(4)全局模型输出简化为 ${Y_j} = \sum\limits_{k = 1}^N {{\mu _k} \cdot y_{_{jk}}^*} $

1.3 逻辑的矩阵表示

矩阵的半张量积是中科院系统所程代展教授在文献[12]中提出的一种新的矩阵乘法,即设 ${{A}} \in {{\cal{M}}_{m \times n}},\;{{{ B}}} \in {{\cal{M}}_{p \times q}}$

1)如果 $n = p$ ,则称AB满足等维数关系;

2)果 $n = tp$ (记为 ${{A}}{ \succ _t}{{B}}$ ),或者 $nt = p$ (记 ${{A}}{ \prec _t}{{B}}$ ),则称 ${{A}}$ ${{B}}$ 满足倍维数关系,否则称一般维数关系。矩阵乘积在倍维数关系下的一种推广如定义1。

定义1[13]  设 ${{X}}$ $n = pq$ 维行向量, ${{Y}}$ $p$ 维列向量。将X等分成 ${{X}} = ({{{X}}^1},\;{{{X}}^2},\; \cdots ,\;{\kern 1pt} {{{X}}^p})$ ,这里 ${{{X}}^i} \in {{\bf R}^q},$ $i = 1,2, \cdots ,p$ 。那么, ${{X}}$ ${{Y}}$ 的半张量积记作 ${{X}} \triangleright {{Y}}$ ,定义为一个行向量:

${{X}} \triangleright {{Y}} = \sum\limits_{i = 1}^p {{{{X}}^i}{{{y}}_i}} \in {{\bf R}^q}$ (5)

类似的,

${{{Y}}^{\rm T}} \triangleright {{{X}}^{\rm T}} = \sum\limits_{i = 1}^p {{{{y}}_i}({{\bf X}^i}} {)^{\rm T}} \in {{\bf R}^q}$ (6)

为一列向量。

普通矩阵乘法是半张量积的特殊形式,普通矩阵乘法具有的性质,对于半张量积几乎都成立,此外还具有一些特有的性质。

定理1[13]  设 ${{f}}({x_1},{x_2}, \cdots ,{x_n})$ 为一个逻辑函数,在向量形式下 $f:{\varDelta _{{2^n}}} \to \varDelta $ 则存在唯一的逻辑矩阵 ${{{M}}_f} \in {{\cal{L}}_{2 \times {2^n}}}$ ,称为 $f$ 的结构矩阵,使得

${{f}}({{{x}}_1},{{{x}}_2}, \cdots ,{{{x}}_n}) = {{{M}}_f} \triangleright {{x}}$ (7)
$式中 ${{x}} = \triangleright _{i = 1}^n{{{x}}_i}$ 。常用的逻辑算子及其结构矩阵分别为 $\begin{array}{c}{{{M}}_ \wedge } = {\delta _2}[1\;\;2\;\;2\;\;2],\quad{{{M}}_ \vee } = {\delta _2}[1\;\;1\;\;1\;\;2],\\{{{M}}_ \to } = {\delta _2}[1\;\;2\;\;1\;\;1],\quad{{{M}}_{\overline \vee }} = {\delta _2}[2\;\;1\;\;1\;\;2].\end{array}$$
1.4 布尔(控制)网络

定义2[13]  1)布尔网络的动态方程为

$\left\{ \begin{gathered} {x_1}(t + 1) = {f_1}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t)) \\ {x_2}(t + 1) = {f_2}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t)) \\ \vdots \\ {x_n}(t + 1) = {f_n}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t)) \\ \end{gathered} \right.{\kern 1pt} $ (8)

式中: ${f_i}:{{\cal{D}}^{n+m}} \to {\cal{D}},i=1,2, \cdots ,n;$ 为逻辑函数; ${x_i}(t) \in $ $ {\cal{D}},$ $i = 1,$ $2, \cdots ,n$ 为状态变量。

2)布尔控制网络是指一个含有输入输出的布尔网络,其动态方程为

$\left\{ \begin{gathered} {x_1}(t + 1) = {f_1}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t),{u_1}(t),\; \cdots ,{u_m}(t)) \\ {x_2}(t + 1) = {f_2}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t),{u_1}(t),\; \cdots ,{u_m}(t)) \\ \vdots \\ {x_n}(t + 1) = {f_n}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t),{u_1}(t),\; \cdots ,{u_m}(t)) \\\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{y_j}(t) = {h_j}({x_1}(t),{x_2}(t), \cdots ,{x_n}(t) \\ \end{gathered} \right.$ (9)

式中: ${x_i}(t) \in \cal{D},$ $i = 1,2, \cdots ,n$ 为状态变量; ${u_i}(t) \in $ $ {\cal{D}},i =1,2,\cdots ,m$ 为控制变量; ${y_i}(t) \in {\cal{D}},$ $i = 1,2, \cdots ,p$ 为输出变量; ${f_i}:{{\cal{D}}^{n + m}} \to {\cal{D}},i = 1,2, \cdots ,n;\;{\kern 1pt} {h_i}:{{\cal{D}}^n} \to {\cal{D}},i = $ $ 1,2, \cdots ,p$ 为逻辑函数。

在向量表达式下令 ${{x}} = \triangleright _{i = 1}^n{{{x}}_i}{\text{、}}$ ${{u}} = \triangleright _{i = 1}^m{{{u}}_i}$ ${{y}} = \triangleright _{i = 1}^p{{{y}}_i}$ ,利用定理1可知:

定理2[13]  利用向量表达式

1)布尔网络的动态方程式(8)可表示为

${{X}}(t + 1) = {{LX}}(t),\;\;\;\;{{L}} \in {{\cal{L}}_{{2^n} \times {2^n}}}$ (10)

2)布尔控制网络的动态方程式(9)可表示为

$\left\{ \begin{aligned}& {{X}}(t + 1) = {{LU}}(t){{X}}(t) \\&{{Y}}(t) = {{HX}}(t) \\ \end{aligned} \right.$ (11)

式中: ${{L}} \in {{\cal{L}}_{{2^n} \times {2^n}}},$ ${{H}} \in {{\cal{L}}_{{2^p} \times {2^n}}}$ , ${{L}}{\text{、}}\!\!\!{{H}}$ 称为结构矩阵。

2 非线性布尔网络系统的基本概念

定义3 1)一个多输入多输出的非线性布尔网络系统可以表示成模糊动态布尔网络模型(fuzzy dynamic Boolean network model,FDBNM),即

$\begin{aligned}& {R^k}:\;{\rm{if}}\;{z_1}\;{\rm{is}}\;F_1^k\;{\rm{and}}\; \cdots \;{z_i}\;{\rm{is}}\;F_i^k\; \cdots \;{z_n}\;{\rm{is}}\;F_n^k\\& {\rm{Then}}\;\left\{ \begin{aligned}& {{X}}(t + 1) = {{{L}}_k}{{U}}(t){{X}}(t)\\& {{Y}}(t) = {{{H}}_k}{{X}}(t)\end{aligned} \right.,\quad k = 1,\;2, \cdots ,{\kern 1pt} \;n\end{aligned}$ (12)

简记为

$\begin{aligned}& {R^k}:\;{\rm{if}}\;{z_1}\;{\rm{is}}\;F_1^k\;{\rm{and}}\; \cdots \;{z_i}\;{\rm{is}}\;F_i^k\; \cdots \;{z_n}\;{\rm{is}}\;F_n^k\\& {\rm{Then}}\;\;{\rm{L}}{{\rm{M}}_k} = [\;{\mu _k}(z),{\kern 1pt} ({{{L}}_k},\;{{{H}}_k})],\quad k = 1,\;2, \cdots \;,\;{\kern 1pt} n\end{aligned}$ (13)

2)第 $k$ 个FDBNM为

${\rm{FDBNML}}{{\rm{M}}_k} = [{\;\mu _k}(z),({{{L}}_k},{\kern 1pt} {\kern 1pt} {{{H}}_k})]$ (14)

式中: ${\mu _k}\left( z \right)$ 是模糊推理集合 ${F^k}$ 的隶属函数; $ {F^k} =$ $\prod _{i = 1}^nF_i^k,\;{\mu _k}\left( z \right) = {\mu _{F_1^k}}(z_{_1}^{}) \wedge {\mu _{F_2^k}}(z_{_2}^{}) \;\wedge\; \cdots \;\wedge\; {\mu _{F_n^k}}(z_{_m}^{}),$ $\displaystyle\sum\limits_{k = 1}^N {{\mu _k}\left( z \right)} $ $({{{L}}_k},{{{H}}_k})$ 为局部模型的结构矩阵, $({{{L}}_k},{{{{ {{H}}}}}_k})$ 也称为FDBNM的第 $k$ 个子系统。

这里, ${R_k},k = 1,2, \cdots ,N$ 为系统的第 $k$ 条模糊规则,也称为第 $k$ 个模糊子系统,N为总的模糊规则数; ${z_1},{z_2}, \cdots {z_n}$ 为规则前件语言变量; ${F_{i}^k}$ 为模糊集,其隶属度函数设为三角形函数,记作 ${\mu _{{F_{i}^k}}}$ $({{X}}({{t + 1}}),{{Y}}(t))$ 是系统的输出; ${{U}}(t)$ 是输出部分布尔控制网络的控制变量; $({{{L}}_k},{{{{ H}}}_k})$ 为系统的结构矩阵。

对给定的输入信号 ${z_1} = z_1^*,{\kern 1pt} {\kern 1pt} {z_2} = z_2^*,{\kern 1pt} \cdots ,{\kern 1pt} {z_n} = z_n^*,{\kern 1pt} $ 利用三角形隶属函数将其模糊化,对每条规则的激活度 ${\mu _k}(z)$ 采用max-min方法:

$\left\{ \begin{gathered} {\mu _1}\left( z \right) = {\mu _{F_1^1}}(z_{_1}^*) \wedge {\mu _{F_2^1}}(z_{_2}^*) \wedge \cdots \wedge {\mu _{F_n^1}}(z_{_m}^*) \\ \vdots \\ {\mu _k}\left( z \right) = {\mu _{F_1^k}}(z_{_1}^*) \wedge {\mu _{F_2^k}}(z_{_2}^*) \wedge \cdots \wedge {\mu _{F_n^k}}(z_{_m}^*) \\ \vdots \\ {\mu _N}\left( z \right) = {\mu _{F_1^N}}(z_{_1}^*) \wedge {\mu _{F_2^N}}(z_{_2}^*) \wedge \cdots \wedge {\mu _{F_n^N}}(z_{_m}^*) \\ \end{gathered} \right.$

3)使用加权平均法解模糊,可得FDBNM的全局模型为

$\left\{ \begin{aligned} & {{X}}(t + 1) = {{LU}}(t){{X}}(t) \\ &{{Y}}(t) = {{HX}}(t) \end{aligned} \right.$ (15)

简记为

$GM = [\;\mu (z),({{L}},{\kern 1pt} {\kern 1pt} {{H}})]$

式中: ${{L}} = \displaystyle\sum\limits_{k = 1}^N {{\mu _k}(z){{{L}}_k}} ,{\kern 1pt} {\kern 1pt} {{H}} = \sum\limits_{k = 1}^N {{\mu _k}(z){{{H}}_k}} ,{{L}},{{H}}$ 为全局模型的结构矩阵。 $ \mu (z) = {\mu _1}(z) \vee \cdots $ $\vee {\mu _k}(z) \vee \cdots \vee {\mu _N}(z)$

定义4 1)对于模糊动态布尔网络控制系统的局部模型式(14), $G$ 固定,如果存在控制变量 ${U_0}$ ,能使方程式(14)从初始状态 $X(U,0) = {X_0}$ 到达终端状态 $X(U,t) = {X_d}$ ,则称 ${X_d}$ ${X_0}$ 经过 $t$ 步是能控的。如果存在控制变量 $U$ ,使式(14)能从任意初始状态 ${X_0}$ 到达 $X(T) = {X_d}$ ,则称模糊动态布尔网络控制系统的局部模型是能控的。

2)对于模糊动态布尔网络控制系统的全局模型式(15),如果存在控制变量 ${U_0}$ ,能使方程式(15)从初始状态 $X(U,0) = {X_0}$ 到达终端状态 $ X(U,t) =$ $ {X_d}$ ,则称 ${X_d}$ ${X_0}$ 经过 $t$ 步是能控的。如果式(15)能从任意初始状态 ${X_0}$ 到达 $X(T) = {X_d}$ ,则称模糊动态布尔网络控制系统的全局模型是能控的。

定义5 1)对于模糊动态布尔网络控制系统的局部模型式(14),对任意给定的初始状态如果至少存在一个布尔控制序列,使初始状态能由输出序列唯一地确定,则称局部模型是状态能观测的。

2)对于模糊动态布尔网络控制系统的全局模型式(15),对任意给定的初始状态如果至少存在一个布尔控制序列,使初始状态能由输出序列唯一地确定,则称全部模型是状态能观测的。

定义6 1)对于模糊动态布尔网络控制系统的局部模型式(14),如果经过固定步数 ${T_0}$ ,存在一个不动点 ${X_e}$ ,使得对于局部模型的任意的初始状态 $X(0) = \left( {x_1^0,x_2^0, \cdots ,x_n^0} \right)$ ,都有 $X(t) = X(e),t \geqslant {T_0}$ ,则称系统的局部模型是能稳定的。

2)对于模糊动态布尔网络控制系统的全局模型式(15),如果经过固定步数 ${T_0}$ ,存在一个不动点 ${X_e}$ ,使得对于全局模型的任意的初始状态 $X(0) = $ $ \left( {x_1^0, x_2^0,\cdots ,x_n^0} \right)$ ,都有 $X(t) = X(e),t \geqslant {T_0}$ , 则称系统的全局模型是能稳定的。

3 非线性布尔网络系统动态性能分析

系统局部模型的布尔控制网络系统(14)的结构矩阵可以等分为 ${2^m}$ [13],即

$\begin{aligned}{{{L}}_k} = & [{\rm{Bl}}{{\rm{k}}_1}({{{L}}_k}){\rm{Bl}}{{\rm{k}}_2}({{{L}}_k}) \cdots {\rm{Bl}}{{\rm{k}}_{{2^m}}}({{{L}}_k})] =\\& [{B_{k1}},{B_{k2}}, \cdots ,{B_{k{2^m}}}]\end{aligned}$ (16)

式中: ${\rm{Bl}}{{\rm{k}}_i}({{{L}}_k})$ 是矩阵 ${{{L}}_k}$ 的第 $i$ $n \times n$ 的块; ${B_i} \in {{\cal{L}}_{{2^n} \times {2^n}}}$ $i = 1,2, \cdots ,{2^m}$ 。令

${{{M}}_k} = \sum\limits_{i = 1}^{{2^m}} {{\rm{Bl}}{{\rm{k}}_i}} ({{{L}}_k})$ (17)
3.1 能控性

首先引用布尔矩阵的布尔乘法及布尔幂的定义[13]

1)设 $\alpha ,\beta ,{\alpha _i} \in {\cal{D}},i = 1,2, \cdots ,n$ ,则布尔加法定义为

$\left\{ \begin{aligned} & \alpha { + _{\cal{B}}}\beta = \alpha \vee \beta \\& \sum\limits_{i = {1_{\cal{B}}}}^n {{\alpha _i} = {\alpha _1} \vee {\alpha _2}\vee \cdots \vee {\alpha _n}} \end{aligned} \right.$ (18)

2)设 ${{A}} = ({a_{ij}}) \in {{\cal{B}}_{m \times n}},{{B}} = ({b_{ij}}) \in {{\cal{B}}_{n \times p}}$ ,则布尔乘法定义为

${{A}}{ \triangleright _{\cal{B}}}{{B}} = {{C}} \in {{\cal{B}}_{m \times p}}$ (19)

式中: ${c_{ij}} = \displaystyle\sum\limits_{k = {1_{\cal{B}}}}^n {{a_{ik}}{b_{kj}}} ,i = 1,2, \cdots ,m,j = 1,2, \cdots, p$

3)设 ${{A}} \triangleright {{A}}$ 有定义,则布尔幂定义为

${{{A}}^{(k)}} = \underbrace {{{A}}{ \triangleright _{\cal{B}}}{{A}}{ \triangleright _{\cal{B}}} \cdots { \triangleright _{\cal{B}}}{{A}}}_k$ (20)

对于模糊动态布尔网络控制系统的局部模型式(14),定义能控性矩阵为[13]

${{{C}}_k} = \sum\limits_{j = 1}^{{2^{(m + n)}}} {_{\cal{B}}\sum\limits_{i = 1}^{{2^m}} {_{\cal{B}}} } {\rm{Bl}}{{\rm{i}}_i}({{L}}_{_k}^{(j)}) = \sum\limits_{j = 1}^{{2^{(m + n)}}} {_{\cal{B}}} M_{_k}^{(j)} \in {{\cal{B}}_{{2^m} \times {2^n}}}$ (21)

定理3 1)当且仅当矩阵 ${{{C}}_k} > 0$ 时,局部模型式(14)是能控的。

2)当且仅当矩阵 ${{C}} > 0$ 时,全局模型式(15)是能控的,即 ${{C}} = \displaystyle\sum\limits_{k = 1}^N {{\mu _k}(z){{{C}}_k}} $

证明 通过数学归纳法来证明。

$j = 1$ 时,由式(17)知,当 ${M_k} > 0$ 时,存在一个控制序列使状态 ${X_0}$ ${X_d}$ ,显然局部模型式(14)是能控的;假设当 $j = k'$ 时, ${{{C}}_k} = \displaystyle\sum\limits_{k' = 1}^{{2^{m + n}}} {M_k^{(k')}} > 0$ ,式(14)能控,则当 $j = k' + 1$ 时, ${{{C}}_k} = \displaystyle\sum\limits_{k' = 0}^{{2^{m + n}}} {M_k^{(k' + !)}} = \displaystyle\sum\limits_{k' = 1}^{{2^{m + n}}} {M_k^{(k')}} $ $ + {M_k}$ ,因为 $\displaystyle\sum\limits_{k' = 1}^{{2^{m + n}}} {M_k^{(k')}} > 0,$ ${M_k} > 0,$ 可知 ${{{C}}_k}>0$ ,且存在控制序列使局部模型式(14)能控。

同理可证,当矩阵 ${{C}} > 0$ 时,全局模型式(15)是能控的。

3.2 能观性

能观性所表示的是输出 $y(t)$ 反映状态变量 $x(t)$ 的能力,与控制作用没有直接关系。

为了找到能观性矩阵,对于第k个局部模型,定义一组矩阵集合[13]: ${{{\varOmega}} _{ki}} \in {{\cal{L}}_{{2^p} \times {2^n}}},\; i = 0,\;1,\; 2 \cdots $ 其表达式为

$\left\{ \begin{aligned} & {{{\varOmega}} _{k0}} = \{ {{{H}}_k}\} \\ & {{{\varOmega}} _{k1}} = \{ {{{H}}_k}{{{B}}_i}\left| {i = 1,2, \cdots ,{2^m}} \right.\} \\ & \quad\quad\quad \vdots \\ & {{{\varOmega}} _{ks}} = \{ {{{H}}_k}{{{B}}_{{i_1}}}{{{B}}_{{i_2}}} \cdots {{{B}}_{{i_s}}}\left| {{i_1},{i_2}, \cdots ,{i_s} = 1,2, \cdots ,{2^m}} \right.\} \end{aligned} \right.$ (22)

式中 ${{{\varOmega}} _{k\left( {s + 1} \right)}} \subset \bigcup\limits_{i = 1}^s {{{{\varOmega}} _i}} $ ,记 ${s^*}$ 为使之成立的最小正整数。从 ${{{\varOmega}} _{ki}}$ 中选取矩阵 ${{{\varGamma}} _{ki}}$

$\left\{ {\begin{aligned}& {{{{\varGamma}} _{k0}} = {{{H}}_k}}\\& {{{{\varGamma}} _{k1}} = {{\left[ {{{{H}}_k}{{{B}}_1}\;\;{{{H}}_k}{{{B}}_2} \cdots {{{H}}_k}{{{B}}_{{2^m}}}} \right]}^{\rm{T}}}}\\& {{{{\varGamma}} _{k2}} = {{\left[ {{{{H}}_k}{{{B}}_1}\;\;{{{B}}_1}{{{H}}_k}{{{B}}_1}{{{B}}_2} \cdots {{{H}}_k}{{{B}}_{{2^m}}}{{{B}}_{{2^m}}}} \right]}^{\rm{T}}}}\end{aligned}} \right.$ (23)

${{{O}}_k} = \left[ {{{{\varGamma}} _{k0}}\;} \right.{{{\varGamma}} _{k1}}\; \cdots \;{\left. {{{{\varGamma}} _{ks*}}} \right]^{\rm T}}$ (24)

式中 ${{{O}}_k}$ 称为系统局部模型式(14)的能观测矩阵。

定理4 1)设系统局部模型式(14)是能控的,那么局部模型式(14)是能观的,当且仅当

${\rm{Rank}}({{{O}}_k}) = {2^n}$ (25)

2)设系统全局模型式(14)是能控的,那么全局模型式(17)是能观的,当且仅当

${\rm{Rank}}({{O}}) = {2^n}$ (26)

式中 ${{O}} = \displaystyle\sum\limits_{k = 1}^N {{\mu _k}(z){{{O}}_k}} $

证明[16]  对给定的状态 ${X_0}$ ,可以观测 ${{HX}}_0$ ,因为系统是能控的,所以使用不同的控制序列 ${U_i}$ ,可以观测 ${{{HLU}}_i}$ ,故 ${{{HLU}}_{{i_1}}}{{{LU}}_{{i_2}}} \cdots {{{LU}}_{{i_s}}}{{{X}}_0}$ 是可观测的。因为 $s \geqslant {k_0}$ ,没有增加之前集合的线性独立的行数,线性独立行数对初始状态的辨识是无用的,只有当 ${{O}}$ 包含所有不同列,即 ${\rm{Rank}}({{O}}) = {2^n}$ 时,初始状态才能辨识,全局模型是能观测的。

3.3 稳定性

对于布尔系统式(8)和布尔控制系统式(9)记 $\chi = {{\cal{D}}^{n}}$ 为它们的状态空间。点 $X \in \chi $ 可以表示为 $X = {[{x_1},{x_2}, \cdots ,{x_n}]^{\rm T}}$ ,逻辑映射 $F:\chi \to \chi $ ,逻辑映射形式为[13]

$\left\{ \begin{gathered} {z_1} = {f_1}({x_1},{x_2}, \cdots ,{x_n}) \\ \;\;\;\;\;\;{\kern 1pt} {\kern 1pt} \vdots \\ {z_n} = {f_n}({x_1}, {x_2},\cdots ,{x_n}) \\ \end{gathered} \right.$ (27)

简记为 $Z = F(X)$ , $X,Z \in \chi $

定义7[13]  逻辑映射 ${{F}}$ 的关联矩阵 ${{{I}}}(\cal F)$ =(bij)是一个 $n \times n$ 矩阵,定义为

${b_{ij}} = \left\{ {\begin{aligned}& {1,\quad {x_j}\left( {t + 1} \right)\text{依赖于}{x_i}\left( t \right)}\\& {0,\quad \text{其他}}\end{aligned}} \right.$ (28)

${{X}} = {[{x_1},{x_2}, \cdots ,{x_n}]^{\rm T}},{{F}}= {[{f_1},{f_2} \cdots ,{f_n}]^{\rm T}}$ ,则式(14)对应的布尔网络的逻辑映射可简记为

$\left\{\begin{aligned} &{{X}}(t + 1) = F({{X}}(t),{{U}}(t)) \\ &{{X}}(t) \in {{\cal{D}}^n},{{U}}(t) \in {{\cal{D}}^m}\end{aligned}\right. $ (29)

同理式(15)对应的布尔控制网络的逻辑映射可简记为

${{X}}(t + 1) = {F_k}({{X}}(t),{{U}}(t)){\kern 1pt} $ (30)

式中 ${{X}}(t) \in {{\cal{D}}^n},{{U}}(t) \in {{\cal{D}}^m}$

定理5[13]  1)设 $\xi $ 是式(14)的一个不动点,则 $X(k)\bar \vee \xi \leqslant {{\cal I}}{({F_k})^j} \times \left( {X(0)\bar \vee \xi } \right)$ ,如果存在 $j > 0$ ,使得 ${[{\cal{I}}({F_k})]^{(j)}} = 0$ ,则称局部模型是能稳定的。

2)设 $\xi $ 是式(15)的一个不动点,则 $X(k)\bar \vee \xi \leqslant $ ${\cal I}{(F)^j} \times \left( {X(0)\bar \vee \xi } \right)$ ,如果存在 $j > 0$ ,使得 $[{\cal{I}}(F)]$ (j)=0,则称全部模型是能稳定的。其中 $F = \displaystyle\sum\limits_{k = 1}^N {{\mu _k}(z){F_k}} $

证明[14]  ①必要性:如果系统的局部模型是稳定的,即系统是收敛的,则 $T$ 步后( $T$ 为极短的时间段),所有的状态收敛到 $\xi $ ,所以当 $j > T$ 时成立。

②充分性:假设存在 $j > 0$ ,使 ${[{\cal{I}}({F_k})]^{(j)}} = 0$ 成立,那么对于任意的 $X$ $F_k^j(X) = \xi $ $X(t) \in {{\cal{D}}^n}$ 。故对任意步数 $t \geqslant j, F_k^t(X) = F_k^j(F_k^{t - j}(X)) = \xi $ ,得证。

同理可证全局模型时系统的稳定性。

4 实验仿真

根据第3节介绍的模糊动态布尔网络控制系统,本节选取多输入多输出模糊模型进行非线性模糊建模及其能控性、能观性分析。

$\begin{aligned} & {R^1}:{\rm{if}}\;{z_1}\;{\rm{is}}\;{F_1}, {z_2}\;{\rm{is}}\;{F_2}\\ & {\rm{Then}}\left\{ \begin{aligned} & {x_1}(t + 1) = {x_2}(t) \wedge {u_1}(t) \\ & {x_2}(t + 1) = {x_3}(t) \vee {u_2}(t) \\ & {x_3}(t + 1) = {x_1}(t) \\ & {y_1}(t) = {x_1}(t) \\ & {y_2}(t) = \neg {x_2}(t) \\ \end{aligned} \right. \\ & {R^2}:{\rm{if}}\;{z_1}\;{\rm{is}}\;{F_2}, {z_2}\;{\rm{is}}\;{F_1} \\ &{\rm{Then}}\left\{ \begin{aligned} & {x_1}(t + 1) = {x_3}(t) \wedge {u_1}(t) \\ & {x_2}(t + 1) = \neg {u_2}(t) \\ & {x_3}(t + 1) = {x_1}(t) \vee {x_2}(t) \\ & {y_1}(t) = {x_1}(t) \\ & {y_2}(t) = {x_2}(t) \vee {x_3}(t) \\ \end{aligned} \right. \\ \end{aligned} $ (31)

其代数表达式为

$\begin{aligned} & {R^1}:{\rm{if}}\;{z_1}\;{\rm{is}}\;{F_1}, {z_2}\;{\rm{is}}\;{F_2}\\ & {\rm{Then}}\left\{ \begin{aligned} & x(t + 1) = {L_1}u(t)x(t) \\ & y(t) = {H_1}x(t) \\ \end{aligned} \right. \\ & {R^2}:{\rm{if}}\;{z_1}\;{\rm{is}}\;{F_2},{z_2}\;{\rm{is}}\;{F_1} \\& {\rm{Then}}\left\{ \begin{aligned} & x(t + 1) = {L_2}u(t)x(t) \\ & y(t) = {H_2}x(t) \\ \end{aligned} \right. \end{aligned} $ (32)

则整个系统模糊状态方程可表示为式(15),其中

$\begin{aligned}&{{{L}}_1} = {\delta _8}[1\,1\,5\,5\,2\,2\,6\,6\,1\,3\,5\,7\,2\,4\,6\,8\,5\,5\,5\,5\,6\,6\,6\,6\,5\,7\,5\,7\,6\,8\,6\,8]\end{aligned}$
$\begin{aligned}& {{{L}}_2} = {\delta _8}[3\,1\,7\,5\,3\,1\,7\,5\,7\,5\,7\,5\,7\,5\,7\,5\,3\,1\,7\,5\,3\,1\,7\,5\,8\,6\,8\,6\,8\,6\,8\,6]\end{aligned} $
${{{H}}_1} = {\delta _4}[2{\kern 1pt}\;\; 1{\kern 1pt}\;\; 2{\kern 1pt}\;\; 1{\kern 1pt}\;\; 4{\kern 1pt}\;\; 3{\kern 1pt}\;\; 4{\kern 1pt}\;\; 3],\quad{{{H}}_2} = {\delta _4}[1{\kern 1pt}\;\; 1{\kern 1pt}\;\; 1{\kern 1pt}\;\; 2{\kern 1pt}\;\; 3{\kern 1pt}\;\; 3{\kern 1pt}\;\; 3{\kern 1pt}\;\; 4]$
${{L}} = {\mu _1}(z){{{L}}_1} + {\mu _2}(z){{{L}}_2},\quad{{H}} = {\mu _1}(z){{{H}}_1} + {\mu _2}(z){{{H}}_2},$
${\mu _1} = {\mu _{{F_{11}}}}({z_1}) \wedge {\mu _{{F_{21}}}}({z_2}),\quad{\mu _2} = {\mu _{{F_{12}}}}({z_1}) \wedge {\mu _{{F_{22}}}}({z_2}){\text{。}}$

由选取的模糊模型,规则1和规则2时的状态变量、输入变量和输出变量的初始状态可分别设为

$\left\{ \begin{aligned} & {x_1}(0) = 0,\quad {x_2}(0) = 1,\quad{x_3}(0) = 0 \\ & {u_1}(0) = 1,\quad{u_2}(0) = 0 \\ & {y_1}(0) = 0,\quad{y_2}(0) = 0 \\ \end{aligned} \right.$
4.1 能控性

由式(16),将规则1时的结构矩阵等分为4块:

$\begin{aligned}{{{L}}_1} & = [{\rm{Bl}}{{\rm{k}}_1}({{{L}}_1})\;\;{\rm{Bl}}{{\rm{k}}_2}({{{L}}_1})\;\;{\rm{Bl}}{{\rm{k}}_3}({{{L}}_1})\;\;{\rm{Bl}}{{\rm{k}}_4}({{{L}}_1})]=\\& [{{{B}}_{11}}\;{{{B}}_{12}}\;{{{B}}_{13}}\;{{{B}}_{14}}]\end{aligned}$
${{{B}}_{11}} = {\delta _8}[1\;1\;5\;5\;2\;2\;6\;6],\quad {{{B}}_{12}} = {\delta _8}[1\;3\;5\;7\;2\;4\;6\;8],$
${{{B}}_{13}} = {\delta _8}[5\;5\;5\;5\;6\;6\;6\;6],\quad {{{B}}_{14}} = {\delta _8}[5\;7\;5\;7\;6\;8\;6\;8]$

$\begin{array}{c} {{{M}}_1} = \sum\limits_{i = 1}^{{2^m}} {_{\cal{B}}{\rm{Bl}}{{\rm{k}}_i}({{{L}}_1})} = \sum\limits_{i = 1}^4 {_{\cal{B}}{\rm{Bl}}{{\rm{k}}_i}({{{L}}_1})} = \\ \left[ {\begin{array}{*{20}{l}}1&1&0&0&0&0&0&0\\0&0&0&0&1&1&0&0\\0&1&0&0&0&0&0&0\\0&0&0&0&0&1&0&0\\1&1&1&1&0&0&0&0\\0&0&0&0&1&1&1&1\\0&1&0&1&0&0&0&0\\0&0&0&0&0&1&0&1\end{array}} \right]\end{array}$

其能控性矩阵为

${{{C}}_1} = \sum\limits_{j = 1}^{{2^{m + n}}} {_{\cal{B}}\sum\limits_{i = 1}^m {_{\cal{B}}{B_{1i}}^{(j)}} } = \sum\limits_{j = 1}^{{2^5}} {_{\cal{B}}{M_1}^{(j)}} = {\rm{10^{12}}} \times \\\left[ \begin{gathered} {\rm{1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;0}}{\rm{.810\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0}} \\ {\rm{2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;5}}{\rm{.555\;6\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5}} \\ {\rm{0}}{\rm{.810\;5\;\;1}}{\rm{.311\;5\;\;0}}{\rm{.500\;9\;\;0}}{\rm{.810\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;0}}{\rm{.810\;5\;\;1}}{\rm{.311\;5}} \\ {\rm{1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;0}}{\rm{.810\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0}} \\ {\rm{2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;5}}{\rm{.555\;6\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5}} \\ {\rm{3}}{\rm{.433\;5\;\;5}}{\rm{.555\;6\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;5}}{\rm{.555\;6\;\;8}}{\rm{.989\;1\;\;3}}{\rm{.433\;5\;\;5}}{\rm{.555\;6}} \\ {\rm{1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;0}}{\rm{.810\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0}} \\ {\rm{2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;1}}{\rm{.311\;5\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5\;\;5}}{\rm{.555\;6\;\;2}}{\rm{.122\;0\;\;3}}{\rm{.433\;5}} \end{gathered} \right] > 0$

故规则1时的局部模型是可控的。

由MATLAB仿真知,规则1时控制变量和状态变量的关系如图1(a)所示。

Download:
图 1 控制变量和状态变量的关系 Fig. 1 Schematic diagram of relationship between the control variables and state variables

同理,规则2时有:

${{{B}}_{21}} = {\delta _8}[3\;1\;7\;5\;3\;1\;7\;5],\quad {{{B}}_{22}} = {\delta _8}[7\;5\;7\;5\;7\;5\;7\;5],$
${{{B}}_{23}} = {\delta _8}[3\;1\;7\;5\;3\;1\;7\;5],\quad{{{B}}_{24}} = {\delta _8}[8\;6\;8\;6\;8\;6\;8\;6]$
$\begin{aligned} {{{M}}_2} & = \displaystyle\sum\limits_{i = 1}^{{2^m}} {_{\cal{B}}{\rm{Bl}}{{\rm{k}}_i}({{{L}}_2})} = \displaystyle\sum\limits_{i = 1}^4 {_{\cal{B}}{\rm{Bl}}{{\rm{k}}_i}({{{L}}_2})}=\\[-4pt]&\left[ {\begin{array}{*{20}{l}}1&1&0&0&0&0&0&0\\[-4pt]0&0&0&0&1&1&0&0\\[-4pt]0&1&0&0&0&0&0&0\\[-4pt]0&0&0&0&0&1&0&0\\[-4pt]1&1&1&1&0&0&0&0\\[-4pt]0&0&0&0&1&1&1&1\\[-4pt]0&1&0&1&0&0&0&0\\[-4pt]0&0&0&0&0&1&0&1\end{array}} \right]\end{aligned}$

其能控性矩阵有

${{{C}}_2} = \sum\limits_{j = 1}^{{2^{m + n}}} {_{\cal{B}}\sum\limits_{i = 1}^m {_{\cal{B}}{B_{2i}}^{(j)}} } = \sum\limits_{j = 1}^{{2^5}} {_{\cal{B}}{M_2}^{(j)}} = {\rm{10^{11}}} \times \\\left[ \begin{gathered} {\rm{1}}{\rm{.298\;6\;\;1}}{\rm{.836\;5\;\;0}}{\rm{.918\;2\;\;1}}{\rm{.298\;6\;\;1}}{\rm{.298\;6\;\;1}}{\rm{.836\;5\;\;0}}{\rm{.918\;2\;\;1}}{\rm{.298\;6}} \\ {\rm{0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0}} \\ {\rm{1}}{\rm{.836\;5\;\;2}}{\rm{.597\;2\;\;1}}{\rm{.298\;6\;\;1}}{\rm{.836\;5\;\;1}}{\rm{.836\;5\;\;2}}{\rm{.597\;2\;\;1}}{\rm{.298\;6\;\;1}}{\rm{.836\;5}} \\ {\rm{0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0 \quad\quad\quad 0}} \\ {\rm{3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7\;\;2}}{\rm{.216\;8\;\;3}}{\rm{.135\;1\;\;3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7\;\;2}}{\rm{.216\;8\;\;3}}{\rm{.135\;1}} \\ {\rm{3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7\;\;2}}{\rm{.216\;8\;\;3}}{\rm{.135\;1\;\;3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7\;\;2}}{\rm{.216\;8\;\;3}}{\rm{.135\;1}} \\ {\rm{4}}{\rm{.433\;7\;\;6}}{\rm{.270\;1\;\;3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7\;\;4}}{\rm{.433\;7\;\;6}}{\rm{.270\;1\;\;3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7}} \\ {\rm{4}}{\rm{.433\;7\;\;6}}{\rm{.270\;1\;\;3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7\;\;4}}{\rm{.433\;7\;\;6}}{\rm{.270\;1\;\;3}}{\rm{.135\;1\;\;4}}{\rm{.433\;7}} \end{gathered} \right] = 0$

故规则2时的局部模型是不可控的。

图1可知,规则1时的局部模型是可控的,规则2时的局部模型是不可控的,与实验计算结果是一致的。

4.2 能观性

${{{\varOmega}} _{k\left( {s + 1} \right)}} \subset \bigcup\limits_{i = 1}^s {{{{\varOmega}} _i}} ,$ $s = 3$ ,则式(25)中 ${{{\varGamma}} _{10}}\text{、}{{{\varGamma}} _{11}}$ ${{{\varGamma}} _{12}}$ ${{{\varGamma}} _{13}} $ 分别为 ${{{H}}_1} = {\delta _4}[2\;1\;2\;1\;4\;3\;4\;3]$

${{{\varGamma}} _{10}} = {{{H}}_1} = \left[ {\begin{array}{*{20}{l}}0&1&0&1&0&0&0&0\\1&0&1&0&0&0&0&0\\0&1&0&0&0&1&0&1\\0&0&0&0&1&0&1&0\end{array}} \right]$
${{{\varGamma}} _{11}}=\left[ {\begin{aligned}{{{H}}{}_1{{{B}}_{11}}}\\{{{H}}{}_1{{{B}}_{12}}}\\{{{H}}{}_1{{{B}}_{13}}}\\{{{H}}{}_1{{{B}}_{14}}}\end{aligned}} \right]=\\\left[ {\begin{array}{*{20}{l}}0&1&0&0&0&1&0&0&0&0&0&1&0&0&0&1\\0&1&0&0&0&1&0&0&0&0&0&1&0&0&0&1\\0&0&0&1&0&0&0&1&0&0&0&1&0&0&0&1\\0&1&0&0&0&1&0&0&0&0&0&1&0&0&0&1\\1&0&0&0&1&0&0&0&0&0&1&0&0&0&1&0\\1&0&0&0&1&0&0&0&0&0&1&0&0&0&1&0\\0&0&1&0&0&0&1&0&0&0&1&0&0&0&1&0\\0&0&1&0&0&0&1&0&0&0&1&0&0&0&1&0\end{array}} \right]$
${{{\varGamma}} _{12}} = \left[ \begin{array}{c}H{}_1{B_{11}}{B_{11}}\\ \vdots \\{H_1}{B_{12}}{B_{11}}\\ \vdots \\{H_1}{B_{13}}{B_{11}}\\ \vdots \\{H_1}{B_{14}}{B_{14}}\end{array} \right] = \left[\begin{array}{*{20}{c}}0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad\cdots\quad0\quad0\quad1\quad0\\1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\\1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad\cdots\quad0\quad0\quad1\quad0\\0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad\cdots\quad0\quad0\quad1\quad0\\1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\\1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad\cdots\quad0\quad0\quad1\quad0\end{array} \right]_{8 \times 64}^{\rm T}$

Γ12为64×8维的布尔矩阵,由于篇幅的限制,上式中只列出了一部分。

${{{\varGamma}} _{13}} = \left[ \begin{array}{c}{H_1}{B_{11}}{B_{11}}{B_{11}}\\{H_1}{B_{11}}{B_{11}}{B_{12}}\\ \vdots \\{H_1}{B_{14}}{B_{14}}{B_{13}}\\{H_1}{B_{14}}{B_{14}}{B_{14}}\end{array} \right] = \left[\begin{array}{*{20}{c}}0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\\0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\end{array} \right]_{8 \times 256}^{\rm T}$

Γ13为256×8维的布尔矩阵,由于篇幅的限制,上式中只列出了一部分。

${{{O}}_{1}} = \left[ \begin{array}{l}{{{\varGamma}}_{10}}\\{{{\varGamma}}_{11}}\\{{{\varGamma}}_{12}}\\{{{\varGamma}}_{13}}\end{array} \right] = \left[\begin{array}{*{20}{c}}0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\1\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\0\quad1\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\1\quad0\quad0\quad0\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad\cdots\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\\0\quad0\quad0\quad1\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\\0\quad0\quad1\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\\0\quad0\quad0\quad1\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\\0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\quad\cdots\quad0\quad0\quad1\quad0\quad0\quad0\quad1\quad0\end{array} \right]_{8 \times 340}^{\rm T}$

O1为340×8维的布尔矩阵,由于篇幅的限制,文中只列出了一部分。

可得,Rank (O1)-7<23=8。

由MATLAB仿真知,规则1时输出变量和状态变量的关系如图2所示。由图2知,虽然规则1时的局部模型是能控的,但是是不能观测的;规则2时的局部模型是不可控的,故也是不能观测的,与实验计算结果是一致的。

Download:
图 2 输出变量和状态变量的关系 Fig. 2 Schematic diagram of relationship between the output variables and state variables under rule 1
4.3 稳定性

规则1时,式(33)对应的布尔网络的逻辑映射为 $X(t + 1) = {F_1}(X(t)),X \in {{\cal{D}}^3}$

关联矩阵为

${{\cal{I}}}({F_1}) = \left[ {\begin{array}{*{20}{c}} 0&1&0 \\ 0&0&1 \\ 1&0&0 \end{array}} \right]$

${[{\cal{I}}({F_1})]^{(5)}} = {\cal{I}}({F_1}) \ne 0$ ,即 ${[{\cal{I}}({F_1})]^{(k)}} \ne 0$ ,则规则1下的局部模型是不稳定的。

规则2时,式(31)对应的布尔网络的逻辑映射为 $X(t + 1) = {F_2}(X(t)),X \in {{\cal{D}}^3}$

关联矩阵为

${\cal{I}}({F_2}) = \left[ {\begin{array}{*{20}{c}} 0&0&1 \\ 0&0&0 \\ 1&1&0 \end{array}} \right]$

${[{\cal{I}}({F_2})]^{(3)}} = 0$ ,则规则2下的局部模型是稳定的。

综上分析,整个系统模糊状态方程可表示为式(15),其中,取 ${\mu _1} = 1,{\mu _2} = 0,$ ${{L}} = {{{L}}_1}, {{H}} = {{{H}}_1},$ $ {{F}} =$ $ {{{F}}_1},$ ${{C}} = {{{C}}_1} > 0,$ ${\rm{Rank}}({{O}}) = {\rm{Rank}}({{{O}}_1}) = 7 < {2^3} = 8$ , 因此系统的全局模型是能控不能观且不稳定的。

5 结束语

通过具体实例仿真分析可知,对于布尔网络系统,可以实现非线性模糊建模和动态性能分析。利用模糊动态模型的非线性特点,将模糊动态模型和布尔控制网络相结合,分别建立了模糊动态布尔网络控制系统的局部模型和全局模型,并且对其能控性、能观性和稳定性进行了分析。

参考文献
[1] ZADEH L A. Fuzzy sets[J]. Information and control, 1965, 8(3): 338-353. DOI:10.1016/S0019-9958(65)90241-X (0)
[2] MAMDANI E H. Application of fuzzy algorithms for control of simple dynamic plant[J]. Proceedings of the institution of electrical engineers, 1974, 121(12): 1585-1588. DOI:10.1049/piee.1974.0328 (0)
[3] TAKAGI T, SUGENO M. Fuzzy identification of systems and its applications to modeling and control[J]. IEEE transactions on systems man and cybernetics, 1985, 15(1): 116-132. (0)
[4] ABONYI J, NAGY L, SZEIFERT F. Fuzzy model-based predictive control by instantaneous linearization[J]. Fuzzy sets and systems, 2001, 120(1): 109-122. (0)
[5] FENG G. A survey on analysis and design of model-based fuzzy control systems[J]. IEEE transactions on fuzzy systems, 2006, 14(5): 676-697. DOI:10.1109/TFUZZ.2006.883415 (0)
[6] CAO S G, REES N W, FENG G. Analysis and design for a class of complex control systems part I: fuzzy modelling and identification[J]. Automatica, 1997, 33(6): 1017-1028. DOI:10.1016/S0005-1098(97)00010-1 (0)
[7] CAO S G, REES N W, FENG G. Analysis and design for a class of complex control systems part II: fuzzy controller design[J]. Automatica, 1997, 33(6): 1029-1039. DOI:10.1016/S0005-1098(97)00011-3 (0)
[8] SILVA L F P, LEITE V J S, CASTELAN E B, et al. Stabilization of time-delay nonlinear discrete-time systems with saturating actuators through T-S models[J]. IFAC proceedings volumes, 2014, 47(3): 11000-11005. DOI:10.3182/20140824-6-ZA-1003.01502 (0)
[9] ZHANG Changzhu, FENG Gang, Qiu Jianbin, et al. T-S fuzzy-model-based piecewise H output feedback controller design for networked nonlinear systems with medium access constraint [J]. Fuzzy sets and systems, 2014, 248: 86-105. DOI:10.1016/j.fss.2013.09.006 (0)
[10] CAO Shuguang, REES N W, FENG Gang. Analysis and design of fuzzy control systems using dynamic fuzzy-state space models[J]. IEEE transactions on fuzzy systems, 1999, 7(2): 192-200. DOI:10.1109/91.755400 (0)
[11] KAUFFMAN S A. Metabolic stability and epigenesis in randomly constructed genetic nets[J]. Journal of theoretical biology, 1969, 22(3): 437-467. (0)
[12] Cheng D Z, Qi H S, Semi-tensor Product of Matrices Theory and Application[M]. Beijing: Science Press, 2007. (0)
[13] 程代展, 齐洪胜, 赵寅. 布尔网络的分析与控制——矩阵半张量积方法[J]. 自动化学报, 2011, 37(5): 529-540.
CHENG Daizhan, QI Hongsheng, ZHAO Yan. Analysis and control of boolean networks: a semi-tensor product approach[J]. Acta automatica sinica, 2011, 37(5): 529-540. (0)
[14] CHENG Daizhan, QI Hongsheng, LI Zhiqiang, et al. Stability and stabilization of Boolean networks[J]. International journal of robust and nonlinear control, 2011, 21(2): 134-156. (0)
[15] CHENG Daizhan, LI Zhiqiang, QI Hongsheng. Realization of boolean control networks[J]. Automatica, 2010, 46(1): 62-69. DOI:10.1016/j.automatica.2009.10.036 (0)
[16] CHENG Daizhan, QI Hongsheng. Controllability and observability of boolean control networks[J]. Automatica, 2009, 45(7): 1659-1667. DOI:10.1016/j.automatica.2009.03.006 (0)
[17] CHENG Daizhan. Disturbance decoupling of boolean control networks[J]. IEEE transactions on automatic control, 2011, 56(1): 2-10. DOI:10.1109/TAC.2010.2050161 (0)
[18] CHENG Daizhan, QI Hongsheng. A linear representation of dynamics of boolean networks[J]. IEEE transactions on automatic control, 2010, 55(10): 2251-2258. DOI:10.1109/TAC.2010.2043294 (0)
[19] CHENG Daizhan, QI Hongsheng, LI Zhiqiang. Model construction of boolean network via observed data[J]. IEEE transactions on neural networks, 2011, 22(4): 525-536. DOI:10.1109/TNN.2011.2106512 (0)
[20] CHENG Danzhan, QI Hongsheng. State–space analysis of boolean networks[J]. IEEE transactions on neural networks, 2010, 21(4): 584-594. DOI:10.1109/TNN.2009.2039802 (0)
[21] CHENG Daizhan, QI Hongsheng, LI Zhiqiang. Feedback decomposition of boolean control networks[M]//Analysis and Control of Boolean Networks. London: Springer, 2011: 297–311. (0)
[22] 宋金利, 肖会敏, 李志强. 布尔控制网络的部分变量能控性[J]. 中国科学: 信息科学, 2016, 46(3): 338-349.
SONG Jinli, XIAO Huimin, LI Zhiqiang. Partial variables controllability of Boolean control networks[J]. Scientia sinica informationis, 2016, 46(3): 338-349. (0)
[23] 李志强, 宋金利. 布尔控制网络的能控性与能观性[J]. 控制理论与应用, 2013, 30(6): 760-764.
LI Zhiqiang, SONG Jinli. Controllability and observability of boolean control networks[J]. Control theory and applications, 2013, 30(6): 760-764. (0)