发明名称 基于偏置继电反馈的空调系统模型辨识方法
摘要 一种基于偏置继电反馈的空调系统模型辨识方法,涉及空调控制技术领域,所解决的是缩短测试时间的技术问题。该方法先在空调系统处于闭环模式正常运行状态时,将空调系统控制回路中的PID控制器切换为带滞环的偏置继电器,使得空调系统转换为非线性的闭环系统,再实施偏置继电器反馈实验,再根据偏置继电器反馈实验过程中产生的极限环特性数据,获取被控对象的临界频率特性;再根据被控对象的临界频率特性,得到被控对象的一阶加纯滞后模型。本发明提供的方法,适用于空调系统模型的辨识。
申请公布号 CN104950679A 申请公布日期 2015.09.30
申请号 CN201510337213.8 申请日期 2015.06.17
申请人 上海建坤信息技术有限责任公司;上海市建筑科学研究院(集团)有限公司 发明人 李建;周慎;陈烈
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 上海申汇专利代理有限公司 31001 代理人 林炜
主权项 一种基于偏置继电反馈的空调系统模型辨识方法,涉及配置有PID控制器的空调系统,该方法的具体步骤如下:1)实施偏置继电器反馈实验,实验方法为:先在空调系统处于闭环模式正常运行状态时,将空调系统控制回路中的PID控制器切换为带滞环的偏置继电器,使得空调系统转换为非线性的闭环系统,并将偏置继电器的初始输出值设置为CO‑μ+μ<sub>0</sub>;其中,CO为PID控制器在当前时刻的输出值,μ为偏置继电器的跳变幅值,μ<sub>0</sub>为偏置继电器中的偏差值,且|μ<sub>0</sub>|<μ;然后监测被控对象一段时间,并在监测被控对象过程中根据被控对象的e值来设置偏置继电器的输出值;如果e>val,则在滞后ε时间后将偏置继电器的输出值设置为CO+μ+μ<sub>0</sub>;如果e<‑val,则在滞后ε时间后将偏置继电器的输出值设置为CO‑μ+μ<sub>0</sub>;其中,e为被控对象的输入设定值与pv输出值的差值,val为预先设定的误差阈值,val按被控对象的pv输出值的百分比取值,val的典型值是被控对象的pv输出值的1.2%,ε为预先设定的时间常数;2)根据偏置继电器反馈实验过程中产生的极限环特性数据,获取被控对象的临界频率特性;设PID控制器切换为偏置继电器的时刻为0时刻,偏置继电器的输出值从0时刻起,第i‑1次发生跳变的时刻至第i次发生跳变的时刻为t<sub>i</sub>时间段;根据偏置继电器的极限环特性,有:t<sub>2</sub>=t<sub>4</sub>=t<sub>6</sub>......=T<sub>u1</sub>t<sub>3</sub>=t<sub>5</sub>=t<sub>7</sub>......=T<sub>u2</sub>根据偏置继电器反馈特性,得出被控对象在t时刻的pv输出值为:y(t)=y(0)e<sup>‑(t‑L)/T</sup>+(μ<sub>0</sub>‑μ)K(1‑e<sup>‑(t‑L)/T</sup>)=(y(0)e<sup>L/T</sup>‑μ<sub>0</sub>Ke<sup>L/T</sup>)e<sup>‑t/T</sup>+μ<sub>0</sub>K‑μK(1‑e<sup>‑(t‑L)/T</sup>)=y'(0)e<sup>‑t/T</sup>+μ<sub>0</sub>K‑μK(1‑e<sup>‑(t‑L)/T</sup>)其中,y(t)为被控对象在t时刻的pv输出值,y(0)为被控对象在0时刻的pv输出值,K为被控对象的稳态增益,T为时间常量,L为被控对象的纯滞后时间;得出被控对象的临界频率特性为:A<sub>u</sub>=(μ<sub>0</sub>+μ)K(1‑e<sup>‑L/T</sup>)+εe<sup>‑L/T</sup>A<sub>d</sub>=(μ<sub>0</sub>‑μ)K(1‑e<sup>‑L/T</sup>)‑εe<sup>‑L/T</sup><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>=</mo><mi>T</mi><mi>ln</mi><mfrac><mrow><mn>2</mn><mi>&mu;</mi><msup><mi>Ke</mi><mrow><mi>L</mi><mo>/</mo><mi>T</mi></mrow></msup><mo>+</mo><msub><mi>&mu;</mi><mn>0</mn></msub><mi>K</mi><mo>-</mo><mi>&mu;K</mi><mo>+</mo><mi>&epsiv;</mi></mrow><mrow><mi>&mu;K</mi><mo>+</mo><msub><mi>&mu;</mi><mn>0</mn></msub><mi>K</mi><mo>-</mo><mi>&epsiv;</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000740422960000021.GIF" wi="819" he="166" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>=</mo><mi>T</mi><mi>ln</mi><mfrac><mrow><mn>2</mn><mi>&mu;</mi><msup><mi>Ke</mi><mrow><mi>L</mi><mo>/</mo><mi>T</mi></mrow></msup><mo>-</mo><msub><mi>&mu;</mi><mn>0</mn></msub><mi>K</mi><mo>-</mo><mi>&mu;K</mi><mo>+</mo><mi>&epsiv;</mi></mrow><mrow><mi>&mu;K</mi><mo>-</mo><msub><mi>&mu;</mi><mn>0</mn></msub><mi>K</mi><mo>-</mo><mi>&epsiv;</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000740422960000022.GIF" wi="829" he="170" /></maths>其中,A<sub>u</sub>为被控对象的pv输出峰值,A<sub>d</sub>为被控对象的pv输出谷值;3)根据被控对象的临界频率特性,得到被控对象的一阶加纯滞后模型为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>K</mi><mrow><mi>Ts</mi><mo>+</mo><mn>1</mn></mrow></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>Ls</mi></mrow></msup></mrow>]]></math><img file="FDA0000740422960000023.GIF" wi="407" he="140" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfrac><mrow><msub><mi>A</mi><mi>d</mi></msub><mo>+</mo><msub><mi>A</mi><mi>u</mi></msub></mrow><mrow><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><mi>o</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>L</mi><mo>/</mo><mi>T</mi></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&epsiv;</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>d</mi></msub><mo>+</mo><msub><mi>A</mi><mi>u</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><mi>o</mi></msub><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>0</mn></msub><mo>-</mo><mi>&mu;</mi></mrow><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><mn>0</mn></msub></mfrac><msub><mi>A</mi><mi>u</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>0</mn></msub><mo>+</mo><mi>&mu;</mi></mrow><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><mn>0</mn></msub></mfrac><msub><mi>A</mi><mi>d</mi></msub><mo>-</mo><mi>&epsiv;</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000740422960000024.GIF" wi="1214" he="267" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>u</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mrow><mi>u</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><msup><mi>ln</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mrow><msup><mrow><mo>(</mo><mn>2</mn><mi>&mu;</mi><mo>&CenterDot;</mo><mi>K</mi><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>L</mi><mo>/</mo><mi>T</mi></mrow></msup><mo>-</mo><mi>&mu;</mi><mo>&CenterDot;</mo><mi>K</mi><mo>+</mo><mi>&epsiv;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>0</mn></msub><mo>&CenterDot;</mo><mi>K</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><msup><mrow><mo>(</mo><mi>&mu;</mi><mo>&CenterDot;</mo><mi>K</mi><mo>-</mo><mi>&epsiv;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><mi>&mu;</mi><mo>&CenterDot;</mo><mi>K</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>]</mo></mrow>]]></math><img file="FDA0000740422960000025.GIF" wi="1337" he="273" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mo>-</mo><mi>T</mi><mi>ln</mi><mo>[</mo><mfrac><mn>1</mn><mi>&epsiv;</mi></mfrac><mrow><mo>(</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>0</mn></msub><mo>-</mo><mi>&mu;</mi></mrow><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><mn>0</mn></msub></mfrac><msub><mi>A</mi><mi>u</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>&mu;</mi><mn>0</mn></msub><mo>+</mo><mi>&mu;</mi></mrow><msub><mrow><mn>2</mn><mi>&mu;</mi></mrow><mn>0</mn></msub></mfrac><msub><mi>A</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>.</mo></mrow>]]></math><img file="FDA0000740422960000026.GIF" wi="888" he="195" /></maths>其中,G(s)为被控对象的一阶加纯滞后传递函数。
地址 200032 上海市徐汇区宛平南路75号1号楼22层