• <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蜜桃精品久久| 久久久成人免费电影| 国产精品人妻久久久影院| 超碰97精品在线观看| 国产精品国产三级国产专区5o| 日韩一区二区三区影片| 国产精品久久久久成人av| 国产高清有码在线观看视频| 午夜福利影视在线免费观看| 国产乱来视频区| 中文字幕制服av| 欧美一级a爱片免费观看看| 一级黄片播放器| 亚洲国产精品国产精品| a级一级毛片免费在线观看| 男人添女人高潮全过程视频| 国产一区二区三区综合在线观看 | 久久av网站| 亚洲第一av免费看| 99久久精品一区二区三区| 精品视频人人做人人爽| 亚洲人成网站在线观看播放| 高清视频免费观看一区二区| 少妇高潮的动态图| 网址你懂的国产日韩在线| 色网站视频免费| 天堂中文最新版在线下载| 国产免费福利视频在线观看| 亚洲激情五月婷婷啪啪| 欧美xxⅹ黑人| 高清视频免费观看一区二区| 极品少妇高潮喷水抽搐| 国产69精品久久久久777片| 最后的刺客免费高清国语| 日本一二三区视频观看| 日韩 亚洲 欧美在线| 人妻系列 视频| 一级黄片播放器| 免费看av在线观看网站| 最近最新中文字幕大全电影3| av在线老鸭窝| 不卡视频在线观看欧美| 国产成人精品久久久久久| 国产高潮美女av| 国产高潮美女av| 男女啪啪激烈高潮av片| 日本av免费视频播放| 亚洲成人中文字幕在线播放| 爱豆传媒免费全集在线观看| 亚洲国产欧美人成| 高清不卡的av网站| 国产人妻一区二区三区在| 51国产日韩欧美| 青春草国产在线视频| 免费黄色在线免费观看| 日韩一区二区三区影片| 久久精品国产亚洲av涩爱| 日日啪夜夜撸| 嫩草影院新地址| 热re99久久精品国产66热6| 18+在线观看网站| 日韩一本色道免费dvd| 高清av免费在线| 久久国内精品自在自线图片| 亚洲丝袜综合中文字幕| 亚洲精品亚洲一区二区| 青青草视频在线视频观看| 天堂俺去俺来也www色官网| 国产又色又爽无遮挡免| 亚洲精品日韩av片在线观看| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 日本免费在线观看一区| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品成人久久小说| 人人妻人人看人人澡| 国内少妇人妻偷人精品xxx网站| 亚洲美女视频黄频| 国国产精品蜜臀av免费| 久久久久久久久大av| 七月丁香在线播放| 偷拍熟女少妇极品色| 嫩草影院入口| 91精品国产九色| 国产精品国产av在线观看| 国产在线视频一区二区| 国产黄色视频一区二区在线观看| 久久韩国三级中文字幕| 下体分泌物呈黄色| 国产深夜福利视频在线观看| 免费黄频网站在线观看国产| 久久久a久久爽久久v久久| 日本黄色日本黄色录像| 欧美激情极品国产一区二区三区 | 国产中年淑女户外野战色| 国产成人aa在线观看| 嘟嘟电影网在线观看| 久热这里只有精品99| 久久亚洲国产成人精品v| av专区在线播放| 国产成人a∨麻豆精品| 能在线免费看毛片的网站| 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 国产爽快片一区二区三区| 伊人久久国产一区二区| 国产精品国产av在线观看| 99视频精品全部免费 在线| av国产免费在线观看| 国产精品爽爽va在线观看网站| 一本久久精品| 亚洲精品成人av观看孕妇| 极品少妇高潮喷水抽搐| 欧美+日韩+精品| 欧美亚洲 丝袜 人妻 在线| 久久久国产一区二区| 春色校园在线视频观看| 国产欧美日韩一区二区三区在线 | 中文字幕av成人在线电影| 大香蕉97超碰在线| 国产免费一区二区三区四区乱码| 狂野欧美白嫩少妇大欣赏| 亚洲国产毛片av蜜桃av| 久热这里只有精品99| 欧美xxⅹ黑人| 国产精品一及| 观看美女的网站| 在线 av 中文字幕| 国产精品伦人一区二区| 日本欧美视频一区| 亚洲久久久国产精品| 1000部很黄的大片| 国产国拍精品亚洲av在线观看| 欧美zozozo另类| 欧美丝袜亚洲另类| 99热这里只有是精品50| 美女脱内裤让男人舔精品视频| 亚洲av欧美aⅴ国产| 国产男女内射视频| 青青草视频在线视频观看| 亚洲精品日本国产第一区| 日日啪夜夜爽| 九草在线视频观看| 成人毛片a级毛片在线播放| 亚洲伊人久久精品综合| 亚洲欧美日韩卡通动漫| 日本黄色片子视频| 国产午夜精品一二区理论片| 你懂的网址亚洲精品在线观看| 久久精品国产自在天天线| 青春草亚洲视频在线观看| 亚洲人与动物交配视频| 国产在线男女| 欧美精品人与动牲交sv欧美| 亚洲av国产av综合av卡| 国产一区亚洲一区在线观看| 嫩草影院入口| 欧美成人精品欧美一级黄| 精品视频人人做人人爽| 国产男女超爽视频在线观看| 草草在线视频免费看| 久久久久久久精品精品| 欧美丝袜亚洲另类| 欧美日本视频| 亚洲精品日韩在线中文字幕| 国国产精品蜜臀av免费| 水蜜桃什么品种好| 99久久精品一区二区三区| 在线天堂最新版资源| av国产免费在线观看| 日韩强制内射视频| 日本欧美国产在线视频| 纯流量卡能插随身wifi吗| 女人久久www免费人成看片| 国产真实伦视频高清在线观看| 99热网站在线观看| 观看美女的网站| 国产无遮挡羞羞视频在线观看| 国产精品人妻久久久久久| 午夜福利在线观看免费完整高清在| 国产一区二区在线观看日韩| 女人十人毛片免费观看3o分钟| 日韩一区二区三区影片| 蜜桃久久精品国产亚洲av| 午夜视频国产福利| 亚洲av二区三区四区| 日韩国内少妇激情av| 久久鲁丝午夜福利片| 亚洲精品一二三| 联通29元200g的流量卡| a级一级毛片免费在线观看| 亚洲欧美成人精品一区二区| 人人妻人人爽人人添夜夜欢视频 | 啦啦啦中文免费视频观看日本| 久久国内精品自在自线图片| 超碰97精品在线观看| 妹子高潮喷水视频| 精品久久久噜噜| 欧美老熟妇乱子伦牲交| 大陆偷拍与自拍| 只有这里有精品99| 色婷婷av一区二区三区视频| 亚洲国产精品一区三区| 热99国产精品久久久久久7| 中文在线观看免费www的网站| 午夜福利高清视频| 99九九线精品视频在线观看视频| 国产欧美日韩一区二区三区在线 | 乱码一卡2卡4卡精品| 亚洲av欧美aⅴ国产| 国产综合精华液| 亚洲图色成人| 午夜福利网站1000一区二区三区| 亚洲伊人久久精品综合| 美女cb高潮喷水在线观看| 在线观看免费视频网站a站| 亚洲一区二区三区欧美精品| 免费大片18禁| 国产在线一区二区三区精| 两个人的视频大全免费| 国产视频首页在线观看| 少妇人妻精品综合一区二区| 国产精品人妻久久久久久| 欧美日韩视频精品一区| 一区二区三区四区激情视频| 久久精品国产亚洲av天美| 一个人看的www免费观看视频| 久久久成人免费电影| 一本一本综合久久| 欧美日本视频| 成人黄色视频免费在线看| 午夜福利视频精品| 亚洲人成网站高清观看| 男的添女的下面高潮视频| 一级a做视频免费观看| 狂野欧美激情性bbbbbb| 日日摸夜夜添夜夜添av毛片| 国产午夜精品久久久久久一区二区三区| 久久97久久精品| 热re99久久精品国产66热6| 91狼人影院| 精华霜和精华液先用哪个| 午夜老司机福利剧场| 国产精品麻豆人妻色哟哟久久| 午夜福利网站1000一区二区三区| 男女国产视频网站| 欧美丝袜亚洲另类| 内射极品少妇av片p| 亚洲欧美成人综合另类久久久| 亚洲精品国产av成人精品| av专区在线播放| av免费观看日本| 免费大片18禁| 热99国产精品久久久久久7| 欧美日韩亚洲高清精品| 日韩 亚洲 欧美在线| 天堂俺去俺来也www色官网| 一级爰片在线观看| 亚洲,欧美,日韩| 亚洲av二区三区四区| 熟女电影av网| 街头女战士在线观看网站| 亚洲av男天堂| 日本猛色少妇xxxxx猛交久久| 在线观看三级黄色| 日本猛色少妇xxxxx猛交久久| 日韩av免费高清视频| 日韩中文字幕视频在线看片 | 久久精品久久久久久久性| 免费看日本二区| 纯流量卡能插随身wifi吗| 九九爱精品视频在线观看| 国产在线一区二区三区精| 精品久久久噜噜| 久久久精品免费免费高清| 一个人看的www免费观看视频| 国产精品久久久久久久电影| av在线蜜桃| 亚洲欧美精品专区久久| 国产成人精品婷婷| 黄色视频在线播放观看不卡| a级一级毛片免费在线观看| 亚洲av国产av综合av卡| 99九九线精品视频在线观看视频| 国产精品久久久久久精品古装| 久久久欧美国产精品| 久久99蜜桃精品久久| 国产黄片视频在线免费观看| 亚洲精品乱码久久久v下载方式| av天堂中文字幕网| av免费在线看不卡| 男人和女人高潮做爰伦理| 日韩电影二区| 婷婷色麻豆天堂久久| 精品少妇久久久久久888优播| 性色av一级| 国产一区二区三区综合在线观看 | 一本久久精品| 国产精品99久久久久久久久| 精品午夜福利在线看| 中文字幕制服av| 男女边摸边吃奶| 久久国产亚洲av麻豆专区| 精品人妻一区二区三区麻豆| 一区在线观看完整版| 狂野欧美激情性bbbbbb| 欧美xxⅹ黑人| 97热精品久久久久久| 亚洲欧美日韩东京热| 舔av片在线| 国产欧美日韩精品一区二区| 97超视频在线观看视频| 蜜桃久久精品国产亚洲av| 日本av免费视频播放| 97精品久久久久久久久久精品| 久久人妻熟女aⅴ| 午夜福利视频精品| 久久精品国产亚洲av涩爱| 日韩强制内射视频| 青春草视频在线免费观看| 久久久久久伊人网av| 久久久久久久国产电影| 国产黄频视频在线观看| 欧美精品亚洲一区二区| 亚洲av日韩在线播放| 国模一区二区三区四区视频| 黄色配什么色好看| 欧美日韩视频高清一区二区三区二| 国产伦精品一区二区三区四那| 国产一区亚洲一区在线观看| 日日啪夜夜爽| 久久精品久久久久久噜噜老黄| 妹子高潮喷水视频| 国产极品天堂在线| 国产精品一区www在线观看| 亚洲欧美一区二区三区黑人 | 久久久久久久国产电影| 久久国产乱子免费精品| 建设人人有责人人尽责人人享有的 | 国产永久视频网站| 大香蕉97超碰在线| 在线天堂最新版资源| 国产免费一区二区三区四区乱码| 女人久久www免费人成看片| 日本爱情动作片www.在线观看| kizo精华| av一本久久久久| 人妻 亚洲 视频| 精品久久久噜噜| 七月丁香在线播放| 亚洲精品一区蜜桃| 成人国产麻豆网| 女的被弄到高潮叫床怎么办| 日本av免费视频播放| 天堂俺去俺来也www色官网| 亚洲精品日韩在线中文字幕| 国产美女午夜福利| 交换朋友夫妻互换小说| 在线 av 中文字幕| 身体一侧抽搐| 国产亚洲精品久久久com| 国产有黄有色有爽视频| 丰满少妇做爰视频| 亚洲色图av天堂| 久久精品久久精品一区二区三区| 欧美人与善性xxx| 亚洲国产精品一区三区| 尾随美女入室| av一本久久久久| 另类亚洲欧美激情| 欧美国产精品一级二级三级 | 男女免费视频国产| 我要看黄色一级片免费的| 国产精品久久久久久久电影| 国精品久久久久久国模美| 一个人看视频在线观看www免费| 视频中文字幕在线观看| 91狼人影院| 青青草视频在线视频观看| 亚洲性久久影院| 国产精品一区二区在线不卡| 亚洲精品一二三| 久久影院123| 日日摸夜夜添夜夜添av毛片| 丰满迷人的少妇在线观看| 国产精品国产av在线观看| 伦理电影免费视频| 另类亚洲欧美激情| 在线观看国产h片| 国产av一区二区精品久久 | 亚洲人成网站高清观看| 大码成人一级视频| 美女中出高潮动态图| 中文精品一卡2卡3卡4更新| 国产成人精品一,二区| 亚洲一级一片aⅴ在线观看| videos熟女内射| 成人国产av品久久久| 我要看黄色一级片免费的| 国产熟女欧美一区二区| 18禁裸乳无遮挡免费网站照片| 人妻制服诱惑在线中文字幕| 欧美zozozo另类| 久久鲁丝午夜福利片| 免费观看性生交大片5| 亚洲最大成人中文| 91精品伊人久久大香线蕉| 国产黄片美女视频| 王馨瑶露胸无遮挡在线观看| 国产熟女欧美一区二区| 日本欧美国产在线视频| 亚洲av日韩在线播放| 老司机影院成人| 人妻夜夜爽99麻豆av| 男人爽女人下面视频在线观看| 久久人妻熟女aⅴ| 日本欧美国产在线视频| 中文天堂在线官网| 国产成人aa在线观看| 亚洲精品久久久久久婷婷小说| 建设人人有责人人尽责人人享有的 | 亚洲欧洲日产国产| 成年女人在线观看亚洲视频| 激情五月婷婷亚洲| 国产成人一区二区在线| 免费观看性生交大片5| 成人免费观看视频高清| 六月丁香七月| 亚洲国产高清在线一区二区三| av不卡在线播放| 欧美精品一区二区大全| 韩国av在线不卡| 久久国产乱子免费精品| 欧美成人午夜免费资源| 丰满人妻一区二区三区视频av| 亚洲人成网站高清观看| 亚洲欧洲日产国产| 亚洲国产精品999| av国产精品久久久久影院| 三级经典国产精品| 午夜福利影视在线免费观看| 两个人的视频大全免费| 国产中年淑女户外野战色| 国产精品偷伦视频观看了| 蜜桃亚洲精品一区二区三区| 亚洲国产av新网站| tube8黄色片| 一级片'在线观看视频| 国产在线男女| 午夜免费鲁丝| 免费av不卡在线播放| 久久国内精品自在自线图片| 国产女主播在线喷水免费视频网站| 麻豆成人午夜福利视频| 免费高清在线观看视频在线观看| 精品99又大又爽又粗少妇毛片| 99久久中文字幕三级久久日本| 国产淫语在线视频| 亚洲av电影在线观看一区二区三区| 久久鲁丝午夜福利片| 国模一区二区三区四区视频| 午夜福利高清视频| 极品教师在线视频| 色吧在线观看| av线在线观看网站| 最黄视频免费看| 亚洲精品456在线播放app| 欧美精品一区二区免费开放| 韩国高清视频一区二区三区| 亚洲精品日韩在线中文字幕| 一区二区三区乱码不卡18| 久久久欧美国产精品| 一级毛片 在线播放| 欧美丝袜亚洲另类| 国产成人精品久久久久久| 久久久久人妻精品一区果冻| 亚洲精品一区蜜桃| 亚洲欧美精品专区久久| 国内精品宾馆在线| 国产探花极品一区二区| 精品久久久久久电影网| 少妇 在线观看| 夫妻性生交免费视频一级片| 伦理电影免费视频| 国产欧美日韩一区二区三区在线 | 国产男女内射视频| 国产视频内射| 日本-黄色视频高清免费观看| 精品久久国产蜜桃| 熟妇人妻不卡中文字幕| 国产永久视频网站| 毛片一级片免费看久久久久| 亚洲丝袜综合中文字幕| 亚洲美女视频黄频| 99热6这里只有精品| 大陆偷拍与自拍| 在线观看国产h片| 免费少妇av软件| 观看av在线不卡| 中文乱码字字幕精品一区二区三区| 亚洲在久久综合| 国产亚洲91精品色在线| 久久久a久久爽久久v久久| 国产欧美另类精品又又久久亚洲欧美| 亚洲av在线观看美女高潮| 日韩中字成人| 香蕉精品网在线| 久久国产乱子免费精品| 欧美成人精品欧美一级黄| 亚洲熟女精品中文字幕| 性色av一级| 中文欧美无线码| 国产无遮挡羞羞视频在线观看| 亚洲精品,欧美精品| 成年免费大片在线观看| 天美传媒精品一区二区| 18+在线观看网站| 国产精品爽爽va在线观看网站| 欧美高清成人免费视频www| 男女国产视频网站| 国产精品不卡视频一区二区| 国产精品女同一区二区软件| 欧美人与善性xxx| 国产高清有码在线观看视频| 亚洲中文av在线| 欧美zozozo另类| 波野结衣二区三区在线| videossex国产| 精品国产一区二区三区久久久樱花 | 少妇人妻精品综合一区二区| 国产精品秋霞免费鲁丝片| 人妻 亚洲 视频| 欧美成人a在线观看| 国产永久视频网站| 国产久久久一区二区三区| 一级毛片电影观看| 欧美最新免费一区二区三区| 少妇 在线观看| 亚洲国产精品专区欧美| 久久亚洲国产成人精品v| 欧美成人一区二区免费高清观看| 亚洲精品乱码久久久久久按摩| 亚洲人成网站在线播| 精品久久久久久久久av| 三级经典国产精品| 久久人人爽人人片av| 日日撸夜夜添| 高清视频免费观看一区二区| 久久国内精品自在自线图片| 亚洲天堂av无毛| 男女无遮挡免费网站观看| 美女xxoo啪啪120秒动态图| 99视频精品全部免费 在线| 在线免费观看不下载黄p国产| 少妇人妻久久综合中文| 免费看不卡的av| 国语对白做爰xxxⅹ性视频网站| 一区二区三区精品91| 久久久a久久爽久久v久久| 观看美女的网站| 亚洲精品aⅴ在线观看| 99热这里只有精品一区| 成人特级av手机在线观看| 黑丝袜美女国产一区| 成年免费大片在线观看| 777米奇影视久久| 少妇高潮的动态图| 一级黄片播放器| 亚洲一级一片aⅴ在线观看| 99re6热这里在线精品视频| 在线 av 中文字幕| 日韩精品有码人妻一区| 18+在线观看网站| 美女cb高潮喷水在线观看| 哪个播放器可以免费观看大片| 久久精品熟女亚洲av麻豆精品| 卡戴珊不雅视频在线播放| 久久久午夜欧美精品| 国产色爽女视频免费观看| av国产精品久久久久影院| 在线观看av片永久免费下载| 老熟女久久久| 青青草视频在线视频观看| 国产女主播在线喷水免费视频网站| 九九爱精品视频在线观看| 91久久精品国产一区二区三区| 国产有黄有色有爽视频| 久久精品久久久久久久性| 日韩大片免费观看网站| 国产av一区二区精品久久 | 亚洲av在线观看美女高潮| 国产男人的电影天堂91| 国产精品人妻久久久久久| 免费观看无遮挡的男女| 成年美女黄网站色视频大全免费 | 中国三级夫妇交换| 国产欧美日韩精品一区二区| 日本av手机在线免费观看| 人体艺术视频欧美日本| 亚洲人与动物交配视频| 最黄视频免费看| 在线观看一区二区三区激情| 狂野欧美激情性xxxx在线观看| 黄色配什么色好看| 国产精品久久久久成人av| 国产亚洲5aaaaa淫片| 高清毛片免费看| 国产精品免费大片| 伊人久久精品亚洲午夜| 免费观看的影片在线观看| av在线观看视频网站免费| av线在线观看网站| 国产精品欧美亚洲77777| 国产在线视频一区二区|