[发明专利]一种基于PSO‑PID算法的燃气自动控制阀门智能调控方法在审

专利信息
申请号: 201611192138.1 申请日: 2016-12-21
公开(公告)号: CN106773649A 公开(公告)日: 2017-05-31
发明(设计)人: 周鑫浩;张海军 申请(专利权)人: 成都千嘉科技有限公司
主分类号: G05B11/42 分类号: G05B11/42
代理公司: 成都行之专利代理事务所(普通合伙)51220 代理人: 王记明
地址: 610000 四川省成都*** 国省代码: 四川;51
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 基于 pso pid 算法 燃气 自动控制 阀门 智能 调控 方法
【权利要求书】:

1.一种基于PSO-PID算法的燃气自动控制阀门智能调控方法,其特征在于,包括以下步骤:

(1)建立燃气阀门被控对象的传递函数模型,所述传递函数模型为

(2)在传递函数模型的基础上,进行PID算法控制器的设置;

(3)在进行PID算法控制器设置的同时,初始化PSO算法中的粒子群参数;

(4)利用粒子群算法对PID控制器的参数进行优化,比较计算结果是否达到精度或迭代次数,记录最佳Kp、Ki、Kd参数值;

(5)若计算结果达到性能目标,则结束计算,输出Kp、Ki、Kd参数值作为最佳控制参数;若计算结果达不到性能要求,则重复第(3)步、第(4)步,在K时步的结果基础上更新粒子群的速度和位置,输出最佳Kp、Ki、Kd参数值。

2.根据权利要求1所述的一种基于PSO-PID算法的燃气自动控制阀门智能调控方法,其特征在于,所述步骤(1)中燃气阀门被控对象的传递函数模型建模是依据如下步骤完成的:

首先,设定电动机的初始转速为ω,减速后的转速为ω’,将电磁惯性和机械惯性均设置为0,其减速输出:

ω’=k1k2Ur 公式(1)

其中,k1为电机转换系数,k2为减速比,Ur为电动机的工作电压;

其次,根据阀门联轴器的传递作用,使电机减速后的速度与阀门转轴的速度相等,电机减速后的转速为ω’与阀门转角的关系为:

接着,将公式(1)和公式(2)联立,积分可得:

其中,t0为电磁惯性和机械惯性;

最后,对公式(3)进行拉普拉斯变换可得到传递函数G(s)为:

其中,s为复数。

3.根据权利要求1所述的一种基于PSO-PID算法的燃气自动控制阀门智能调控方法,其特征在于,所述步骤(2)中的PID算法控制器的设置,包括如下步骤:

首先,根据给定的目标值r(t)与实际输出值y(t)构成的偏差e(t),将偏差的Kp、Ki和Kd通过线性组合构成控制项,对被控对象进行控制,其控制规律为:

e(t)=r(t)-y(t) 公式(5)

然后,对公式(6)进行离散化处理,即得到燃气自动控制阀增量型PID算法函数Δu(k):

Δu(k)=Kp[e(k)-e(k-1)]+Kie(k)+Kd[e(k)+e(k-2)-2e(k-1)公式(8)

其中,比例系数积分系数微分系数

4.根据权利要求1所述的一种基于PSO-PID算法的燃气自动控制阀门智能调控方法,其特征在于,所述步骤(4)中利用粒子群算法对PID控制器的参数进行优化,包括如下步骤:

首次,将初始化的粒子群参数赋值到PID控制器的Kp、Ki、Kd参数数组中;

其次,根据系统偏差与时间的关系,在此以时间绝对偏差积分ITAE的倒数作为适应度函数,适应度函数变化过程如下:

对于增量型的PID控制器,将适应度函数进行离散化处理,处理后适应度函数如下:

然后,根据适应度函数公式(10)和公式(11)计算每一粒子的适应度值,并找出粒子群中适应度最佳个体的Kp、Ki、Kd参数值。

5.根据权利要求4所述的一种基于PSO-PID算法的燃气自动控制阀门智能调控方法,其特征在于,所述的更新粒子群的速度和位置的步骤如下:

根据公式(10)和公式(11)计算粒子新位置及更新后的适应度值,然后根据公式(12)进行适应度比对:

若当前粒子适应度优于粒子本身前一时刻的最优适应度,则把当前粒子位置作为自身最优位置

若粒子当前适应度优于整个粒子群前一时刻的最优适应度,则把当前位置粒子群作为全局最优

其中,

<mrow><mi>&omega;</mi><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>s</mi><mi>t</mi><mi>a</mi><mi>r</mi><mi>t</mi></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>&omega;</mi><mrow><mi>s</mi><mi>t</mi><mi>a</mi><mi>r</mi><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>e</mi><mi>n</mi><mi>d</mi></mrow></msub></mrow><msub><mi>t</mi><mi>max</mi></msub></mfrac><mi>t</mi><mo>;</mo></mrow>

<mrow><msubsup><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>d</mi></mrow><mi>t</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mi>d</mi></mrow><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>;</mo></mrow>

在公式(12)中,r1和r2是在(0,(1)之间均匀分布的随机数,ωstart和ωend分别为起始权重和终止权重,tmax是最大迭代次数,t是当前迭代次数。

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于成都千嘉科技有限公司,未经成都千嘉科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201611192138.1/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top