首页> 中国专利> 一种飞机推出模糊滑模控制方法

一种飞机推出模糊滑模控制方法

摘要

一种飞机推出模糊滑模控制方法,本发明涉及飞机推出模糊滑模控制方法。本发明是要解决机场地面交通拥挤,机场地面运行效率低,飞机地面油耗高和环境污染的问题而提出的一种飞机推出模糊滑模控制方法。该方法是通过1、明确飞机离港过程的约束条件,并提出飞机离港过程建模的假设条件;2、运用排队论对飞机离港过程进行建模;3、将机场地面信息的模拟信号转化为数字信号;4、计算得到单位时间内飞机平均排队长度Q(t),平均轮候时间T(t),排队长度误差E(t)以及排队长度误差变化量EC(t);5、推出飞机数量U(t);6、得到不会对系统造成波动的飞机推出数量的U(t)或U′(t)统称为N(t);7、经转换器将N(t)的模拟信号反馈给执行机构等步骤实现的。本发明应用于飞机推出模糊滑模控制领域。

著录项

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2017-06-23

    授权

    授权

  • 2015-10-07

    实质审查的生效 IPC(主分类):G05B13/04 申请日:20150515

    实质审查的生效

  • 2015-09-09

    公开

    公开

说明书

技术领域

本发明涉及飞机推出模糊滑模方法,特别涉及一种飞机推出模糊滑模控制方法。

背景技术

随着我国航空运输的快速发展,民航业务需求量持续上升。2013年,全国民航旅客吞 吐量达到3.54亿人,且近几年的年均增长率均在10%以上。然而,我国机场基础设施的 建设却因受到各方面的限制和影响而难以与日益增长的旅客客流相匹配,造成机场“容量” 与客流“需求”之间的失衡。一方面,作为飞机离港过程中的瓶颈,机场跑道本身具有物理 限制(如修建周期长,过程复杂及投资大),短期内难以得到缓解。另一方面,我国机场业 务不均衡问题突出,航空运力集中于少数几个大型枢纽机场的现象加剧了我国机场现有运 行资源的紧张局面。当机场容量小于航班流量时,就会发生交通拥挤和延误现象,给民航 运输业的发展带来挑战。

机场地面交通拥堵所导致的直接后果就是航班延误加剧、飞机轮候时间延长。据统计, 2010年,我国总体航正点率仅为75.8%。解决机场地面交通拥堵的问题有以下两种方法: 一是,加强机场基础设施建设,增加机场容量;二是,通过智能控制飞机离港过程,提高 现有机场基础设施的使用效率。就当前情况来看,后者才是解决机场地面交通拥堵问题的 最直接、快速、有效途径。

发明内容

本发明的目的是为了解决机场地面交通拥挤,机场地面运行效率低,飞机地面油耗高 和环境污染的问题而提出的一种飞机推出模糊滑模控制方法。

上述的发明目的是通过以下技术方案实现的:

步骤一:分析飞机离港过程,明确飞机离港过程的约束条件,并提出飞机离港过程建 模的假设条件;

其中,飞机离港过程的约束条件如下:一条跑道仅容纳一架飞机,跑道上飞机与飞机 之间的间隔需大于尾流间隔;对飞机离港过程建模做如下假设条件为:飞机都取相同的最 小尾流安全距离;

步骤二:运用排队论对飞机离港过程进行建模;即将整个飞机离港过程看作一个串联 排队系统,并计算出飞机离港过程的串联排队系统的平均排队长度和飞机离港过程串 联排队系统的平均轮候时间

步骤三:采用机场现有电子设备以单位时间间隔采集机场地面信息,并将机场地面信 息反馈给塔台控制中心,利用塔台控制中心的转换器将机场地面信息的模拟信号转化为数 字信号;其中,采集机场地面信息包括:飞机滑行时间、飞机在跑道上加速至起飞所用时 间、单位时间内飞机进入滑行道和跑道的到达率λ1和λ2、单位时间内滑行道和跑道的服 务率μ1和μ2

步骤四:根据步骤三中的数字信号利用和计算得到单位时间内飞机平均排队长 度Q(t),平均轮候时间T(t),排队长度误差E(t)以及排队长度误差变化量EC(t);

步骤五:将单位时间内飞机排队长度误差E(t)及飞机排队长度变化量EC(t)作为输入 量,输入模糊控制器即模糊逻辑控制器进行模糊控制器运算得到输出量即推出飞机数量 U(t);

步骤六:将模糊控制器的输出量即推出飞机数量U(t)输入到滑模控制器中;如果切 换增益能消除飞机推出数量滑模控制系统干扰项,将U(t)直接输出;如果切换增益不能 消除飞机推出数量滑模控制系统干扰项,对U(t)进行纠正得到可靠值U′(t),并将纠正后 的U′(t)输出给执行结构,其中,将不会对系统造成波动的飞机推出数量的U(t)或U′(t)统 称为N(t);

步骤七:模糊滑模控制器验证后的可靠值N(t),经转换器将N(t)的数字信号转化成 模拟信号,并将模拟信号反馈给执行机构;执行机构向推出飞机数量信息的电子设备发送 推出指令;即完成了一种飞机推出模糊滑模控制方法。

发明效果

为了缓解机场地面交通拥挤,提高机场地面运行效率,减少飞机地面油耗和环境污染, 实现飞机离港过程的智能控制,本发明公开了一种飞机推出模糊滑模控制方法,包括飞机 离港过程分析,飞机离港过程建模,定时数据采集及处理,模糊控制器运算,滑模控制器 验证,输出执行指令。通过分析飞机离港过程,将该过程视为一个串联排队系统,运用排 队论构建飞机离港队列模型,得到飞机离港队列的平均排队长度和平均轮候时间的计算公 式;采集单位时间内飞机滑行时间及在跑道加速至起飞所用时间等数据,求得单位时间内 飞机离港队列的平均排队长度,平均轮候时间,队列长度误差及其变化量;以排队长度误 差及其变化量为输入、推出飞机数量为输出,设计一个双输入单输出结构的模糊控制器, 根据模糊控制规则和推理运算确定推出飞机数量,将其输入以状态空间、切换函数、控制 函数构成的滑模控制器中,以判断该输出量(推出飞机数量)是否可靠,若可靠,则直接输 出;若不可靠,则经纠正后输出;输出量经转换器由数字信号转化成模拟信号,反馈给执 行机构,下达飞机推出指令。本发明以模糊控制器为主控运算,辅以滑模控制器的验证纠 错功能,实现飞机推出的双保险控制,可有效地缓解机场地面交通拥堵,提高机场运行效 率,还可最大程度地将飞机在滑行道和跑道队列中的轮候时间转化为停机坪的等待时间, 还可减少飞机待飞油耗,降低环境污染,具有广阔的应用前景。

本发明对西安咸阳国际机场T3航站楼单跑道的飞机推出进行控制,并用软件进行仿 真,仿真结果显示采用基于模糊滑模控制方案后,196架次航班的平均轮候时间从仿真前 的6.6min下降到了5.6min,同比下降了15%。而平均排队长度则从仿真前的5架次下降 到4架次,同比下降了20%。

附图说明

图1为具体实施方式一提出的本发明飞机离港过程示意图;

图2(a)为具体实施方式一提出的模糊控制器输入输出的误差E的隶属函数;

图2(b)为具体实施方式一提出的模糊控制器输入输出的误差变化量EC的隶属函数;

图2(c)为具体实施方式一提出的模糊控制器输入输出的可推出飞机数量U的隶属函 数;

图3为具体实施方式一提出的基于模糊滑模理论的飞机推出总体控制方案流程图。

具体实施方式

具体实施方式一:结合图1、图2(a)~(c)和图3说明本实施方式的一种飞机推出模糊 滑模控制方法,具体是按照以下步骤制备的:

一种飞机推出模糊滑模控制方法包括飞机离港过程分析与建模,定时数据采集及处 理,模糊控制器运算,滑模控制器验证,输出执行指令,

步骤一:飞机离港过程分析;分析飞机离港过程,明确飞机离港过程的约束条件,并 提出飞机离港过程建模的假设条件;

其中,飞机离港过程(见附图1)为:空管中心发出起飞指令,由牵引车将飞机从停机 坪推出至滑行道入口,飞机脱离牵引车并打开辅助动力系统进入滑行道队列,滑行至滑行 道末端即进入跑道队列,待前方飞机起飞后便加速起飞,正常离港;本发明对飞机离港过 程设定一些约束条件并做适当简化,以保证建模准确性和合理性;飞机离港过程的约束条 件如下:一条跑道仅容纳一架飞机,跑道上飞机与飞机之间的间隔需大于尾流间隔;对飞 机离港过程建模做如下假设条件为:忽略入港飞机的影响,暂不考虑飞机机型,飞机都取 相同的最小尾流安全距离;

步骤二:飞机离港过程建模;运用排队论对飞机离港过程进行建模;即将整个飞机离 港过程看作一个串联排队系统,并计算出飞机离港过程的串联排队系统的平均排队长度 和飞机离港过程串联排队系统的平均轮候时间

步骤三:定时数据采集及处理;采用机场现有电子设备以单位时间间隔采集机场地面 信息,并将机场地面信息反馈给塔台控制中心,利用塔台控制中心的转换器将机场地面信 息的模拟信号转化为数字信号;其中,采集机场地面信息包括:飞机滑行时间、飞机在跑 道上加速至起飞所用时间、单位时间内飞机进入滑行道和跑道的到达率λ1和λ2、单位时 间内滑行道和跑道的服务率μ1和μ2

步骤四:根据步骤三中的数字信号利用和计算得到单位时间内飞机平均排队长 度Q(t),平均轮候时间T(t),排队长度误差E(t)以及排队长度误差变化量EC(t);

步骤五:模糊控制器运算;将机场地面串联排队系统视为受控对象,塔台人员指令视 为执行机构,进行模糊控制器运算;将单位时间内飞机排队长度误差E(t)及飞机排队长 度变化量EC(t)作为输入量,输入模糊控制器即模糊逻辑控制器进行模糊控制器运算得到 输出量即推出飞机数量U(t);

步骤六:滑模控制器验证;将模糊控制器的输出量即推出飞机数量U(t)输入到滑模 控制器中,判断模糊控制器输出的推出飞机数量U(t)值是否可靠即U(t)是否对系统造成 波动;如果U(t)值不会对系统造成波动即切换增益能消除飞机推出数量滑模控制系统干 扰项,将U(t)直接输出;如果U(t)值会对系统造成波动即切换增益不能消除飞机推出数 量滑模控制系统干扰项,对U(t)进行纠正得到可靠值U′(t)(对U(t)进行纠正得到可靠值 U′(t)是利用滑模控制器来进行纠正的),并将纠正后的U′(t)输出给执行结构,其中,将不 会对系统造成波动的飞机推出数量的U(t)或U′(t)统称为N(t);

步骤七:输出执行指令;模糊滑模控制器验证后的可靠值N(t),经转换器将N(t)的数 字信号转化成模拟信号,并将模拟信号反馈给执行机构;执行机构向推出飞机数量信息的 电子设备发送推出指令;空管中心工作人员将电子设备的推出指令下达给机场牵引车和在 停机坪等候的飞机;即完成了一种飞机推出模糊滑模控制方法。

本实施方式效果:

为了缓解机场地面交通拥挤,提高机场地面运行效率,减少飞机地面油耗和环境污染, 实现飞机离港过程的智能控制,本实施方式公开了一种飞机推出模糊滑模控制方法,包括 飞机离港过程分析,飞机离港过程建模,定时数据采集及处理,模糊控制器运算,滑模控 制器验证,输出执行指令。通过分析飞机离港过程,将该过程视为一个串联排队系统,运 用排队论构建飞机离港队列模型,得到飞机离港队列的平均排队长度和平均轮候时间的计 算公式;采集单位时间内飞机滑行时间及在跑道加速至起飞所用时间等数据,求得单位时 间内飞机离港队列的平均排队长度,平均轮候时间,队列长度误差及其变化量;以排队长 度误差及其变化量为输入、推出飞机数量为输出,设计一个双输入单输出结构的模糊控制 器,根据模糊控制规则和推理运算确定推出飞机数量,将其输入以状态空间、切换函数、 控制函数构成的滑模控制器中,以判断该输出量(推出飞机数量)是否可靠,若可靠,则直 接输出;若不可靠,则经纠正后输出;输出量经转换器由数字信号转化成模拟信号,反馈 给执行机构,下达飞机推出指令。本实施方式以模糊控制器为主控运算,辅以滑模控制器 的验证纠错功能,实现飞机推出的双保险控制,可有效地缓解机场地面交通拥堵,提高机 场运行效率,还可最大程度地将飞机在滑行道和跑道队列中的轮候时间转化为停机坪的等 待时间,还可减少飞机待飞油耗,降低环境污染,具有广阔的应用前景。

本实施方式对西安咸阳国际机场T3航站楼单跑道的飞机推出进行控制,并用软件进 行仿真,仿真结果显示采用基于模糊滑模控制方案后,196架次航班的平均轮候时间从仿 真前的6.6min下降到了5.6min,同比下降了15%。而平均排队长度则从仿真前的5架次 下降到4架次,同比下降了20%。

具体实施方式二:本实施方式与具体实施方式一不同的是:步骤二中运用排队论对飞 机离港过程进行建模具体过程为:

(1)应用排队论对飞机离港过程进行建模;将整个飞机离港过程作为一个串联排队系 统,飞机进入滑行道为第一个排队系统,飞机离开滑行道进入跑道为第二个排队系统;

(2)假设飞机到达滑行道和跑道均服从泊松分布,到达率分别为λ1和λ2;滑行道和 跑道的输出过程均服从负指数分布,滑行道服务率为μ1和跑道服务率为μ2,且μ1和 μ2相互独立;

则飞机离港过程串联排队系统的平均排队长度为:

L=ρ1+ρ2-2ρ1ρ2(1-ρ1)(1-ρ2)

式中,(ρi<1,i=1,2)

λ1为飞机进入滑行道的到达率;

λ2为飞机进入跑道的到达率;

μ1为滑行道服务率,等于滑行道服务时间的倒数,滑行道服务时间即飞机滑行时间;

μ2为跑道服务率,等于跑道服务时间的倒数,跑道服务时间即飞机在跑道上加速 至起飞所用时间;

(3)飞机离港过程串联排队系统的平均轮候时间为:

其它步骤及参数与具体实施方式一相同。

具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤四中根据步骤三 中的数字信号利用和计算得到单位时间内飞机平均排队长度Q(t),平均轮候时间 T(t),排队长度误差E(t)以及排队长度误差变化量EC(t)具体过程为:

(1)单位时间内飞机平均排队长度Q(t),即为飞机离港过程串联排队系统的平均排队长 度L:

Q(t)=L=ρ1+ρ2-2ρ1ρ2(1-ρ1)(1-ρ2)

(2)单位时间内飞机平均轮候时间T(t),即为飞机离港过程串联排队系统的平均轮候时 间

T(t)=W=ρ1ρ2μ1(1-ρ1)+ρ2μ2(1-ρ2)

(3)单位时间内飞机排队长度误差E(t),即为单位时间内飞机平均排队长度Q(t)与期 望排队长度Q之差:

E(t)=Q(t)-Q

(4)单位时间内飞机排队长度误差变化量EC(t),即为本单位时间内飞机排队长度误差 E(t)与前一单位时间内飞机排队长度误差E(t-1)之差;

EC(t)=E(t)-E(t-1)。其它步骤及参数与具体实施方式一或二相同。

具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:步骤五中模糊控 制器的模糊控制的总体规则为:

模糊控制器运算包括确定该控制的模糊控制总体规则,确定控制器的结构,确定输入 输出的语言变量,确定各语言隶属度函数,建立模糊控制规则,确定模糊关系,模糊推理 运算,解模糊;

(1)当排队系统的实际排队长度Q(t)大于期望排队长度Q时,控制系统不予以推出指 令,要求飞机在停机坪持续等待;

(2)当排队系统的实际排队长度Q(t)等于期望排队长度Q时,控制系统允许飞机推 出,并按照逐个推出原则发出推出指令;

(3)当排队系统的实际排队长度Q(t)小于期望排队长度Q时,控制系统允许飞机推出, 推出并按照逐个推出原则发出推出指令;

推出不是指从排队系统推出而是从停机坪推出飞机的推出是指飞机在将其发动机引 擎关闭的状况下,被牵引车从停机坪牵引到滑行道入口位置的过程;整个过程是飞机从停 机坪退出后才进入排队系统,如果排队系统实际排队长度大于期望排队长度了就不要再从 停机坪推出飞机进入排队系统了,如果相等还可以推出,如果小于的话说明就可以推出让 其进入排队系统。其它步骤及参数与具体实施方式一至三之一相同。

具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:步骤五中模糊控 制器运算;将机场地面串联排队系统视为受控对象,塔台人员指令视为执行机构,进行模 糊控制器运算;将单位时间内飞机排队长度误差E(t)及飞机排队长度变化量EC(t)作为输 入量,输入模糊控制器即模糊逻辑控制器进行模糊控制器运算得到输出量即推出飞机数量 U(t)具体过程为:

(1)确定模糊控制器的结构;

模糊逻辑控制器的结构确定为:以单位时间内飞机排队误差E(t)和误差变化量EC(t)为输 入、以推出飞机数量U(t)为输出的双输入单输出二维结构;

(2)确定模糊控制器的输入输出语言变量;

输入参数和输出参数的论域、相应语言值、模糊子集及所代表的意义,如表1、表2、 表3所示;

(3)确定模糊控制器的输入输出各语言变量隶属度函数为三角形函数得到模糊集合 A、B和C的隶属度;

该模糊控制器各语言值模糊化隶属度函数为便于计算的三角形函数,且三个语言值的 隶属度函数图像如附图2所示;为后期计算方便,将模糊集隶属度函数以矢量表示如表4、 表5、表6所示;

(4)建立模糊控制器的19条模糊控制规则;

由飞机推出实际情况确定模糊控制规则,使用的模糊规则是:若输入量E(t)和EC(t)分 别为A与B,则输出量U(t)为C;

1)If E=PB and EC=NB or NS or ZE,then U=NB

2)If E=PB and EC=PB,then U=NM

3)If E=PM and EC=NB,then U=NB

4)If E=PM and EC=PB,then U=NS

5)If E=PM and EC=PS,then U=NM

6)If E=PS and EC=NB,then U=NB

7)If E=PS and EC=NS or ZE,then U=NS

8)If E=PS and EC=PB or PS,then U=ZE

9)If E=ZE and EC=NB or NS,then U=NS

10)If E=ZE and EC=PB or PS,then U=PS

11)If E=ZE and EC=ZE,then U=ZE

12)If E=NS and EC=NB or NS,then U=ZE

13)If E=NS and EC=ZE or PS,then U=PS

14)If E=NS and EC=PB,then U=PM

15)If E=NM and EC=NB,then U=PS

16)If E=NM and EC=ZE or NS,then U=PM

17)If E=NB and EC=NS or NB,then U=PM

18)If E=NB and EC=NS or NB,then U=PB

19)If E=NB and EC=ZE or PS or PB,then U=PB

由19条模糊控制规则可得到模糊控制规则表,如表7所示;

(5)根据建立的模糊控制规则确定模糊关系:

采用第一类推理方式,由每一条模糊控制规则给出一个模糊关系Rl: Rl=RAl∩RBl(l=1,2,...,19),将19条模糊控制规则表示为:

R1=RA1RB1=(A1×C7)(j=35Bj×C7)

R2=RA2∩RB2=(A1×C6)∩((B1∪B2)×C6)

R3=RA3∩RB3=(A2×C7)∩((B4∪B5)×C7)

R4=RA4∩RB4=(A2×C5)∩(B1×C5)

R5=RA5∩RB5=(A2×C6)∩((B2∪B3)×C6)

R6=RA6∩RB6=(A3×C6)∩(B5×C6)

R7=RA7∩RB7=(A3×C5)∩((B3∪B4)×C5)

R8=RA8∩RB8=(A3×C4)∩((B1∪B2)×C4)

R9=RA9∩RB9=(A4×C5)∩((B4∪B5)×C5)

R10=RA10∩RB10=(A4×C3)∩((B1∪B5)×C3)

R11=RA11∩RB11=(A4×C4)∩(B3×C4)

R12=RA12∩RB12=(A5×C4)∩((B4∪B5)×C4)

R13=RA13∩RB13=(A5×C3)∩((B2∪B3)×C3)

R14=RA14∩RB14=(A5×C2)∩(B1×C2)

R15=RA15∩RB15=(A6×C3)∩(B5×C3)

R16=RA16∩RB16=(A6×C2)∩((B3∪B4)×C2)

R17=RA17∩RB17=(A7×C2)∩((B4∪B5)×C2)

R18=RA18∩RB18=(A6×C1)∩((B1∪B2)×C1)

R19=RA19RB19=(A7×C1)(j=13Bj×C1)

其中,R1~R19就是19条模糊控制规则对应的19个模糊关系;每一条模糊控制规则对 应的矩阵运算中得到矩阵RAl和RBl,A1~A7 B1~B5 C1~C7就是模糊控制器定义的输入输出 参数的模糊子集如表1、表2、表3;

根据模糊集合A、B和C的隶属度即表4、表5和表6进行矩阵计算,得到RAl和RBl

根据RAl和RBl计算求得总的模糊关系R=RA∩RB

(6)根据(5)确定的模糊关系R1~R19得到的矩阵RA和RB进行模糊推理运算;

求得模糊控制器的模糊关系R1~R19即得到RA、RB的矩阵结构后,根据合成推理方 法进行模糊推理,得到输出量即推出飞机数量U(t)的模糊矢量值Ul;针对规则库(5)建 立模糊控制规则中的19条控制规则中所设定的任何一条规则l,其输出量的模糊矢量值 Ul为:

代表模糊控制原理里面的一种计算-合成运算合成算子“ο”用来代表两个模糊矩阵 的相乘,与线性代数中的矩阵乘极为相似只是将普通矩阵乘运算中对应元素间相乘用取小 运算“∧”来代替,而运算间相加用取大“∨”来代替;

控制器总的输出控制量的模糊矢量值U为:

(7)解模糊

采用加权平均法进行模糊控制器隶属度解模糊,根据步骤(3)确定的各语言变量隶属 度如表4~表6、模糊控制规则如表7和模糊关系R1~R19利用加权平均法公式得:

xavr=Σi=1nxiμC(xi)Σi=1nμC(xi)

其中:xi——输出量:推出飞机数量U(t)论域上的元素;

μC(xi)——推出飞机数量U(t)的模糊矢量值U;

xavr——推出飞机数量U(t)值。其它步骤及参数与具体实施方式一至四之一相同。

具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:步骤六中对U(t) 进行纠正得到可靠值U′(t)的具体过程:

(1)滑模控制器系统的状态空间为当前串联排队系统中飞机平均排队长度的导数值,表示如下:

Q·(t)=Q(t)-Q(t-1)T=μ2+U(t)+d(t)T

其中:Q(t)——当前串联排队系统中飞机平均排队长度;

Q(t-1)——单位时长前串联排队系统中飞机平均排队长度;

T——单位时长;

μ2——跑道服务率;

U(t)——推出飞机数量;

d(t)——未知干扰,即推出飞机数量的误差;

(2)确定滑模控制器切换函数和机场地面串联排队系统的广义误差E(t)为实际排队长 度与期望排队长度的误差得到滑模控制器的切换函数S(t)为:

S(t)=cE(t)+E·(t)

式中:E(t)——当前各飞机排队长度误差;

——当前各飞机排队长度误差的导数,

系数c大于0,且满足Hurwitz条件;

(3)确定滑模控制器控制函数U(t);Lyapunov函数为:

V=12S2

其中,S为切换函数即S(t);

且,

S·=cE·(t)+E··(t)=cQ(t)-Q(t-1)T+E·(t)-E·(t-1)T=cQ(t)-Q(t-1)T+Q(t)-Q(t-1-[Q(t-1)-Q(t-2)])T2=cQ·(t)+Q(t)-Q(t-2)T2=c[-μ2+U(t)+d(t)T]+Q(t)-Q(t-2)T2

则,

V.=SS.=S{c[-μ2+U(t)+d(t)T]+Q(t)-Q(t-2)T2}

(4)滑模控制器的控制函数U(t)即滑模控制律(推出飞机数量)为:

U(t)=μ2T-Q(t)-Q(t-2)cT+K(t)sgn(S)

K(t)=max|d(t)|+η,η>0,η为误差项;

(5)若S>0,S为切换函数即S(t),则sgn(S)=1;若S=0,则sgn(S)=0;若S<0,则 sgn(S)=-1;

将控制律U(t)带入中,得

V·=S{c[K(t)sgn(S)+d(t)T]}=cT(K(t)|S|+d(t)S)cTη|S|---(1)

在滑模控制律中,切换增益K(t)值是造成抖振的原因,K(t)用于补偿不确定项d(t),以

保证滑模存在条件得到满足;如果d(t)时变,则为了降低抖振,K(t)也时变;可采用 模糊控制实现K(t)的变化;

滑模控制中采用模糊规则,并根据切换增益K(t)的变化,以消除系统干扰项,其中, 模糊规则如下:

确定模糊规则设计;滑模存在条件为当系统达到滑模面后,将会保持在滑模 面上;K(t)必须为保证系统运动得以到达滑模面的增益,K(t)值必须足以消除不确定项的影 响,才能保证滑模存在条件成立;由此可得模糊规则如下:

若则增大K(t);

若则减小K(t);从而确定K(t)的变化量以保证滑模存在条件成立;

(6)采用二维模糊控制器计算△K(t);

a、确定模糊控制器结构,模糊控制器的输入为S和输出为△K(t),即推出飞机数 量的变化量△U(t);

b、确定输入参数S的论域、输入参数的论域、输出参数△K(t)的论域、输入参数S的 相应语言值、输入参数的相应语言值、输出参数△K(t)的相应语言值、输入参数S的模 糊子集、输入参数的模糊子集和输出参数△K(t)的模糊子集如表9、表10、表11所示;

c、确定输入参数S的隶属度、输入参数的隶属度和输出参数△K(t)的隶属度;

滑模控制器中的三个语言变量S、和△K(t)均采用三角形函数,为便于后期计算,将 隶属度函数以矢量表示,如表12、表13、表14所示;

d、根据K(t)的变化量确定滑模控制器的模糊控制规则;根据模糊控制规则按照模糊控 制原理进行计算得到△K(t),即为推出飞机数量变化量△U(t);

满足滑模存在条件成立下设计U(t),控制器的模糊控制规则如表15所示,使用 模糊规则:若输入量S和的取值分别为M和N,则输出量△K(t)为L;

△K(t)为控制量的变化量,即推出飞机数量的变化量△U(t),也可以理解为K(t)需要修正 的值;

(7)对U(t)进行纠正得到可靠值U′(t)具体为:U′(t)=U(t)+△K(t)。其它步骤及参数与 具体实施方式一至五之一相同。

采用以下实施例验证本发明的有益效果:

实施例一:

本实施例一种飞机推出模糊滑模控制方法,具体是按照以下步骤制备的:

一种飞机推出模糊滑模控制方法包括飞机离港过程分析与建模,定时数据采集及处 理,模糊控制器运算,滑模控制器验证,输出执行指令,

步骤一:飞机离港过程分析;分析飞机离港过程,明确飞机离港过程的约束条件,并 提出飞机离港过程建模的假设条件;

其中,飞机离港过程(见附图1)为:空管中心发出起飞指令,由牵引车将飞机从停机 坪推出至滑行道入口,飞机脱离牵引车并打开辅助动力系统进入滑行道队列,滑行至滑行 道末端即进入跑道队列,待前方飞机起飞后便加速起飞,正常离港;本发明对飞机离港过 程设定一些约束条件并做适当简化,以保证建模准确性和合理性;飞机离港过程的约束条 件如下:一条跑道仅容纳一架飞机,跑道上飞机与飞机之间的间隔需大于尾流间隔;对飞 机离港过程建模做如下假设条件为:忽略入港飞机的影响,暂不考虑飞机机型,飞机都取 相同的最小尾流安全距离;

步骤二:飞机离港过程建模;运用排队论对飞机离港过程进行建模;即将整个飞机离 港过程看作一个串联排队系统,并计算出飞机离港过程的串联排队系统的平均排队长度 和飞机离港过程串联排队系统的平均轮候时间

运用排队论对飞机离港过程进行建模具体过程为:

(1)应用排队论对飞机离港过程进行建模;将整个飞机离港过程作为一个串联排队系 统,飞机进入滑行道为第一个排队系统,飞机离开滑行道进入跑道为第二个排队系统;

(2)假设飞机到达滑行道和跑道均服从泊松分布,到达率分别为λ1和λ2;滑行道和 跑道的输出过程均服从负指数分布,滑行道服务率为μ1和跑道服务率为μ2,且μ1和 μ2相互独立;

则飞机离港过程串联排队系统的平均排队长度为:

L=ρ1+ρ2-2ρ1ρ2(1-ρ1)(1-ρ2)

式中,(ρi<1,i=1,2)

λ1为飞机进入滑行道的到达率;

λ2为飞机进入跑道的到达率;

μ1为滑行道服务率,等于滑行道服务时间的倒数,滑行道服务时间即飞机滑行时间;

μ2为跑道服务率,等于跑道服务时间的倒数,跑道服务时间即飞机在跑道上加速 至起飞所用时间;

(3)飞机离港过程串联排队系统的平均轮候时间为:

W=ρ1ρ2μ1(1-ρ1)+ρ2μ2(1-ρ2).

步骤三:定时数据采集及处理;采用机场现有电子设备以单位时间间隔(每30min) 采集机场地面信息,并将机场地面信息反馈给塔台控制中心,利用塔台控制中心的转换器 将机场地面信息的模拟信号转化为数字信号;其中,采集机场地面信息包括:飞机滑行时 间、飞机在跑道上加速至起飞所用时间、单位时间内飞机进入滑行道和跑道的到达率λ1和 λ2、单位时间内滑行道和跑道的服务率μ1和μ2

步骤四:根据步骤三中的数字信号利用和计算得到单位时间(每30min)内飞 机平均排队长度Q(t),平均轮候时间T(t),排队长度误差E(t)以及排队长度误差变化量 EC(t);

(1)单位时间内飞机平均排队长度Q(t),即为飞机离港过程串联排队系统的平均排队长 度

Q(t)=L=ρ1+ρ2-2ρ1ρ2(1-ρ1)(1-ρ2)

(2)单位时间内飞机平均轮候时间T(t),即为飞机离港过程串联排队系统的平均轮候时 间

T(t)=W=ρ1ρ2μ1(1-ρ1)+ρ2μ2(1-ρ2)

(3)单位时间内飞机排队长度误差E(t),即为单位时间内飞机平均排队长度Q(t)与期 望排队长度Q之差:

E(t)=Q(t)-Q

(4)单位时间内飞机排队长度误差变化量EC(t),即为本单位时间内飞机排队长度误差 E(t)与前一单位时间内飞机排队长度误差E(t-1)之差;

EC(t)=E(t)-E(t-1)。

步骤五:模糊控制器运算;将机场地面串联排队系统视为受控对象,塔台人员指令视 为执行机构,进行模糊控制器运算;将单位时间内(每30min)飞机排队长度误差E(t)及 飞机排队长度变化量EC(t)作为输入量,输入模糊控制器即模糊逻辑控制器进行模糊控制 器运算得到输出量即推出飞机数量U(t);

模糊控制器的模糊控制的总体规则为:

模糊控制器运算包括确定该控制的模糊控制总体规则,确定控制器的结构,确定输入 输出的语言变量,确定各语言隶属度函数,建立模糊控制规则,确定模糊关系,模糊推理 运算,解模糊;

(1)当排队系统的实际排队长度Q(t)大于期望排队长度Q时,控制系统不予以推出指 令,要求飞机在停机坪持续等待;

(2)当排队系统的实际排队长度Q(t)等于期望排队长度Q时,控制系统允许飞机推 出,并按照逐个推出原则发出推出指令;

(3)当排队系统的实际排队长度Q(t)小于期望排队长度Q时,控制系统允许飞机推出, 推出并按照逐个推出原则发出推出指令;

推出不是指从排队系统推出而是从停机坪推出飞机的推出是指飞机在将其发动机引 擎关闭的状况下,被牵引车从停机坪牵引到滑行道入口位置的过程;整个过程是飞机从停 机坪退出后才进入排队系统,如果排队系统实际排队长度大于期望排队长度了就不要再从 停机坪推出飞机进入排队系统了,如果相等还可以推出,如果小于的话说明就可以推出让 其进入排队系统。

(1)确定模糊控制器的结构;

模糊逻辑控制器的结构确定为:以单位时间内飞机排队误差E(t)和误差变化量EC(t)为输 入、以推出飞机数量U(t)为输出的双输入单输出二维结构;

(2)确定模糊控制器的输入输出语言变量;

输入参数和输出参数的论域、相应语言值、模糊子集及所代表的意义,如表1、表2、 表3所示;

表1误差E(t)的语言变量设定值

表2误差变化量EC(t)的语言变量设定值

表3输出控制量U(t)的语言变量设定值

(3)确定模糊控制器的输入输出各语言变量隶属度函数为三角形函数得到模糊集合 A、B和C的隶属度;

该模糊控制器各语言值模糊化隶属度函数为便于计算的三角形函数,且三个语言值的 隶属度函数图像如附图2所示;为后期计算方便,将模糊集隶属度函数以矢量表示如表4、 表5、表6所示;

表4模糊集合A的隶属度

uA(x)是指模糊集合A在x=-3\-2\-1\0\1\2\3时的隶属度,这种表达方式是模糊控制原理 里常用的与U(t)完全没有关系;

表5模糊集合B的隶属度

表6模糊集合C的隶属度

(4)建立模糊控制器的19条模糊控制规则;

由飞机推出实际情况确定模糊控制规则,使用的模糊规则是:若输入量E(t)和EC(t)分 别为A与B,则输出量U(t)为C;

1)If E=PB and EC=NB or NS or ZE,then U=NB

2)If E=PB and EC=PB,then U=NM

3)If E=PM and EC=NB,then U=NB

4)If E=PM and EC=PB,then U=NS

5)If E=PM and EC=PS,then U=NM

6)If E=PS and EC=NB,then U=NB

7)If E=PS and EC=NS or ZE,then U=NS

8)If E=PS and EC=PB or PS,then U=ZE

9)If E=ZE and EC=NB or NS,then U=NS

10)If E=ZE and EC=PB or PS,then U=PS

11)If E=ZE and EC=ZE,then U=ZE

12)If E=NS and EC=NB or NS,then U=ZE

13)If E=NS and EC=ZE or PS,then U=PS

14)If E=NS and EC=PB,then U=PM

15)If E=NM and EC=NB,then U=PS

16)If E=NM and EC=ZE or NS,then U=PM

17)If E=NB and EC=NS or NB,then U=PM

18)If E=NB and EC=NS or NB,then U=PB

19)If E=NB and EC=ZE or PS or PB,then U=PB

由19条模糊控制规则可得到模糊控制规则表,如表7所示;

表7模糊控制规则

(5)根据建立的模糊控制规则确定模糊关系:

采用第一类推理方式,由每一条模糊控制规则给出一个模糊关系Rl: Rl=RAl∩RBl(l=1,2,...,19),将19条模糊控制规则表示为:

R1=RA1RB1=(A1×C7)(j=35Bj×C7)

R2=RA2∩RB2=(A1×C6)∩((B1∪B2)×C6)

R3=RA3∩RB3=(A2×C7)∩((B4∪B5)×C7)

R4=RA4∩RB4=(A2×C5)∩(B1×C5)

R5=RA5∩RB5=(A2×C6)∩((B2∪B3)×C6)

R6=RA6∩RB6=(A3×C6)∩(B5×C6)

R7=RA7∩RB7=(A3×C5)∩((B3∪B4)×C5)

R8=RA8∩RB8=(A3×C4)∩((B1∪B2)×C4)

R9=RA9∩RB9=(A4×C5)∩((B4∪B5)×C5)

R10=RA10∩RB10=(A4×C3)∩((B1∪B5)×C3)

R11=RA11∩RB11=(A4×C4)∩(B3×C4)

R12=RA12∩RB12=(A5×C4)∩((B4∪B5)×C4)

R13=RA13∩RB13=(A5×C3)∩((B2∪B3)×C3)

R14=RA14∩RB14=(A5×C2)∩(B1×C2)

R15=RA15∩RB15=(A6×C3)∩(B5×C3)

R16=RA16∩RB16=(A6×C2)∩((B3∪B4)×C2)

R17=RA17∩RB17=(A7×C2)∩((B4∪B5)×C2)

R18=RA18∩RB18=(A6×C1)∩((B1∪B2)×C1)

R1=RA1RB1=(A1×C7)(j=35Bj×C7)

其中,R1~R19就是19条模糊控制规则对应的19个模糊关系;每一条模糊控制规则对 应的矩阵运算中得到矩阵RAl和RBl,A1~A7 B1~B5 C1~C7就是模糊控制器定义的输入输出 参数的模糊子集如表1、表2、表3;

根据模糊集合A、B和C的隶属度即表4、表5和表6进行矩阵计算,得到RAl和RBl

求RB1之前需要先求并集

其中,Bj为第j个模糊子集B

同样地,求得RAl和RBl(l=2,3,…,19);

然后分别根据矩阵RAl和RBl,求得RA和RB

根据RAl和RBl计算求得总的模糊关系R=RA∩RB

(6)根据(5)确定的模糊关系R1~R19得到的矩阵RA和RB进行模糊推理运算;

求得模糊控制器的模糊关系R1~R19即得到RA、RB的矩阵结构后,根据合成推理方 法进行模糊推理,得到输出量即推出飞机数量U(t)的模糊矢量值Ul;针对规则库(5)建 立模糊控制规则中的19条控制规则中所设定的任何一条规则l,其输出量的模糊矢量值 Ul为:

代表模糊控制原理里面的一种计算-合成运算合成算子“ο”用来代表两个模糊矩阵 的相乘,与线性代数中的矩阵乘极为相似只是将普通矩阵乘运算中对应元素间相乘用取小 运算“∧”来代替,而运算间相加用取大“∨”来代替;

控制器总的输出控制量的模糊矢量值U为:

(7)解模糊

采用加权平均法进行模糊控制器隶属度解模糊,根据步骤(3)确定的各语言变量隶属 度如表4~表6、模糊控制规则如表7和模糊关系R1~R19利用加权平均法公式得:

xavr=Σi=1nxiμC(xi)Σi=1nμC(xi)

其中:xi——输出量:推出飞机数量U(t)论域上的元素;

μC(xi)——推出飞机数量U(t)的模糊矢量值U;

xavr——推出飞机数量U(t)值;

(9)把(7)得到的值带进去(8)的公式里计算求得在不同排队误差和误差变化量情况 下的解模糊值E(t)、EC(t)、U(t),如表8所示;

表8采用加权平均法得到的解模糊值

步骤六:滑模控制器验证;将模糊控制器的输出量即推出飞机数量U(t)输入到滑模 控制器中,判断模糊控制器输出的推出飞机数量U(t)值是否可靠即U(t)是否对系统造成 波动;如果U(t)值不会对系统造成波动即切换增益能消除飞机推出数量滑模控制系统干 扰项,将U(t)直接输出;如果U(t)值会对系统造成波动即切换增益不能消除飞机推出数 量滑模控制系统干扰项,对U(t)进行纠正得到可靠值U′(t)(对U(t)进行纠正得到可靠值 U′(t)是利用滑模控制器来进行纠正的),并将纠正后的U′(t)输出给执行结构,其中,将不 会对系统造成波动的飞机推出数量的U(t)或U′(t)统称为N(t);

对U(t)进行纠正得到可靠值U′(t)的具体过程:

(1)滑模控制器系统的状态空间为当前串联排队系统中飞机平均排队长度的导数值,表示如下:

Q.(t)=Q(t)-Q(t-1)T=-μ2+U(t)+d(t)T

其中:Q(t)——当前串联排队系统中飞机平均排队长度;

Q(t-1)——单位时长前串联排队系统中飞机平均排队长度;

T——单位时长;

μ2——跑道服务率;

U(t)——推出飞机数量;

d(t)——未知干扰,即推出飞机数量的误差;

(2)确定滑模控制器切换函数和机场地面串联排队系统的广义误差E(t)为实际排队长 度与期望排队长度的误差得到滑模控制器的切换函数S(t)为:

S(t)=cE(t)+E·(t)

式中:E(t)——当前各飞机排队长度误差;

——当前各飞机排队长度误差的导数,E·(t)=E(t)-E(t-1)T=Q(t)-Q(t-1)T;

系数c大于0,且满足Hurwitz条件;

(3)确定滑模控制器控制函数U(t);Lyapunov函数为:

V=12S2

其中,S为切换函数即S(t);

且,

S·=cE·(t)+E··(t)=cQ(t)-Q(t-1)T+E·(t)-E·(t-1)T=cQ(t)-Q(t-1)T+Q(t)-Q(t-1-[Q(t-1)-Q(t-2)])T2=cQ·(t)+Q(t)-Q(t-2)T2=c[-μ2+U(t)+d(t)T]+Q(t)-Q(t-2)T2

则,

V.=SS.=S{c[-μ2+U(t)+d(t)T]+Q(t)-Q(t-2)T2}

(4)滑模控制器的控制函数U(t)即滑模控制律(推出飞机数量)为:

U(t)=μ2T-Q(t)-Q(t-2)cT+K(t)sgn(S)

K(t)=max|d(t)|+η,η>0,η为误差项;

(5)若S>0,S为切换函数即S(t),则sgn(S)=1;若S=0,则sgn(S)=0;若S<0,则 sgn(S)=-1;

将控制律U(t)带入中,得

V·=S{c[K(t)sgn(S)+d(t)T]}=cT(K(t)|S|+d(t)S)cTη|S|---(1)

在滑模控制律中,切换增益K(t)值是造成抖振的原因,K(t)用于补偿不确定项d(t),以 保证滑模存在条件得到满足;如果d(t)时变,则为了降低抖振,K(t)也时变;可采用 模糊控制实现K(t)的变化;

滑模控制中采用模糊规则,并根据切换增益K(t)的变化,以消除系统干扰项,其中, 模糊规则如下:

确定模糊规则设计;滑模存在条件为当系统达到滑模面后,将会保持在滑模 面上;K(t)必须为保证系统运动得以到达滑模面的增益,K(t)值必须足以消除不确定项的影 响,才能保证滑模存在条件成立;由此可得模糊规则如下:

若则增大K(t);

若则减小K(t);从而确定K(t)的变化量以保证滑模存在条件成立;

(6)采用二维模糊控制器计算△K(t);

a、确定模糊控制器结构,模糊控制器的输入为S和输出为△K(t),即推出飞机数 量的变化量△U(t);

b、确定输入参数S的论域、输入参数的论域、输出参数△K(t)的论域、输入参数S的 相应语言值、输入参数的相应语言值、输出参数△K(t)的相应语言值、输入参数S的模 糊子集、输入参数的模糊子集和输出参数△K(t)的模糊子集如表9、表10、表11所示;

表9滑模控制器的输入量S的模糊集

表10滑模控制器的输入量的模糊集

表11滑模控制器的输出量△K(t)的模糊集

c、确定输入参数S的隶属度、输入参数的隶属度和输出参数△K(t)的隶属度;

滑模控制器中的三个语言变量S、和△K(t)均采用三角形函数,为便于后期计算,将 隶属度函数以矢量表示,如表12、表13、表14所示;

表12滑模控制器的模糊集合M的隶属度

S是滑模变量uM(x)是模糊集合M的隶属度

表13滑模控制器的模糊集合N的隶属度

为滑模变量的导数uM(x)是模糊集合N的隶属度,

表14滑模控制器的模糊集合L的隶属度

d、根据K(t)的变化量确定滑模控制器的模糊控制规则;根据模糊控制规则按照模糊控 制原理进行计算得到△K(t),即为推出飞机数量变化量△U(t);

满足滑模存在条件成立下设计U(t),控制器的模糊控制规则如表15所示,使用 模糊规则:若输入量S和的取值分别为M和N,则输出量△K(t)为L;

△K(t)为控制量的变化量,即推出飞机数量的变化量△U(t),也可以理解为K(t)需要修正 的值;

表15滑模控制器的模糊控制规则

(7)对U(t)进行纠正得到可靠值U′(t)具体为:U′(t)=U(t)+△K(t)。

步骤七:输出执行指令;模糊滑模控制器验证后的可靠值N(t),经转换器将N(t)的 数字信号转化成模拟信号,并将模拟信号反馈给执行机构;执行机构向推出飞机数量信息 的电子设备发送推出指令;空管中心工作人员将电子设备的推出指令下达给机场牵引车和 在停机坪等候的飞机;即完成了一种飞机推出模糊滑模控制方法。

本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术 人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发 明所附的权利要求的保护范围。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号