发明名称 基于预测函数控制的分子量输出PDF控制方法
摘要 本发明公开了属于随机系统输出分布控制领域,特别涉及一种基于预测函数控制的分子量输出PDF控制方法。采用线性B样条函数逼近方法,将分子量的输出概率密度函数表示为预先确定的基函数的线性组合,在确定所有基函数后,通过控制权值向量就可以逼近系统输出概率密度函数,构建聚合材料的分子量输出PDF模型;预测函数控制算法实现分子量输出分布的跟踪控制,采用quadprog函数来求解在限定范围内的最优解。本发明进一步完善分子量分布控制的研究,对比于广义预测控制算法,计算量小,缩短了控制算法的时间,提高了控制的实时性;同时预测函数控制具有规律性,使得系统能够更快的收敛到期望的输出,提高了系统的控制性能和稳定性。
申请公布号 CN105676637A 申请公布日期 2016.06.15
申请号 CN201610015582.X 申请日期 2016.01.11
申请人 华北电力大学 发明人 张金芳;郭萍;赵建勋;李进
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 张文宝
主权项 一种基于预测函数控制的分子量输出PDF控制方法,其特征在于,包括以下步骤:步骤1:首先考虑在聚合反应过程中,聚合物的分子量分布形状的复杂性,分子量分布形状是多峰的或不对称的,不能简单的用集中参数,平均分子量进行描述,而是需要直接考虑完整的分布;采用广义预测控制方法控制分子量输出分布,采用线性B样条函数逼近方法,构建聚合材料的分子量输出PDF模型;为描述模型,采用线性B样条函数逼近方法,分子量输出PDF用一系列基函数的线性组合来表示,这些基函数是预先确定的,因此,控制权值向量就能够控制分子量输出PDF;假定u(k)为k时刻控制输入向量,γ(y,u(k))表示概率密度函数,其中y为函数变量,那么在u(k)的控制下,概率密度函数可以用线性B样条函数的加权来表示,具体如下:<maths num="0001"><math><![CDATA[<mrow><mi>&gamma;</mi><mrow><mo>(</mo><mi>y</mi><mo>,</mo><mi>u</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&omega;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><msub><mi>B</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000904281150000011.GIF" wi="921" he="134" /></maths>其中ωi(u(k))是与u(k)有关的权值,B<sub>i</sub>(y)≥0(i=1,2,...,n)是基函数,n为基函数的个数,上述B样条表示形式能够逼近任意形状的分子量输出PDF;定义<maths num="0002"><math><![CDATA[<mrow><msub><mi>b</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>c</mi><mn>2</mn></msub></munderover><msub><mi>B</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mi>d</mi><mi>y</mi><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>...</mo><mi>n</mi><mo>,</mo></mrow>]]></math><img file="FDA0000904281150000012.GIF" wi="495" he="152" /></maths><maths num="0003"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>B</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><msub><mi>b</mi><mrow><mn>0</mn><mi>n</mi></mrow></msub></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000904281150000013.GIF" wi="306" he="128" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>c</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>B</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mfrac><msub><mi>b</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><msub><mi>b</mi><mrow><mn>0</mn><mi>n</mi></mrow></msub></mfrac><msub><mi>B</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>...</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000904281150000014.GIF" wi="741" he="134" /></maths>C(y)=[c<sub>1</sub>(y),c<sub>2</sub>(y),...,c<sub>n‑1</sub>(y)],V<sub>k</sub>=[ω<sub>1</sub>(u(k)),ω<sub>2</sub>(u(k)),...,ω<sub>n‑1</sub>(u(k))]<sup>T</sup>,则分子量输出PDF可简化表示成如下式子:γ(y,u(k))=C(y)V<sub>k</sub>+L(y)  (2)其中C(y)和L(y)只与B样条函数相关,由此可知,当基函数选定之后,C(y)和L(y)是已知量,为进一步简化分子量输出PDF表达式,令f<sub>m</sub>(y,k)=γ(y,u(k))‑L(y)=C(y)V<sub>k</sub>  (3)步骤2:对于动态系统而言,根据权值和控制输入的动态关系,将分子量输出PDF的模型写成滑动平均模型形式;权值与控制输入的动态关系表示为:V<sub>k+1</sub>=AV<sub>k</sub>+Bu(k),  (4)其中A,B是适当维数的模型参数矩阵,V<sub>k</sub>是由独立权值组成的向量,将式(4)引入位移算子z<sup>‑1</sup>并带入式(3)得到:f<sub>m</sub>(y,k)=C(y)(I‑z<sup>‑1</sup>A)<sup>‑1</sup>Bu(k‑1)  (5)其中I为(n‑1)维单位阵,u(k‑1)是k‑1时刻的输入控制作用,C(y)与B样条函数相关;将(I‑z<sup>‑1</sup>A)<sup>‑1</sup>B根据矩阵理论展开,则式(5)表示为模型输出PDF表达式:<maths num="0005"><math><![CDATA[<mrow><msub><mi>f</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>a</mi><mi>i</mi></msub><msub><mi>f</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>,</mo><mi>k</mi><mo>-</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></munderover><mi>C</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>D</mi><mi>j</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000904281150000021.GIF" wi="1347" he="144" /></maths>其中a<sub>i</sub>和D<sub>j</sub>是(I‑z<sup>‑1</sup>A)<sup>‑1</sup>B的展开式系数;以上所有的参数均由最小二乘法辨识方法得到;步骤3:在步骤2的基础上,采用预测函数控制分子量输出。
地址 102206 北京市昌平区朱辛庄北农路2号
您可能感兴趣的专利