2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于CPS框架的微粉生产过程多模型自适应控制

李晓理 王康 于秀明 苏伟

褚睿鸿, 王红军, 杨燕, 李天瑞. 基于密度峰值的聚类集成. 自动化学报, 2016, 42(9): 1401-1412. doi: 10.16383/j.aas.2016.c150864
引用本文: 李晓理, 王康, 于秀明, 苏伟. 基于CPS框架的微粉生产过程多模型自适应控制. 自动化学报, 2019, 45(7): 1354-1365. doi: 10.16383/j.aas.2018.c180387
CHU Rui-Hong, WANG Hong-Jun, YANG Yan, LI Tian-Rui. Clustering Ensemble Based on Density Peaks. ACTA AUTOMATICA SINICA, 2016, 42(9): 1401-1412. doi: 10.16383/j.aas.2016.c150864
Citation: LI Xiao-Li, WANG Kang, YU Xiu-Ming, SU Wei. CPS-based Multiple Model Adaptive Control of GGBS Production Process. ACTA AUTOMATICA SINICA, 2019, 45(7): 1354-1365. doi: 10.16383/j.aas.2018.c180387

基于CPS框架的微粉生产过程多模型自适应控制

doi: 10.16383/j.aas.2018.c180387
基金项目: 

北京市科技重大专项 Z181100003118012

国家重点研发计划项目 2018YFB1702704

国家重点研发计划项目 2018YFC1602704

国家自然科学基金 61873006

北京市科技新星交叉学科项目 Z161100004916041

国家自然科学基金 61673053

国家自然科学基金 61473034

详细信息
    作者简介:

    王康   北京工业大学信息学部博士后.分别于2012年和2018年获得北京科技大学学士学位和博士学位.主要研究方向为最优控制, 智能控制.E-mail:wangkang@bjut.edu.cn

    于秀明   中国电子技术标准化研究院软件工程与评估中心副主任.2004年获北京航空航天大学学士学位, 2013年获对外经济贸易大学硕士学位.主要研究方向为信息物理系统, 工业互联网以及两化融合.E-mail:yuxiuming@cesi.cn

    苏伟   中国电子技术标准化研究院助理工程师.分别于2014年和2016年获得延边大学学士学位和硕士学位.主要研究方向为信息物理系统, 智能制造, 两化融合, 工业互联网平台.E-mail:suwei@cesi.cn

    通讯作者:

    李晓理   北京工业大学信息学部教授.2000年获得东北大学博士学位.主要研究方向为复杂系统的建模优化与控制, 智能控制.本文通信作者.E-mail:lixiaolibjut@bjut.edu.cn

CPS-based Multiple Model Adaptive Control of GGBS Production Process

Funds: 

Beijing Major Science and Technology Special Projects Z181100003118012

National Key Research and Development Project 2018YFB1702704

National Key Research and Development Project 2018YFC1602704

National Natural Science Foundation of China 61873006

Beijing Nova Programme Interdisciplinary Cooperation Project Z161100004916041

National Natural Science Foundation of China 61673053

National Natural Science Foundation of China 61473034

More Information
    Author Bio:

       Postdoctor at Faculty of Information Technology, Beijing University of Technology. He received his bachelor degree and Ph. D. degree from the University of Science and Technology Beijing in 2012 and 2018, respectively. His research interest covers optimal control and intelligent control

       Deputy director of Software Engineering and Appraisal Center, China Electronic Standardization Institute. She received her bachelor degree from Beijing University of Aeronautics and Astronautics in 2004, and master degree from University of International Business and Economics in 2013. Her research interest covers cyber-physical system, industrial internet, integration of informatization and industrialization

       Assistant engineer at China Electronic Standardization Institute. He received his bachelor degree and master degree from Yanbian University in 2014 and 2016, respectively. His research interest covers cyber-physical system, intelligent manufacturing, integration of informatization and industrialization, industrial internet platform

    Corresponding author: LI Xiao-Li   Professor at the Faculty of Information Technology, Beijing University of Technology. He received his Ph. D degree from Northeastern University in 2000. His research interest covers modeling, control and optimization of complex system, and intelligent control. Corresponding author of this paper
  • 摘要: 针对矿渣微粉(Ground granulated blast-furnace slag,GGBS)生产这一多变量、强耦合、多工况的复杂非线性过程,本文根据大量生产数据,提炼出矿渣微粉生产过程的三个典型工况.求解多工况多目标优化问题以求得最优设定值.建立多工况下的递归神经网数据驱动模型,并采用自适应动态规划方法,建立多个控制器,结合加权多模型控制,实现矿渣微粉生产过程在多工况切换情况下的自适应控制.通过过程运行优化、跟踪控制优化、通讯、工业以太网等信息资源与矿渣微粉生产物理资源之间的融合,构建基于信息物理系统(Cyber-physical system,CPS)的矿渣微粉生产优化控制系统.实验分析表明,本文提出的基于CPS的多模型自适应控制器,能够有效实现多工况条件下矿渣微粉生产过程的自适应控制,减小超调量,提高控制品质.
  • 视频信号在捕捉、记录和传输的过程中都可能引入噪声.引入的噪声严重降低视频画面质量, 影响观众的视觉体验.而视频去噪是将数字视频图像的特点与现有的信号处理技术相结合, 尽可能地降低视频图像中噪声干扰的一种多媒体信息处理技术.目前存在的大多数视频去噪算法的性能都不同程度地依赖于含噪视频中的噪声参数.虽然通过各种各样的算法, 可以达到很理想的去噪效果, 但大多数算法都有个假设前提, 即噪声强度是预先知道的[1-3].人为给定噪声参数或者噪声参数不准确, 都会导致去噪效果不理想.因此对于含噪视频的噪声参数估计是视频去噪研究中一个关键性问题.近年来, 图像去噪算法在空域和变换域等取得了较大发展.其中较为优秀的是Knaus等的基于空域和频域的双域滤波图像去噪算法[4], 去噪后的图像细节信息丰富, 但是需要设置与噪声有关的经验参数. Pierazzo等[5]利用非局部贝叶斯去噪替换了文献[4]中的引导层, 构造图像快速去噪算法, 但图片的适用性不高.为了实现视频噪声去除, 肖进胜等[6]将文献[4]和三维块匹配算法相结合将图像去噪拓展到了视频去噪领域, 主客观效果较好, 但对未知噪声水平时鲁棒性较差.另外Knaus等[7]基于双域滤波引入了鲁棒噪声估计, 部分解决了人为设定噪声值的问题. Dabov等[8]对视频图像进行3D稀疏变换后再滤波(Video block-matching and 3-D filtering, VBM3D), 取得了较优秀的去噪性能.但该算法耗时高, 且去噪视频存在块效应. Matteo等提出了VBM4D[9]算法, 该算法将VBM3D拓展到4维结构, 能更好地保存视频的细节.上述去噪算法均没有对噪声进行有效的估计, 对于未知噪声的视频序列不能获得其最佳的去噪效果.

    高斯白噪声是视频图像采集中最常见的一种噪声, 针对该噪声的估计主要分为空域、时域和空时域[10]三种.空域法通常分为基于块[11-12], 基于滤波[13]和变换域三种方法[14]. Amer等[15]采用自适应平均的方法对所有块的方差进行选取和平均, 基于块进行噪声估计, 但该算法对平滑块的数量有严格的要求.基于图像块的算法估计结果受图像内容和噪声强度影响很大, 而Pyatykh等[16]提出的基于主成分分析(Principal component analysis, PCA)的方法则对含噪图像没有严格的要求, 且估计结果较精确.柳薇[17]利用PCA的思想对图像块进行噪声估计, 图像块的协方差矩阵最小特征值作为噪声方差的估计值, 该方法无需图像含有许多同种类区域. Aditya等[18]是一种基于奇异值分解的比较稳定的噪声估计方法, 用奇异值的尾部数据进行噪声强度估计, 降低图像信息对噪声估计的干扰.而时域的方法主要考虑帧与帧之间的关系, 对视频的整体运动较难把握, 因此需要进行运动检测或者运动补偿. Yin等[19]提出基于运动估计的视频噪声估计算法, 考虑到了视频编码, 该算法能对各类型的视频信号进行准确的估计.目前主要有算法[10, 20-21]利用了时空域的信息对视频噪声估计, Zlokolica等[20]主要用小波变换系数对空时域进行分析, 然而该方法计算复杂性较高. Ghazal等[10]利用5个域来探索空时域的信息, 每个域的局部相似性主要利用了高斯拉普拉斯算子, 该算法估计效果较好. Yang等[21]利用Sobel梯度作为同种类块的衡量标准, 利用了3个域进行噪声估计, 计算复杂性降低.总体上说, 噪声估计算法的效果都有待提高.

    考虑到PCA对含噪图像和噪声类型的鲁棒性, 本文提出一种基于PCA的分块视频噪声估计.本文所提出的方法有如下创新: 1)首先通过前后帧块匹配寻找相似块, 充分利用了视频序列的相关性, 并进行前后帧的差分运算以消除视频运动的影响. 2)使用正态分布函数作为选择弱纹理块的阈值函数, 使得计算复杂度降低, 同时简化了算法模型. 3)设置了明确的迭代指标, 使得最终结果更加精确.通过理论分析和实验结果表明本文的视频噪声估计算法具有较大的应用范围, 成功的运用于盲视频去噪.

    文献[14]提出的基于PCA的图像噪声水平估计算法, 对于纹理信息丰富的噪声图, 基于PCA的噪声估计会高估其噪声水平.若先选择噪声图中的弱纹理块, PCA则可以精确地估计出噪声水平, 因此本文算法首先选择含噪图中的弱纹理块, 再基于弱纹理块估计图像的噪声水平.含噪图像块模型为

    $ \begin{equation} \label{eq1} {y_i}={z_i}+{n_i} \end{equation} $

    (1)

    其中, $z_i$是原始图像第$i$个矢量块, $y_i$是观测到的矢量块, $n_i$是零均值高斯噪声.图像块可以认为是欧氏空间的数据, 利用PCA计算最小方差向量.协方差矩阵定义如下:

    $ \begin{equation} \label{eq2} \Sigma_y=\frac{1}{M}\sum\limits_{i=1}^{N}(y_i-\mu)(y_i-\mu )^{\rm T} \end{equation} $

    (2)

    其中, $M$是数据数目, $\mu$是数据集$y_i$的平均数.投影到最小方差方向数据的方差与协方差矩阵最小特征值相等, 可得到下述公式:

    $ \begin{equation} \label{eq3} \lambda_{\min}(\Sigma_y)=\lambda_{\min}(\Sigma_z)+\sigma_n^2 \end{equation} $

    (3)

    $\Sigma_y$和$\Sigma_z$分别指噪声块y和不含噪块z的协方差矩阵, $\lambda_{\min}$表示最小特征值.对于一般的图像块, $\lambda_{\min}(\Sigma_z)$是未知的.弱纹理块只适用于低维子空间, 它的协方差矩阵的最小特征值约为零, 对于弱纹理块噪声水平$\hat{\sigma }_{n}^{2}$可以简化为

    $ \begin{equation} \label{eq4} \hat{\sigma }_{n}^{2}={{\lambda}_{\min }}({{\Sigma }_{y'}}) \end{equation} $

    (4)

    ${{\sum }_{y'}}$为噪声图中所有弱纹理块的协方差矩阵.因此, 只要从噪声图像中选择出弱纹理块就能估计出噪声水平.梯度协方差矩阵能反应出图像的纹理信息, $N\times N$的含噪图像块$n$表示成列向量形式为${{N}^{2}}\times 1$.若${{D}_{h}}$和${{D}_{v}}$是由3阶滤波算子构造的水平和垂直方向的算子, 均为${{N}^{2}}\times {{N}^{2}}$的常对角矩阵, 从而梯度矩阵为

    $ \begin{equation} \label{eq5} {{G}_{n}}=[{{D}_{h}}n\quad {{D}_{v}}n] \end{equation} $

    (5)

    梯度协方差矩阵$C_n=G_n^{\rm T} G_n$的期望为

    $ \begin{align} \label{eq6} {\rm E}({{C}_{n}})=\left[\begin{matrix} {\rm E}({{n}^{\rm T}}D_{h}^{\rm T}{{D}_{h}}n)&0 \\ 0&{\rm E}({{n}^{\rm T}}D_{v}^{\rm T}{{D}_{v}}n) \\ \end{matrix} \right] \end{align} $

    (6)

    含噪块的梯度矩阵如下(其中${z_f}$是不含噪的平坦块)

    $ \begin{align} \label{eq7} {{G}_{y}}=\, &[{{D}_{h}}({{z}_{f}}+n)\ {{D}_{v}}({{z}_{f}}+n)]=\nonumber\\ &[D_h n \ D_v n]=G_n \end{align} $

    (7)

    对角线的元素有相同的特性, 设

    $ \begin{equation} \label{eq8} \varepsilon(n) ={{n}^{\rm T}}D_{h}^{\rm T}{{D}_{h}}n+{{n}^{\rm T}}D_{v}^{\rm T}{{D}_{v}}n \end{equation} $

    (8)

    $\varepsilon (n)$的生成函数决定了其分布[15], 且$\varepsilon(n)$的生成函数与伽马(Gamma)分布的生成函数形式一致, 因此对应的Gamma分布的形状参数$\alpha$和尺度参数$\beta $, 分别为

    $ \begin{equation} \label{eq9} \begin{cases} \alpha =\displaystyle\frac{N^2}{2} \\[2mm] \beta =\displaystyle\frac{2}{N^2} \sigma_n^2{\rm tr}(D_h^{\rm T} D_h+D_v^{\rm T} D_v) = 2\sigma_n^2 \\ \end{cases} \end{equation} $

    (9)

    该算法主要通过Gamma函数来逼近从而求得最终的阈值函数, 但Gamma函数形式较复杂, 不利于广泛的应用.同时原算法文献[14]中设置经验的迭代次数作为迭代停止条件, 虽然效果不错, 但始终存在适应性[22-23]的隐患.

    本文视频噪声估计算法首先通过前后帧块匹配寻找相似块, 对匹配效果最佳的块(弱纹理块)进行噪声估计, 并进行前后帧的差分运算以消除视频运动的影响.其次修改了选择弱纹理块的阈值函数, 使得判断标准更加精确, 减少图像纹理信息对估计结果的干扰.

    视频图像序列前后两帧图像在时间上具有很强的相关性.原本拥有很多纹理和细节信息的块, 如果它在前后两帧具有较强的相关性, 那么差分的结果仍然会产生一个平滑块, 降低了图像信息对噪声估计的影响, 同时降低了噪声水平对估计结果的影响.算法首先对当前帧图像进行块划分, 然后对于每一图像块在前一帧图像内进行块匹配寻找相似块, 并根据最小代价得到最佳匹配块, 所有块与最佳匹配块的差值块合并得到帧间差分图像, 接下来将差分图像作为原始含噪图像进行噪声估计.较朱磊等[24]对整幅图像进行差分运算, 本文算法的效率和精确度得到了提高.

    假设观察到的视频图像为

    $ \begin{equation} \label{eq10} I(i, j, n)=S(i, j, n)+N(i, j, n) \end{equation} $

    (10)

    其中, $S(i, j, n)$为原始不含噪的视频帧, $N(i, j, n)$为噪声信号, $n$为帧编号, $i$, $j$是像素坐标.经过块匹配, 我们选择的是匹配效果最佳的块(弱纹理块)进行噪声估计.通过选择弱纹理块, 降低了图像纹理对噪声估计的干扰, 对于最佳匹配块有:

    $ \begin{equation} \label{eq11} S(i, j, n)=S({{i}^\prime}, {{j}^\prime}, n+1) \end{equation} $

    (11)

    此时的差分图像可以表示为:

    $ \begin{align} \label{eq12} D(i, j, n)=\, &I(i, j, n)-I({{i}^\prime}, {{j}^\prime}, n+1)= \nonumber\\ &N(i, j, n)-N({{i}^\prime}, {{j}^\prime}, n+1) \end{align} $

    (12)

    假如原始图像的实际噪声方差为$\sigma _N^2$, 根据概率计算式(12)中差分图像的噪声方差$\sigma _D^2$有:

    $ \begin{equation} \label{eq13} \sigma _{D}^{2}=2\sigma _{N}^{2} \end{equation} $

    (13)

    因此, 如果原始图像中前后两帧的最佳匹配块足够接近, 那么差分图像块将会变为平滑块, 该平滑块的方差可以认为是噪声造成的, 可以用式(13)进行噪声方差的估计.

    对于伽马(Gamma)分布, 通过如下图 1所示可以得到:当形状参数$\alpha $越大其峰值越远离$Y$轴, 同时越逼近正态分布.当$\alpha$为正整数时, 分布可看作$\alpha $个独立的指数分布之和, 当$\alpha $趋向于较大数值时, 分布近似于正态分布.普通参数的Gamma分布可以看作多个独立的Gamma分布之和.按照中心极限定理, 独立同分布的随机变量之和趋于正态分布[25].

    图 1  Gamma概率密度图
    Fig. 1  Gamma probability density function

    图 1, $\alpha $, $\beta$的变化也可看出当$\alpha$足够大时Gamma分布近似正态分布.通过计算可以得到:

    $ \begin{eqnarray*} \label{eq14} \underset{\alpha \to \infty }{\mathop{\lim}}\Gamma(\alpha, \beta ) =N(\mu, {{\sigma }^{2}}), \mu=\alpha \beta, \sigma ^2=\alpha \beta ^2 \end{eqnarray*} $

    (14)

    由于梯度算子的大小只与块的大小有关, 不受含噪图像的影响, 因此当用3阶滤波算子时, 水平和垂直梯度算子为常矩阵, $D_{h}^{\rm T}{{D}_{h}}+D_{v}^{\rm T}{{D}_{v}}$也为常矩阵.本文中, $\alpha$值为${N}^{2}/{2}$, 块的大小为5, 可以用正态分布近似Gamma分布, 于是由式(9)和(14)可得到正态分布的均值和方差为

    $ \begin{equation} \label{eq15} \mu =\sigma_{n}^{2}{{N}^{2}}, {{\sigma }^{2}}=2{{N}^{2}}\sigma _{n}^{4} \end{equation} $

    (15)

    可得到本文的阈值为

    $ \begin{equation} \label{eq16} \tau =\sigma _n^2{{F}^{-1}}(\delta, {{N}^{2}}, \sqrt{2}N) \end{equation} $

    (16)

    其中, $F^{-1}(\delta, \alpha, \beta)$是正态分布的逆累加分布函数, $\delta$是人为给定的显著性水平(本文为0.99).当显著性水平和噪声水平给定时, 阈值随之确定.这里经过推导简化了Gamma分布的形状参数$\alpha $和尺度参数$\beta$, 同时用正态分布替换了Gamma分布, 对算法的运算函数也进行了简化.而当图像块梯度协方差矩阵的最大特征值小于上述阈值时, 即为选定的弱纹理块, 再对图像中的弱纹理区域利用PCA进行噪声估计.实验结果表明当块大小小于5时, 同样可以达到较好的估计效果.改进的正态分布函数和Gamma分布的噪声估计准确度(单位: dB)对比如下表所示:

    表 1  不同分布函数噪声估计对比
    Table 1  Comparison of noise estimation for different function
    Noise level (dB) Lena Akiyo Bus Coastguard
    ${{\sigma }_{n}}=10$ Liu等[14] 9.68 9.88 9.69 9.79
    ${{\sigma }_{n}}=10$ Proposed 9.86 9.79 9.97 9.93
    ${{\sigma }_{n}}=20$ Liu等[14] 19.56 19.61 19.64 19.67
    ${{\sigma }_{n}}=20$ Proposed 19.72 19.65 19.78 19.74
    ${{\sigma }_{n}}=30$ Liu等[14] 29.54 29.54 29.13 29.59
    ${{\sigma }_{n}}=30$ Proposed 29.65 29.81 29.50 29.67
    ${{\sigma }_{n}}=40$ Liu等[14] 38.43 39.34 38.99 39.37
    ${{\sigma }_{n}}=40$ Proposed 39.50 39.38 39.64 39.60
    下载: 导出CSV 
    | 显示表格

    由上表可以看出, 替换后的函数估计效果更好.因此噪声估计采用的分布函数改为正态分布完全可行且取得更好的估计效果.大量实验结果表明迭代次数达到6次时, 噪声水平可基本达到稳定.但是为了进一步提高噪声估计算法的稳定性, 对于迭代次数的设置, 本文采用经验次数和前后两次噪声估计水平的差值比例也可确定来决定迭代停止条件.下面将讨论一下针对不同噪声水平, 改进算法相对于其他算法噪声估计的效果.当噪声强度小于10 dB时, 大部分图像算法的有效性都会降低, 但Shin等[26]算法例外, 它仅在低噪环境下工作稳定可靠.由于Shin等提出的自适应高斯滤波的算法只适用于低噪环境, 因此仅在$\sigma<15$的情况下比较Liu等[14]、Pyatykh等[16]、Shin等[26]和本文算法的实验结果.在低噪声强度下噪声估计对比结果如图 2所示, 噪声估计的误差定义为$\delta (\sigma )=|\hat{\sigma }-\sigma |$ (dB).分别选取了Lena (含有细节较少)和Baboon (含有细节较多)等图片进行测试.由图可以看出本文算法比Pyatykh等[16]、Aditya等[18]和Shin等[26]误差较小.对于Liu等[14], 本文算法在$\sigma <10$时也明显误差较小.当图像含有复杂的纹理结构, 图像块的最小特征值大于0, 估计误差较大.由图 2 (a)(b)对比看出, 细节较少的图片噪声估计更为精确.本文算法的噪声强度估计结果对细节较少和细节较多的图片均较对比算法精确, 说明本文的噪声估计算法在低噪条件下有较高的准确性.

    图 2  在低噪声强度下估计效果对比
    Fig. 2  Comparisons in low noise case
    图 3  在高噪声强度下估计效果对比
    Fig. 3  Comparisons in high noise case

    在较高噪声水平情况下, Pyatykh等[16]是目前广为应用的基于主成分分析的噪声估计方法, 有较高的精确性和较快的速度.由曲线图可以看出, Pyatykh等[16]估计误差随着噪声水平增加, 当误差大于40 dB时逐渐减小, Aditya等[18]估计误差相对较稳定, Liu等[14]算法估计误差较小, 本文算法估计误差比这3种算法小.因此, 本文算法在高噪声水平下对噪声的估计结果较其他算法精确, 受噪声强度的影响较小.通常一段视频噪声的分布是均匀的, 每帧的噪声水平是类似的.考虑到视频噪声的突变情况, 本文在对视频进行噪声估计时, 利用前后帧信息和PCA得到当前帧的噪声水平, 因此当下一帧出现噪声突变时, 对估计结果没有影响.

    本文视频噪声估计算法首先通过前后帧块匹配寻找相似块, 充分利用了视频序列的相关性, 并进行前后帧的差分运算以消除视频运动的影响.其次修改了选择弱纹理块的阈值函数, 用正态分布函数简化了计算复杂度降低.最后设置了明确的迭代指标, 使得估计结果更加精确.本文的视频噪声水平估计过程如下图 4所示:

    图 4  本文算法迭代噪声水平估计流程图
    Fig. 4  Flowchart of the iterative noise level estimation for proposed algorithm

    为验证本文算法对视频进行噪声估计的效果, 使用两个CIF格式的标准测试视频, 两个测试视频为: Flowergarden和Football.算法运行环境为Windows XP, CPU-Intel Core i5-2500 K, 主频3.30 GHz, 内存3 GB, 32 bit.由于Shin等[26]数据取自原文, 因此本文对这两个测试视频同样加入均值为0, 标准差为20, 30, 40的高斯白噪声与算法Amer等[15]、Pyatykh等[16]、Liu等[14]进行对比. 图 5~图 7是本文算法和其他算法在视频序列前40帧不同噪声水平的估计误差, 对于不同的视频序列和噪声水平, 本文算法的估计误差较低.

    图 5  加噪20 dB视频序列(Flower, Football)估计误差
    Fig. 5  Noise estimation error for 20 dB noisy sequences (Flower, Football)
    图 6  加噪30 dB视频序列(Flower, Football)估计误差
    Fig. 6  Noise estimation error for 30 dB noisy sequences (Flower, Football)
    图 7  加噪40 dB视频序列(Flower, Football)估计误差
    Fig. 7  Noise estimation error for 40 dB noisy sequences (Flower, Football)

    PSNR和SSIM是两个比较常用的评价去噪效果的客观指标, 本文选用这两个指标来对VBM3D算法、PID算法、VBM4D算法和文献[6]算法加噪声估计进行对比.由于VBM3D和VBM4D未加入噪声估计, 因此测试程序中设定随机取真实值附近的噪声水平进行视频去噪.不同噪声水平下各算法的去噪效果的PSNR和SSIM对比如表 2所示:

    表 2  VBM3D、PID、VBM4D和本文算法的PSNR和SSIM对比
    Table 2  The comparisons of PSNR and SSIM results of VBM3D, PID, VBM4D and proposed algorithm
    Noise level (dB) Algorithm Akiyo PSNR/SSIM Mobile PSNR/SSIM Flowergarden PSNR/SSIM Foreman PSNR/SSIM Football PSNR/SSIM
    ${{\sigma }_{n}}=10$ VBM3D 35.488/0.877 32.374/0.954 34.250/0.984 34.313/0.902 33.048/0.951
    ${{\sigma }_{n}}=10$ PID 31.396/0.763 29.538/0.917 31.276/0.962 31.094/0.844 30.048/0.921
    ${{\sigma }_{n}}=10$ VBM4D 30.290/0.730 29.727/0.915 29.972/0.842 30.078/0.813 29.833/0.923
    ${{\sigma }_{n}}=10$ Proposed 37.842/0.944 32.454/0.972 33.675/0.983 36.295/0.938 33.488/0.960
    ${{\sigma }_{n}}=20$ VBM3D 30.239 /0.715 27.653/0.897 29.444/0.961 29.587/0.786 28.171/0.868
    ${{\sigma }_{n}}=20$ PID 26.011/0.836 24.572/0.833 25.945/0.896 25.937/0.675 24.889/0.793
    ${{\sigma }_{n}}=20$ VBM4D 29.594/0.767 27.249/0.898 27.599/0.852 28.904/0.797 27.727/0.885
    ${{\sigma }_{n}}=20$ Proposed 34.986/0.925 28.322/0.936 29.386/0.957 33.098/0.878 29.890/0.894
    ${{\sigma }_{n}}=30$ VBM3D 27.463/0.609 25.025/0.853 26.611/0.932 26.934/0.698 25.406/0.783
    ${{\sigma }_{n}}=30$ PID 23.051/0.815 21.803/0.769 22.895/0.825 23.017/0.565 21.929/0.676
    ${{\sigma }_{n}}=30$ VBM4D 29.026/0.780 25.884/0.883 26.185/0.848 28.208/0.782 26.557/0.842
    ${{\sigma }_{n}}=30$ Proposed 32.579/0.866 26.045/0.898 27.098/0.932 31.379/0.835 27.995/0.825
    下载: 导出CSV 
    | 显示表格

    表 2中的对比结果可以看出, 本文算法的PSNR和SSIM比VBM3D算法在大多数情况下高出很多, 主要是因为VBM3D处理过后像素值减小, 图像的亮度降低, 导致PSNR与SSIM均较低.对于Flowergarden视频序列, 本文客观效果在噪声水平较低情况下比VBM3D略低, 主要因为在对视频去噪时, 本文指定的是真实噪声水平利用VBM3D进行测试, 而实际情况中并不可能准确获取图像的真实噪声水平. PID算法中加入了噪声估计, 但综合而言, 本文的客观指标在各种噪声水平均比PID略高.上述表格表明, 本文算法与VBM3D、PID、VBM4D进行对比时, 客观效果优势较明显.因此说明本文算法不仅去除了图像中的噪声, 而且较好地保持了图像本身结构的信息, 去噪效果较为优秀, 同时加入了噪声估计应用范围更广.

    本文提出了一种新颖的视频噪声估计算法, 充分利用了视频序列的相关性, 利用帧间进行相似块的搜索, 基于最小代价准则获得帧间的差分图像, 消除了视频运动的影响, 得到初步的弱纹理差分图像:引入基于块的噪声估计, 能自适应地获取噪声水平参数, 并提出正态分布函数作为文中选择弱纹理块的阈值函数, 降低了计算复杂度, 另设置了明确的迭代指标使得估计的噪声水平更加精确.最后考虑像素会受到噪声饱和效应的影响, 避免了在高噪声水平下的低估现象.通过理论分析和实验结果表明本文提出的视频噪声估计算法估计精确, 可运用到盲视频去噪领域, 具有广阔的应用前景.


  • 本文责任编委 乔俊飞
  • 图  1  信息物理系统

    Fig.  1  Cyber-physical system

    图  2  矿渣微粉生产系统流程图

    Fig.  2  Flow chart of GGBS production process

    图  3  优化控制结构图

    Fig.  3  Structure of optimal contol

    图  4  矿渣微粉生产最优控制系统CPS硬件结构

    Fig.  4  The CPS hardware structure of GGBS production optimal control system

    图  5  微粉参数作用机理

    Fig.  5  Interaction among parameters of GGBS production process

    图  6  多目标优化设定值流程

    Fig.  6  Flow chart of set-point optimization using multi-objective optimization algorithm

    图  7  概率加权多模型ADP结构

    Fig.  7  Structure of weighted multiple model ADP

    图  8  矿渣微粉生产过程CPS框架

    Fig.  8  The CPS structure of GGBS production process

    图  9  Pareto最优解集

    Fig.  9  Obtained Pareto solutions

    图  10  采用控制器3时的质量曲线

    Fig.  10  The quality curve using controller 3

    图  11  多模型ADP微粉质量曲线

    Fig.  11  The quality curve using multiple model ADP

    图  12  多模型ADP控制输入曲线

    Fig.  12  The curve of control using multiple model ADP

    图  13  多模型ADP工况切换识别

    Fig.  13  Working condition identification using the multiple model ADP algorithm

    图  14  矿渣微粉粉磨系统运行数据

    Fig.  14  Operation data of GGBS production system

    表  1  各控制变量允许变化范围

    Table  1  Permitted range for each variable

    名称 变量 最小值 最大值 单位
    喂料量 $u_1$ 75 115 $10^3 {\rm kg/h}$
    选粉机转速 $u_2$ 850 1 250 ${\rm r/min}$
    入磨风温 $u_3$ 190 300
    冷风阀开度 $u_4$ 30 95 $\%$
    下载: 导出CSV

    表  2  微粉厂3号矿渣微粉生产线生产工况1运行数据

    Table  2  Process data for GGBS production line 3 in condition 1

    编号 喂料量 电机转速 入磨风温 冷风阀开度 比表面积 产品产量
    ($10^3$ kg/h) (r/min) (℃) (%) ($\rm{m}^2$/kg) ($10^3$ kg/h)
    1 101.76 1 090.31 240.71 63.46 451.75 95.12
    2 103.44 1 089.70 241.10 64.96 436.13 102.17
    3 108.15 1 099.65 265.08 57.08 419.30 106.16
    $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$
    198 108.95 1 089.32 249.54 65.82 435.96 102.49
    199 103.84 1 089.06 244.98 66.06 431.47 102.40
    200 102.12 1 119.46 266.82 61.03 427.28 92.61
    下载: 导出CSV

    表  3  微粉厂3号矿渣微粉生产线生产工况2运行数据

    Table  3  Process data for GGBS production line 3 in condition 2

    编号 喂料量 电机转速 入磨风温 冷风阀开度 比表面积 产品产量
    ($10^3$ kg/h) (r/min) (℃) (%) ($\rm{m}^2$/kg) ($10^3$ kg/h)
    1 84.24 1 249.09 234.88 63.80 439.75 79.33
    2 86.89 1 251.25 231.13 66.85 431.42 82.05
    3 82.54 1 249.28 241.51 63.24 428.76 76.55
    $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$
    198 84.94 1 159.96 228.21 65.72 427.26 80.49
    199 85.64 1 239.06 242.95 62.55 439.88 78.54
    200 85.42 1 248.68 233.01 69.12 424.21 78.75
    下载: 导出CSV

    表  4  微粉厂3号矿渣微粉生产线生产工况3运行数据

    Table  4  Process data for GGBS production line 3 in condition 3

    编号 喂料量 电机转速 入磨风温 冷风阀开度 比表面积 产品产量
    ($10^3$ kg/h) (r/min) (℃) (%) ($\rm{m}^2$/kg) ($10^3$ kg/h)
    1 104.09 1 015.28 216.60 61.44 435.62 96.91
    2 104.66 998.94 250.78 57.66 423.72 95.04
    3 102.15 1 000.03 237.65 59.70 445.55 94.23
    $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$ $\vdots$
    198 92.77 1 020.70 263.63 55.51 433.80 83.13
    199 106.54 1 011.66 236.26 62.65 443.47 97.64
    200 93.11 1 009.24 209.42 56.66 426.29 87.44
    下载: 导出CSV
  • [1] Saranya P, Nagarajan P, Shashikala A P. Eco-friendly GGBS concrete:a state-of-the-art review. IOP Conference Series:Materials Science and Engineering, 2018, 330(1):012057
    [2] Li X L, Jia C, Liu D X, Ding D W. Nonlinear adaptive control using multiple models and dynamic neural networks. Neurocomputing, 2014, 136:190-200 doi: 10.1016/j.neucom.2014.01.013
    [3] Li X L, Jia C, Wang K, Wang J. Trajectory tracking of nonlinear system using multiple series-parallel dynamic neural networks. Neurocomputing, 2015, 168:1-12 doi: 10.1016/j.neucom.2015.06.024
    [4] Wei Q L, Song R Z, Yan P F. Data-driven zero-sum neuro-optimal control for a class of continuous-time unknown nonlinear systems with disturbance using ADP. IEEE Transactions on Neural Networks and Learning Systems, 2016, 27(2):444-458 doi: 10.1109/TNNLS.2015.2464080
    [5] 王康, 李晓理, 贾超, 宋桂芝.基于自适应动态规划的矿渣微粉生产过程跟踪控制.自动化学报, 2016, 42(10):1542-1551 http://www.aas.net.cn/CN/abstract/abstract18941.shtml

    Wang Kang, Li Xiao-Li, Jia Chao, Song Gui-Zhi. Optimal tracking control for slag grinding process based on adaptive dynamic programming. Acta Automatica Sinica, 2016, 42(10):1542-1551 http://www.aas.net.cn/CN/abstract/abstract18941.shtml
    [6] 信息物理系统白皮书(2017).中国电子技术标准化研究院, 2017

    Cyber-physical systems white paper (2017). China Electronics Standardization Institute, 2017
    [7] Zhao H Z, Sun D H, Yue H, Zhao M, Cheng S L. Using CSTPNs to model traffic control CPS. IET Software, 2017, 11(3):116-125 doi: 10.1049/iet-sen.2016.0119
    [8] 马大中, 胡旭光, 孙秋野.基于大维数据驱动的油气管网泄漏监控模糊决策方法.自动化学报, 2017, 43(8):1370-1382 http://www.aas.net.cn/CN/abstract/abstract19111.shtml

    Ma Da-Zhong, Hu Xu-Guang, Sun Qiu-Ye. A large dimensional data-driven fuzzy detection method for oil-gas pipeline network leakage. Acta Automatica Sinica, 2017, 43(8):1370-1382 http://www.aas.net.cn/CN/abstract/abstract19111.shtml
    [9] Li D, Zhan M Y, Zhang X Z, Fang Z P, Liu H Q. ISAR imaging of nonuniformly rotating target based on the multicomponent CPS model under low SNR environment. IEEE Transactions on Aerospace and Electronic Systems, 2017, 53(3):1119-1135 doi: 10.1109/TAES.2017.2667538
    [10] Zhang Y, Qiu M K, Tsai C W, Hassan M M, Alamri A. Health-CPS:healthcare cyber-physical system assisted by cloud and big data. IEEE Systems Journal, 2017, 11(1):88-95 doi: 10.1109/JSYST.2015.2460747
    [11] Khan M U, Li S, Wang Q X, Shao Z L. CPS oriented control design for networked surveillance robots with multiple physical constraints. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 2016, 35(5):778-791 doi: 10.1109/TCAD.2016.2524653
    [12] Higuera-Toledano M T, Risco-Martin J L, Arroba P, Ayala J L. Green adaptation of real-time web services for industrial CPS within a cloud environment. IEEE Transactions on Industrial Informatics, 2017, 13(3):1249-1256 doi: 10.1109/TII.2017.2693365
    [13] Chai T Y, Wu Z W, Wang H. A CPS based optimal operational control system for fused magnesium furnace. IFAC-PapersOnLine, 2017, 50(1):14992-14999 doi: 10.1016/j.ifacol.2017.08.2566
    [14] Pal R, Prasanna V. The STREAM mechanism for CPS security:the case of the smart grid. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 2017, 36(4):537-550 doi: 10.1109/TCAD.2016.2565201
    [15] Wang K, Li X L. Multiple set-points tracking control based on online ADP. In: Proceedings of the 2016 Chinese Control and Decision Conference (CCDC). Yinchuan, China: IEEE, 2016. 1214-1219
    [16] Zhang X, Zhang H G, Sun Q Y, Luo Y H. Adaptive dynamic programming-based optimal control of unknown nonaffine nonlinear discrete-time systems with proof of convergence. Neurocomputing, 2012, 91:48-55 doi: 10.1016/j.neucom.2012.01.025
    [17] Wang K, Li X L, Jia C, Yang S X, Li M Q, Li Y. Multiobjective optimization of the production process for ground granulated blast furnace slags. Soft Computing, DOI: 10.1007/s00500-017-2761-x
  • 期刊类型引用(10)

    1. 段现银,邵宇轩,彭芳瑜,向峰,杨岑岑,周林. 基于数字孪生的复杂件多机器人智能加工方法. 华中科技大学学报(自然科学版). 2024(06): 1-9 . 百度学术
    2. 马洪浩,刘钊,王孝红,李凡军. 水泥粉磨细度控制研究进展. 济南大学学报(自然科学版). 2024(05): 634-643 . 百度学术
    3. 毕素环,蒋一翔,于树松,丁香乾,牟亮亮,王彬. 基于强化学习的减少烘丝过程中烟丝“干头”量的方法. 自动化学报. 2023(08): 1679-1687 . 本站查看
    4. 赵耕,柳军,孙文权,张哲,刘向国,潘健华,郭瑞春. 人工智能在钢铁能源管控中的应用. 金属世界. 2022(03): 29-35 . 百度学术
    5. 张岩,吴鲲魁,孙瑞琪,刘洪宇. 基于信息物理系统的轧制过程能源负荷时间序列预测. 冶金自动化. 2022(03): 42-47 . 百度学术
    6. 于子航,王改云. 基于近似Q-学习算法的数据驱动控制仿真. 计算机仿真. 2022(05): 344-347+379 . 百度学术
    7. 徐亦卿,陆海澎. 基于扩张状态观测器的发动机转速双闭环自适应控制系统设计. 计算机测量与控制. 2022(07): 123-128 . 百度学术
    8. 解治宇,马连成,张岩,王怀远,钟浩. 人工智能技术在金属露天矿能源管控中的应用探讨. 冶金与材料. 2022(04): 100-104+107 . 百度学术
    9. 曹朝阳,吴庆涛. 信息数据融合技术支持下的自动化制造管理系统设计. 制造业自动化. 2020(05): 125-128 . 百度学术
    10. 朱飞燕. 多模型自适应加权的无线网络安全风险评估. 国外电子测量技术. 2020(06): 106-109 . 百度学术

    其他类型引用(4)

  • 加载中
  • 图(14) / 表(4)
    计量
    • 文章访问数:  2554
    • HTML全文浏览量:  323
    • PDF下载量:  489
    • 被引次数: 14
    出版历程
    • 收稿日期:  2018-05-31
    • 录用日期:  2018-09-12
    • 刊出日期:  2019-07-20

    目录

    /

    返回文章
    返回