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

    基于粒子群優(yōu)化算法和變步長擾動(dòng)觀察法的局部陰影情況下MPPT控制

    2016-10-12 05:17:44張慶超
    電源學(xué)報(bào) 2016年3期
    關(guān)鍵詞:觀察法輸出功率陰影

    馬 昊,張慶超

    (天津大學(xué)電氣與自動(dòng)化工程學(xué)院,天津300072)

    基于粒子群優(yōu)化算法和變步長擾動(dòng)觀察法的局部陰影情況下MPPT控制

    馬昊,張慶超

    (天津大學(xué)電氣與自動(dòng)化工程學(xué)院,天津300072)

    針對(duì)光伏陣列的輸出特性在局部陰影情況下具有高度非線性、時(shí)變性以及多個(gè)局部功率極值點(diǎn)等特點(diǎn),并導(dǎo)致傳統(tǒng)MPPT(maximum power point tracking)算法失效的問題,提出一種基于粒子群優(yōu)化算法和變步長擾動(dòng)觀察法的改進(jìn)MPPT算法。其中粒子群優(yōu)化算法用于系統(tǒng)啟動(dòng)和光照情況發(fā)生突變后迅速定位近似最大功率點(diǎn),變步長擾動(dòng)觀察法則根據(jù)實(shí)際狀況使光伏陣列精確穩(wěn)定在最大功率點(diǎn),以克服使用數(shù)學(xué)模型與實(shí)際輸出特性偏差或微小擾動(dòng)所導(dǎo)致的功率損失。通過建立Matlab/Simulink模型進(jìn)行仿真實(shí)驗(yàn),實(shí)驗(yàn)結(jié)果表明所提算法使光伏陣列在不同陰影情況下以及發(fā)生光照強(qiáng)度突變時(shí)都具有迅速精確的跟蹤能力。

    局部陰影;MPPT;粒子群算法;變步長擾動(dòng)觀察法

    引言

    由于太陽能具有取之不盡、用之不竭、無污染、無噪音等優(yōu)點(diǎn),光伏發(fā)電在新能源應(yīng)用領(lǐng)域受到越來越廣泛的重視并取得快速進(jìn)步,但發(fā)展中也存在許多問題亟需解決,其中太陽能轉(zhuǎn)換效率不高是主要問題[1]。工程上通常采用最大功率點(diǎn)跟蹤MPPT (maximum power poiot tracking)算法使光伏陣列輸出功率始終維持在最大功率點(diǎn)。在局部陰影等復(fù)雜光照條件下,光伏陣列P(UL)特性存在多個(gè)功率極值點(diǎn),導(dǎo)致傳統(tǒng)的擾動(dòng)觀察法、電導(dǎo)增量法等無法跟蹤全局功率最大值[2]。常見的發(fā)生局部遮擋時(shí)MPPT方法如下:復(fù)合MPPT算法是先把光伏陣列的工作點(diǎn)設(shè)在最大功率點(diǎn)的附近范圍,再利用擾動(dòng)觀察發(fā)或者增量電導(dǎo)法等常規(guī)算法進(jìn)行MPP的定位。因此在搜索局部最優(yōu)時(shí),在極值點(diǎn)處仍會(huì)有振蕩,造成功率損失;Fibonacci搜索法基于Fibonacci序列,該算法搜索精度和適應(yīng)環(huán)境突變能力較強(qiáng),但是收斂速度較慢,實(shí)用性不強(qiáng);短路電流脈沖法是利用光伏陣列的最佳輸出電流與短路電流存在的比例關(guān)系,周期性引入短路電流脈沖獲得光伏陣列輸出特性和短路電流來計(jì)算最佳輸出電流,但引入脈沖對(duì)系統(tǒng)的控制性能造成影響,在此期間光伏陣列無法輸出功率,因此造成能量的損失。

    PSO是一種新型優(yōu)化算法,在多峰函數(shù)優(yōu)化和多變量系統(tǒng)尋優(yōu)方面具有良好的性能,并且收斂快、精度高、不易陷入局部最優(yōu)點(diǎn),已經(jīng)有效地應(yīng)用到解決實(shí)際問題當(dāng)中[3]。但將其用在基于理論光伏電池?cái)?shù)學(xué)模型的局部陰影條件下MPPT控制時(shí),由于該模型與實(shí)際光伏特性有出入,導(dǎo)致輸出功率與真正功率最大值有一定差距,且當(dāng)系統(tǒng)有微小擾動(dòng)時(shí),此方法不能根據(jù)具體情況迅速跟蹤最大功率點(diǎn)[4];再由于算法初始化的隨機(jī)性,對(duì)收斂速度有一定影響,可能影響MPPT控制效果。變步長擾動(dòng)觀察法和常規(guī)擾動(dòng)觀察法一樣,在復(fù)雜光照條件下容易陷入局部最優(yōu)點(diǎn),但在最大功率點(diǎn)附近的跟蹤精度、發(fā)生微小擾動(dòng)后響應(yīng)速度以及穩(wěn)態(tài)功率損失等方面體現(xiàn)出更好的適用性[5]。

    因此,本文結(jié)合光伏陣列在局部陰影條件下的工作特征,提出基于PSO和變步長擾動(dòng)觀察法的復(fù)合MPPT控制方法,利用PSO快速定位近似最大功率點(diǎn),然后通過變步長擾動(dòng)觀察法尋找準(zhǔn)確值。搭建Matlab/Simulink仿真模型進(jìn)行實(shí)驗(yàn),證明了該方法在局部陰影光照條件下以及陰影情況發(fā)生突變時(shí)都能夠快速準(zhǔn)確地定位、跟蹤最大功率點(diǎn),始終保持良好的控制性能。

    1 局部陰影條件下光伏陣列的等效模型

    1.1光伏電池?cái)?shù)學(xué)模型

    利用PN結(jié)接受太陽光照產(chǎn)生光生伏特效應(yīng)制成的光伏電池可以視作一個(gè)恒定電流源和一支正向二極管的并聯(lián)回路,其等效電路如圖1所示。其中Iph為光生電流,其值與光伏電池的受光面積、入射光的光照強(qiáng)度以及電池本體溫度成正比;Id為光伏電池內(nèi)部暗電流;Ish為通過旁路電阻Rsh的電流;UL、IL為負(fù)載RL上的電壓和電流;Rs為串聯(lián)電阻。由等效電路可得光伏電池的ILUL特性方程為

    式中:I0為太陽能板的逆向飽和電流;q為電子電荷量;k為波爾茲曼常數(shù);A為PN結(jié)曲線常數(shù);T為溫度。理想情況下Rsh→∞,式(1)中最后一項(xiàng)即可忽略。

    圖1 光伏電池等效電路Fig.1 Equivalent circuit of PV cell

    實(shí)際應(yīng)用中通過串、并聯(lián)電池單元來構(gòu)成光伏陣列以獲得期望的輸出電壓、電流和功率。在均勻光照條件下光伏陣列的ILUL特性方程[6]為

    式中,Ns、Np分別為串、并聯(lián)電池單元的個(gè)數(shù)。此時(shí)光伏陣列的輸出IL-UL、P-UL特性曲線如圖2、圖3所示。

    圖3 光伏電池輸出特性隨溫度的變化曲線(s=1 000 W/m2)Fig.3 Relationship of output characteristics between PV cell and temperture(s=1 000 W/m2)

    由圖可見,在均勻光照條件下PUL曲線具有單一峰值特性,并且光伏電池輸出功率隨著溫度升高和光照強(qiáng)度降低而減小,從而導(dǎo)致發(fā)電效率下降。

    1.2局部陰影條件下的光伏陣列模型

    一旦光伏陣列部分組件被遮擋,導(dǎo)致光照不均勻,則被遮擋的組件變成負(fù)載,可能產(chǎn)生熱斑效應(yīng)造成永久性損壞[7]。工程中通過并聯(lián)旁路二極管以避免熱斑效應(yīng),提高光伏陣列在局部陰影情況下的輸出功率。當(dāng)光伏陣列處于均勻光照情況下,各電池單元都正常工作,旁路二極管處于反向截止?fàn)顟B(tài);當(dāng)有電池單元被遮蔽時(shí),其光生電流減小,電壓變?yōu)樨?fù)壓,此時(shí)旁路二極管導(dǎo)通將其短路,防止被反向雪崩電流擊穿,提高輸出功率。并列旁路二極管的光伏陣列結(jié)構(gòu)如圖4所示。

    圖4 并聯(lián)旁路二極管的光伏陣列Fig.4 PV array with bypass diode

    此時(shí)式(2)不再適用,光伏陣列等效模型需重新建立。本文引入光線遮擋率用來表示對(duì)光線的遮擋程度,用E表示為

    式中:Sbehind為光伏電池單元被陰影遮擋后實(shí)際的光照強(qiáng)度;Sbefore為被遮擋前的光照強(qiáng)度。為便于比較,本文中Sbefore均為參考值l kW/m2。

    光生電流變?yōu)?/p>

    式中,Iph0為光伏電池在參考光照強(qiáng)度和參考溫度下的光生電流。

    2個(gè)光伏電池單元串聯(lián)的單串陣列結(jié)構(gòu)如圖5所示,圖中a單元陰影遮擋,因此Ipha<Iphb,Ia<Ib,具體工作模式有如下2種情況:

    (1)當(dāng)RL較大,負(fù)載電流IL≤Ia時(shí),電池單元a、b均能正常工作,串聯(lián)陣列的輸出電流由a、b共同提供,輸出電壓、功率為兩者輸出電壓、功率之和;

    (2)當(dāng)RL較小,負(fù)載電流IL>Ia時(shí),此時(shí)電池單元a不能輸出相應(yīng)的電流,被旁路二極管Da短路,于是串聯(lián)陣列的輸出電流由電池單元b提供,輸出功率也只由單元b提供。

    圖5 2個(gè)光伏單元串聯(lián)等效電路Fig.5 Equivalent circuit of two series-connected photovoltaic cells

    對(duì)應(yīng)于這種情況下的數(shù)學(xué)模型可表示為

    由式(5)建立的光伏電池等效模型的輸出特性如圖6所示。

    圖6 局部陰影條件下的輸出特性曲線Fig.6 Output characteristics curves of two series-connected PV cells under partial shading

    推廣到Ns個(gè)遮擋面積逐漸減小的電池單元串聯(lián)構(gòu)成的光伏陣列,各單元輸出電流分別為Ia、Ib、Ic…INs,且Ipha<Iphb<Iphc<…<IphNs,則其輸出特性為

    式中,Ix、Ux為單串陣列的輸出電流、電壓。

    2 基于粒子群算法和變步長擾動(dòng)觀察法的MPPT控制

    2.1粒子群算法

    PSO算法從鳥類捕食的模型中得到啟示并用于優(yōu)化問題的求解。在PSO算法中,每個(gè)粒子都有一個(gè)由被優(yōu)化的目標(biāo)函數(shù)f(x)決定的適應(yīng)度值,優(yōu)化問題的解為搜索空間中具有最優(yōu)適應(yīng)度值的某個(gè)粒子。種群內(nèi)的所有粒子通過不斷搜索最終收斂到該粒子周圍,粒子的搜索過程通過位置和速度來描述。在MPPT控制中定義光伏陣列輸出功率為目標(biāo)函數(shù),以光伏電池輸出電壓值來表示粒子的位置,粒子的速度是電壓的更新步長,最大功率點(diǎn)所對(duì)應(yīng)的電壓為全局最優(yōu)解。每個(gè)粒子的飛行速度和位置更新方程為[3]

    式中:i為第i個(gè)粒子;v(t)為粒子t時(shí)刻的速度;x(t)為粒子t時(shí)刻的位置,搜索位置范圍為[xmin,xmax];pbest(t)為個(gè)體最優(yōu)位置,表示個(gè)體所經(jīng)歷的位置中對(duì)應(yīng)的適用度值最好的位置;gbest(t)為全局最優(yōu)位置,表示當(dāng)前所有的粒子中適應(yīng)度值最好的位置;c1、c2為學(xué)習(xí)因子,一般取2;r1、r2為介于(0,1)之間的隨機(jī)數(shù),用來調(diào)節(jié)粒子向pbest(t)、gbest(t)的飛行步長;w為慣性因子,表示粒子歷史速度對(duì)當(dāng)前速度的影響。則個(gè)體最優(yōu)、全局最優(yōu)位置分別為

    式中,N為粒子的個(gè)數(shù)。

    PSO算法需要設(shè)置的參數(shù)包括種群粒子數(shù),各個(gè)粒子的初始位置以及粒子位置的搜索范圍和慣性因子。N的取值大小影響算法的收斂速度以及精度。對(duì)于標(biāo)準(zhǔn)PSO算法,種群通過在搜索空間內(nèi)隨機(jī)生成N個(gè)粒子來初始化,其隨機(jī)性可能帶來收斂速度損失。本文結(jié)合局部陰影條件下光伏陣列工作的特殊性對(duì)其改進(jìn):對(duì)于{m×n}的光伏陣列,最多對(duì)應(yīng)n+1個(gè)極值點(diǎn)[8],因此令N=n+1,其初始位置分別為U1,U2,…,UN。為了更快速地收斂到全局最優(yōu)解,首先令 Boost開關(guān)管占空比 D=0,得到Uoc_series,然后按式(11)設(shè)置初始電壓,N取5。粒子位置搜索范圍為[0,Uoc_series]。

    在算法開始階段,w取較大值可使算法不易陷入局部最優(yōu);到算法的后期,w取較小值可使收斂速度加快,使收斂更平穩(wěn),不至于出現(xiàn)震蕩現(xiàn)象。本文選擇對(duì)w進(jìn)行凹函數(shù)衰減處理[4],初始值為0.8。

    當(dāng)達(dá)到最大迭代次數(shù)時(shí)停止迭代,此時(shí)所有粒子中對(duì)應(yīng)功率最大的判定為最大功率Pm,其位置即為Um。而當(dāng)陰影情況或太陽光輻照度發(fā)生突變時(shí),光伏陣列輸出特性曲線也隨之變化,為了提高太陽能利用率需要重新啟動(dòng)粒子群算法求解此時(shí)最大功率點(diǎn)。因此系統(tǒng)需要實(shí)時(shí)檢測(cè)功率變化量ΔP,即

    式中:P實(shí)時(shí)檢測(cè)到的功率;Pm為之前判定的功率最大值。當(dāng)ΔP>0.1時(shí)重啟PSO算法。

    2.2變步長擾動(dòng)觀察法

    對(duì)于定步長擾動(dòng)觀察法,快速性和精確性存在不可調(diào)和的矛盾。擾動(dòng)步長取較小值時(shí),能夠穩(wěn)態(tài)在最大功率點(diǎn),但是跟蹤速度慢、效率低;而取較大值時(shí)會(huì)縮短跟蹤的時(shí)間,但會(huì)造成穩(wěn)態(tài)時(shí)在最大功率點(diǎn)附近的震蕩幅度增大,造成功率損失。因此本文采用步長逐次逼近的方法來搜索最大功率點(diǎn),具體擾動(dòng)觀察過程為:在PSO得到初始工作點(diǎn)的基礎(chǔ)上,以較小步長ΔUconst進(jìn)行擾動(dòng),如果在第i+1次擾動(dòng)之前,始終有Pn+1>Pn(n<i),說明當(dāng)前工作點(diǎn)在最大功率點(diǎn)的左側(cè),則繼續(xù)以該步長搜索;當(dāng)出現(xiàn)Pi+1<Pi時(shí),說明當(dāng)前工作點(diǎn)變化到最大功率點(diǎn)的右側(cè),由此推斷出最大功率點(diǎn)的范圍應(yīng)在當(dāng)前工作點(diǎn)的2ΔUconst范圍內(nèi)。接著以ΔUconst/2為步長進(jìn)行反方向擾動(dòng),直到出現(xiàn)功率方向的下一次改變,此時(shí)最大功率點(diǎn)的范圍在ΔUconst范圍以內(nèi),即跟蹤精度提高1倍。重復(fù)上述過程,每當(dāng)發(fā)生功率方向改變則步長減半,并且改變擾動(dòng)方向繼續(xù)搜索,直到最大功率點(diǎn)落在當(dāng)前工作點(diǎn)的給定精度范圍內(nèi)時(shí),判定為搜索到了最大功率點(diǎn)。

    變步長擾動(dòng)觀察法的步長計(jì)算公式為

    顯然,逐步逼近法的精度在擾動(dòng)搜索過程中以指數(shù)形式提高,較好地解決了定步長擾動(dòng)法跟蹤速度和精度間的矛盾,使系統(tǒng)盡快穩(wěn)定在最大功率點(diǎn)。與此同時(shí),當(dāng)光照強(qiáng)度緩慢持續(xù)變化,ΔP較小,不滿足PSO算法突變重啟條件,但光伏陣列的工作效率也會(huì)降低,此時(shí)依據(jù)檢測(cè)到的前后兩次功率差的大小判定是否啟動(dòng)逐次逼近法進(jìn)行最大功率跟蹤,及時(shí)調(diào)整系統(tǒng)輸出功率,使太陽能利用率最大化。

    2.3控制策略

    本文提出的復(fù)合MPPT控制算法由2部分組成,第1部分是粒子群尋優(yōu)過程,第2部分是變步長擾動(dòng)觀察法跟蹤過程??刂葡到y(tǒng)流程如圖7所示。各部分的具體步驟如下。

    粒子群尋優(yōu)步驟:①利用式(11)求得各個(gè)粒子的初始電壓;②檢測(cè)各個(gè)電壓對(duì)應(yīng)輸出功率,即為粒子的初始pbest,所有粒子對(duì)應(yīng)的最大功率即為初始gbest;③按式(8)更新各個(gè)粒子對(duì)應(yīng)電壓,將此時(shí)的輸出功率與該粒子所經(jīng)歷的最大輸出功率進(jìn)行比較,如果當(dāng)前電壓對(duì)應(yīng)的功率較大,則更新pbest;④將所有粒子對(duì)應(yīng)的功率與當(dāng)前gbest比較,若存在新粒子的輸出功率大于gbest,則更新gbest;⑤重復(fù)步驟③和④,當(dāng)滿足最大迭代次數(shù)時(shí)輸出PPSO以及對(duì)應(yīng)的UPSO。

    變步長擾動(dòng)算法步驟:①將粒子群尋優(yōu)得到的電壓UPSO作為初始電壓;②以ΔUconst為步長進(jìn)行擾動(dòng),檢測(cè)功率變化ΔP的大小、方向,判斷是否需要調(diào)整步長;③利用新得到的步長重復(fù)上述擾動(dòng)過程,直到擾動(dòng)步長減小到某一特定值ε1時(shí),判定系統(tǒng)已穩(wěn)定在最大功率點(diǎn);而當(dāng)ΔP>0.1時(shí)判斷為外界環(huán)境突變,此時(shí)重啟PSO算法尋優(yōu)過程。

    圖7 基于復(fù)合算法的MPPT控制流程Fig.7 Flow chart of MPPT based on PSO and the improved P&O

    3 仿真實(shí)驗(yàn)

    利用Matlab/Simulink建立4個(gè)光伏電池單元組成的串聯(lián)陣列模型,其參考數(shù)據(jù)采用BP585U太陽能電池組件的標(biāo)稱參數(shù):在參考溫度25℃,參考光照1 000 W/m2的情況下,光伏電池單元Pm=85 W,Uoc=22.1 V,Isc=5 A,Um=18 V,Im=4.72 A[7]。在基于本文提出的復(fù)合算法的MPPT控制系統(tǒng)基礎(chǔ)上,分別進(jìn)行局部遮擋情況下MPPT尋優(yōu)以及光照強(qiáng)度突變后尋優(yōu)試驗(yàn),以及系統(tǒng)達(dá)到穩(wěn)態(tài)光照強(qiáng)度發(fā)生小幅擾動(dòng)后功率跟蹤試驗(yàn)來測(cè)驗(yàn)算法性能。在第1個(gè)實(shí)驗(yàn)中,陰影情況1對(duì)應(yīng)光線遮擋率E1=[0.0,0.0,0.2,0.4],陰影情況2對(duì)應(yīng)E2=[0.2,0.2,0.6,0.6]。首先在陰影情況1觀察算法的收斂速度和精度,檢測(cè)算法的靜態(tài)性能;而當(dāng)這種情況下系統(tǒng)達(dá)到穩(wěn)態(tài),光伏陣列輸出功率維持在最大功率點(diǎn)后陰影情況發(fā)生突變,觀察系統(tǒng)重新達(dá)到穩(wěn)態(tài)的時(shí)間和振蕩過程的強(qiáng)度以檢測(cè)算法的動(dòng)態(tài)性能??偡抡鏁r(shí)長為0.2 s,在t=0.1 s時(shí)由E1突變?yōu)镋2,對(duì)應(yīng)的光伏特性曲線如圖8所示;在t=0.15 s時(shí)發(fā)生輻射度小幅擾動(dòng),進(jìn)行穩(wěn)態(tài)擾動(dòng)實(shí)驗(yàn)。

    圖8 2種陰影情況下的光伏特性曲線Fig.8 Output characteristics of two partial shading modes

    在2種陰影情況下的PSO尋優(yōu)過程初始電壓設(shè)置如表1所示,其迭代搜索過程以及采用遺傳算法進(jìn)行對(duì)比實(shí)驗(yàn)的結(jié)果如圖9所示。遺傳算法采用二進(jìn)制編碼,種群個(gè)體數(shù)與PSO算法相同,變異率取0.06,交叉率取0.6。每一代進(jìn)化過程中種群中適應(yīng)度值最優(yōu)的個(gè)體不參加交叉變異操作,以保留優(yōu)良特性,圖9(b)描述的是遺傳算法獲得的最有適應(yīng)度值的變化過程。

    表1 2種陰影情況下各粒子的初始位置Tab.1 Initial position of agents in two partial shoding modes V

    圖9 PSO算法與遺傳算法尋優(yōu)對(duì)比實(shí)驗(yàn)Fig.9 Optimization experiment between PSO and GA

    由圖可見,PSO算法利用式(11)的初始化機(jī)制后能夠在不同陰影條件下,經(jīng)過較少的迭代次數(shù)即收斂到全局最優(yōu)點(diǎn):E1下經(jīng)過19次迭代收斂到全局最優(yōu)解,(Um1,Pm1)=(52.83 V,196.82 W);E2下經(jīng)過 12次迭代收斂到全局最優(yōu)解,(Um2,Pm2)=(35.36 V,131.27 W)。利用pbest、gbest指導(dǎo)進(jìn)化,這種記憶特性使得迭代過程始終向有利進(jìn)化的方向進(jìn)行。與之相比,遺傳算法不僅尋優(yōu)速度不及PSO算法,而且E1下在全局最優(yōu)點(diǎn)周圍搜索緩慢,有限的迭代次數(shù)內(nèi)無法穩(wěn)定在全局最優(yōu)點(diǎn);在陰影情況2下收斂到了非全局最優(yōu)點(diǎn),算法失效。

    在PSO輸出電壓參考值的基礎(chǔ)上,變步長擾動(dòng)觀察法跟蹤過程以及定步長擾動(dòng)觀察法對(duì)比實(shí)驗(yàn)結(jié)果如圖10所示。由圖10(a)可見,E1下經(jīng)過0.021 s即可穩(wěn)定在最大功率點(diǎn),系統(tǒng)超調(diào)和輸出功率振蕩都比定步長擾動(dòng)觀察法小,而且穩(wěn)態(tài)下輸出功率波動(dòng)小,光伏陣列損失功率少,效率較高;光照強(qiáng)度突變到E2后,又經(jīng)過0.008 s成功穩(wěn)定在新的最大功率點(diǎn)。當(dāng)t=0.15 s發(fā)生小幅擾動(dòng)后,先以較大步長進(jìn)行跟蹤,然后步長迅速減小,經(jīng)過0.006 s即回復(fù)到原始狀態(tài),功率振蕩時(shí)間短,證明了算法具有很好的快速性和精確性。

    由圖10(b)可見,在E1下功率穩(wěn)定在167.5 W附近,即第1個(gè)局部功率極大點(diǎn)而非全局最大點(diǎn),并且算法迭代過程中功率振蕩嚴(yán)重,超調(diào)大,對(duì)系統(tǒng)產(chǎn)生巨大沖擊;突變后,經(jīng)過0.014 s才能達(dá)到穩(wěn)定,并且突變前后穩(wěn)態(tài)始終存在振蕩,系統(tǒng)功率損失較大。發(fā)生小幅擾動(dòng)后,該算法雖然能夠識(shí)別擾動(dòng),但是經(jīng)過0.02 s才能穩(wěn)定,依然存在穩(wěn)態(tài)誤差。證明其無法在局部陰影條件下有效地使光伏陣列輸出功率跟蹤全局最大功率點(diǎn)以及在穩(wěn)態(tài)情況下無法使光伏發(fā)電系統(tǒng)效率最大化。

    圖10 復(fù)合控制算法與定步長擾動(dòng)觀察法尋優(yōu)對(duì)比實(shí)驗(yàn)Fig.10 Optimization experiment between the recommended MPPT and P&O

    4 結(jié)語

    本文針對(duì)局部陰影情況下光伏陣列輸出特性具有多個(gè)局部極值點(diǎn)的特性,提出一種基于PSO和變步長擾動(dòng)觀察法相結(jié)合的復(fù)合MPPT算法,并利用Matlab/Simukink搭建模型,與傳統(tǒng)的擾動(dòng)觀察法進(jìn)行對(duì)比仿真實(shí)驗(yàn)。由仿真結(jié)果可看出,采用PSO可快速定位到近似全局最大功率點(diǎn),響應(yīng)時(shí)間幾乎不受搜索范圍大小和極值點(diǎn)個(gè)數(shù)和分布位置的限制,尋優(yōu)效果優(yōu)于遺傳算法。變步長擾動(dòng)觀察法可以精確穩(wěn)定在最大功率點(diǎn),并且在陰影情況突變或者光照強(qiáng)度緩慢變化的情況下都能迅速有效地調(diào)整光伏陣列輸出功率,具有良好的穩(wěn)態(tài)和動(dòng)態(tài)性能,提高了太陽能利用率;而擾動(dòng)觀察法在局部遮擋條件下無法保證有效性。復(fù)合算法實(shí)現(xiàn)簡單,運(yùn)行可靠,減少了MPPT控制系統(tǒng)成本,提高光伏發(fā)電的經(jīng)濟(jì)性,體現(xiàn)出對(duì)環(huán)境的變化良好的適應(yīng)能力,具備工程應(yīng)用的條件。

    [1]趙爭(zhēng)鳴,劉建政,孫曉瑛,等.太陽能光伏發(fā)電及其應(yīng)用[M].北京∶科學(xué)出版社,2006:1-18.

    [2]劉艷莉,周航,程澤.基于粒子群優(yōu)化的光伏系統(tǒng)MPPT控制方法[J].計(jì)算機(jī)工程,2010,36(15)∶265-267.Liu Yanli,Zhou Hang,Cheng Ze.MPPT control method of PV system based on PSO[J].Computer Engineering,2010,36(15)∶265-267(in Chinese).

    [3]徐剛,瞿金平,楊智韜,等.一種改進(jìn)的自適應(yīng)粒子群優(yōu)化算法[J].華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,36(9)∶6-10.Xu Gang,Qu Jinping,Yang Zhitao,et al.An improved adaptive particle swarm optimization algorithm[J].Journal of South China University of Technology(Natural Science Edition),2008,36(9)∶6-10(in Chinese).

    [4]王雨,胡仁杰.基于粒子群優(yōu)化和爬山法的MPPT算法[J].太陽能學(xué)報(bào),2014,35(1)∶149-153.Wang Yu,Hu Renjie.MPPT algorithm based on particle swam optimization with hill climbing method[J].Acta En-ergiae Solaris Sinica,2014,35(1)∶149-153(in Chinese).

    [5]張興.太陽能光伏并網(wǎng)發(fā)電及其逆變控制[M].北京∶機(jī)械工業(yè)出版社,2011.

    [6]錢念書,劉闊,郭建業(yè),等.光伏電池建模及其輸出特性研究[J].電源學(xué)報(bào),2012,10(5)∶78-82.Qian Nianshu,Liu Kuo,Guo Jianye,et al.Modeling and output features simulation of photovoltaic cells[J].Journal of Power Supply,2012,10(5)∶78-82(in Chinese).

    [7]劉曉艷,祁新梅,鄭壽森,等.局部陰影條件下光伏陣列的建模與分析[J].電網(wǎng)技術(shù),2010,11∶192-197.Liu Xiaoyan,Qi Xinmei,Zheng Shousen,et al.Model and analysis of photovoltaic array under partial shading[J].Power System Technology,2010,11∶192-197(in Chinese).

    [8]朱艷偉,石新春,但揚(yáng)清,等.粒子群優(yōu)化算法在光伏陣

    MPPT Control Under Partial Shading Condition Based on PSO and Variable Step Size P&O

    MA Hao,ZHANG Qingchao
    (School of Electrical Engineering and Automation,Tianjin University,Tianjin 300072,China)

    For the photovoltaic array in the case of partial shading has high nonlinearity,time-varying uncertainty and several local peak power characteristics,which resulting in conventional MPPT(maximum power point tracking)algorithms failing to work,an algorithm was proposed based on particle swarm optimization(PSO)and variable step size perturbation and observation(P&O)method.PSO was used to quickly locate approximate maximum power point(MPP)when the MPPT system started or illumination intensity motution occurred,and variable step size P&O method could find precise MPP based on the actual status to overcome power loss caused by the deviation between mathematical model and the actual output or small pertubarions.Matlab/Simulink simulation model was established.The results show that this method has a high precision and rapidity tracking capabilities in the case of partial shading as well as illumination intensity motution.

    partial shading;maximum power point tracking(MPPT);particle swarm optimization;variable step size perturbation and observation method

    10.13234/j.issn.2095-2805.2016.3.94

    TM 615

    A

    2015-09-06

    猜你喜歡
    觀察法輸出功率陰影
    你來了,草就沒有了陰影
    文苑(2020年11期)2020-11-19 11:45:11
    觀察法求解至少有一個(gè)有理根的一元高次方程
    中國詩歌(2019年6期)2019-11-15 00:26:47
    一種改進(jìn)的變步長擾動(dòng)觀察法在光伏MPPT中的應(yīng)用
    讓光“驅(qū)走”陰影
    基于改進(jìn)擾動(dòng)觀察法的光伏MPPT研究
    陰影魔怪
    一種改進(jìn)的擾動(dòng)觀察法在光伏MPPT中的應(yīng)用
    適用于智能電網(wǎng)的任意波形輸出功率源
    基于雙層BP神經(jīng)網(wǎng)絡(luò)的光伏電站輸出功率預(yù)測(cè)
    avwww免费| 欧美日韩瑟瑟在线播放| 桃色一区二区三区在线观看| 国产精华一区二区三区| 国产精华一区二区三区| 久久久久久久午夜电影| 国产精品一及| 欧美午夜高清在线| 黄色女人牲交| 草草在线视频免费看| 亚洲国产精品999在线| 悠悠久久av| 日韩国内少妇激情av| 内射极品少妇av片p| 欧美+日韩+精品| 亚洲国产日韩欧美精品在线观看 | av福利片在线观看| 欧美激情久久久久久爽电影| 成年版毛片免费区| 国产精品98久久久久久宅男小说| 少妇高潮的动态图| 欧美乱色亚洲激情| 国产免费男女视频| 中文字幕久久专区| 天堂网av新在线| 两个人看的免费小视频| 成人18禁在线播放| 校园春色视频在线观看| 亚洲国产精品久久男人天堂| 欧美成人免费av一区二区三区| 色视频www国产| 老司机午夜十八禁免费视频| 亚洲成人久久爱视频| 热99re8久久精品国产| 可以在线观看的亚洲视频| 久99久视频精品免费| 好男人在线观看高清免费视频| 桃红色精品国产亚洲av| 99在线视频只有这里精品首页| 成人av在线播放网站| 成人性生交大片免费视频hd| 久9热在线精品视频| 日本一二三区视频观看| 给我免费播放毛片高清在线观看| 亚洲av免费高清在线观看| 久久午夜亚洲精品久久| 久久性视频一级片| 亚洲最大成人手机在线| 国产亚洲精品久久久com| 黄片大片在线免费观看| tocl精华| 男人和女人高潮做爰伦理| 久久久色成人| 啪啪无遮挡十八禁网站| av视频在线观看入口| 欧美黑人巨大hd| 岛国在线免费视频观看| 亚洲欧美一区二区三区黑人| 观看美女的网站| 香蕉丝袜av| 国产蜜桃级精品一区二区三区| 热99re8久久精品国产| 免费在线观看亚洲国产| 天美传媒精品一区二区| 国产aⅴ精品一区二区三区波| 国产欧美日韩一区二区精品| 午夜老司机福利剧场| 久久久久久久精品吃奶| 99国产精品一区二区蜜桃av| 久久6这里有精品| www.www免费av| 成年人黄色毛片网站| 真人一进一出gif抽搐免费| 别揉我奶头~嗯~啊~动态视频| 日韩有码中文字幕| 非洲黑人性xxxx精品又粗又长| 老司机午夜十八禁免费视频| x7x7x7水蜜桃| 99久国产av精品| 亚洲在线自拍视频| 午夜日韩欧美国产| 欧美绝顶高潮抽搐喷水| 久久99热这里只有精品18| 3wmmmm亚洲av在线观看| 香蕉av资源在线| 成人特级av手机在线观看| 嫩草影院精品99| 午夜福利免费观看在线| 琪琪午夜伦伦电影理论片6080| 深爱激情五月婷婷| 乱人视频在线观看| 亚洲欧美一区二区三区黑人| 亚洲aⅴ乱码一区二区在线播放| 免费搜索国产男女视频| 好男人在线观看高清免费视频| 亚洲成人久久爱视频| 女人被狂操c到高潮| 特大巨黑吊av在线直播| 日本三级黄在线观看| 精品一区二区三区av网在线观看| 欧美最黄视频在线播放免费| 内射极品少妇av片p| 午夜视频国产福利| 波多野结衣巨乳人妻| 香蕉久久夜色| 白带黄色成豆腐渣| 亚洲在线自拍视频| 免费人成在线观看视频色| 国产亚洲欧美98| 亚洲第一欧美日韩一区二区三区| 成人18禁在线播放| 午夜精品一区二区三区免费看| 成人欧美大片| 香蕉久久夜色| 免费无遮挡裸体视频| 亚洲色图av天堂| 欧美不卡视频在线免费观看| 手机成人av网站| 欧美日韩黄片免| 波多野结衣高清无吗| 搡老岳熟女国产| 久久精品国产清高在天天线| 亚洲av成人av| 久久久久久久午夜电影| 免费大片18禁| 日本 欧美在线| 欧美日韩福利视频一区二区| 国产视频内射| 婷婷丁香在线五月| 国产高清videossex| 一级毛片女人18水好多| 欧美日韩国产亚洲二区| 欧美精品啪啪一区二区三区| 婷婷精品国产亚洲av| 亚洲精品在线美女| 此物有八面人人有两片| 亚洲国产精品999在线| 国产伦精品一区二区三区视频9 | 叶爱在线成人免费视频播放| 69av精品久久久久久| 男人和女人高潮做爰伦理| 国产麻豆成人av免费视频| 亚洲不卡免费看| 国产又黄又爽又无遮挡在线| 中文字幕高清在线视频| 亚洲国产精品合色在线| 精品国产亚洲在线| 亚洲精品亚洲一区二区| 亚洲成人中文字幕在线播放| 久久久成人免费电影| 五月玫瑰六月丁香| 亚洲性夜色夜夜综合| 18美女黄网站色大片免费观看| 国产淫片久久久久久久久 | 九九热线精品视视频播放| 一夜夜www| 免费看十八禁软件| 真人一进一出gif抽搐免费| 又爽又黄无遮挡网站| 蜜桃久久精品国产亚洲av| 99精品在免费线老司机午夜| 色综合婷婷激情| 黄色成人免费大全| 亚洲真实伦在线观看| 亚洲片人在线观看| 欧美一区二区国产精品久久精品| 亚洲av熟女| 麻豆成人午夜福利视频| 国产精品 欧美亚洲| 中文字幕高清在线视频| h日本视频在线播放| 成人国产一区最新在线观看| 九九久久精品国产亚洲av麻豆| a级毛片a级免费在线| 亚洲五月婷婷丁香| 国产69精品久久久久777片| 99久久成人亚洲精品观看| 国产熟女xx| 日本黄大片高清| 五月伊人婷婷丁香| 国产欧美日韩精品亚洲av| 免费观看人在逋| 午夜老司机福利剧场| 在线免费观看不下载黄p国产 | 色噜噜av男人的天堂激情| 变态另类丝袜制服| 夜夜爽天天搞| 国产午夜精品论理片| 成人特级黄色片久久久久久久| 美女免费视频网站| 90打野战视频偷拍视频| 国模一区二区三区四区视频| 色综合欧美亚洲国产小说| 波多野结衣巨乳人妻| 欧美+日韩+精品| 怎么达到女性高潮| 国产精品99久久久久久久久| 成人国产一区最新在线观看| 亚洲熟妇熟女久久| 欧美成人一区二区免费高清观看| 欧美黄色淫秽网站| 日本在线视频免费播放| 久久久久久久午夜电影| 俄罗斯特黄特色一大片| 欧美乱妇无乱码| 法律面前人人平等表现在哪些方面| 成人永久免费在线观看视频| 欧美一区二区精品小视频在线| 99国产极品粉嫩在线观看| 男女做爰动态图高潮gif福利片| 淫秽高清视频在线观看| 日韩欧美精品免费久久 | 国产高清有码在线观看视频| 国产亚洲av嫩草精品影院| 校园春色视频在线观看| 国产精品99久久久久久久久| 国产高潮美女av| 国产精品久久久人人做人人爽| 我要搜黄色片| 欧美色视频一区免费| 久久久久精品国产欧美久久久| 国产亚洲精品一区二区www| 成人av一区二区三区在线看| 少妇丰满av| a级毛片a级免费在线| 男插女下体视频免费在线播放| 亚洲人成伊人成综合网2020| 人妻丰满熟妇av一区二区三区| 老司机在亚洲福利影院| 国产精品久久久久久亚洲av鲁大| 国产主播在线观看一区二区| 美女黄网站色视频| 一个人免费在线观看电影| 国产精品美女特级片免费视频播放器| 高清在线国产一区| 亚洲一区二区三区不卡视频| av在线蜜桃| 久久这里只有精品中国| 欧美黑人欧美精品刺激| 久久天躁狠狠躁夜夜2o2o| 国产视频内射| 波多野结衣巨乳人妻| 三级国产精品欧美在线观看| 国产精品一及| 国产高清有码在线观看视频| 亚洲成a人片在线一区二区| 人妻丰满熟妇av一区二区三区| 亚洲成人免费电影在线观看| 男女做爰动态图高潮gif福利片| 亚洲av熟女| e午夜精品久久久久久久| 美女被艹到高潮喷水动态| 小说图片视频综合网站| 久久精品国产自在天天线| 俄罗斯特黄特色一大片| 在线观看日韩欧美| 人妻丰满熟妇av一区二区三区| 麻豆成人av在线观看| 五月伊人婷婷丁香| 偷拍熟女少妇极品色| 成人性生交大片免费视频hd| 观看免费一级毛片| 亚洲av二区三区四区| 啦啦啦免费观看视频1| 国产黄a三级三级三级人| 午夜精品在线福利| 蜜桃久久精品国产亚洲av| 欧美乱色亚洲激情| 亚洲真实伦在线观看| 有码 亚洲区| 国产真实乱freesex| 怎么达到女性高潮| 日韩免费av在线播放| АⅤ资源中文在线天堂| a在线观看视频网站| 午夜福利成人在线免费观看| 18+在线观看网站| 一进一出抽搐gif免费好疼| 欧美黄色片欧美黄色片| 欧美xxxx黑人xx丫x性爽| 精品久久久久久成人av| 国产精品99久久久久久久久| 日韩欧美国产一区二区入口| 琪琪午夜伦伦电影理论片6080| 精品一区二区三区视频在线 | 2021天堂中文幕一二区在线观| 一进一出好大好爽视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一欧美日韩一区二区三区| 在线观看舔阴道视频| 最好的美女福利视频网| 少妇的逼好多水| av女优亚洲男人天堂| 美女被艹到高潮喷水动态| 久久久久久人人人人人| 国产国拍精品亚洲av在线观看 | 久久伊人香网站| 国产精品99久久久久久久久| 久久国产乱子伦精品免费另类| 天天添夜夜摸| 无人区码免费观看不卡| 精品久久久久久久久久免费视频| 少妇高潮的动态图| 91字幕亚洲| 亚洲国产欧美人成| 午夜福利视频1000在线观看| 国产精品美女特级片免费视频播放器| 成人欧美大片| 日韩欧美国产一区二区入口| 露出奶头的视频| 最近在线观看免费完整版| 国产精品一区二区免费欧美| 亚洲中文字幕一区二区三区有码在线看| 亚洲av美国av| 亚洲成av人片免费观看| 亚洲天堂国产精品一区在线| 黄色日韩在线| 国产欧美日韩一区二区三| 国产成人影院久久av| 亚洲国产中文字幕在线视频| av福利片在线观看| 免费av毛片视频| 欧美成狂野欧美在线观看| 欧美最黄视频在线播放免费| 久久久久久九九精品二区国产| 久久久久久大精品| 日韩有码中文字幕| 非洲黑人性xxxx精品又粗又长| 国产成人av激情在线播放| 国内精品美女久久久久久| 中文字幕人妻丝袜一区二区| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频| 午夜福利高清视频| 丰满乱子伦码专区| 精品无人区乱码1区二区| 99久久99久久久精品蜜桃| 日韩欧美精品v在线| 国产精品日韩av在线免费观看| 美女 人体艺术 gogo| 香蕉久久夜色| 亚洲片人在线观看| 最近最新免费中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 国产视频一区二区在线看| 五月玫瑰六月丁香| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 亚洲成人中文字幕在线播放| 色精品久久人妻99蜜桃| 国产成人a区在线观看| 精品人妻1区二区| 国产高清有码在线观看视频| 国产av一区在线观看免费| 成人特级黄色片久久久久久久| 国产欧美日韩精品亚洲av| 51国产日韩欧美| 欧美乱色亚洲激情| 男女下面进入的视频免费午夜| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久久久免费视频| 日本免费一区二区三区高清不卡| 丰满乱子伦码专区| 91麻豆av在线| 又紧又爽又黄一区二区| 九色成人免费人妻av| 亚洲欧美激情综合另类| 内射极品少妇av片p| 久久这里只有精品中国| 中国美女看黄片| 观看免费一级毛片| 亚洲av成人不卡在线观看播放网| 欧美性猛交╳xxx乱大交人| 99久久综合精品五月天人人| 成人特级黄色片久久久久久久| 久久久国产精品麻豆| 成年免费大片在线观看| 亚洲欧美日韩东京热| 久久亚洲精品不卡| 毛片女人毛片| 亚洲成av人片在线播放无| eeuss影院久久| 日本一二三区视频观看| 国产精品野战在线观看| 无遮挡黄片免费观看| 母亲3免费完整高清在线观看| 亚洲国产精品久久男人天堂| 麻豆久久精品国产亚洲av| 亚洲av第一区精品v没综合| 精品一区二区三区av网在线观看| 亚洲精品在线观看二区| 亚洲欧美日韩卡通动漫| 欧美一级a爱片免费观看看| 欧美一区二区精品小视频在线| 成人午夜高清在线视频| 淫妇啪啪啪对白视频| 国产一区在线观看成人免费| 在线观看午夜福利视频| 亚洲在线自拍视频| 午夜福利18| 欧美成人性av电影在线观看| 天堂影院成人在线观看| 91字幕亚洲| 最后的刺客免费高清国语| 一级毛片高清免费大全| 亚洲美女视频黄频| tocl精华| 欧美+日韩+精品| 亚洲av美国av| 两人在一起打扑克的视频| 久久国产精品影院| 深爱激情五月婷婷| 中文在线观看免费www的网站| 亚洲专区中文字幕在线| 亚洲欧美一区二区三区黑人| 51国产日韩欧美| 国产精品1区2区在线观看.| 在线免费观看的www视频| 97超级碰碰碰精品色视频在线观看| 老司机在亚洲福利影院| 女生性感内裤真人,穿戴方法视频| 中出人妻视频一区二区| 精品一区二区三区视频在线观看免费| 90打野战视频偷拍视频| 色尼玛亚洲综合影院| 日韩人妻高清精品专区| 日日干狠狠操夜夜爽| 免费观看精品视频网站| 美女大奶头视频| 毛片女人毛片| 亚洲中文日韩欧美视频| 日韩欧美一区二区三区在线观看| 亚洲欧美日韩高清专用| 国产精品国产高清国产av| 久久婷婷人人爽人人干人人爱| 欧美色视频一区免费| 一本综合久久免费| 亚洲美女视频黄频| 美女 人体艺术 gogo| 日本黄色视频三级网站网址| 舔av片在线| 欧美大码av| 午夜福利欧美成人| 亚洲成人久久爱视频| 亚洲国产精品合色在线| 日本与韩国留学比较| 中亚洲国语对白在线视频| 午夜免费男女啪啪视频观看 | 久久人人精品亚洲av| 亚洲精品久久国产高清桃花| 少妇的逼水好多| 午夜免费观看网址| 欧美3d第一页| 亚洲成人免费电影在线观看| 丰满乱子伦码专区| 亚洲精品乱码久久久v下载方式 | 亚洲精品久久国产高清桃花| 国产亚洲精品av在线| 精品乱码久久久久久99久播| 国产97色在线日韩免费| 亚洲 欧美 日韩 在线 免费| 免费看美女性在线毛片视频| 别揉我奶头~嗯~啊~动态视频| 最近在线观看免费完整版| 岛国在线观看网站| 亚洲国产欧美人成| 国产 一区 欧美 日韩| 国产成人系列免费观看| 久9热在线精品视频| 国语自产精品视频在线第100页| 亚洲国产欧美人成| 国产色爽女视频免费观看| 日韩欧美精品v在线| 三级毛片av免费| 三级男女做爰猛烈吃奶摸视频| or卡值多少钱| 丰满人妻熟妇乱又伦精品不卡| 午夜影院日韩av| 成人av在线播放网站| 久久香蕉精品热| 亚洲aⅴ乱码一区二区在线播放| 真人做人爱边吃奶动态| 色吧在线观看| 欧美3d第一页| 在线观看av片永久免费下载| 亚洲黑人精品在线| 老司机深夜福利视频在线观看| 亚洲国产欧美人成| 脱女人内裤的视频| 丰满人妻一区二区三区视频av | 天堂av国产一区二区熟女人妻| 深夜精品福利| 亚洲国产中文字幕在线视频| 亚洲精品日韩av片在线观看 | 一个人看视频在线观看www免费 | 午夜免费男女啪啪视频观看 | a在线观看视频网站| 亚洲色图av天堂| 69人妻影院| 欧美日韩国产亚洲二区| 噜噜噜噜噜久久久久久91| a级一级毛片免费在线观看| 18禁国产床啪视频网站| 亚洲avbb在线观看| 欧美在线黄色| 婷婷丁香在线五月| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 亚洲乱码一区二区免费版| 99热6这里只有精品| 一进一出抽搐gif免费好疼| 亚洲最大成人手机在线| 在线a可以看的网站| 亚洲精品影视一区二区三区av| 成人高潮视频无遮挡免费网站| 99国产综合亚洲精品| 欧美乱码精品一区二区三区| 欧美性猛交黑人性爽| 一级毛片女人18水好多| 九九久久精品国产亚洲av麻豆| 欧洲精品卡2卡3卡4卡5卡区| av国产免费在线观看| 夜夜夜夜夜久久久久| 此物有八面人人有两片| 亚洲欧美日韩高清专用| 国产99白浆流出| 久久久色成人| 精品国产亚洲在线| 免费观看的影片在线观看| 国产毛片a区久久久久| 天天添夜夜摸| 精品久久久久久久久久久久久| 老汉色∧v一级毛片| 男女之事视频高清在线观看| 国产探花在线观看一区二区| www日本黄色视频网| 亚洲成av人片免费观看| 99久久精品一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲国产欧洲综合997久久,| 亚洲人成伊人成综合网2020| 精品不卡国产一区二区三区| 91在线精品国自产拍蜜月 | 国产一区在线观看成人免费| 婷婷丁香在线五月| 亚洲精品成人久久久久久| 色综合欧美亚洲国产小说| eeuss影院久久| 亚洲国产中文字幕在线视频| 男女做爰动态图高潮gif福利片| 国产av一区在线观看免费| 他把我摸到了高潮在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲中文日韩欧美视频| 亚洲av日韩精品久久久久久密| 一级毛片高清免费大全| 亚洲狠狠婷婷综合久久图片| 一级毛片高清免费大全| 十八禁网站免费在线| 国产av麻豆久久久久久久| 色综合亚洲欧美另类图片| 精品日产1卡2卡| 国产精品久久电影中文字幕| 蜜桃久久精品国产亚洲av| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 一个人看视频在线观看www免费 | 国产精品亚洲av一区麻豆| 亚洲精品成人久久久久久| 色综合欧美亚洲国产小说| 国产高清视频在线播放一区| 国产精品嫩草影院av在线观看 | 久久久久久九九精品二区国产| 中国美女看黄片| 亚洲人成网站在线播放欧美日韩| 有码 亚洲区| 小蜜桃在线观看免费完整版高清| 18禁国产床啪视频网站| 两个人看的免费小视频| 一区福利在线观看| 精品久久久久久成人av| 久久国产乱子伦精品免费另类| 两人在一起打扑克的视频| 亚洲国产欧美网| 一本精品99久久精品77| 欧美在线黄色| 在线播放国产精品三级| 他把我摸到了高潮在线观看| 在线观看日韩欧美| av天堂在线播放| 精品无人区乱码1区二区| 欧美日韩综合久久久久久 | 桃红色精品国产亚洲av| 国产亚洲精品久久久久久毛片| 亚洲aⅴ乱码一区二区在线播放| 国产精品99久久99久久久不卡| 亚洲av一区综合| 亚洲av五月六月丁香网| 狠狠狠狠99中文字幕| 国产黄片美女视频| www.www免费av| 亚洲中文字幕日韩| 国产真实乱freesex| 日本黄色视频三级网站网址| 久9热在线精品视频| 伊人久久大香线蕉亚洲五| 亚洲五月婷婷丁香| 99久久久亚洲精品蜜臀av| 国产淫片久久久久久久久 | 亚洲人成网站在线播| ponron亚洲| 国产亚洲欧美98| 亚洲精品一区av在线观看| 99热这里只有精品一区| 日韩欧美国产一区二区入口|