專利名稱：一種結晶器armax模型辨識方法
技術領域：
本發明涉及鋼鐵冶金技術中連鑄機結晶器控制系統，尤其涉及ー種結晶器ARMAX(Auto-Regressive and Moving Average Model)模型辨識方法。
背景技術：
結晶器振動對鑄坯脫模及表面質量有著直接、重要的影響，在板坯連鑄實際澆鑄過程中，拉速通常是隨著工況條件(如澆鑄溫度)的變化而發生變化的，為確保獲得良好的鑄坯脫模效果和鑄坯表而質量，應在保證振動エ藝參數基本穩定的前提下，適當地調整頻率、振幅等振動基本參數。然而，要獲得良好的頻率、振幅控制效果，必須設計合理的結晶器控制系統以快速、準確跟蹤頻率、振幅給定值，而優秀的控制系統 是以模型為基礎進行系統分析和設計的，鑒于目前結晶器控制系統基于經驗的PID控制器設計方法，有必要首先對結晶器進行模型辨識，在合理模型基礎上再進行控制系統設計以獲得良好的控制效果。
發明內容
有鑒于此,本發明的主要目的在于提供一種結晶器ARMAX (Auto-Regressive andMoving Average Model)模型辨識方法，其以結晶器油缸閥開度為輸入數據，以結晶器位置為輸出數據，在采樣數據基礎上利用預測誤差法對結晶器ARMAX模型參數進行離線辨識，通過迭代捜索能夠逼近ARMAX模型未知參數全局最優解，能夠獲得帶有干擾模型的結晶器模型，為工程應用中結晶器控制系統設計提供科學的模型基礎。為達到上述目的，本發明的技術方案是這樣實現的一種結晶器ARMAX模型辨識方法，該方法主要包括A、采集輸入、輸出數據，以結晶器油缸閥開度u(t)為輸入數據，以結晶器位置I (t)為輸出數據，采集N對數據樣本Zn ;N為自然數；B、構建結晶器 ARMAX 模型 A(q)y(t) = B (q) u (t)+C (q) e (t);其中A(q)為第一結晶器輸出多項式；B(q)為第二結晶器輸出多項式；C(q)為第三結晶器輸出多項式；e(t)為高斯白噪聲；C、設ARMAX模型待辨識參數為Θ、令基于參數Θ的模型輸出預測值為タひI 0)、令辨識過程的目標函數為V ( Θ，Zn);D、針對所述目標函數，利用預測誤差法進行迭代捜索以逼近Θ全局最優解；E、設預測誤差法捜索停止條件為目標函數值V( Θ , ZN) < tol，其中tol正實數；F、利用隨機賦值的方法為預測誤差法初始化搜索起始點Q1 ;G、計算目標函數ν(θ，ZN)關于Θ的梯度；H :計算第i+Ι代Θ估計值Θ i+1 ；I :若V( Θ i+1,ZN) < tol，則停止運算，輸出參數估計值Θ i+1，否則運行步驟G，繼續進行迭代計算。其中步驟B中所述第一結晶器輸出多項式A(q) = l+a^+ai2+...+anaq-na ;第二結晶器輸出多項式B(q) =M-Wb2q-V^bnbCfb;第三結晶器輸出多項式C(q)=l+c^+c；^2+···+cnecfne ；q_1為后向移動算子；na為結晶器輸出多項式階次，nb為結晶器輸入多項式階次。步驟C所述ARMAX模型待辨識參數為Θ = [&1 a2 …ana Id1 b2 …bnb C1 C2 …cnc]；所述基于參數Θ的模型輸出預測值為
權利要求
1.一種結晶器ARMAX模型辨識方法，其特征在于，該方法主要包括 A、采集輸入、輸出數據，以結晶器油缸閥開度u(t)為輸入數據，以結晶器位置y(t)為輸出數據，采集N對數據樣本Zn ;N為自然數； B、構建結晶器ARMAX模型A(q)y(t)= B(q)u(t)+C(q) e (t);其中A(q)為第一結晶器輸出多項式;B(q)為第二結晶器輸出多項式;C(q)為第三結晶器輸出多項式;e(t)為高斯白噪聲； C、設ARMAX模型待辨識參數為Θ、令基于參數Θ的模型輸出預測值為バ/|の、令辨識過程的目標函數為V ( Θ，Zn)； D、針對所述目標函數，利用預測誤差法進行迭代捜索以逼近Θ全局最優解； E、設預測誤差法捜索停止條件為目標函數值V(Θ , ZN) < tol，其中tol正實數； F、利用隨機賦值的方法為預測誤差法初始化搜索起始點Θi ； G、計算目標函數ν(θ，ZN)關于Θ的梯度； H :計算第i+1代Θ估計值Θ i+1 ； I :若V( Θ i+1,ZN) < tol，則停止運算，輸出參數估計值Θ i+1，否則運行步驟G，繼續進行迭代計算。
2.根據權利要求I所述的結晶器ARMAX模型辨識方法，其特征在于，步驟B中所述第一結晶器輸出多項式A(q) = l+a^+a^+.u+aj^qi1 ;第二結晶器輸出多項式B(q)=b1q_1+b2q_2+---+bnbq_nb ;第三結晶器輸出多項式 C(q) = l+c^+c；^2+...+cMcfn。；q_1 為后向移動算子；na為結晶器輸出多項式階次，nb為結晶器輸入多項式階次。
3.根據權利要求I所述的結晶器ARMAX模型辨識方法，其特征在于，步驟C所述ARMAX模型待辨識參數為Θ = La1 a2 ... ana b1 b2 ... bnb C1 C2 ... cnc]； 所述基于參數Θ的模型輸出預測值為沖Iの=·^^)+y{t)； 所述辨識過程的目標函數為バ仏ZN) = -￡(y(i)-y(t\0))2
4.根據權利要求I或3所述的結晶器ARMAX模型辨識方法，其特征在干，步驟G所述目標函數ν(θ，ΖΝ)關于Θ的梯度其中T=一點が)，f=—士今た)，茅か-叫s{t-k,G) = y(t)-y(t\0) ο
5.根據權利要求4所述的結晶器ARMAX模型辨識方法，其特征在干， QV 步驟H所述第i+1代Θ估計值^+1其中u為迭代步長，u的選擇滿足V( Θ i+1, ZN) < V(OjjZn) ο
全文摘要
本發明公開了一種結晶器ARMAX模型辨識方法，包括A、采集輸入、輸出數據，以結晶器油缸閥開度u(t)為輸入，以結晶器位置y(t)為輸出，采集N對數據樣本ZN；B、構建結晶器ARMAX模型；C、設所述模型待辨識參數為θ、令基于參數θ的模型輸出預測值為令辨識過程的目標函數為V(θ，ZN)；D、針對目標函數，進行迭代搜索以逼近θ全局最優解；E、設預測誤差法搜索停止條件為目標函數值V(θ，ZN)＜tol；F、利用隨機賦值的方法為預測誤差法初始化搜索起始點θ1；G、計算目標函數V(θ，ZN)關于θ的梯度；H計算第i+1代θ估計值θi+1；I若V(θi+1，ZN)＜tol，則停止運算，輸出參數估計值θi+1，否則運行步驟G，繼續進行迭代計算。
文檔編號B22D11/16GK102672128SQ201210130088
公開日2012年9月19日 申請日期2012年4月28日 優先權日2012年4月28日
發明者張華軍 申請人:中冶南方工程技術有限公司