发明名称 一种空调的冷却水回水温度预测控制方法
摘要 一种空调的冷却水回水温度预测控制方法,包括如下步骤:(1)记录当前冷却水回水的温度和冷水水泵的频率,设定冷却水的回水目标温度为37℃;(2)建立预测模型,采用两层的BP神经网络模型;(3)神经网络预测模型的学;(4)反馈校正;(5)参考曲线选取;(6)优化计算:采用RBF神经网络,设<img file="DDA00003340941400011.GIF" wi="578" he="130" />通梯度下降法学网络的各个参数,计算出下一步最优控制量u(k+1),以此作为冷却水泵的下一步频率。本发明提出一种动态性能良好、节能效果明显的空调的冷却水回水温度预测控制方法。
申请公布号 CN103322646B 申请公布日期 2016.01.20
申请号 CN201310236529.9 申请日期 2013.06.13
申请人 浙江工业大学 发明人 董辉;李晓宇;吴祥;邢科新;蒋文成;高阳;罗立锋;仲晓帆
分类号 F24F11/00(2006.01)I 主分类号 F24F11/00(2006.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种中央空调的冷却水回水温度预测控制方法,其特征在于:所述预测控制方法包括如下步骤:(1)记录当前冷却水回水温度和或冷却水泵的频率,设定冷却水回水温度的目标温度为37℃;(2)建立预测模型,采用两层的BP神经网络模型,设温度预测值有如下表达式:y<sub>m</sub>(k)=f[u(k‑1),u(k‑2),...,u(k‑m),y(k‑1),y(k‑2),..,y(k‑n)]y<sub>m</sub>(k+1)=f[u(k),u(k‑1),...,u(k+1‑m),y<sub>m</sub>(k),y(k‑1),..,y(k+1‑n)]y<sub>m</sub>(k+p)=f[u(k+p‑1),u(k+p‑2),...,u(k+p‑m),y<sub>m</sub>(k+p‑1),...,y<sub>m</sub>(k),y(k‑1),..,y(k+p‑n)]其中,y<sub>m</sub>(k)为第k次模型温度输出值,u(k‑1)为第k‑1次冷却水泵的频率值,y(k‑1)为第k‑1次实际温度值,m、n分别为输入和输出的维度,p为预测步长;设在以后的j步中,u(k+j)=u(k+j‑1)=...=u(k+1)=u(k)则有:x(k+j)=[u(k),u(k),...,u(k‑1),...,u(k+j‑m),y<sub>m</sub>(k+j‑1),...,y<sub>m</sub>(k),y(k‑1),..,y(k+j‑n)]<sup>T</sup>;w=[w<sub>1</sub>,w<sub>2</sub>,...,w<sub>m+n</sub>]<sup>T</sup>则:y<sub>m</sub>(k+j)=g[w*x(k+j)];j=1,2,...,p其中,x(k+j)为第k+j次模型的输入变量,p为预测步长;g(x)取单极性sigmoid函数,<img file="FDA0000758271260000011.GIF" wi="439" he="162" />(3)神经网络预测模型的学习,过程如下:采集各个控制对象的阶跃响应d<sub>p</sub>(k),k=1,2,…,n,p表示样本个数m中的每一个,然后根据预测模型的对应阶跃响应输出y<sub>p</sub>(k),k=1,2,…,n;采用梯度下降法来修正参数w;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>&lsqb;</mo><msub><mi>y</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>d</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA0000758271260000012.GIF" wi="668" he="175" /></maths>j(k)对w求导得:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><mo>&dtri;</mo><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>&Delta;</mi><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>m</mi><mo>+</mo><mi>n</mi></mrow></munderover><mfrac><mrow><mo>&part;</mo><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><msub><mi>w</mi><mi>j</mi></msub></mrow></mfrac><msub><mi>&Delta;w</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000758271260000013.GIF" wi="968" he="136" /></maths>若按下式选择Δw<sub>j</sub>(k),<img file="FDA0000758271260000021.GIF" wi="395" he="143" />其中j=1,2,...,m+n;η>0;(4‑1)则Δj(k)≤0,所以只要保证Δw<sub>j</sub>(k)按式(4‑1)来调整,则j(k)就会随着这一调整一直减小;当其误差达到设定范围内后便可停止学习,保存好学习好的w参数;(4)反馈校正:设预测模型输出与实际输出之间的误差为:e(k)=y(k)‑y<sub>m</sub>(k),修正后的反馈为:y<sub>c</sub>(k+p)=y<sub>m</sub>(k+p)+he(k),其中,h为修正因子;(5)参考曲线选取:参考曲线选为:y<sub>r</sub>(k+i)=α<sup>i</sup>y(k)+(1‑α<sup>i</sup>)y<sub>r</sub>,i=1,2,…,p,其中<img file="FDA0000758271260000022.GIF" wi="202" he="102" />T为采样周期,t为时间常数,y(k)为现时刻实际输出,y<sub>r</sub>为设定值;(6)优化计算:采用非线性RBF神经网络来实现滚动优化,设x(k)=[y<sub>r</sub>(k+p);y<sub>c</sub>(k+p),...,y<sub>c</sub>(k+p‑n+1);u(k),...,u(k‑m)]<sup>T</sup><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>h</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><mo>|</mo><mo>|</mo><mi>x</mi><mo>-</mo><msub><mi>c</mi><mi>j</mi></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow><msubsup><mi>&sigma;</mi><mi>j</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000758271260000023.GIF" wi="515" he="209" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>v</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>h</mi><mi>j</mi></msub><mo>&lsqb;</mo><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000758271260000024.GIF" wi="555" he="140" /></maths>其中y<sub>r</sub>为参考轨迹;x为N维输入向量;c<sub>j</sub>为RBF的中心向量;σ<sub>j</sub>为隐节点的宽度;m为隐层节点个数;控制指标为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>E</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>&lsqb;</mo><msub><mi>y</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>j</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA0000758271260000025.GIF" wi="639" he="127" /></maths>通过使控制指标E达到最小,按照梯度下降法调节,来调节连接权值系数、函数的中心向量和隐节点的宽度;根据学习好的非线性RBF神经网络得到下一步预测输入值u(k+1),以此作为冷却水泵的下一步频率。
地址 310014 浙江省杭州市下城区朝晖六区潮王路18号