Recursive Bilinear Subspace Modeling and Model-free Adaptive Control of Wastewater Treatment
-
摘要: 污水处理过程中, 生化反应硝态氮浓度和溶解氧浓度是决定出水水质好坏的两个最关键变量, 难以采用常规基于模型的方法进行有效控制. 本文基于数据驱动建模与控制技术, 提出一种污水处理过程递推双线性子空间辨识(Recursive bilinear subspace identification, RBLSI)建模和无模型自适应控制方法. 首先, 针对污水处理过程的非线性时变动态特性, 采用最小二乘递推双线性子空间辨识方法建立污水处理生化反应过程具有参数自适应能力的递推双线性模型; 其次, 基于建立的数据驱动模型, 采用基于多参数灵敏度分析(Multi-parameter sensitivity analysis, MPSA)和遗传粒子群优化(Genetic algorithm-particle swarm optimization, GA-PSO)算法的无模型自适应控制(Model-free adaptive control, MFAC)方法对硝态氮和溶解氧浓度进行直接数据驱动控制; 最后, 数据实验及其比较分析表明了所提方法的有效性和优越性.
-
关键词:
- 污水处理 /
- 递推双线性子空间辨识 /
- 无模型自适应控制 /
- 多参数灵敏度分析
Abstract: Conventional model-based approaches are unable to control the nitrate nitrogen concentration and dissolved oxygen concentration of biochemistry reaction effectively, which are the two most critical variables that determine the quality of effluent in wastewater treatment process. In this paper, a recursive bilinear subspace identification (RBLSI) modeling and model-free adaptive control method for wastewater treatment based on data-driven modeling and control technology was proposed. Firstly, according to the nonlinear time-varying dynamic characteristics of wastewater treatment process, a recursive bilinear model with parameter adaptability for biochemical reaction process of wastewater treatment was established by using the least square recursive bilinear subspace identification method. Secondly, the model-free adaptive control (MFAC) method based on the multi-parameter sensitivity analysis (MPSA) and genetic algorithm−particle swarm optimization (GA-PSO) algorithm was used to directly control the nitrate concentration and dissolved oxygen concentration in a data-driven mode based on the established data-driven model. Finally, data experiments and comparative analysis show the effectiveness and superiority of the proposed method. -
发票是财务系统的重要组成部分, 发票的识别和处理是财务人员的首要工作. 人工智能已被广泛应用在人脸识别、车辆识别、物体检测等各个领域, 而基于人工智能的发票识别, 可为财务人员节省大量的时间. 目前发票自动识别根据使用场景不同, 可分为移动客户端、云端和本地客户端. 移动端客户端发票识别是通过手机端的离线SDK, 集成到公司财务系统的APP内, 自动检测发票的边线并拍照识别, 提取发票上的关键信息; 云端发票识别是通过调用云平台的API接口, 将拍照或者扫描的发票图像传送到云平台上进行识别; 而本地客户端识别是识别软件与扫描仪硬件配合使用进行识别. 三种方式本质都是基于识别软件进行的. 由于发票种类繁多, 格式不固定, 有的发票存在手写的内容, 且不同人手写字的风格不一, 因此基于传统数字图像处理的识别软件, 相应识别准确率较低. 随着卷积神经网络的兴起, 包括表格线定位、手写文字识别等问题得到很好的解决. 文献[1]采用卷积神经网络实现增值税发票批量识别, 通过对发票图像进行采集、预处理、字符切割, 基于CNN神经网络进行字符识别, 同时融入人工修改, 提高识别率. 但是发票上的印章对发票识别准确率依旧具有一定影响, 且印章是每张发票必不可少的, 因此如何消除印章也是现在OCR (Optical character recognition)领域的研究热点. 目前大多数研究者对于发票印章的研究主要是印章识别与验证[2-4]. 针对发票印章的去除问题, 若通过印章定位, 而后直接挖除印章则会丢失印章下的发票内容; 若采用过滤颜色的方式来消除, 会导致发票中与印章颜色相近的文字变得模糊不清, 而且不同发票的印章颜色不同. 文献[5]提出一种基于彩色图像的二次分割算法来去除印章, 但是对于发票的要求较高, 需要保证印章的颜色为红色或者蓝色, 票据不能出现明显的扭曲和变形, 字体还需要求是黑色或灰色, 其实用性较差. 文献[6]针对印章覆盖、发票折痕等干扰因素影响发票号码分割的问题, 采用基于数字结构特征的识别算法, 通过对噪声粘连区域进行修复, 解决干扰因素对数字分割的影响, 实现发票号码识别. 但是如何判断噪声区域最小连通面积的阈值是算法关键, 当发票数量增多, 印章颜色深度和位置不统一时, 消除印章干扰则很困难.
各种图像处理软件也可用在印章消除上, 如PhotoShop, 但需要基于软件进行手动处理, 大量发票的印章消除则会新增大量工作量, 反倒得不偿失. 因此研究如何自动消除发票印章, 对提高发票识别准确率具有重要意义. 生成式对抗网络GAN (Generative adversarial network)是Goodfellow等在2014年提出的一种生成式网络[7]. 在GAN提出之前, 常见的生成式网络有最大似然估计法、近似法[8]、马尔科夫链法[9]等. 这些方法的基本步骤是, 先对样本数据进行分布假设, 然后利用数据样本直接对假设分布的参数进行拟合, 这会导致学习到的生成式模型受到限制[10]. 而GAN不同于上述的生成网络, 该方法采用对抗学习方式, 先通过判别器学习生成分布和真实分布之间的差异, 再驱使生成器去缩小差异. GAN相比于其他的生成网络, 对数据的分布不做显性的限制, 从而避免了人工设计网络分布[11]. GAN目前常用于实现图像的风格迁移以及超分辨图像的生成. 本文基于生成对抗网络提出用于消除印章的SealGAN网络. SealGAN借鉴了CycleGAN网络的循环结构, 采用两个独立的分类器来取代判别网络, 并针对印章的特点去设计生成网络, 实现发票图像的带印章到不带印章的风格迁移, 从而达到消除发票印章的效果.
1. 生成式对抗网络和CycleGAN网络
1.1 生成式对抗网络
生成式对抗网络结构如图1所示, 包括一个生成网络(G)和一个判别网络(D). 生成网络用于生成目标图像, 而判别网络用于判断图像是生成的虚假图像还是真实目标图像. 在训练时, 生成网络会在损失函数的牵引下让生成的图像尽可能地逼近目标图像, 判别网络则将生成图像和真实目标图像区分开, 由此构成对抗博弈的训练方式. 通过竞争方式, 生成网络生成的图像会越来越真实, 而判别网络的判别能力也会越来越强[7].
图1中
$ {{L}}_{G} $ 是更新生成网络G的损失,$ {{L}}_{D} $ 是更新判别网络D的损失, 如下式所示:$$ {{L}}_{G}=-\dfrac{1}{m}\displaystyle\sum _{i}^{m}{\rm{log}}\left(D\right(G\left({z}^{i}\right)\left)\right) $$ (1) $$ \begin{split} {{L}}_{D}=\,&-\frac{1}{m}\sum _{i}^{m}[{\rm{log}}\left(D\left({x}^{i}\right)\right)+ \\ &{\rm{log}}\left(1-D\left(G\left({z}^{i}\right)\right)\right)] \end{split}$$ (2) 式中
$ {z}^{i} $ 表示的是第i个输入,$ {x}^{i} $ 表示对应的真实图像. 生成网络和判别网络的损失函数均为交叉熵损失, 只是判别网络的目标是让真实图像的判别结果趋近于1, 生成图像结果趋近于0, 而生成网络则希望生成图像的判别结果趋近于1[7].1.2 CycleGAN网络
生成式对抗网络的概念提出后, 便成为学术界研究的热点. 在此基础上, 一些改进网络结构相继被提出, 如引入卷积神经网络的DCGAN[12], 用于生成超分辨图像的SRGAN、LapGAN[13-14], 用于风格迁移的pix2pix、cVAE-GAN[15-16]. 但是大部分生成式对抗网络都要求待处理图像和目标图像一一匹配. 针对带印章的发票, 则很难获取印章图像与其对应的无章图像. 2017年, 伯克利大学的研究团队在ICCV上提出了CycleGAN网络, 其最大特点是不需要待处理图像和目标图像一一匹配, 只需要两个分布不同的数据集即可实现数据集之间的相互风格迁移.
CycleGAN的基本思想是采用两个基础的生成式对抗网络, 采用循环的方式训练, 其网络结构如图2所示[17].
假设有两个不同分布的数据集A与B, CycleGAN中包含两个基础GAN网络, 一个负责把A类数据转化成B类数据, 记作
$ {GAN}_{A\; to\; B} $ , 另一个负责把B类数据转化成A类数据, 记作$ {GAN}_{B\; to\; A} .$ 设原始的A类图像为$ {X}_{A} $ , 通过$ {GAN}_{A\; to\; B} $ , 可以得到原始图像对应的B类风格图像$ {X'}_{B} $ . 如果采用传统图像对模式, 则存在$ {X}_{A} $ 对应的真实图像$ {X}_{B} $ , 这时只需将$ {X'}_{B} $ 和$ {X}_{B} $ 做均方差之类的损失计算, 通过反向传播更新生成网络参数即可驱使生成式网络生成内容相匹配的图像. 当无法获取图像对时, CycleGAN则把生成的$ {X'}_{B}$ 输入$ {GAN}_{B\;to\;A} $ , 生成相应的A类图像$ {X'}_{A} $ , 换言之就是试图把生成的$ {X'}_{B} $ 还原成原图像$ {X}_{A} $ , 此时计算$ {X}_{A} $ 和$ {X'}_{A} $ 之间的损失即可反向传播去更新两个生成对抗网络. 这种循环的训练模式巧妙地实现了不存在图像对时的两类数据集之间相互转化.2. 基于CycleGAN改进的SealGAN网络
2.1 SealGAN网络结构
SealGAN网络结构如图3所示,
$ {GEN}_{A \;to \;B} $ 和$ {GEN}_{B\; to\; A} $ 是生成式对抗网络中抽取出来的生成网络, 采用循环训练的方式. 与CycleGAN网络相比, 本文所提出的SealGAN中, 每个生成网络并没有配备相应的专属判别网络, 而是采用原始图像和所有生成图像单独训练了两个分类器, 其中$ {Classifier}_{AB} $ 用于区分A类和B类数据,$ {Classifier}_{RF} $ 用于区分真实图像和生成图像, 然后将两个分类器所计算出的损失加入到总损失中来更新整体网络的参数.在CycleGAN网络中, 采用
$ {GAN}_{A\; to\; B} $ 和$ {GAN}_{B\; to\; A} $ 两个基础的生成式对抗网络进行循环训练, 且有各自的判别网络. 以$ {GAN}_{A\; to\; B} $ 为例, 它的判别网络既需要判断输入图像是属于A类还是B类, 还需要判断输入图像是生成图像还是原始真实图像. 其实, 生成对抗网络中的判别网络是一个分类器, 对于分类器而言, 分类要求越高则分类效果越差, 而且两个判别网络虽然相互独立, 但是分类要求相似, 采用的训练样本交集较大. 因此本文所提出的SealGAN网络采用两个独立的分类器取代两个判别网络, 一个分类器用于区分A类和B类数据, 另一个分类器用于区分生成图像和真实图像. 相当于用两个独立的分类器来取代原本的判别网络, 网络整体参数数量没有太大变化, 但是对于单个分类器其分类要求降低, 分类效果更好.2.2 SealGAN网络中的生成网络
CycleGAN中有两种生成网络结构, 残差网络ResNet结构和UNet结构. ResNet结构用于解决层数较深的网络难以训练的问题. 在ResNet结构提出之前, 研究人员在提升网络性能上往往是单纯的叠加网络层数. 尽管网络层数的增加能够提高网络的特征容量, 理论上能更好地对样本进行拟合. 但是实际训练中, 当网络层数过大时, 参数的训练难度也大幅提升, 过分引入的卷积层和非线性激活不仅无法提升网络性能, 而且还难以维持浅层网络所达到的效果, 从而出现了模型退化问题[18]. ResNet结构的核心是采用了跨层连接思想, 引入了捷径连接的概念[19]. 残差网络结构如图4所示, 由于捷径连接的存在, 在训练过程中, 拟合的不再是目标分布, 而是目标分布和上一层网络输出之间的残差. 当浅层网络能够达到较好效果时, 残差趋近于零. 尽管在拟合残差时也引入了卷积和非线性激活, 但是要想拟合零只需要把所有参数都趋向零即可, 网络训练易于实现. 因此, ResNet结构的引入能够在加深网络层数, 提高网络拟合能力的同时, 也能保证达到浅层网络的效果, 从而消除网络退化问题. 图4中
$ F\left(x\right) $ 为待拟合的残差函数,$ H\left(x\right) $ 为待拟合的完整目标函数.UNet结构常用于图像分割领域, 其网络结构如图5所示. UNet结构包括下采样和上采样的过程. 在上采样过程中, 先把下采样和上采样中尺度相同的层拼接在一起, 再进行卷积和反卷积操作. 上采样层和下采样层之间的连接是UNet结构的亮点, 这些连接使得在上采样过程中可保留了更多的原始图像的纹理信息[20].
本文提出的SealGAN网络中, 生成网络结合ResNet和UNet两种结构, 搭建了一个下采样−精炼−上采样的网络结构, 如图6所示. 由于在消除印章的同时, 也希望能尽可能地保留印章下的内容, 而UNet结构的优势正是能更多地保留原始图像纹理信息, 因此在SealGAN网络的生成网络中, 上下采样部分采用了UNet结构, 而精炼部分采用了ResNet结构以保证网络的效果.
图 6 SealGAN生成网络结构示意图Fig. 6 Schematic diagram of SealGAN generative networks structure3. 实验结果与讨论
3.1 网络参数设置
输入图像的尺寸统一缩放为
$ 256\times 256 $ , 并把图像像素归一化到[−1, 1]区间. 为了保证消除印章后能够尽可能保留印章下的内容, 下采样通道数应多于上采样通道数, 以保证在上采样时能够有足够多的原始图像信息. 多次实验发现, 下采样的通道数与上采样的通道数的比例设置为1.5 : 1时效果最佳. 在生成网络中, 除了输出层, 其余层的激活函数均采用LeakyReLU.$$ {\rm{\sigma }}_{LeakyReLU}\left(x\right)=\left\{\begin{aligned}&x,\qquad\;\;\;{\rm{if}} \;x>0\\ &0.2x,\;\;\;\;\;\;{\rm{else}}\end{aligned}\right. $$ (3) 输出层采用tanh激活函数, 保证输出结果在[−1, 1]区间内.
$$ {\rm{\sigma }}_{{\rm{tanh}}}\left(x\right)=\dfrac{{\rm{sinh}}x}{{\rm{cosh}}x}=\dfrac{{{\rm{e}}}^{x}-{{\rm{e}}}^{-x}}{{\rm{{e}}}^{x}+{{\rm{e}}}^{-x}} $$ (4) 分类器网络前四层的激活函数采用LeakyReLU, 最后一层没有激活. 生成网络以及分类器的相关参数如表1所示.
表 1 生成网络和分类器参数表Table 1 Parameters of the generative network and classifier生成网络 分类器 下采样 精炼 上采样 7×7 conv, 96 Residual_block(3×3, 384) ×7 4×4 deconv, 256, ×2 4×4 conv, 64, /2 4×4 conv, 192, /2 4×4 deconv, 256, ×2 4×4 conv, 128, /2 4×4 conv, 384, /2 4×4 deconv, 256, ×2 4×4 conv, 256, /2 4×4 conv, 384, /2 4×4 deconv, 256, ×2 4×4 conv, 512, /2 4×4 conv, 384, /2 4×4 deconv, 128, ×2 4×4 conv, 1 4×4 conv, 384, /2 4×4 deconv, 64, ×2 4×4 conv, 384, /2 7×7 conv, 3 3.2 损失函数以及权重
损失函数包括两个分类器的损失
${Loss}_{AB}$ 、${Loss}_{RF} $ 和一个循环回来的图像之间的像素损失$ {Loss}_{G} $ . 设z为原始图像, G为生成网络, D是分类器, A为带印章图像, B对应的无印章图像. 原始的生成式对抗网络所使用的损失均为交叉熵损失, 在最近的研究中发现, 采用交叉熵损失存在梯度消失和模型坍塌的风险[21]. 相比于交叉熵损失, 最小二乘损失效果更佳[22]. 因此, SealGAN采用了最小二乘损失替换了交叉熵损失. 除此以外, 在计算生成图像和原始图像的内容损失时, 也不再使用常见的L2损失, 而是采用L1损失, 因为L2损失驱动下的生成图像更为模糊. 在训练分类器时, 损失函数如式(5)、(6)所示. 当训练生成网络时, 损失函数如式(7)、(8)、(9)所示. 其中,$ {Loss}_{AB} $ 和$ {Loss}_{RF} $ 的权重为1,$ {Loss}_{G} $ 的权重为10时, 生成效果最佳.$$ {Loss}_{AB}=\sum {\{\left(D\left({z}_{A}\right)\right)}^{2}+{\left(D\left({z}_{B}\right)-1\right)}^{2}\} $$ (5) $$ \begin{split} {Loss}_{RF}=&\sum \{{\left(D\left({z}_{A}\right)\right)}^{2}+{\left(D\left({z}_{B}\right)\right)}^{2}+ \\ &{\left(D\left(G\left({z}_{A}\right)\right)-1\right)}^{2}+{\left(D\left(G\left({z}_{B}\right)\right)-1\right)}^{2}\} \end{split}$$ (6) $$ {Loss}_{AB}=\sum \{{\left(D\left(G\left({z}_{A}\right)-1\right)\right)}^{2}+{\left(D\left(G\left({z}_{B}\right)\right)\right)}^{2}\} $$ (7) $$ {Loss}_{RF}=\sum \{{\left(D\left(G\left({z}_{A}\right)\right)\right)}^{2}+{\left(D\left(G\left({z}_{B}\right)\right)\right)}^{2}\} $$ (8) $$\begin{split} {Loss}_{G}=\,&\sum \{\left|{z}_{A}-{G}_{BtoA}\left({G}_{AtoB}\left({z}_{A}\right)\right)\right|+ \\ &\left|{z}_{B}-{G}_{AtoB}\left({G}_{BtoA}\left({z}_{B}\right)\right)\right|\} \end{split} $$ (9) 3.3 网络性能综合评价
为了能够更好地评价网络的性能, 本文提出针对印章消除的综合评价指标Evaluation-score (ES), 包含风格评价指标Class-score (CS1)和内容评价指标Content-score (CS2), 如下式所示.
$$ ES={\rm{\alpha }}\times{CS}_{1}+{\rm{\beta }}\times{CS}_{2} $$ (10) 式中的
$ {\rm{\alpha }} $ 和$ {\rm{\beta }} $ 是超参数, 可基于需求重要性进行调整, 本文实验中$ {\rm{\alpha }} $ 和$ {\rm{\beta }} $ 均为1.风格评价指标需要搭建一个分类网络, 采用的训练样本是真实的发票样本, 包含带印章和不带印章的样本. 本论文所搭建的分类网络借鉴GoogleNet网络结构[23], 损失函数采用余弦损失函数, 目标是让带印章样本的输出趋向于0, 不带印章样本的输出趋向于1. 采用余弦损失函数能够更好地增大类间距离, 公式如式(11)所示[24]. 将生成的图像通过训练好的分类网络, 所得到的结果即为CS1, 得分越高, 则消除印章的效果越好.
$$ {\rm{L}}_{cos}=\frac{1}{N}\sum _{i}-\log\frac{{{\rm{e}}}^{s\left({\rm{cos}}\left({\theta }_{i,i}\right)-m\right)}}{{{\rm{e}}}^{s\left({\rm{cos}}\left({\theta }_{i,i}\right)-m\right)}+\displaystyle\sum _{j\ne i}{{\rm{e}}}^{s{\rm{cos}}\left({\theta }_{i,j}\right)}} $$ (11) 内容评价指标则是利用生成图像和原带章图像之间的均方根误差计算的, 均方根误差越小, 意味着保留的内容越多. 定义的内容评价指标如下:
$$ {CS}_{2}=\frac{1}{{\rm{log}}\left(1+\sqrt{\dfrac{\displaystyle\sum {\left({x}_{i}-{y}_{i}\right)}^{2}}{N}}\right)} $$ (12) 3.4 结果与讨论
训练图像一共有六千余张发票图, 其中带印章样本和不带印章样本的比例约为1 : 1, 印章样本包含了出租车票印章、增值税发票印章、高速公路收费票据印章、定额发票印章以及各类企业单位加盖的公章. 在训练网络过程中, 为选择最佳的训练集和测试集比例, 对CycleGAN-ResNet网络和CycleGAN-UNet网络, 以及本文所提的SealGAN网络等三种网络, 采用不同的训练集和测试集划分比例, 包括7 : 3、8 : 2和9 : 1三种, 进行对比实验, 并基于网络性能综合评价指标ES进行评价. 三种网络不同数据集划分比例的综合评价指标如图7所示.
图 7 三种网络在不同数据集划分比例下的性能指标Fig. 7 Performance indices of three networks under different division proportion of data set由图7可知, 随训练集和测试集比例增大, 网络综合性能评价指标也在提升. 为了让每种网络效果达到最佳, 因此在论文实验中采用训练集和测试集比例为9 : 1分布进行实验比对.
本文分别使用了文献[5]的二次分割方法、采用ResNet和UNet生成网络的CycleGAN网络, 以及SealGAN网络, 进行发票印章消除实验. 训练模型时, 针对三种常见的优化器Momentum、RMSProp和Adam进行实验对比, 最终选择稳定性最佳的Adam优化器. 每种网络均训练了20个Epoch, Batch Size均为4. 学习率的初始值是0.0002, 每训练一个Epoch, 学习率下降10 %. 在CycleGAN模型中, 针对生成器和判别器不同的损失权重比, 包括1 : 1, 5 : 1, 10 : 1和15 : 1, 通过网络综合评价指标ES进行测试对比, 可知10:1为最佳损失权重比. 基于上述的参数设定进行实验, 图8是四种方法对发票印章消除的效果.
图 8 基于二次分割、CycleGAN-ResNet、CycleGAN-UNet和SealGAN的印章消除对比Fig. 8 Comparsion of effect of the seal elimination based on re-segmentation, CycleGAN-ResNet, CycleGAN-UNet and SealGAN由图中可知, 采用二次分割的方法效果最差, 由于方法的核心思想是对颜色进行阈值分割来获取印章部分再进行修改, 所以当发票颜色复杂时, 对印章部分的定位不完整, 导致印章消除不全, 而且对于获取到的印章部分采用直接填充背景方法, 对发票内容造成影响. 采用ResNet生成网络的CycleGAN, 生成的结果比较模糊, 色彩与原图相比, 色差巨大; 除此以外, 内容丢失也很严重. 采用UNet生成网络的CycleGAN, 由于UNet在进行上采样时融合了下采样层的特征, 保证了上采样过程中保留更多的原始信息, 因此色彩问题得到解决, 但是依旧有明显的内容丢失, 特别是当内容的颜色与印章的颜色相近时. 基于SealGAN网络的印章消除, 采用了下采样−精炼−上采样的生成结构, 而且上采样过程也融入了下采样层的原始信息, 因此色差问题基本不存在, 生成的结果相比于前两个网络的实验结果要清晰很多. 除此以外, 由于SealGAN并不是简单地用两个基础GAN网络搭建的循环结构, 而是用两个专门的分类器来替代判别器, 其中一个分类器专门区分真实样本和生成样本, 另一个专门区分带章样本和不带章样本. 这种替换使得网络能够更好地学习到哪些部分属于印章, 哪些部分不属于印章, 因此SealGAN网络可更好地保留非章内容.
三种不同网络消除发票印章的评价如表2所示, 可知SealGAN具有最高得分.
表 2 三种网络性能评价指标Table 2 Performance evaluation indices of three kinds of network网络类型 CS1 CS2 $ ES $ 二次分割 0.340 0.991 1.331 CycleGAN-ResNet 0.678 0.699 1.377 CycleGAN-UNet 0.703 0.679 1.382 SealGAN 0.699 0.740 1.439 4. 结语
本文针对发票上消除印章的问题, 分析了CycleGAN网络的优缺点, 并基于此进行相应改进, 提出SealGAN网络结构, 设计了包括风格评价和内容评价的网络综合评价指标. SealGAN网络用两个独立的分类器来取代原本的判别网络, 可降低单个分类器的分类要求, 提高网络的学习效果, 同时生成网络结合ResNet和Unet两种结构, 实现下采样−精炼−上采样的过程, 能够生成更加清晰的图像. 基于发票印章消除实验, 与原始的CycleGAN进行对比可看出, 本文所提出的SealGAN无论是在模糊程度、生成色彩等图像细节上, 还是内容的保留程度上都要明显优于原始的CycleGAN网络, 且综合评价指标高于CycleGAN网络, 实验结果表明了所提出的SealGAN网络的有效性. 未来可进一步提高非章内容的保留效果, 尤其是当印章与印章下的发票内容颜色很相近时, 并可针对发票复印件进行印章消除研究.
-
图 2 不同算法
$S_{\rm{NO},2}$ 和$D_{\rm{O},5}$ 的模型预测效果Fig. 2 Model prediction effects of
$S_{\rm{NO},2 }$ and$D_{\rm{O},5}$ with different algorithms
图 3 不同算法
$S_{\rm{NO},2}$ 和$D_{\rm{O},5}$ 的预测误差PDF曲线Fig. 3 Error PDF shapes of
$S_{\rm{NO},2 }$ and$D_{\rm{O},5}$ prediction with different algorithms
图 5 控制器参数累计频率分布曲线
Fig. 5 Cumulative frequency distribution curve of controller parameters
图 6 CFDL-MFAC控制器灵敏参数优化收敛过程
Fig. 6 The optimization convergence process of sensitive parameter based on CFDL-MFAC controller
表 1 不同算法的RMSE比较
Table 1 Comparison of RMSE based on different algorithms
算法 RMSE $(S_{\rm{NO},2 })$ RMSE $(D_{\rm{O},5})$ RBLSI 0.0165 0.0046 RLSI 0.0211 0.0085 
下载: 导出CSV
表 2 CFDL-MFAC控制器参数灵敏度分析结果
Table 2 Sensitivity analysis results of CFDL-MFAC controller parameters
序号 参数 含义 DS 不灵敏参数固定值 1 $\lambda$ 输入权重因子 0.9981 0.5 2 $\mu $ PJM权重因子 0.9985 0.6 3 $\eta $ PJM步长因子 0.9996 0.5 4 $\rho $ 输入步长因子 0.9378 1.0 5 $\alpha $ PJM重置限定参数 0.998 1.5 6 $b_1 $ PJM重置限定参数 0.9991 0.55 7 $b_2 $ PJM重置限定参数 0.9995 0.8 8 $\phi _{11}(0)$ PJM初值 0.5551 − 9 $\phi _{12}(0)$ PJM初值 −0.6339 − 10 $\phi _{21}(0)$ PJM初值 0.7289 − 11 $\phi _{22}(0)$ PJM初值 0.5779 −
下载: 导出CSV
表 3 GA-PSO算法参数
Table 3 GA-PSO algorithm parameters
前期 GA 参数 后期 PSO 参数 种群规模 $M=40$ 种群规模 $M=40$ 指定遗传代数 60 最大迭代次数 240 基因重组概率 0.7 加速系数 ${c_1} = 4,\;{c_2} = 2$ 变异概率 0.25 最大惯性
权重系数$W_{\max } =0.9$ 子代选择系数 T = 1,
${\alpha _T} = 0.5$最小惯性
权重系数$W_{\min } =0.1$
下载: 导出CSV
表 4 不同算法控制性能对比
Table 4 Comparison of control performance based on different algorithms
算法 RBL-MPC CFDL-MFAC 控制量平均求解时间 (s) 0.0987 0.00003528 方波扰动 RMSE $S_{\rm{NO},2} $ 0.0865 0.0297 $D_{\mathrm{O},5} $ 0.0158 0.0107 正弦扰动 RMSE $S_{\mathrm{NO},2} $ 0.0859 0.0249 $D_{\mathrm{O},5} $ 0.0099 0.0089
下载: 导出CSV
-
[1] 魏伟, 左敏, 李伟, 王小艺, 刘载文. 污水处理过程溶解氧浓度的自抗扰控制. 控制理论与应用, 2018, 35(1): 24−30 doi: 10.7641/CTA.2017.70229Wei Wei, Zuo Min, Li Wei, Wang Xiao-Yi, Liu Zai-Wen. Control of dissolved oxygen for a wastewater treatment process by active disturbance rejection control approach. Control Theory and Applications, 2018, 35(1): 24−30 doi: 10.7641/CTA.2017.70229 [2] Phuc B D H, You S S, Hung B M, Kim H S. Robust control synthesis for the activated sludge process. Environmental Science Water Research and Technology, 2018, 4: 992−1001 doi: 10.1039/C8EW00032H [3] 栗三一, 乔俊飞, 李文静, 顾锞. 污水处理决策优化控制. 自动化学报, 2018, 44(12): 2198−2209Li San-Yi, Qiao Jun-Fei, Li Wen-Jing, Gu Ke. Advanced decision and optimization control for wastewater treatment plants. Acta Automatica Sinica, 2018, 44(12): 2198−2209 [4] 宋贺达, 周平, 王宏, 柴天佑. 高炉炼铁过程多元铁水质量非线性子空间建模及应用. 自动化学报, 2016, 42(11): 1664−1679Song He-Da, Zhou Ping, Wang Hong, Chai Tian-You. Nonlinear subspace modeling of multivariate molten iron quality in blast furnace ironmaking and its application. Acta Automatica Sinica, 2016, 42(11): 1664−1679 [5] Virote B, Tharach J. Data-driven stochastic subspace identification of flutter derivatives of bridge decks. Journal of Wind Engineering and Industrial Aerodynamics, 2010, 98(12): 784−799 doi: 10.1016/j.jweia.2010.07.003 [6] Wahab N A Katebi R, Balderud J, Rahmat M F. Datadriven adaptive model-based predictive control with application in wastewater systems. IET Control Theory and Applications, 2011, 5(6): 803−812 doi: 10.1049/iet-cta.2010.0068 [7] Zhou P Dai P Song H D. Data-driven recursive subspace identification based online modelling for prediction and control of molten iron quality in blast furnace ironmaking. IET Control Theory and Applications, 2017, 11(14): 2343−2351 doi: 10.1049/iet-cta.2016.1474 [8] Zhou P, Song He D, Wang H, Chai T Y. Data-driven nonlinear subspace modeling for prediction and control of molten iron quality indices in blast furnace ironmaking. IEEE Transactions on Control Systems Technology, 2017, 25(5): 1761−1774 doi: 10.1109/TCST.2016.2631124 [9] Gijsvander V, Janwillem W , Bergamasco M, Lovera M, Verhaegen M. Closed-loop subspace identification methods: an overview. IET Control Theory and Applications, 2013, 7(10): 1339−1358 doi: 10.1049/iet-cta.2012.0653 [10] Qin S. J. An overview of subspace identification. Computers and Chemical Engineering, 2006, 30(2006): 1502−1513 [11] Samsudin S I, Rahmat M F, Norhaliza A W, Zulfatman H, Siti N S M, Mashitah C R. Two-time scales matrix decomposition for wastewater treatment plant. In: Proceedings of the 8th IEEE Eighth International Colloquium on Signal Processing and Its Applications. Melaka, Malaysia: IEEE, 2012. 347−351 [12] Dong N, Wu C H, Gao Z K, Chen Z Q, Ip W H. Data-driven control based on simultaneous perturbation stochastic approximation with adaptive weighted gradient estimation. IET Control Theory and Applications, 2016, 10(2): 201−209 doi: 10.1049/iet-cta.2015.0636 [13] Han H G, Zhang L, Qiao J F, Data-based predictive control for wastewater treatment process. IEEE Access, 2018, 6: 1498−1512 doi: 10.1109/ACCESS.2017.2779175 [14] Qiao J F, Zhou H B. Modeling of energy consumption and effluent quality using density peaks-based adaptive fuzzy neural network. IEEE/CAA Journal of Automatica Sinica, 2018, 5(5): 968−976 doi: 10.1109/JAS.2018.7511168 [15] Han H G, Qiao J F, Nonlinear model-predictive control for industrial processes: an application to wastewater treatment process. IEEE Transactions on Industrial Electronics, 2014, 61(4): 1970−1982 doi: 10.1109/TIE.2013.2266086 [16] Han H G, Zhang L, Hou Y, Qiao J F. Nonlinear model predictive control based on a self-organizing recurrent neural network. IEEE Transactions on Neural Networks and Learning Systems, 2016, 27(2): 402−415 doi: 10.1109/TNNLS.2015.2465174 [17] Wang D, Ha M M, Qiao J F. Self-Learning optimal regulation for discrete-time nonlinear systems under eventdriven formulation. IEEE Transactions on Automatic Control, 2019, 65(3): 1272−1279 [18] Wang D, He H B, Liu D R. Adaptive critic nonlinear robust control: a survey. IEEE Transactions on Cybernetics, 2017, 47(10): 3429−3451 doi: 10.1109/TCYB.2017.2712188 [19] 侯忠生. 非线性系统参数辨识、自适应控制和无模型学习自适应控制 [博士学位论文]. 东北大学, 中国, 1994.Hou Zhong-Sheng. The Parameter Identification, Adaptive Control and Model Free Learning Adaptive Control for Nonlinear Systems [Ph.D. dissertation], Northeastern University, China, 1994. [20] Hou Z S, Jin S T, Data-driven model-free adaptive control for a class of MIMO nonlinear discrete-time systems. IEEE Transactions on Neural Networks, 2011, 22(12): 2173−2188 doi: 10.1109/TNN.2011.2176141 [21] Hou Z S, Jin S T. A novel data-driven control approach for a class of discrete-time nonlinear systems. IEEE Transactions on Control Systems Technology, 2011, 19(6): 1549−1558 doi: 10.1109/TCST.2010.2093136 [22] Hou Z S, Bu X H. Model-free adaptive control with data dropouts. Express Systems with Applications, 2011, 38(8): 10709−10717 doi: 10.1016/j.eswa.2011.01.158 [23] Zhu Y M, Hou Z S. Controller dynamic linearisation-based model-free adaptive control framework for a class of nonlinear system. IET Control Theory and Applications, 2015, 9(7): 1162−1172 doi: 10.1049/iet-cta.2014.0743 [24] Bu X H, Qiao Y X, Hou Z S, Yang J Q. Model free adaptive control for a class of nonlinear systems using quantized information. Asian Journal of Control, 2018, 20(2): 962−968 doi: 10.1002/asjc.1610 [25] 侯忠生, 金尚泰. 无模型自适应控制: 理论与应用. 北京: 科学出版社, 2013. 34−122Hou Zhong-Sheng, Jin Shang-Tai. Model Free Adaptive Control: Theory and Application. Beijing: Science Press, 2013. 34−122 [26] Hou Z S, Chi R H, Gao H J. An overview of dynamiclinearization-based data-driven control and applications. IEEE Transactions on Industral Electronics, 2017, 64(5): 4076−4090 doi: 10.1109/TIE.2016.2636126 [27] 乔俊飞, 薄迎春, 韩广. 基于ESN 的多指标DHP控制策略在污水处理过程中的应用. 自动化学报, 2013, 39(7): 1146−1151Qiao Jun-Fei, Bo Ying-Chun, Han Guang. Application of ESN-based multi indices dual heuristic dynamic programming on wastewater treatment process. Acta Automatica Sinica, 2013, 39(7): 1146−1151 [28] Vincent V, Verhaegenz M. Identification of multivariable bilinear state space systems based on subspace techniques and separable least squares optimization. International Journal of Control, 2010, 74(18): 1824−1836 [29] Francesco V, Minh Q P, Raimondo B, Longman R W. Linear state representations for identification of bilinear discretetime models by interaction matrices. Nonlinear Dynamics, 2014, 77(4): 1561−1576 doi: 10.1007/s11071-014-1399-9 [30] Van B H. Theory of adjoint structures. AIAA Journal, 1976, 14(7): 977−979 doi: 10.2514/3.7174 [31] Wang D Y, Zeng P, Li L W, Wang K. A novel particle swarm optimization algorithm. In: Processings of the 2010 IEEE International Conference on Software Engineering and Service Sciences. Beijing, China: IEEE, 2010. 408−411 [32] Eberhart R C, Shi Y H. Particle swarm optimization: developments, applications and resources. In: Proceedings of the 2001 Congress on Evolutionary Computation. Seoul, South Korea: IEEE, 2001. 81−86 [33] 吉根林. 遗传算法研究综述. 计算机应用与软件, 2004, 21(2): 69−73 doi: 10.3969/j.issn.1000-386X.2004.02.032Ji Gen-Lin. Survey on genetic algorithm. Computer Applications and Software, 2004, 21(2): 69−73 doi: 10.3969/j.issn.1000-386X.2004.02.032 [34] Tang K S, Man K F, Kwong S, He Q H. Genetic algorithms and their applications. IEEE Signal Processing Magazine, 1996, 13(6): 22−37 doi: 10.1109/79.543973 -
下载:
计量
- 文章访问数: 810
- HTML全文浏览量: 147
- PDF下载量: 211
- 被引次数: 0
