• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    隨機(jī)激勵下非線性水輪機(jī)系統(tǒng)的控制設(shè)計(jì)*

    2017-07-03 15:08:08朱晨烜丁云飛
    關(guān)鍵詞:哈密頓概率密度水輪機(jī)

    朱晨烜 丁云飛

    (上海電機(jī)學(xué)院,上海 200240)

    隨機(jī)激勵下非線性水輪機(jī)系統(tǒng)的控制設(shè)計(jì)*

    朱晨烜 丁云飛?

    (上海電機(jī)學(xué)院,上海 200240)

    針對一類隨機(jī)非線性哈密頓系統(tǒng)提出了一種全新的反饋跟蹤控制方法.該控制策略可以準(zhǔn)確地控制系統(tǒng)輸出的概率密度分布特性.閉環(huán)系統(tǒng)的穩(wěn)定性也通過李雅普諾夫函數(shù)法得到嚴(yán)格的數(shù)學(xué)證明.最后,以隨機(jī)非線性水輪機(jī)系統(tǒng)為例,詳細(xì)演示了控制設(shè)計(jì)過程及其有效性.仿真結(jié)果表明,新的反饋控制策略可以使水輪機(jī)系統(tǒng)的輸出滿足預(yù)先指定的平穩(wěn)概率密度函數(shù).

    隨機(jī)非線性系統(tǒng), 反饋?zhàn)粉櫩刂? 概率密度函數(shù), 李雅普諾夫函數(shù)

    引言

    非線性系統(tǒng)控制是最近幾十年的一個(gè)重要研究領(lǐng)域.常規(guī)非線性系統(tǒng)控制設(shè)計(jì)是基于確定性非線性系統(tǒng),例如Lipschitz非線性系統(tǒng)[1,2],單調(diào)非線性[3]和高增益觀察器[4]等.然而,這些方法忽略了系統(tǒng)測量和過程中的隨機(jī)噪聲(如隨機(jī)振動或沖擊).幾乎所有的控制過程都不可避免地受到噪聲或隨機(jī)擾動的影響,相應(yīng)地噪音可以顯著改變動力學(xué)系統(tǒng)的確定性.換句話說,它可以使隨機(jī)擾動下的非線性系具有新的穩(wěn)定狀態(tài),它甚至可以增強(qiáng)一個(gè)非線性系統(tǒng)對外部信號的響應(yīng)[5-7].因此,噪聲干擾下的非線性過程應(yīng)該更確切地被定義為隨機(jī)非線性系統(tǒng)的輸出.

    傳統(tǒng)意義上,隨機(jī)非線性控制的設(shè)計(jì)多基于線性化方法.許多行之有效的控制方法只關(guān)心兩個(gè)被控量,即均值和方差或協(xié)方差[8-11].這類控制基于所述系統(tǒng)輸出符合高斯分布的假定.然而,對于具有本質(zhì)非線性特性的許多實(shí)際過程,如化工,煉鋼或造紙的粒度分布的控制等,這些系統(tǒng)不能被線性化.此時(shí),均值和方差或協(xié)方差不足以表征輸出過程,因?yàn)橐粋€(gè)隨機(jī)非線性系統(tǒng)的輸出通常為非高斯.因此,在最近十幾年,一些基于隨機(jī)非線性系統(tǒng)輸出的概率密度函數(shù)(PDF)的反饋控制方法研究得到了越來越多的科研人員的關(guān)注.

    在過去的幾十年里,一些有效的控制策略被開發(fā)出來以達(dá)到控制輸出PDF的形狀的目的[12-23].其中,伊藤[12]和Fokker-Planck(FP)方程[13]可以用來描述SPDF.在由一個(gè)多項(xiàng)式反饋部分和一個(gè)處理不確定性參數(shù)的切換項(xiàng)所組成的控制設(shè)計(jì)的基礎(chǔ)上,含有不確定參數(shù)的伊藤微分方程被用來表示一類連續(xù)時(shí)間隨機(jī)過程.最近,為了使輸出PDF追隨目標(biāo)PDF,一種遞歸PDF控制策略被用到具有信道時(shí)間延遲和非高斯噪聲的非線性ARMAX系統(tǒng)中,一種多步前饋非線性累計(jì)成本函數(shù)也被用于改善其追蹤性能[14].另一種可行的解決方案是引入平方根B樣條擴(kuò)展神經(jīng)網(wǎng)絡(luò),通過跟蹤動態(tài)權(quán)重逼近輸出PDF[15].同樣,在復(fù)雜的動態(tài)非高斯隨機(jī)過程中,應(yīng)用廣義迭代神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)算法也能估算輸出PDF,隨后輸出PDF的跟蹤控制可以簡化為兩個(gè)相鄰的重復(fù)過程之間的參數(shù)調(diào)整問題[16].一類開關(guān)線性控制器被用來近似一維系統(tǒng)的PDF形狀[17].Herzallah和Karny[18]為了避開明確評估最優(yōu)值函數(shù)而提出了基于自適應(yīng)評判方法的一種完全概率控制算法.Annunziato和Brozi[19]根據(jù)FP方程提出一個(gè)有效的框架來解決多維隨機(jī)過程PDF的最優(yōu)控制.其相應(yīng)的優(yōu)化問題轉(zhuǎn)變?yōu)樵跐L動時(shí)域控制策略中的開環(huán)最優(yōu)系統(tǒng)序列.Abharian和Fadaei[20]討論了PDF控制系統(tǒng)的性能評估問題,評估對象是一個(gè)復(fù)雜核反應(yīng)堆系統(tǒng),并且將評估原則用于分析核研究反應(yīng)堆的能量PDF的性能.該控制器的設(shè)計(jì)使得能量PDF能夠跟蹤給定形狀.一類無界隨機(jī)不確定非線性系統(tǒng)的模型預(yù)測控制方法是通過采用FP方程追蹤PDF, 同時(shí)根據(jù)隨機(jī)控制李雅普洛夫函數(shù)來設(shè)計(jì)穩(wěn)定性約束條件[21].其他隨機(jī)模型預(yù)測控制相關(guān)的最新的文獻(xiàn)綜述可以參見文獻(xiàn)[22].

    應(yīng)當(dāng)注意的是,上述所有提出的設(shè)計(jì)方法大多是基于近似,它們在實(shí)際應(yīng)用中還存在著理論和實(shí)際上的問題[23].這些問題包括三點(diǎn):首先,如果忽略了權(quán)重向量的約束,輸出的PDF的跟蹤可能偏離,這樣會導(dǎo)致分析穩(wěn)定性和性能的困難, 從而不能保證閉環(huán)系統(tǒng)是穩(wěn)定的.其次,大部分方法都只關(guān)注低維隨機(jī)過程模型或一些相對簡單的物理模型.此外,隨著系統(tǒng)維數(shù)的增加,用于閉環(huán)實(shí)現(xiàn)的計(jì)算量過大.

    為了克服這些問題,本文針對一類含隨機(jī)擾動的非線性系統(tǒng)設(shè)計(jì)了一種新的反饋控制方法.該方法利用獲得高維非線性隨機(jī)系統(tǒng)的精確穩(wěn)態(tài)解的技術(shù)來追蹤預(yù)先設(shè)定的SPDF.并引入李雅普洛夫函數(shù)方法保證閉環(huán)系統(tǒng)的輸出收斂于一個(gè)預(yù)先設(shè)定的SPDF,以此來確保閉環(huán)穩(wěn)定性.隨后,通過工程應(yīng)用中的一個(gè)隨機(jī)非線性水輪機(jī)系統(tǒng)的仿真實(shí)驗(yàn)來驗(yàn)證該方法的有效性.本文的結(jié)構(gòu)安排如下:在第1節(jié)中簡要回顧由實(shí)際物理系統(tǒng)轉(zhuǎn)換為相應(yīng)的哈密頓模型的過程;在第2節(jié)中給出一個(gè)跟蹤預(yù)先設(shè)定的SPDF的隨機(jī)反饋控制律的推導(dǎo)過程,同時(shí)還分析了系統(tǒng)的穩(wěn)定性,以保證控制系統(tǒng)的收斂性,詳細(xì)的論證過程將會給出.同時(shí)為了驗(yàn)證所提出的理論,該方法將被用于一個(gè)在彈性水擊和隨機(jī)擾動共同作用下的非線性水輪機(jī)系統(tǒng),所得到的精確最優(yōu)控制律和結(jié)果在第3節(jié)給出.最后一部分給出全文的結(jié)論.

    1 隨機(jī)哈密頓建模

    哈密頓理論被廣泛用于描述實(shí)際系統(tǒng)內(nèi)部的能量耗散及系統(tǒng)與外界的能量交換.直接法作為哈密頓建模理論之一,是由實(shí)際系統(tǒng)的基本動態(tài)特性經(jīng)數(shù)學(xué)演繹而得.由于直接對實(shí)際系統(tǒng)建模十分困難,有時(shí)甚至無法實(shí)現(xiàn),于是基于微分方程與哈密頓模型轉(zhuǎn)換方法的哈密頓實(shí)現(xiàn)理論得到了發(fā)展[25].下面將簡要介紹轉(zhuǎn)換方法和約束條件:

    考慮一個(gè)控制動力系統(tǒng):

    (1)

    其中M是一個(gè)γ維流形,f(x)是向量,g(x)=(g1(x),…,gn(x)),u=(u1,…,un)T,稱系統(tǒng)(1)有一個(gè)狀態(tài)反饋哈密頓實(shí)現(xiàn),如果存在一個(gè)合適的坐標(biāo)卡,一個(gè)哈密頓函數(shù)H和一個(gè)合適的狀態(tài)反饋控制律u=α(x)+u′,使得系統(tǒng)(1)可以表示為:

    (2)

    如果結(jié)構(gòu)矩陣T(x)可以表示為:T(x)=J(x)+Z(x),其中反對稱陣J(x)為:

    (3)

    對稱陣Z(x)為:

    (4)

    其中<.,.>表示內(nèi)積.稱式(2)為反饋哈密頓耗散實(shí)現(xiàn).Z(x)可以進(jìn)一步分解為一個(gè)(半)正定陣R(x)和一個(gè)對稱陣S(x):

    Z(X)=-R(x)+S(x)

    (5)

    其中,

    可選擇如下控制器:

    (6)

    此時(shí),隨機(jī)激勵下的系統(tǒng)(2)可以表示為:

    (7)

    其中ζ(x)代表隨機(jī)激勵幅值;W(t)為高斯白噪聲,均值為零,相關(guān)函數(shù)為E[Wk(t)Wl(t+τ)]=2Dklδ(τ).

    2 追蹤控制設(shè)計(jì)

    2.1 問題的表述

    考慮如下隨機(jī)動力學(xué)系統(tǒng):

    (8)

    (9)

    其中,Bk(t)是標(biāo)準(zhǔn)維納過程;σik=(ζL)ik且有ζ=[ζik],LLT=2D,D=[Dkl].注意愛因斯坦慣例(即:重復(fù)的標(biāo)注表示對該指標(biāo)的求和)被應(yīng)用于全文.式(9)等價(jià)于如下伊藤隨機(jī)微分方程:

    (10)

    其中,

    (11)

    式(11)右側(cè)的第二項(xiàng)為Wong-Zakai修正項(xiàng).與伊藤方程(10)相對應(yīng)的簡化FPK方程(?ρ/?t=0)為:

    (12)

    其中,ρ=ρ(x)是系統(tǒng)(10)的穩(wěn)態(tài)概率密度;bij=(σσT)ij=2Dklζikζjl;文獻(xiàn)[26-27]中指出式(12)的解取決于哈密頓系統(tǒng)的可積性和共振型.對于不可積系統(tǒng),式(12)的解為:

    (13)

    其中N是歸一化常數(shù).λ(H)為概率勢并由以下n個(gè)一階線性常微分方程求得:

    (14)

    這一求解技巧將被用于多自由度隨機(jī)非線性系統(tǒng)對給定平穩(wěn)概率密度函數(shù)的追蹤控制設(shè)計(jì)中.

    2.2 被控過程

    考慮受控的隨機(jī)非線性系統(tǒng).假設(shè)可以由如下受控的耗散的哈密頓系統(tǒng)描述:

    (15)

    (16)

    (17)

    2.3 輸出過程的收斂性

    系統(tǒng)(10)可以簡化為:

    (18)

    定理1

    對任一高斯白噪聲激勵系統(tǒng):

    (19)

    (20)

    此定理的證明詳見文獻(xiàn)[28].

    3 算例

    考慮如下單機(jī)單管彈性水擊下,非線性水輪機(jī)模型:

    (21)

    其中,

    參數(shù)x1和x2是系統(tǒng)內(nèi)部變量,x3=q是水輪機(jī)流量,x4=g是導(dǎo)葉開度,g0,gr分別是導(dǎo)葉開度初始值和額定值,h0是靜態(tài)水頭(每單位),k為引水系統(tǒng)水頭損失參數(shù),Rn是壓力管的特性阻抗(每單位),Te是彈性時(shí)間常量(秒),TW是伺服時(shí)間常量(秒).

    考慮到內(nèi)部能量損失,水輪機(jī)出力為:

    Δpz0[e-m(qz-x3)2-e-m(qz-qnl)2]

    (22)

    其中,pm為輸出機(jī)械功率(每單位),γ=9.81,qnl是每單位水輪機(jī)空載流量,qz為在每單位水頭最高效率下的流量,At是增益系數(shù),Kh是水輪機(jī)管道的水頭損失系數(shù),Δpz0是沖擊損失的系數(shù),n1是每單位的沖擊損失.

    水輪機(jī)系統(tǒng)(21)滿足式(1)的形式,因此該系統(tǒng)的哈密頓模型可由第二節(jié)中給出的哈密頓建模方法得到.可選擇如下哈密頓函數(shù):

    e-n1(qz-qnl)2]x4+x1x2

    (23)

    式子的最后一部分表示內(nèi)部變量在轉(zhuǎn)移過程中產(chǎn)生的能量.變量x3和x4分別是流量和導(dǎo)葉開度,當(dāng)水輪機(jī)開啟并與電力系統(tǒng)連接后,x3>qnl>0,x1,x2,x4>0,因此所選擇的哈密頓函數(shù)為正.

    哈密頓系統(tǒng)的自然輸出為:

    y=gTH=-pm

    (24)

    根據(jù)公式(4)和(5),Z(x)可以分解為:

    =-r(x)+s(x)

    (25)

    其中,

    (26)

    (27)

    代入公式(6)得到控制律u為:

    (28)

    系統(tǒng)(21)可以重建為:

    (29)

    (30)

    (31)

    表1 系統(tǒng)參數(shù)

    為了證明系統(tǒng)的轉(zhuǎn)移概率密度的確隨時(shí)間變化收斂到目標(biāo)概率密度,設(shè)如下Lyapunov函數(shù):

    (32)

    根據(jù)公式(20)可知V的導(dǎo)數(shù)為:

    (33)

    顯然有V(x)≥0,V(x)→∞,當(dāng)|x|→∞且在區(qū)間R4-Ω中L*V<0,其中,

    由此可見,受控系統(tǒng)(29)的轉(zhuǎn)移概率密度會隨著時(shí)間t→∞逐漸逼近到目標(biāo)穩(wěn)態(tài)概率密度.

    然后,通過一個(gè)閉環(huán)系統(tǒng)仿真實(shí)驗(yàn)來評估該方法的有效性.圖1所示是系統(tǒng)能量H在控制過程中的轉(zhuǎn)移概率密度隨時(shí)間的變化過程.當(dāng)t>4s的時(shí)候,H的轉(zhuǎn)移PDF的均值和方差分別是0.28和0.01,這個(gè)結(jié)果和其參考SPDF值N(0.28;0.01)2一致,可見轉(zhuǎn)移概率密度能夠逼近目標(biāo)SPDF.因此,由公式(28)和(31)確定的控制設(shè)計(jì)的確可以使系統(tǒng)(29)追蹤到公式(30)給出的期望SPDF值.圖2展示了輸出功率的響應(yīng)曲線.從中看出,未控系統(tǒng)的輸出為實(shí)線所示,在隨機(jī)激勵下水力輸出響應(yīng)有較大波動;而受控系統(tǒng)的輸出pm(虛線)可以滿足在短時(shí)間內(nèi),大約3~4秒后,經(jīng)過小幅震蕩后穩(wěn)定在一個(gè)特定區(qū)間.

    圖1 轉(zhuǎn)移概率密度ρ(H,t)隨時(shí)間的變化Fig.1 Evolution of ρ(H,t) with the time

    圖2 水力輸出響應(yīng)Fig.2 Hydraulic output response

    4 結(jié)論

    雖然非線性系統(tǒng)的控制設(shè)計(jì)在很多工程領(lǐng)域應(yīng)用中取得了一定進(jìn)展,但是對隨機(jī)擾動下非線性系統(tǒng)的控制研究的不足可能會導(dǎo)致一些問題,如針對更復(fù)雜的隨機(jī)非線性系統(tǒng)而言,可能面臨不準(zhǔn)確的追蹤性能,穩(wěn)定性分析不足或較高的計(jì)算成本等.由于輸出PDF的近似值仍不能準(zhǔn)確地描述系統(tǒng)的特點(diǎn),因此我們研究了一種隨機(jī)非線性系統(tǒng)的反饋?zhàn)粉櫩刂茊栴}.具體來說,該隨機(jī)非線性系統(tǒng)的控制設(shè)計(jì)利用獲得隨機(jī)非線性系統(tǒng)的精確平穩(wěn)解的技術(shù)來追蹤預(yù)先設(shè)定的目標(biāo)SPDF.系統(tǒng)控制方法的閉環(huán)穩(wěn)定性可以通過李雅普洛夫函數(shù)方法證明其收斂性.為了驗(yàn)證理論結(jié)果,所提出的方法被應(yīng)用到彈性水擊下的隨機(jī)非線性水輪機(jī)系統(tǒng)中.系統(tǒng)仿真結(jié)果表明受控系統(tǒng)能夠追蹤任意形狀的參考分布所對應(yīng)的穩(wěn)態(tài)概率分布.將來的工作重點(diǎn)是將這種方法應(yīng)用到更多的工程應(yīng)用領(lǐng)域中去,看看它是否適用于其他種類的隨機(jī)非線性系統(tǒng).

    1Arcak M, Kokotovic P. Nonlinear observers: A circle criterion design and robustness analysis.Automatica, 2001,37(12):1923~1930

    2Chen M S, Chen C C. Robust nonlinear observer for Lipschitz nonlinear systems subject to disturbances.IEEETransactionsonAutomaticControl, 2007,52(12):2365~2369

    3Fan X, Arcak M. Observer design for systems with multivariable monotone nonlinearities.SystemControlLetters, 2003,50(4):319~330

    4Gauthier J P, Hammouri H, Othman S. A simple observer for non-linear systems: Applications to bioreactors.IEEETransactionsonAutomaticControl, 1992,37(6):875~880

    5Jung P. Periodically driven stochastic systems.PhysicsReports, 1993,234(4-5):175~295

    6Lindner B, Garcia-Ojalvo J, Neiman A. Schimansky-Geier L. Effects of noise in excitable systems.PhysicsReports, 2004,392(6):321~424

    7Longtin A, Chialvo D. Stochastic and deterministic resonances in excitable systems.PhysicalReviewLetters, 1998,81(81):4012~4015

    8Astrom K J, Wittenmark B. Self-tuning controllers based on pole-zero placements.IEEProceedings-ControlTheoryandApplications, 1980,127:120 ~130

    9Goodwin G C, SIN K S. Adaptive filtering, prediction and control. New Jersey: Prentice Hall, Englewood Cliffs, 1984

    10 Skelton R E, Iwasaki T, Grigoriadis K M. A unified algebraic approach to linear control design. Bristol, PA:Taylorand Francis, 1998

    11 Lu J, Skelton R R. Covariance control using closed-loop modelling for structures.EarthquakeEngineering&StructuralDynamics, 1998,27(11):1367~1383

    12 Bar-Shalom Y, Li X R. Nonlinear filter design using Fokker-Planck-Kolmogorov probability density evolutions.IEEETransactionsonAerospaceandElectronicSystems, 2000,36(1):309~315

    13 Crespo L G, Sun J Q. Non-linear stochastic control via stationary response design.ProbabilisticEngineeringMechanics, 2003,18(1):73~86

    14 Guo L, Wang H, Wang A P. Optimal probability density function control for NARMAX stochastic systems.Automatica, 2008,44(7):1904~1911

    15 Yang Y, Guo L, Wang H. Constrained PI tracking control for output probability distributions based on two-step neural networks.IEEETransactionsonCircuitsSystemsPartI, 2009,56(7):1416~1426

    16 Yi Y, Sun C Y, Guo L. Probabilistic tracking control for non-Gaussian stochastic process using novel iterative learning algorithms.InternationalJournalofSystemsScience, 2013,44(7):1325~1332

    17 Pigeon B, Perrier M, Srinivasan B. Shaping probability density functions using a switching linear controller.JournalofProcessControl, 2011,21(21):901~908

    18 Herzallah R, Karny. Fully probabilistic control design in an adaptive critic framework.NeuralNetwork, 2011,24(10):1128~1135

    19 Annunziato M, Brozi A. A Fokker-Planck control framework for multidimensional stochastic processes.JournalofComputational&AppliedMathematics, 2013,237(1):487~507

    20 Abharian A E, Fadaei A H. Power probability density function control and performance assessment of a nuclear research reactor.AnnalsofNuclearEnergy, 2014,64:11~20

    21 Buehler E A, Paulson J A, Mesbah A. Lyapunov-based stochastic nonlinear model predictive control: Shaping the state probability density functions. In: Proceedings of the American Control Conference, Boston, To appear, 2016

    22 Mesbah A.Stochastic model predictive control: An overview and perspectives for future research.IEEEControlSystems, 2016,36(6):30~44

    23 Guo L, Wang H. Stochastic distribution control system design: A convex optimization approach. Springer, 2010

    24 Cheng D, Spurgeon S, Xiang J. On the development of generalized Hamiltonian realizations.IEEEConferenceonDecisionandControl, 2000,5(5):5125~5130

    25 Wang Y Z, Li C W, Cheng D Z. New approaches to generalized Hamiltonian realization of autonomous nonlinear systems.ScienceInChinaSeriesF:InformationSciences, 2003,46(5):431~444

    26 Zhu W Q, Cai G Q, Lin Y K. On exact stationary solutions of stochastically perturbed Hamiltonian systems.ProbabilisticEngineeringMechanics, 1990,5(2):84~87

    27 Zhu W Q, Yang Y Q. Exactly stationary solutions of stochastically excited and dissipated integrable Hamiltonian systems.ASMEJournalofAppliedMechanics, 1996,63:493~500

    28 Caughey T K. Nonlinear theory of random vibrations.AdvancesinAppliedMechanics, 1971,11:209~253

    *The project supported by the National Natural Science Foundation of China (11302123), Shanghai Pujiang Program (15PJ1402500), Innovation Program of Shanghai Municipal Education Commission(14YZ163), “Chen Guang” project supported by Shanghai Municipal Education Commission and Shanghai Education Development Function(13CG63)

    ? Corresponding author E-mail:Dingyf@sdju.edu.cn

    2 March 2017,revised 18 April 2017.

    CONTROL DESIGN FOR A NONLINEAR HYDRAULIC TURBINE SYSTEM UNDER STOCHASTIC EXCITATION*

    Zhu Chenxuan Ding Yunfei?

    (ShanghaiDianjiUniversity,Shanghai200240,China)

    A new feedback tracing control method is presented for a class of nonlinear Hamiltonian systems with stochastic disturbances. It aims to accurately control the shape of the output probability density function. Meanwhile, a Lyapunov function method is employed to prove the optimality of the control law and the stability of closed -loop system. Eventually, to evaluate the performance of the proposed approach, a nonlinear stochastic hydraulic turbine system is utilized to be remodeled and to demonstrate the control design procedure and its effectiveness. The system simulation shows that the proposed feedback tracing control strategy for hydro-turbine system can achieve a pre-specified stationary probability density function (SPDF).

    nonlinear stochastic system, feedback tracing control, probability density function, Lyapunov function

    *國家自然科學(xué)基金(11302123)、上海市浦江人才計(jì)劃(15PJ1402500)、上海市教育委員會科研創(chuàng)新項(xiàng)目(14YZ163)、上海市教育委員會和上海市教育發(fā)展基金會“晨光計(jì)劃”(13CG63)資助

    10.6052/1672-6553-2017-031

    2017-3-2收到第1稿,2017-4-18收到修改稿.

    ? 通訊作者 E-mail:Dingyf@sdju.edu.cn

    猜你喜歡
    哈密頓概率密度水輪機(jī)
    水輪機(jī)過流面非金屬材料的修復(fù)及防護(hù)
    大中型水斗式水輪機(jī)的關(guān)鍵技術(shù)
    水輪機(jī)虛擬仿真動畫制作的研究
    連續(xù)型隨機(jī)變量函數(shù)的概率密度公式
    AKNS系統(tǒng)的對稱約束及其哈密頓結(jié)構(gòu)
    一類四階離散哈密頓系統(tǒng)周期解的存在性
    一類新的離散雙哈密頓系統(tǒng)及其二元非線性可積分解
    水輪機(jī)過流部件改造與節(jié)能增效
    Hunt過程在Girsanov變換下的轉(zhuǎn)移概率密度的表示公式
    隨機(jī)變量線性組合的分布的一個(gè)算法
    九九在线视频观看精品| 欧美激情在线99| 在线免费十八禁| 午夜精品国产一区二区电影 | 高清午夜精品一区二区三区 | 91在线精品国自产拍蜜月| 久99久视频精品免费| 国产色爽女视频免费观看| 成人永久免费在线观看视频| 国产午夜精品久久久久久一区二区三区 | 观看美女的网站| 在线免费十八禁| 精品一区二区三区av网在线观看| 在线天堂最新版资源| 最新中文字幕久久久久| 久久精品综合一区二区三区| 别揉我奶头~嗯~啊~动态视频| 观看免费一级毛片| av.在线天堂| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品夜夜夜夜夜久久蜜豆| 18禁在线播放成人免费| 久久99热6这里只有精品| 国产激情偷乱视频一区二区| 亚洲精品亚洲一区二区| 麻豆国产97在线/欧美| 亚洲最大成人手机在线| 国产色婷婷99| 欧美一区二区精品小视频在线| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久中文| 国产女主播在线喷水免费视频网站 | 性欧美人与动物交配| 小蜜桃在线观看免费完整版高清| 神马国产精品三级电影在线观看| 俄罗斯特黄特色一大片| 人妻久久中文字幕网| 特级一级黄色大片| 国产一区二区在线观看日韩| 嫩草影院入口| 乱码一卡2卡4卡精品| 亚洲专区国产一区二区| 丰满人妻一区二区三区视频av| 亚洲五月天丁香| 晚上一个人看的免费电影| 伦精品一区二区三区| av在线观看视频网站免费| 久久久精品欧美日韩精品| 热99re8久久精品国产| 亚洲丝袜综合中文字幕| av在线蜜桃| 亚洲自偷自拍三级| 非洲黑人性xxxx精品又粗又长| 国产伦精品一区二区三区四那| 一区二区三区高清视频在线| 露出奶头的视频| 亚洲av免费高清在线观看| 丝袜喷水一区| 一区二区三区高清视频在线| 偷拍熟女少妇极品色| 久久午夜亚洲精品久久| 天堂网av新在线| 中文亚洲av片在线观看爽| 国产伦精品一区二区三区视频9| 亚洲av一区综合| 国产美女午夜福利| 国产伦一二天堂av在线观看| 2021天堂中文幕一二区在线观| 舔av片在线| aaaaa片日本免费| 国产亚洲精品综合一区在线观看| 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 国产精品伦人一区二区| 亚洲专区国产一区二区| 国产黄a三级三级三级人| 国产片特级美女逼逼视频| 啦啦啦观看免费观看视频高清| 国内精品美女久久久久久| 国产精品一区二区三区四区久久| 三级毛片av免费| 少妇的逼水好多| 神马国产精品三级电影在线观看| 日日摸夜夜添夜夜爱| 欧美另类亚洲清纯唯美| 国产精品久久久久久久电影| 国产一区二区亚洲精品在线观看| 欧美3d第一页| 免费av不卡在线播放| 九九久久精品国产亚洲av麻豆| 国产精品久久视频播放| 99久久精品一区二区三区| 男女啪啪激烈高潮av片| 搡老妇女老女人老熟妇| 又爽又黄a免费视频| 日本黄色视频三级网站网址| 成人特级av手机在线观看| 99热6这里只有精品| 91在线观看av| 美女大奶头视频| av视频在线观看入口| 欧美国产日韩亚洲一区| 悠悠久久av| 成年女人永久免费观看视频| 国产精品久久久久久久久免| 99视频精品全部免费 在线| 国产精品人妻久久久影院| 99在线视频只有这里精品首页| 国产精品电影一区二区三区| 国产三级中文精品| 亚洲av一区综合| 99热这里只有精品一区| 91久久精品电影网| 一个人免费在线观看电影| 亚洲一区高清亚洲精品| 一级黄色大片毛片| 女生性感内裤真人,穿戴方法视频| 狂野欧美白嫩少妇大欣赏| 久久亚洲国产成人精品v| 婷婷色综合大香蕉| 精品一区二区三区视频在线| 熟妇人妻久久中文字幕3abv| 国产一区二区在线观看日韩| 麻豆国产av国片精品| 中文字幕av在线有码专区| 国产男靠女视频免费网站| 少妇裸体淫交视频免费看高清| 99久国产av精品| 黄色日韩在线| 成年女人永久免费观看视频| 如何舔出高潮| 我的老师免费观看完整版| 色综合亚洲欧美另类图片| 最后的刺客免费高清国语| 最近最新中文字幕大全电影3| 欧美一区二区精品小视频在线| 欧美性感艳星| 国产大屁股一区二区在线视频| 国产成人精品久久久久久| 国产精品日韩av在线免费观看| 99热这里只有是精品在线观看| 婷婷亚洲欧美| 国产精品国产三级国产av玫瑰| 99在线视频只有这里精品首页| 日日摸夜夜添夜夜添av毛片| 国产高清激情床上av| 久久久色成人| 搡老熟女国产l中国老女人| 五月伊人婷婷丁香| 国产人妻一区二区三区在| 久久亚洲国产成人精品v| 一a级毛片在线观看| 高清毛片免费观看视频网站| 国产成人一区二区在线| 天堂动漫精品| 精品99又大又爽又粗少妇毛片| 国产极品精品免费视频能看的| 国产精品日韩av在线免费观看| 老司机午夜福利在线观看视频| .国产精品久久| 国产精品一区二区三区四区久久| 午夜福利成人在线免费观看| ponron亚洲| 少妇裸体淫交视频免费看高清| 级片在线观看| 欧美一区二区精品小视频在线| 男女之事视频高清在线观看| 欧美潮喷喷水| 一进一出好大好爽视频| 97超级碰碰碰精品色视频在线观看| 免费看光身美女| 久久久久久久久久黄片| 亚洲精品456在线播放app| 波野结衣二区三区在线| 激情 狠狠 欧美| 国产91av在线免费观看| 日韩人妻高清精品专区| 午夜福利在线观看吧| 国产色婷婷99| 久久婷婷人人爽人人干人人爱| 欧美色视频一区免费| 色5月婷婷丁香| 香蕉av资源在线| 国产精品美女特级片免费视频播放器| 99久久九九国产精品国产免费| 国产亚洲91精品色在线| 99热这里只有精品一区| 最好的美女福利视频网| 一本一本综合久久| 久久久久久久久久成人| 有码 亚洲区| 欧美极品一区二区三区四区| 嫩草影视91久久| 欧美zozozo另类| 99热这里只有是精品50| 亚洲美女视频黄频| 国产视频内射| 日韩欧美 国产精品| 国产黄色小视频在线观看| 国产大屁股一区二区在线视频| 蜜臀久久99精品久久宅男| av卡一久久| 欧美潮喷喷水| 丰满人妻一区二区三区视频av| 麻豆精品久久久久久蜜桃| av天堂中文字幕网| 亚洲av一区综合| 亚洲av成人精品一区久久| 国产三级在线视频| 好男人在线观看高清免费视频| 日本熟妇午夜| 欧美性猛交黑人性爽| 亚洲美女搞黄在线观看 | 亚洲精品456在线播放app| 成年女人看的毛片在线观看| 欧美性感艳星| 精品久久久久久久人妻蜜臀av| 九九久久精品国产亚洲av麻豆| 美女 人体艺术 gogo| h日本视频在线播放| 国产成人福利小说| 中文字幕人妻熟人妻熟丝袜美| 波多野结衣高清无吗| 欧美成人一区二区免费高清观看| 亚洲欧美日韩东京热| 国产黄a三级三级三级人| 欧美高清性xxxxhd video| 欧美+亚洲+日韩+国产| 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| 国产国拍精品亚洲av在线观看| 中文字幕免费在线视频6| 成人av一区二区三区在线看| 一边摸一边抽搐一进一小说| 日本熟妇午夜| 欧美一区二区亚洲| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 成年免费大片在线观看| 老熟妇仑乱视频hdxx| 亚洲国产精品成人综合色| 少妇的逼水好多| 国产精品一及| 99热这里只有精品一区| 97人妻精品一区二区三区麻豆| av女优亚洲男人天堂| 黄色日韩在线| 日日摸夜夜添夜夜爱| 久久久久久久久中文| 99久久中文字幕三级久久日本| 久久精品国产99精品国产亚洲性色| 免费观看在线日韩| 亚洲高清免费不卡视频| 亚洲欧美日韩东京热| 成年女人永久免费观看视频| 国产色婷婷99| 亚洲av免费在线观看| 一级av片app| 国产女主播在线喷水免费视频网站 | av在线亚洲专区| 精品乱码久久久久久99久播| 欧洲精品卡2卡3卡4卡5卡区| 3wmmmm亚洲av在线观看| 色视频www国产| 三级国产精品欧美在线观看| 色播亚洲综合网| 国产亚洲精品av在线| 亚洲av美国av| 在线播放国产精品三级| 村上凉子中文字幕在线| 午夜激情福利司机影院| 日本免费一区二区三区高清不卡| 1024手机看黄色片| 熟女电影av网| 欧美成人一区二区免费高清观看| 网址你懂的国产日韩在线| 我要看日韩黄色一级片| 看黄色毛片网站| 97超视频在线观看视频| 久久久久久久久久久丰满| 欧美日本亚洲视频在线播放| 国产成人aa在线观看| 色综合色国产| 欧美zozozo另类| 久久久久免费精品人妻一区二区| 久久精品综合一区二区三区| 日日摸夜夜添夜夜添小说| 亚洲一级一片aⅴ在线观看| 又爽又黄无遮挡网站| 一区二区三区高清视频在线| 99在线视频只有这里精品首页| 欧美日韩综合久久久久久| 亚洲七黄色美女视频| 赤兔流量卡办理| 99热只有精品国产| 精品久久久久久久久久久久久| 18禁黄网站禁片免费观看直播| 欧美3d第一页| 亚洲自拍偷在线| 热99re8久久精品国产| 亚洲中文字幕一区二区三区有码在线看| 99国产精品一区二区蜜桃av| 国产免费男女视频| 日韩成人av中文字幕在线观看 | 天天躁日日操中文字幕| 亚洲色图av天堂| 搡老妇女老女人老熟妇| 免费看光身美女| 成人漫画全彩无遮挡| 啦啦啦韩国在线观看视频| 99热只有精品国产| 国产精品一区二区免费欧美| 国产精品一及| 老女人水多毛片| 午夜a级毛片| 亚洲一区高清亚洲精品| 国产精品无大码| 在线看三级毛片| 又粗又爽又猛毛片免费看| 国产欧美日韩一区二区精品| 99视频精品全部免费 在线| 男女视频在线观看网站免费| 国产免费男女视频| 偷拍熟女少妇极品色| 中文字幕av在线有码专区| 麻豆国产97在线/欧美| 一本精品99久久精品77| 免费看a级黄色片| 噜噜噜噜噜久久久久久91| 三级毛片av免费| 欧美绝顶高潮抽搐喷水| 亚洲av中文av极速乱| 午夜精品一区二区三区免费看| 成年女人看的毛片在线观看| 欧美性感艳星| 天堂影院成人在线观看| 男女啪啪激烈高潮av片| 99riav亚洲国产免费| 一本精品99久久精品77| 91午夜精品亚洲一区二区三区| 91在线精品国自产拍蜜月| 免费观看精品视频网站| 亚洲精品一区av在线观看| 成人无遮挡网站| 精品日产1卡2卡| 日韩精品青青久久久久久| 毛片一级片免费看久久久久| 亚洲无线观看免费| 亚州av有码| 免费看日本二区| 亚洲精品亚洲一区二区| 男人的好看免费观看在线视频| 久久久久久久久中文| 国产精品免费一区二区三区在线| 美女高潮的动态| 搡老熟女国产l中国老女人| 国产蜜桃级精品一区二区三区| 午夜免费激情av| 日本一本二区三区精品| 欧美一区二区精品小视频在线| 偷拍熟女少妇极品色| 国产 一区 欧美 日韩| 亚洲精品日韩在线中文字幕 | 乱人视频在线观看| 精品人妻一区二区三区麻豆 | а√天堂www在线а√下载| 老女人水多毛片| 亚洲成人精品中文字幕电影| 精品久久久久久久末码| 精品久久久久久久人妻蜜臀av| 国产蜜桃级精品一区二区三区| 桃色一区二区三区在线观看| 在线天堂最新版资源| 精品免费久久久久久久清纯| 成人欧美大片| 免费电影在线观看免费观看| 夜夜夜夜夜久久久久| 99视频精品全部免费 在线| 久久久久国内视频| 国内精品美女久久久久久| 1024手机看黄色片| 国产三级中文精品| 成人二区视频| 国产高清激情床上av| 日本在线视频免费播放| 日韩三级伦理在线观看| 成人一区二区视频在线观看| 一级毛片电影观看 | 99在线人妻在线中文字幕| 天堂√8在线中文| 毛片女人毛片| 国产亚洲精品综合一区在线观看| 搡老岳熟女国产| 亚洲欧美日韩高清专用| 欧美激情国产日韩精品一区| 深夜精品福利| 晚上一个人看的免费电影| 日韩精品青青久久久久久| 日本一二三区视频观看| 亚洲性久久影院| 哪里可以看免费的av片| 99久国产av精品| av在线蜜桃| 性插视频无遮挡在线免费观看| 国产亚洲精品综合一区在线观看| 高清日韩中文字幕在线| 最新在线观看一区二区三区| 精品一区二区三区人妻视频| 亚洲高清免费不卡视频| 国产亚洲精品av在线| 99在线人妻在线中文字幕| 直男gayav资源| 国产精品一区二区性色av| 黄色欧美视频在线观看| 国产午夜福利久久久久久| 国产一区亚洲一区在线观看| 在线观看午夜福利视频| 国产在视频线在精品| 中文字幕av成人在线电影| 性色avwww在线观看| 不卡视频在线观看欧美| 国产精品日韩av在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 搡老岳熟女国产| 精品无人区乱码1区二区| 国内精品美女久久久久久| 精品久久久久久久末码| 久久久久久久久久成人| 亚洲精品粉嫩美女一区| 老女人水多毛片| 在线观看一区二区三区| 免费看av在线观看网站| 一进一出抽搐gif免费好疼| 久久久久国内视频| 亚洲欧美日韩东京热| 我要看日韩黄色一级片| 日本一本二区三区精品| 九色成人免费人妻av| 成人欧美大片| 天美传媒精品一区二区| 国产精品亚洲一级av第二区| 国内精品宾馆在线| 99国产极品粉嫩在线观看| 如何舔出高潮| 日韩欧美三级三区| 久久久久久久久中文| 亚洲av第一区精品v没综合| 白带黄色成豆腐渣| 国产精品嫩草影院av在线观看| ponron亚洲| 欧美3d第一页| 99热网站在线观看| 99精品在免费线老司机午夜| 18禁在线播放成人免费| 赤兔流量卡办理| 亚洲内射少妇av| 成人一区二区视频在线观看| 国产成年人精品一区二区| 精品午夜福利视频在线观看一区| 亚洲精品456在线播放app| 免费观看在线日韩| 69人妻影院| 赤兔流量卡办理| 成人漫画全彩无遮挡| 国产成人一区二区在线| 久久久久久久午夜电影| 精品久久久噜噜| 男女做爰动态图高潮gif福利片| 成人二区视频| 日日啪夜夜撸| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| av视频在线观看入口| 国产精品国产高清国产av| 69av精品久久久久久| 看片在线看免费视频| 久久亚洲精品不卡| 国产美女午夜福利| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩一区二区精品| 91久久精品电影网| 欧美xxxx性猛交bbbb| 免费在线观看影片大全网站| 国产精品一区www在线观看| 成人av一区二区三区在线看| 国产精品一区www在线观看| 一进一出抽搐gif免费好疼| 人妻少妇偷人精品九色| 亚洲激情五月婷婷啪啪| 亚洲人与动物交配视频| 一区二区三区高清视频在线| 亚洲高清免费不卡视频| 国语自产精品视频在线第100页| 亚洲五月天丁香| 精品熟女少妇av免费看| 在线免费十八禁| 亚洲激情五月婷婷啪啪| 99热全是精品| 日韩欧美精品v在线| 精品午夜福利视频在线观看一区| eeuss影院久久| 久久精品国产鲁丝片午夜精品| 国产蜜桃级精品一区二区三区| 亚洲图色成人| 久久精品夜夜夜夜夜久久蜜豆| 国产乱人偷精品视频| av免费在线看不卡| 亚洲在线观看片| 国产男人的电影天堂91| 国产在线男女| or卡值多少钱| 国产成人一区二区在线| 午夜亚洲福利在线播放| 不卡视频在线观看欧美| 两个人视频免费观看高清| 国产精品一区二区三区四区免费观看 | 97超视频在线观看视频| 啦啦啦韩国在线观看视频| 久久中文看片网| 久久精品国产鲁丝片午夜精品| 精华霜和精华液先用哪个| 亚洲成a人片在线一区二区| 精品人妻视频免费看| 如何舔出高潮| 久久欧美精品欧美久久欧美| 99热全是精品| 日日啪夜夜撸| 97超级碰碰碰精品色视频在线观看| 一区二区三区高清视频在线| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 久久久国产成人免费| 不卡视频在线观看欧美| 丰满乱子伦码专区| 国产精品一区二区免费欧美| 国产精品,欧美在线| 又黄又爽又刺激的免费视频.| 天堂√8在线中文| 欧美高清成人免费视频www| 少妇裸体淫交视频免费看高清| 白带黄色成豆腐渣| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 日本黄色片子视频| 午夜a级毛片| 国产在线精品亚洲第一网站| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜添av毛片| 精品免费久久久久久久清纯| 三级男女做爰猛烈吃奶摸视频| 春色校园在线视频观看| 1000部很黄的大片| 成年免费大片在线观看| 国产高清激情床上av| 久久久精品大字幕| 欧美性猛交黑人性爽| 成人永久免费在线观看视频| 日本免费a在线| 国产精品国产高清国产av| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在 | 欧美一区二区国产精品久久精品| 亚洲无线观看免费| 日日摸夜夜添夜夜爱| 国产高清不卡午夜福利| 麻豆国产av国片精品| 婷婷精品国产亚洲av在线| 国产淫片久久久久久久久| 不卡一级毛片| 亚洲国产色片| 日韩精品青青久久久久久| 在线免费观看不下载黄p国产| av中文乱码字幕在线| 久久久国产成人精品二区| 欧美日韩精品成人综合77777| 国产精品永久免费网站| 日韩欧美三级三区| 日韩一区二区视频免费看| 深爱激情五月婷婷| 99久久九九国产精品国产免费| 99riav亚洲国产免费| 亚洲一区二区三区色噜噜| 少妇人妻精品综合一区二区 | 亚洲电影在线观看av| 国产精品一区www在线观看| 久久久午夜欧美精品| 欧美3d第一页| 久久婷婷人人爽人人干人人爱| 久久久久久久亚洲中文字幕| 久久久久久久午夜电影| 久久久久久久久久成人| 天天躁日日操中文字幕| 日本黄大片高清| 在线看三级毛片| 国产精品av视频在线免费观看| 成人综合一区亚洲| 欧美性感艳星| .国产精品久久| 能在线免费观看的黄片| 三级经典国产精品| 久久人人爽人人爽人人片va| 少妇的逼水好多| 亚洲美女视频黄频| 国产欧美日韩精品一区二区| 国产精品免费一区二区三区在线| 赤兔流量卡办理| 在线免费观看的www视频| 国产伦精品一区二区三区四那| 99九九线精品视频在线观看视频| 日本熟妇午夜| 男人的好看免费观看在线视频| 熟女电影av网| 久久久国产成人免费| 美女免费视频网站| 又黄又爽又刺激的免费视频.|