本發(fā)明涉及高速動(dòng)車(chē)組運(yùn)行過(guò)程多個(gè)工況建模與運(yùn)行優(yōu)化控制方法，屬高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制
技術(shù)領(lǐng)域：
。
背景技術(shù)：
：高速動(dòng)車(chē)組以其載客量大，安全舒適、正點(diǎn)快速、節(jié)能環(huán)保和全天候運(yùn)輸?shù)葍?yōu)勢(shì)，為各國(guó)所重視，成為諸多國(guó)家發(fā)展先進(jìn)交通系統(tǒng)的首選。隨著我國(guó)經(jīng)濟(jì)、社會(huì)的發(fā)展，迫切需要規(guī)?；焖侔l(fā)展高速鐵路，建設(shè)以安全、正點(diǎn)、節(jié)能為標(biāo)志的列車(chē)運(yùn)行控制系統(tǒng)。在軌道交通系統(tǒng)中，列車(chē)運(yùn)行控制系統(tǒng)是確保列車(chē)能夠安全運(yùn)行并提高運(yùn)行效率的核心系統(tǒng)，其控制策略的優(yōu)劣直接影響鐵路運(yùn)輸?shù)哪芰?。高速?dòng)車(chē)組運(yùn)行速度高，路程遠(yuǎn)，其運(yùn)行過(guò)程對(duì)相關(guān)影響因素敏感程度高，受到的影響比傳統(tǒng)鐵路更多，更復(fù)雜。現(xiàn)有的高速動(dòng)車(chē)組運(yùn)行控制是駕駛員在列車(chē)自動(dòng)防護(hù)系統(tǒng)(ATP)指導(dǎo)下，基于實(shí)際的運(yùn)營(yíng)工況的人工操縱控制模式，動(dòng)車(chē)組運(yùn)行性能與駕駛員操作經(jīng)驗(yàn)和對(duì)故障反應(yīng)程度密切相關(guān)。因此，為實(shí)時(shí)保障動(dòng)車(chē)組高速安全運(yùn)行，我們需要建立高速動(dòng)車(chē)組運(yùn)行過(guò)程模型并對(duì)其運(yùn)行過(guò)程進(jìn)行實(shí)時(shí)優(yōu)化控制。傳統(tǒng)的高速動(dòng)車(chē)組運(yùn)行過(guò)程的建模，通常采用基于牽引計(jì)算和運(yùn)行阻力經(jīng)驗(yàn)?zāi)Ｐ偷拿枋龇椒ǎ錈o(wú)法完整刻畫(huà)動(dòng)車(chē)組復(fù)雜多變的動(dòng)態(tài)行為。為更準(zhǔn)確地還原高速動(dòng)車(chē)組的動(dòng)態(tài)運(yùn)行過(guò)程，數(shù)據(jù)驅(qū)動(dòng)建模逐漸成為高速動(dòng)車(chē)組建模的研究熱點(diǎn)。針對(duì)列車(chē)運(yùn)行過(guò)程的跟蹤控制，較經(jīng)典的是PID控制方法。相關(guān)學(xué)者設(shè)計(jì)了一種模糊PID方法控制捷運(yùn)列車(chē)跟蹤速度軌跡，由于PID控制不具有自適應(yīng)能力，其只比較適用于環(huán)境較穩(wěn)定，速度較低的地鐵系統(tǒng)。為了解決這個(gè)問(wèn)題，有學(xué)者采用魯棒自適應(yīng)控制方法實(shí)現(xiàn)高速動(dòng)車(chē)組速度、位置跟蹤控制；但自適應(yīng)魯棒控制的參數(shù)分解和控制律設(shè)計(jì)需要大量計(jì)算，不能很好地解決實(shí)際高速動(dòng)車(chē)組運(yùn)行環(huán)境情況?？紤]到廣義預(yù)測(cè)控制方法可有效地克服過(guò)程的不確定性和非線性，并能方便的處理過(guò)程被控變量的各種約束，適用于復(fù)雜不確定系統(tǒng)，目前較多學(xué)者研究高速動(dòng)車(chē)組廣義預(yù)測(cè)控制，實(shí)現(xiàn)了高速動(dòng)車(chē)組運(yùn)行過(guò)程速度、位移高精度控制。但上述控制方法都不具備實(shí)時(shí)優(yōu)化控制性能，難以及時(shí)消除高速動(dòng)車(chē)組運(yùn)行過(guò)程中一些不確定因素帶來(lái)的影響。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的是，針對(duì)高速動(dòng)車(chē)組運(yùn)行環(huán)境復(fù)雜，工況變化頻繁，影響運(yùn)行性能的因素多等特點(diǎn)，為保證其安全、正點(diǎn)、高效自動(dòng)運(yùn)行，本發(fā)明設(shè)計(jì)一種高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制系統(tǒng)，在線調(diào)整模型參數(shù)，提高高速動(dòng)車(chē)組運(yùn)行過(guò)程的安全性和正點(diǎn)性。本發(fā)明的技術(shù)方案是：一種高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制方法，所述方法通過(guò)采集高速動(dòng)車(chē)組實(shí)際運(yùn)營(yíng)數(shù)據(jù)，建立高速動(dòng)車(chē)組ANFIS模型，并基于此模型設(shè)計(jì)相應(yīng)預(yù)測(cè)控制器，實(shí)時(shí)采集預(yù)測(cè)輸出速度和期望速度，分析運(yùn)行速度誤差；當(dāng)高速動(dòng)車(chē)組受到未知環(huán)境或動(dòng)車(chē)組特性改變等不確定因素對(duì)跟蹤控制影響，結(jié)合卡爾曼濾波算法和BP梯度下降法實(shí)時(shí)調(diào)整模型參數(shù)，以適應(yīng)接下來(lái)的高速動(dòng)車(chē)組高精度跟蹤控制，消除特性改變對(duì)動(dòng)車(chē)組運(yùn)行帶來(lái)的影響，從而保證高速動(dòng)車(chē)組的運(yùn)行性能。所述的高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制方法，所述高速動(dòng)車(chē)組離線ANFIS模型為：假設(shè)有m個(gè)高速動(dòng)車(chē)組運(yùn)行過(guò)程數(shù)據(jù)點(diǎn){X1,…,Xi,…,Xm}，其中Xi＝[vi(k-1),ui(k-1),vi(k)]，數(shù)據(jù)點(diǎn)Xi處的密度指標(biāo)定義為式中，δa為設(shè)定的聚類(lèi)中心有效鄰域半徑，是一個(gè)正數(shù)；選擇密度指標(biāo)最高值得到第一個(gè)聚類(lèi)中心則每個(gè)數(shù)據(jù)點(diǎn)Xi的密度指標(biāo)用以下公式修正其中δb是一個(gè)大于δa的正數(shù)；顯然，靠近第一個(gè)聚類(lèi)中心c1的數(shù)據(jù)點(diǎn)的密度指標(biāo)將顯著減小，這樣使得這些點(diǎn)不太可能選為下一個(gè)聚類(lèi)中心；修正了每個(gè)數(shù)據(jù)點(diǎn)的密度指標(biāo)后，選定下一個(gè)聚類(lèi)中心c2，再次修正數(shù)據(jù)點(diǎn)的所有密度指標(biāo)，重復(fù)此過(guò)程，直至得到最后一個(gè)聚類(lèi)中心cn，因此高速動(dòng)車(chē)組運(yùn)行過(guò)程數(shù)據(jù)的聚類(lèi)中心個(gè)數(shù)為n；對(duì)n條規(guī)則后件采用最小方差估計(jì)獲得n個(gè)線性模型；對(duì)這n個(gè)線性模型進(jìn)行融合，獲得以下離線ANFIS模型：其中，是所有規(guī)則適應(yīng)度的歸一化值。所述的高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制方法，所述基于離線ANFIS模型的預(yù)測(cè)控制器為：所述的高速動(dòng)車(chē)組運(yùn)行過(guò)程模型(7)可描述為以下形式式中，和是z-1的多項(xiàng)式，Δ＝1-z-1；其中參數(shù)和由建模過(guò)程獲得。為得到最優(yōu)的控制標(biāo)量，可設(shè)計(jì)性能指標(biāo)函數(shù)為其中，L,H,G和是引入的丟番圖方程參數(shù)矩陣，是加權(quán)系數(shù)矩陣；最小化性能指標(biāo)(即)得到最優(yōu)控制增量為所述的高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制方法，當(dāng)高速動(dòng)車(chē)組受到未知環(huán)境或動(dòng)車(chē)組特性改變等不確定因素對(duì)跟蹤控制影響，啟動(dòng)在線調(diào)整模型策略，以適應(yīng)接下來(lái)的高速動(dòng)車(chē)組高精度跟蹤控制，實(shí)時(shí)優(yōu)化控制策略通過(guò)結(jié)合卡爾曼濾波算法和BP梯度下降法實(shí)時(shí)調(diào)整模型參數(shù)，具體優(yōu)化步驟可表現(xiàn)為：Step1.設(shè)計(jì)基于高速動(dòng)車(chē)組ANFIS模型的廣義預(yù)測(cè)控制器，計(jì)算t時(shí)刻的速度反饋誤差er(t)＝y(tǒng)(t)-yr(t)；Step2.執(zhí)行高速動(dòng)車(chē)組跟蹤運(yùn)行控制；Step3.判斷時(shí)刻t是否到達(dá)總時(shí)刻Tt，如果到達(dá)則結(jié)束進(jìn)程，否則進(jìn)入下一時(shí)刻t+1，轉(zhuǎn)入Step4；Step4.t+1時(shí)刻，將上一時(shí)刻的反饋誤差er(t)與設(shè)定誤差閾值χ比較，若|er(t)|≥χ，則表明模型失配,不適用于當(dāng)時(shí)控制情況，需要Step5對(duì)模型參數(shù)進(jìn)行在線調(diào)整；如果反饋誤差|er(t)|＜χ，則表明跟蹤精度較高，暫不需優(yōu)化模型參數(shù)，直接返回Step1；Step5.結(jié)合卡爾曼濾波算法和BP梯度下降法，在線調(diào)整ANFIS模型前件，后件參數(shù)，優(yōu)化后的模型代入廣義預(yù)測(cè)控制器中，返回Step1。所述的高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制方法，所述Step5.結(jié)合卡爾曼濾波算法和BP梯度下降法，在線調(diào)整ANFIS模型前件，后件參數(shù)的方法為：在已建立的模型(8)的基礎(chǔ)上，首先采用以下卡爾曼濾波算法調(diào)整后件參數(shù)式中，為建模過(guò)程獲得的后件參數(shù)，遺忘因子0＜λKF≤1通常選擇接近于1的正數(shù)(本文中λKF＝0.9995)；PKF(k)＝qKFI∈R2n×2n，qKF是一個(gè)大正數(shù)，通常取為104～1010(本發(fā)明中qKF＝106)；其次，采用BP梯度下降法實(shí)時(shí)優(yōu)化前件參數(shù)cij和σij；計(jì)算誤差指標(biāo)函數(shù)為式中，第k個(gè)數(shù)據(jù)點(diǎn)y(k)和yr(k)分別表示控制部分傳輸過(guò)來(lái)的模型失配t時(shí)刻的實(shí)際輸出速度和期望輸出速度，以此類(lèi)推；前件參數(shù)優(yōu)化算法如下：對(duì)ANFIS模型參數(shù)進(jìn)行校正之后，代入到廣義預(yù)測(cè)控制器中重新計(jì)算，獲得相應(yīng)的控制力對(duì)高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)施實(shí)時(shí)優(yōu)化控制。本發(fā)明與現(xiàn)有技術(shù)比較的有益效果是，高速動(dòng)車(chē)組是一個(gè)運(yùn)行在多變環(huán)境下的復(fù)雜非線性系統(tǒng)。為改善高速動(dòng)車(chē)組運(yùn)行性能，需要設(shè)計(jì)有效的運(yùn)行過(guò)程控制器對(duì)高速動(dòng)車(chē)組進(jìn)行精確的控制。已有的研究學(xué)者設(shè)計(jì)的高速動(dòng)車(chē)組控制方法都不具有實(shí)時(shí)優(yōu)化功能，難以處理動(dòng)車(chē)組特性或環(huán)境變化導(dǎo)致跟蹤性能變差情況，這使得高速動(dòng)車(chē)組運(yùn)行性能得不到保障。本發(fā)明是基于高速動(dòng)車(chē)組運(yùn)行環(huán)境復(fù)雜，工況變化頻繁，影響運(yùn)行性能的因素多等特點(diǎn)，提出的一種新型的跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制方法。首先建立離線的高速動(dòng)車(chē)組運(yùn)行過(guò)程ANFIS模型，并設(shè)計(jì)相應(yīng)的廣義預(yù)測(cè)控制器。當(dāng)動(dòng)車(chē)組特性或環(huán)境變化導(dǎo)致跟蹤性能變差，啟動(dòng)在線調(diào)整策略，采用卡爾曼濾波和BP梯度下降法對(duì)高速動(dòng)車(chē)組運(yùn)行ANFIS模型進(jìn)行在線調(diào)整，從而調(diào)整預(yù)測(cè)控制器的參數(shù)，實(shí)現(xiàn)動(dòng)車(chē)組跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制，改善了其運(yùn)行的安全性和正點(diǎn)性。本發(fā)明適用于高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制。附圖說(shuō)明圖1為高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制系統(tǒng)結(jié)構(gòu)框圖；圖2為高速動(dòng)車(chē)組運(yùn)行過(guò)程受力情況；圖3為高速動(dòng)車(chē)組運(yùn)行過(guò)程實(shí)時(shí)優(yōu)化控制流程；圖4為檢驗(yàn)數(shù)據(jù)的輸出誤差分布曲線；圖5為特性改變速度跟蹤曲線；圖6為特性改變速度跟蹤誤差曲線；圖7為特性改變牽引/制動(dòng)力曲線；圖8為特性改變加速度曲線；圖9為特性改變ANFIS模型的參數(shù)優(yōu)化過(guò)程；具體實(shí)施方式以下結(jié)合具體實(shí)施例，對(duì)本發(fā)明進(jìn)行詳細(xì)說(shuō)明。本發(fā)明采集高速動(dòng)車(chē)組實(shí)際運(yùn)行數(shù)據(jù)，分析動(dòng)車(chē)組跟蹤運(yùn)行控制機(jī)理、結(jié)合其牽引/制動(dòng)特性曲線和實(shí)際運(yùn)行數(shù)據(jù)，建立高速動(dòng)車(chē)組運(yùn)行過(guò)程離線ANFIS模型，并設(shè)計(jì)基于ANFIS模型的廣義預(yù)測(cè)控制算法實(shí)現(xiàn)對(duì)動(dòng)車(chē)組跟蹤運(yùn)行過(guò)程控制；當(dāng)對(duì)象特性或環(huán)境變化導(dǎo)致跟蹤性能變差，啟動(dòng)模型在線調(diào)整策略，采用卡爾曼濾波和BP梯度下降算法實(shí)時(shí)優(yōu)化ANFIS模型從而對(duì)預(yù)測(cè)控制器參數(shù)的在線調(diào)整，實(shí)現(xiàn)高速動(dòng)車(chē)組運(yùn)行過(guò)程的實(shí)時(shí)優(yōu)化控制。本發(fā)明基于ANFIS的高速列車(chē)運(yùn)行過(guò)程建模步驟為：1、高速動(dòng)車(chē)組實(shí)時(shí)優(yōu)化控制原理分析：圖1中闡述了基于ANFIS模型和廣義預(yù)測(cè)控制的動(dòng)車(chē)組運(yùn)行實(shí)時(shí)優(yōu)化控制系統(tǒng)結(jié)構(gòu)?；跀?shù)據(jù)驅(qū)動(dòng)ANFIS建模方法建立動(dòng)車(chē)組運(yùn)行過(guò)程精確模型，并傳輸給廣義預(yù)測(cè)控制器，經(jīng)過(guò)具體計(jì)算獲得并輸出控制量u，控制動(dòng)車(chē)組跟蹤給定的站間運(yùn)行模式曲線(該運(yùn)行模式曲線由實(shí)際ATP限速曲線和最優(yōu)期望速度曲線構(gòu)成。最優(yōu)期望速度曲線是結(jié)合優(yōu)秀駕駛員的動(dòng)車(chē)組操縱經(jīng)驗(yàn)，基于安全、正點(diǎn)和節(jié)能等運(yùn)行指標(biāo)，從大量高速動(dòng)車(chē)組實(shí)際運(yùn)行速度曲線中篩選確定)運(yùn)行。實(shí)時(shí)采集預(yù)測(cè)輸出速度y和期望速度yr，分析運(yùn)行速度誤差。當(dāng)高速動(dòng)車(chē)組受到未知環(huán)境或動(dòng)車(chē)組特性改變等不確定因素對(duì)跟蹤控制影響，致使反饋速度誤差的絕對(duì)值|er|超出設(shè)定誤差閾值χ時(shí)(為保障高速動(dòng)車(chē)組運(yùn)行的正點(diǎn)性和優(yōu)化控制的實(shí)時(shí)性，閾值χ是結(jié)合CTCS-3列控系統(tǒng)可允許的誤差范圍和實(shí)時(shí)控制采樣周期要求來(lái)選擇的)，返回實(shí)際運(yùn)行數(shù)據(jù)，結(jié)合卡爾曼濾波算法和BP梯度下降法實(shí)時(shí)調(diào)整模型參數(shù)，并基于調(diào)整后的ANFIS模型調(diào)整高速動(dòng)車(chē)組的廣義預(yù)測(cè)控制器，以適應(yīng)接下來(lái)的高速動(dòng)車(chē)組高精度跟蹤控制，消除特性改變對(duì)動(dòng)車(chē)組運(yùn)行帶來(lái)的影響。2、高速動(dòng)車(chē)組運(yùn)行過(guò)程離線ANFIS建模：高速動(dòng)車(chē)組運(yùn)行過(guò)程受力情況如圖2所示，將列車(chē)簡(jiǎn)化為單個(gè)剛性質(zhì)點(diǎn)，并把動(dòng)車(chē)組運(yùn)行過(guò)程中所有的受力作用到這個(gè)質(zhì)點(diǎn)上進(jìn)行分析計(jì)算，圖中y是高速列車(chē)運(yùn)行速度，由測(cè)速測(cè)距單元獲得，u為單位控制力(牽引力/制動(dòng)力)，目前是駕駛員在ATP車(chē)載設(shè)備的指導(dǎo)下操縱手柄獲得，從而達(dá)到牽引、恒速、惰行、制動(dòng)的效果，rb為單位基本阻力，rb＝Ar+Bry+Cry2。圖中高速動(dòng)車(chē)組受力情況可用以下數(shù)學(xué)模型進(jìn)行描述。式中，ε是加速度系數(shù)，Ar、Br、Cr是阻力系數(shù)，Cry2代表空氣阻力，隨著列車(chē)運(yùn)行速度的增加，Cry2所占的比例越大，系統(tǒng)非線性特性越明顯。對(duì)公式(1)進(jìn)行差分變換，可描述為關(guān)系式：y(k)＝f{y(k-1),u(k-1)}(2)鑒于ANFIS綜合了神經(jīng)網(wǎng)絡(luò)的自適應(yīng)學(xué)習(xí)特性和T-S模糊模型的非線性建模特性，其模型結(jié)論部分用線性方程代替了一般Mamdani模糊系統(tǒng)中的模糊數(shù)，使系統(tǒng)可用較少的規(guī)則描述一個(gè)復(fù)雜的非線性系統(tǒng)。對(duì)式(2)描述的動(dòng)車(chē)組運(yùn)行過(guò)程，采用以下模糊推理規(guī)則描述Ri表示第i條模糊推理規(guī)則；y(k-1)、u(k-1)是輸入量，y(k)是輸出量；是輸入量的第i個(gè)模糊集；為后件參數(shù)，n是規(guī)則條數(shù)；ξi是常數(shù)項(xiàng)。假設(shè)有m個(gè)高速動(dòng)車(chē)組運(yùn)行過(guò)程數(shù)據(jù)點(diǎn){X1,…,Xi,…,Xm}，其中Xi＝[vi(k-1),ui(k-1),vi(k)]，數(shù)據(jù)點(diǎn)Xi處的密度指標(biāo)定義為式中，δa為設(shè)定的聚類(lèi)中心有效鄰域半徑，是一個(gè)正數(shù)。選擇密度指標(biāo)最高值得到第一個(gè)聚類(lèi)中心則每個(gè)數(shù)據(jù)點(diǎn)Xi的密度指標(biāo)用以下公式修正其中δb是一個(gè)大于δa的正數(shù)。顯然，靠近第一個(gè)聚類(lèi)中心c1的數(shù)據(jù)點(diǎn)的密度指標(biāo)將顯著減小，這樣使得這些點(diǎn)不太可能選為下一個(gè)聚類(lèi)中心。修正了每個(gè)數(shù)據(jù)點(diǎn)的密度指標(biāo)后，選定下一個(gè)聚類(lèi)中心c2，再次修正數(shù)據(jù)點(diǎn)的所有密度指標(biāo)，重復(fù)此過(guò)程，直至得到最后一個(gè)聚類(lèi)中心cn，因此高速動(dòng)車(chē)組運(yùn)行過(guò)程數(shù)據(jù)的聚類(lèi)中心個(gè)數(shù)為n。對(duì)n條規(guī)則后件采用最小方差估計(jì)獲得n個(gè)線性模型。對(duì)這n個(gè)線性模型進(jìn)行融合，獲得以下ANFIS模型：其中，是所有規(guī)則適應(yīng)度的歸一化值。3、高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制基于上述所建立的高速動(dòng)車(chē)組ANFIS模型，相應(yīng)的廣義預(yù)測(cè)控制器設(shè)計(jì)如下：上述得到的高速動(dòng)車(chē)組運(yùn)行過(guò)程模型(7)可描述為以下形式式中，和是z-1的多項(xiàng)式，Δ＝1-z-1。其中參數(shù)和由建模過(guò)程獲得，可表示為和為被控對(duì)象的模型階次。為得到最優(yōu)的控制標(biāo)量，可設(shè)計(jì)性能指標(biāo)函數(shù)為其中，L,H,G和是引入的丟番圖方程參數(shù)矩陣，是加權(quán)系數(shù)矩陣。最小化性能指標(biāo)(即)得到最優(yōu)控制增量為為了消除未建模部分，未知環(huán)境和故障導(dǎo)致的動(dòng)車(chē)組運(yùn)行特性改變對(duì)跟蹤控制帶來(lái)的影響，反饋控制誤差，實(shí)施ANFIS模型在線調(diào)整策略：Step1.設(shè)計(jì)基于高速動(dòng)車(chē)組ANFIS模型的廣義預(yù)測(cè)控制器(如上述所示)，計(jì)算t時(shí)刻的速度反饋誤差er(t)＝y(tǒng)(t)-yr(t)。Step2.執(zhí)行高速動(dòng)車(chē)組跟蹤運(yùn)行控制。Step3.判斷時(shí)刻t是否到達(dá)總時(shí)刻Tt，如果到達(dá)則結(jié)束進(jìn)程，否則進(jìn)入下一時(shí)刻t+1，轉(zhuǎn)入Step4。Step4.t+1時(shí)刻，將上一時(shí)刻的反饋誤差er(t)與設(shè)定誤差閾值χ比較，若|er(t)|≥χ，則表明模型失配,不適用于當(dāng)時(shí)控制情況，需要Step5對(duì)模型參數(shù)進(jìn)行在線調(diào)整；如果反饋誤差|er(t)|<χ，則表明跟蹤精度較高，暫不需優(yōu)化模型參數(shù)，直接返回Step1。Step5.結(jié)合卡爾曼濾波算法和BP梯度下降法，在線調(diào)整ANFIS模型前件，后件參數(shù)(具體如下所示)。優(yōu)化后的模型代入廣義預(yù)測(cè)控制器中，返回Step1。在已建立的模型(8)的基礎(chǔ)上，首先采用以下卡爾曼濾波算法調(diào)整后件參數(shù)式中，為建模過(guò)程獲得的后件參數(shù)，遺忘因子0<λKF≤1通常選擇接近于1的正數(shù)(本文中λKF＝0.9995)；PKF(k)＝qKFI∈R2n×2n，qKF是一個(gè)大正數(shù)，通常取為104～1010(本發(fā)明中qKF＝106)。其次，采用BP梯度下降法實(shí)時(shí)優(yōu)化前件參數(shù)cij和σij。計(jì)算誤差指標(biāo)函數(shù)為式中，第k個(gè)數(shù)據(jù)點(diǎn)y(k)和yr(k)分別表示控制部分傳輸過(guò)來(lái)的模型失配t時(shí)刻的實(shí)際輸出速度和期望輸出速度，以此類(lèi)推。前件參數(shù)優(yōu)化算法如下：對(duì)ANFIS模型參數(shù)進(jìn)行校正之后，代入到廣義預(yù)測(cè)控制器中重新計(jì)算，獲得相應(yīng)的控制力對(duì)高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)施實(shí)時(shí)優(yōu)化控制，具體流程如圖3所示。綜上所述，針對(duì)高速動(dòng)車(chē)組運(yùn)行環(huán)境復(fù)雜，工況變化頻繁，影響運(yùn)行性能的因素多等特點(diǎn)，建立運(yùn)行過(guò)程ANFIS模型，提出高速動(dòng)車(chē)組跟蹤運(yùn)行實(shí)時(shí)優(yōu)化控制方法。當(dāng)動(dòng)車(chē)組特性或環(huán)境變化導(dǎo)致跟蹤性能變差，啟動(dòng)在線調(diào)整策略，實(shí)時(shí)優(yōu)化控制性能，提高動(dòng)車(chē)組的安全性和正點(diǎn)性。本發(fā)明實(shí)施選用CRH380AL型高速動(dòng)車(chē)組為實(shí)驗(yàn)驗(yàn)證對(duì)象。首先，采集該動(dòng)車(chē)組在京滬高鐵的濟(jì)南西到徐州東區(qū)段10天的全程運(yùn)行速度、控制力數(shù)據(jù)，挑選代表牽引、惰行、制動(dòng)所有工況的全程2000組有效數(shù)據(jù)，并全局平均取其中1400組數(shù)據(jù)作為建模數(shù)據(jù)樣本，剩余600組數(shù)據(jù)作為檢驗(yàn)?zāi)Ｐ途鹊臄?shù)據(jù)。首先，根據(jù)1400組建模樣本數(shù)據(jù)，采用ANFIS建模方法，獲得4條最佳模糊規(guī)則(即式(3)中n＝4)，離線ANFIS模型參數(shù)如表1所示。為驗(yàn)證模型有效性，采用剩余600組運(yùn)行數(shù)據(jù)對(duì)建立的模型進(jìn)行檢驗(yàn)，其模型輸出誤差分布曲線如圖4。：表1ANFIS模型參數(shù)圖4中限速曲線是根據(jù)CTCS-3列控系統(tǒng)的定位測(cè)速要求繪制。觀察圖4的模型驗(yàn)證過(guò)程，當(dāng)速度小于30km/h時(shí)，模型輸出誤差范圍：-0.5876～0.5234km/h，速度大于30km/h時(shí)，誤差范圍：-2.1421～1.6899km/h，滿足CTCS-3列控系統(tǒng)的定位測(cè)速要求，表明所建立的ANFIS模型精度高，泛化能力強(qiáng)，有較好預(yù)測(cè)效果?；谒⒌腁NFIS模型，利用廣義預(yù)測(cè)控制，并結(jié)合卡爾曼濾波算法和BP梯度下降法設(shè)計(jì)實(shí)時(shí)優(yōu)化控制策略對(duì)高速動(dòng)車(chē)組在京滬高鐵線路的濟(jì)南西站——徐州東站區(qū)間的跟蹤運(yùn)行實(shí)施實(shí)時(shí)優(yōu)化控制。動(dòng)車(chē)組在實(shí)際運(yùn)行過(guò)程中，若發(fā)生未知環(huán)境引起動(dòng)車(chē)組特性改變，導(dǎo)致高速動(dòng)車(chē)組跟蹤性能變差，以致難以跟蹤上目標(biāo)曲線，速度跟蹤誤差超出設(shè)定閾值(本文設(shè)定閾值χ＝2km/h)。在這種情況下，本發(fā)明實(shí)時(shí)優(yōu)化控制策略及時(shí)優(yōu)化動(dòng)車(chē)組運(yùn)行過(guò)程模型，使得高速動(dòng)車(chē)組快速再次跟蹤上目標(biāo)曲線。本文在動(dòng)車(chē)組運(yùn)行到里程railmileage＝500km和railmileage＝600km時(shí)，加入不確定的干擾因素，導(dǎo)致動(dòng)車(chē)組運(yùn)行特性發(fā)生改變，控制仿真結(jié)果如圖5～8所示，圖9列出了的實(shí)時(shí)優(yōu)化過(guò)程ANFIS模型參數(shù)的優(yōu)化過(guò)程曲線。圖5～8表明動(dòng)車(chē)組發(fā)生特性改變的情況下，速度突然發(fā)生變化，跟蹤不上目標(biāo)曲線。基于ANFIS的實(shí)時(shí)優(yōu)化控制策略能根據(jù)現(xiàn)有的動(dòng)車(chē)組數(shù)據(jù)實(shí)時(shí)優(yōu)化運(yùn)行模型，調(diào)整控制力大小，使得高速動(dòng)車(chē)組快速更正運(yùn)行速度，再一次高精度跟蹤上目標(biāo)曲線。圖9表明，發(fā)生故障后，本發(fā)明方法對(duì)ANFIS模型的前后件參數(shù)cij,σij(i＝1,2,3,4；j＝1,2)進(jìn)行實(shí)時(shí)精確調(diào)整，直至再次獲得最優(yōu)ANFIS模型，實(shí)時(shí)優(yōu)化了高速動(dòng)車(chē)組運(yùn)行模型。表2列出了本發(fā)明方法對(duì)于未知故障發(fā)生后的速度和加速度跟蹤誤差。表2速度和加速度跟蹤誤差特性改變速度(千米/小時(shí))加速度(米/平方秒)最大負(fù)誤差-1.6502-0.4649最大正誤差2.37630.7095均方根誤差0.77930.0203從表2可直觀看出，在動(dòng)車(chē)組運(yùn)行特性發(fā)生變化的情況下，本發(fā)明方法的最大正負(fù)跟蹤誤差和均方根誤差均控制在一定的范圍內(nèi)，滿足CTCS-3列控系統(tǒng)的定位測(cè)速要求，進(jìn)一步定量的表明了本文方法的有效性性。應(yīng)當(dāng)理解的是，對(duì)本領(lǐng)域普通技術(shù)人員來(lái)說(shuō)，可以根據(jù)上述說(shuō)明加以改進(jìn)或變換，而所有這些改進(jìn)和變換都應(yīng)屬于本發(fā)明所附權(quán)利要求的保護(hù)范圍。當(dāng)前第1頁(yè)1 2 3