,
,
,
，
1, 2, 3].如果过程已经有大量的运行数据信息，就可以利用合适的建模方法建立过程的确定性模型来反映过程的输入输出特性，那么利用过程的确定性模型来求解生产设计问题是一种可行方案[3, 4, 5].近年来，生产设计技术受到人们的关注，并取得了一些研究成果[6]，其中，文献[1]提出的零空间概念和潜变量建模方法对求解生产设计问题具有更重要的意义.实际上零空间是模型输入潜变量空间的子空间[7, 8, 9]，它代表了一个过程操作窗口，对应于生产相同等级产品的一系列过程操作条件.在此基础上，文献[10]提出了单变量零空间分析方法，通过每对过程变量间的关系图来确定一个生物过程的运行区域.文献[11, 12]使用同样的方法，利用一个催化裂化(FCC)过程的历史数据确定过程的一个由PCA模型得分所定义的过程运行区域.但是，实际工业过程中的变量数目通常非常大，而且运行数据的内部相关性需要加以考虑，因而这种方法的应用受到限制.目前，潜变量建模及零空间特性的研究取得进展[13, 14]，但用于多变量过程生产设计问题的相关研究鲜见报道.本文从求解生产设计问题的标准回归模型的对比分析入手，深入研究了基于主元回归的潜变量模型的应用策略，并利用实际过程数据仿真分析了具有零空间特性的多变量复杂工业过程的生产设计问题.
1 基于标准回归模型的生产设计
标准回归分析方法是将输入矩阵X和输出矩阵Y经过适当的变换，然后建立操作条件X和输出产品特性Y之间的关系：
Y可以通过最小二乘法(OLS)得到.当X满足m>n时会导致对XTX求逆出现困难，可使用更稳定的修正后的最小二乘估计法，例如岭回归分析[15].
假设已知β的估计值
，建立式(1)的标准回归模型，从期望输出ydesT出发，估计与其对应的过程输入变量
由于式(1)的标准回归模型只定义了X和Y间的外部映射关系，没有考虑输入变量间的内部关联结构，因此不能保证新操作条件xnewT与历史运行策略保持一致[16].
2 基于主元回归模型的生产设计主元回归是一种典型的潜变量建模方法.主元回归方法将n×m维的数据矩阵X线性变换为仅含有A个独立潜变量的得分矩阵 T，每个潜变量是原始输入变量的线性组合：
虽然OLS和主元回归(PCR)方法建立的都是线性回归模型，但是OLS建立的是Y与X列向量之间的映射关系，而PCR方法建立的是Y与T的列向量之间的映射关系：
在求解生产设计问题时，要求新的输出质量变量对应的操作条件xnewT也能具有这种关联结构，为了满足这一要求，新的操作条件xnewT可以表示为
unewT可以通过求解方程式(9)得到：
根据k和A的大小关系，可以分为以下几种情况.
1) 当k>A时，对式(9)求逆可以通过OLS按式(10)求出唯一的最小范数解
：
2) 当k=A时，对式(9)求逆可直接按式(11)求出唯一解：
3) 当k16]，这里只考虑此反应器中的两个区域，其5个可测量质量变量为转化率(Conv)、数均分子量(MWn)、重均分子量(MWw)、长支链(LCB)和短支链(SCB).6个输入变量为进料口温度(Tin)、反应器内压力(P)、两区域的初始进料速率(Fi1，Fi2)和两区域的溶剂浓度(Cs1，Cs2).现有15个不同等级的LDPE数据，前9个等级代表过去生产的产品等级，作为用于建模的历史数据，剩下的6个等级产品的数据用于生产设计结果好坏情况的比较标准.另外，描述生产过程的可测量质量变量通常都是高度相关的.为了选择合适的产品质量变量组成输出矩阵Y，采用主元分析方法对包含所有可测量质量变量的输出矩阵进行预处理，结果表明前3个主成分能够解释全部质量变量的变异信息，解释数据变异信息能力的顺序为MWw，MWn和Conv.
仿真实验分以下两种情况.情况1：Y3=[Conv MWn MWw]，考虑三个质量变量；情况2：Y2=[MWn MWw]，考虑两个质量变量.仿真中建模所用的输入矩阵X包含过去生产9个产品等级中的6个输入变量，两种情况下的主元回归(PCR)模型的输入潜变量空间维数A=3.
3.2 仿真实验情况1 此时A=k=3，不存在零空间.分别采用PCR和OLS求得6个新的期望等级对应过程的操作条件，与实际操作条件比较结果如图 2和图 3所示.
 | 图 2 第一种情况下的Fi1和Cs1 Fig. 2 Fi1 and Cs1 in the first case |
 | 图 3 第一种情况下的P和Tin Fig. 3 P and Tin in the first case |
图 2中除第13个等级外，与OLS模型方法相比，通过PCR模型方法得出的过程操作条件离实际工况更接近.
情况2 输入矩阵X和输出矩阵Y2，此时A－k=1，存在1维零空间.分别采用PCR模型和OLS模型求得6个新的期望等级对应的过程操作条件并与实际工况比较，结果如图 4所示.
 | 图 4 第二情况下的Fi1和Cs1 Fig. 4 Fi1 and Cs1 in the second case |
在图 4中，对应于同一等级的操作条件位于一条直线上，而且对应5个期望等级工况的直线是相互平行的，这也反映了对输出特性没有直接影响的零分量xnullT只能在A-k=1维零空间中移动.图 4中的直线代表一维零空间在Cs1-Fi1平面上的投影.由此可见，通过预测带零空间时的输入变量，可确定零空间操作窗口的变量控制范围，在要求的范围内调整零空间变量，可以帮助系统处于一个稳定状态，为生产一定质量要求的产品提供保障.
4 结论对于具有零空间特性的复杂工业过程的生产设计问题，本文给出的基于隐变量模型的生产设计方法可以确保求得的过程操作条件能够生产出期望等级产品，在技术上具有可行性.这类基于隐变量模型的生产设计方法通过预测带零空间时的输入变量，可确定零空间操作窗口的变量控制范围，在要求的范围内调整零空间变量，为稳定生产和质量控制提供保障.在实际应用中，由生产设计方法得出的过程操作条件可以作为生产过程输入变量调试的初始工作点，有利于减少操作的盲目性，缩短生产调试周期.
| [1] | Jaeckle C M, MacGregor J F.Product design through multivariate statistical analysis of process data[J].AIChE Journal, 1998, 44(5):1105-1118.( 2) |
| [2] | Chen Y D, Shi M H.Application research of knowledge management and engineering for production design organization[J].Electronic & Mechanical Engineering and Information Technology, 2011, 12(8):409-412.(1) |
| [3] | Tomba E, Facco P.Latent variable modeling to assist the implementation of quality-by-design paradigms in pharmaceutical development and manufacturing:a review[J].International Journal of Pharmaceutics, 2013, 74(8):283-297.(2) |
| [4] | Zhang Z H, Tegoshi Y.Representation on tacit knowledge in production design support systems[J].International Conference on Fuzzy Systems and Knowledge Discovery, 2012, 7(10):1546-1549.(1) |
| [5] | Manolakos A, Noam Y.Interference due to null space mismatch in cooperative multipoint MIMO cellular networks[J].Wireless Communications Symposium, 2014, 14(7):5148-5153.(1) |
| [6] | Lakshminarayanan S, Fujii H, Grosman B, et al.New product design via analysis of historical databases[J].Computers & Chemical Engineering, 2000, 24(2):671-676.(1) |
| [7] | Jiang X W, Gao J B.Thin plate spline latent variable models for dimensionality reduction[J]. World Congress on Computational Intelligence, 2012, 47(6):345-352.(1) |
| [8] | Cho M Y, Xu W Y.New algorithms for verifying the null space conditions in compressed sensing[J].Asilomar, 2013, 13(8):1038-1042.(1) |
| [9] | Yun S Y, Kim N J.The group latent variable approach to probity binary classifications[J].Transactions on Neural Networks and Learning Systems, 2014, 25(7):1277-1286.(1) |
| [10] | Zhou Y H, Titchener-Hooker N J.Visualizing integrated bioprocess designs through “windows of operation”[J].Biotechnology and Bioengineering, 1999, 65(5):550-557.(1) |
| [11] | Chen F Z, Wang X Z.Discovery of operational spaces from process data for production of multiple grades of products[J].Industrial & Engineering Chemistry Research, 2000, 39(7):2378-2383.(1) |
| [12] | Chen F Z, Wang X Z.Discovery of operational spaces from process data for production of multiple grades of products[J].Industrial & Engineering Chemistry Research, 2000, 39(7):2378-2383.(1) |
| [13] | Noam Y, Manolakos A.Null space learning with interference feedback for spatial division multiple access[J].Transactions on Wireless Communications, 2014, 13(10):5699-5715.(1) |
| [14] | Bodesheim P, Freytag A.Kernel null space methods for novelty detection[C]// Conference on Computer Vision and Pattern Recognition.Portland, 2013:3374-3381.(1) |
| [15] | Draper N R, Smith H, Pownell E.Applied regression analysis[M].New York:Wiley, 1966.(1) |
| [16] | Kiparissides C, Verros G, Macgregor J F.Mathematical modeling, optimization, and quality control of high-pressure ethylene polymerization reactors[J].Journal of Macromolecular Science, Part C:Polymer Reviews, 1993, 33(4):437-527.(1) |