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

    基于蟻群-粒子群的模糊神經(jīng)網(wǎng)絡(luò)超聲波電動(dòng)機(jī)控制

    2018-04-26 12:45:15喬維德張本法
    微特電機(jī) 2018年4期
    關(guān)鍵詞:螞蟻粒子神經(jīng)網(wǎng)絡(luò)

    喬維德,張本法

    (1.無錫開放大學(xué),無錫 214011;2.常州開放大學(xué),常州 213001)

    0 引 言

    超聲波電動(dòng)機(jī)(以下簡(jiǎn)稱USM)是20世紀(jì)末發(fā)展起來的一種微小型驅(qū)動(dòng)電機(jī)。與傳統(tǒng)的電磁電機(jī)相比,USM結(jié)構(gòu)緊湊、體積小、質(zhì)量小,微位移、低轉(zhuǎn)速、大轉(zhuǎn)矩,運(yùn)行無噪聲,不受電磁干擾,耐低溫這些特性,是一種較為理想的伺服驅(qū)動(dòng)執(zhí)行器件,已廣泛應(yīng)用于航空航天、醫(yī)療器械、精密微動(dòng)機(jī)構(gòu)、辦公自動(dòng)化、機(jī)器人工業(yè)、高檔汽車、軍事及工業(yè)控制等領(lǐng)域。因USM自身參數(shù)的高度非線性、時(shí)變性,特別是電機(jī)工作中伴隨嚴(yán)重的摩擦過程,運(yùn)行特性呈現(xiàn)很多非線性和不確定因素,而且USM受溫度和外加負(fù)載的影響較大,所以難以建立精確的數(shù)學(xué)模型,很難通過傳統(tǒng)控制方法實(shí)現(xiàn)高精度的速度伺服控制[1]。目前很多學(xué)者和科技人員為解決USM的驅(qū)動(dòng)問題,做了許多探索和試驗(yàn),如應(yīng)用PID控制、模糊控制、神經(jīng)網(wǎng)絡(luò)等來優(yōu)化控制USM,然而采用以上方法時(shí),仍然不能有效解決USM轉(zhuǎn)速波動(dòng)較大及運(yùn)行不穩(wěn)定的問題,從而不能保證USM精確的轉(zhuǎn)速跟蹤和優(yōu)良的動(dòng)態(tài)控制性能。綜上所述,本文研究將蟻群算法和粒子群算法組成遞階結(jié)構(gòu)優(yōu)化模糊神經(jīng)網(wǎng)絡(luò)的USM轉(zhuǎn)速控制方法,即應(yīng)用蟻群算法的全局搜索和粒子群算法的局部搜索功能,優(yōu)化模糊神經(jīng)網(wǎng)絡(luò)控制器的結(jié)構(gòu)參數(shù),并且將該控制器引入U(xiǎn)SM的轉(zhuǎn)速控制系統(tǒng)。仿真分析及試驗(yàn)結(jié)果驗(yàn)證了本控制系統(tǒng)具有較強(qiáng)的轉(zhuǎn)速跟蹤與響應(yīng)能力,并且控制精度高,抗干擾能力強(qiáng)。

    1 USM轉(zhuǎn)速控制系統(tǒng)

    圖1是USM轉(zhuǎn)速控制系統(tǒng)原理示意圖。由速度參考模型設(shè)定電機(jī)的轉(zhuǎn)速跟蹤目標(biāo),可以是恒定轉(zhuǎn)速、階躍轉(zhuǎn)速等。在系統(tǒng)轉(zhuǎn)速控制內(nèi)環(huán)中引入電機(jī)定子的振動(dòng)信號(hào)檢測(cè)單元,以增強(qiáng)USM的速度跟蹤準(zhǔn)確度,提高USM的速度響應(yīng)快速性。控制原理結(jié)構(gòu)中的模糊神經(jīng)網(wǎng)絡(luò)(以下簡(jiǎn)稱FNN)控制器的功能是控制USM的轉(zhuǎn)速,其輸入量分別為轉(zhuǎn)速誤差e和誤差e的變化率Δe,e由系統(tǒng)結(jié)構(gòu)中參考模型的給定轉(zhuǎn)速信號(hào)n*與USM轉(zhuǎn)速檢測(cè)部件實(shí)際轉(zhuǎn)速量n比較計(jì)算得出,即e=n*-n。模糊神經(jīng)網(wǎng)絡(luò)控制器結(jié)構(gòu)參數(shù)通過蟻群算法在線優(yōu)化求取,以保證USM控制系統(tǒng)的實(shí)際輸出轉(zhuǎn)速n能及時(shí)準(zhǔn)確跟蹤參考模型的給定轉(zhuǎn)速輸出n*,即e→0。

    圖1 超聲波電動(dòng)機(jī)轉(zhuǎn)速控制原理圖

    2 FNN控制器

    USM參考模型自適應(yīng)轉(zhuǎn)速控制系統(tǒng)中FNN控制器結(jié)構(gòu)如圖2所示。控制器由輸入層(i層)、隸屬函數(shù)層(j層)、規(guī)則層(k層)、輸出層(o層)共4層前向神經(jīng)網(wǎng)絡(luò)組成[2]。

    圖2 FNN控制器結(jié)構(gòu)模型

    在FNN輸入層中,轉(zhuǎn)速誤差e及其變化率Δe由輸入層i節(jié)點(diǎn)傳至第2層,可以表示:

    (1)

    (2)

    FNN的隸屬函數(shù)層(模糊化層)主要負(fù)責(zé)對(duì)輸入變量模糊化,每個(gè)變量各自以模糊子集{N3(負(fù)大)、N2(負(fù)中)、N1(負(fù)小)、Z0(零)、P1(正小)、P2(正中)、P3(正大)}來表征。模糊子集中的隸屬度函數(shù)為該層的輸出量,采用高斯基函數(shù)形式表示隸屬度函數(shù)。該層第j個(gè)節(jié)點(diǎn)的輸出與輸入可表示:

    (3)

    在FNN的規(guī)則層中,每個(gè)節(jié)點(diǎn)對(duì)應(yīng)相應(yīng)的模糊邏輯運(yùn)算規(guī)則。即:

    (4)

    FNN的最后一層是輸出層,它只有1個(gè)節(jié)點(diǎn),其輸出量為規(guī)則層輸出信號(hào)經(jīng)處理后的總和。

    (5)

    3 FNN控制器結(jié)構(gòu)參數(shù)的優(yōu)化

    對(duì)上述FNN控制器的結(jié)構(gòu)參數(shù)優(yōu)化過程,就是通過學(xué)習(xí)訓(xùn)練不斷調(diào)整和更新高斯隸屬度函數(shù)參數(shù)mij和δij,連接權(quán)值ωjk和ωko過程,所以可調(diào)參數(shù)mij,δij,ωjk,ωko對(duì)模糊神經(jīng)網(wǎng)絡(luò)控制器性能產(chǎn)生重要影響。由于本控制系統(tǒng)中e及其變化率Δe2個(gè)輸入變量共有14個(gè)模糊子集,因此需要訓(xùn)練的mij和δij各14個(gè);規(guī)則層神經(jīng)元個(gè)數(shù)為Nj=7×7=49,故模糊神經(jīng)網(wǎng)絡(luò)需要訓(xùn)練的參數(shù)共49+2×14+49=126個(gè)。以往采取傳統(tǒng)的BP梯度下降法,優(yōu)化FNN的結(jié)構(gòu)參數(shù),收斂速度慢,比較容易陷入局部極小值,而且優(yōu)化效率較低。為此本文提出將全局搜索的蟻群算法作為主級(jí),將局部搜索的粒子群算法作為從級(jí),組成主從遞階結(jié)構(gòu)對(duì)FNN參數(shù)進(jìn)行訓(xùn)練優(yōu)化,可以提高FNN的泛化能力和快速全局收斂等。

    3.1 蟻群算法

    蟻群優(yōu)化(以下簡(jiǎn)稱ACO)算法屬于一種隨機(jī)搜索的智能優(yōu)化算法,它是以模仿螞蟻群體間的傳遞信息行為來實(shí)現(xiàn)類似搜索食物源的尋優(yōu)算法。螞蟻在行走路程中,一般會(huì)留下一種特殊的化學(xué)物質(zhì)——信息素,所有螞蟻都是根據(jù)此信息素以及信息素濃度而決定自身的行走路線,并始終向著高濃度信息素的路徑方向運(yùn)動(dòng)。越多的螞蟻通過某條路徑,后面的螞蟻選取這條路徑的機(jī)率就越大。螞蟻個(gè)體間正是以此信息交流方式來選擇最優(yōu)路徑,并實(shí)現(xiàn)最快捷搜索食物源的目的。設(shè)待學(xué)習(xí)優(yōu)化的模糊神經(jīng)網(wǎng)絡(luò)有m個(gè)參數(shù),包括所有的網(wǎng)絡(luò)權(quán)值和閾值(如mij,δij,ωjk,ωko共126個(gè),即m為126),設(shè)定P1,P2,…,Pm,將Pi(1≤i≤m)設(shè)定成N個(gè)非零隨機(jī)數(shù),并形成集合IPi。假設(shè)螞蟻的數(shù)量是S,所有螞蟻均由集合IPi出發(fā)搜索食物,按照集合IPi每個(gè)元素的信息素狀態(tài),從集合IPi隨機(jī)選擇唯一元素,并且對(duì)其信息素進(jìn)行及時(shí)更新[3]。第K(K=1,2,…,S)只螞蟻依此在每個(gè)集合IPi中選擇第j個(gè)元素,直至螞蟻全部到達(dá)食物源。路徑選擇概率規(guī)律為:

    (6)

    通過式(7)和式(8)更新集合IPi元素的信息素:

    τj(IPi)(t+1)=ρτj(IPi)(t)+Δτj(IPi)

    (7)

    (8)

    (9)

    式中:Q是用來控制信息素的調(diào)節(jié)速度,一般設(shè)為常數(shù);eK是蟻群中第K只螞蟻在本次循環(huán)所經(jīng)歷的路徑總長度。為克服ACO算法早熟和局部收斂問題,這里采取自適應(yīng)更新信息素方法。將信息素量值限制在[τmin,τmax],將式(7)變換如下:

    τj(IPi)(t+1)=

    ρ1+φ(m)τj(IPi)(t)+Δτj(IPi)τ>τmax

    (10)

    τj(IPi)(t+1)=

    ρ1-φ(m)τj(IPi)(t)+Δτj(IPi)τ<τmin

    (11)

    式中:φ(m)表示與迭代次數(shù)成正比的函數(shù)。

    3.2 粒子群算法

    粒子群優(yōu)化(以下簡(jiǎn)稱PSO)算法將需要優(yōu)化和解決的問題當(dāng)成某一搜索空間d范圍內(nèi)的鳥,簡(jiǎn)稱“粒子”,粒子所處位置即為待求的潛在解。標(biāo)準(zhǔn)PSO算法根據(jù)式(12)和式(13)進(jìn)行迭代,直至搜尋粒子群中達(dá)到最小誤差要求的最優(yōu)位置[4]。

    Vij(t+1)=ωVij(t)+C1R1[Pj(t)-Xij(t)]+

    C2R2[Gj(t)-Xij(t)]

    (12)

    (13)

    式(12)、式(13)中,j=1,2,…,d;t為迭代次數(shù);Xij(t),Vij(t)分別為第t代粒子位置和速度;C1,C2表示學(xué)習(xí)因子;R1,R2取值范圍為[0,1]。為克服PSO算法的早熟與振蕩現(xiàn)象,讓慣性權(quán)重ω根據(jù)粒子適配值作自動(dòng)調(diào)節(jié),即:

    (14)

    式中:f,fave,fmax分別代表粒子的適應(yīng)度值、平均適應(yīng)度值以及最大適應(yīng)度值。

    3.3 蟻群-粒子群算法優(yōu)化網(wǎng)絡(luò)參數(shù)過程

    ACO比較容易產(chǎn)生停滯現(xiàn)象,搜索時(shí)間較長。PSO又具有較快的尋優(yōu)求解速度,因此將ACO和PSO組成ACO-PSO混合算法,以優(yōu)化模糊神經(jīng)網(wǎng)絡(luò)參數(shù)。ACO將網(wǎng)絡(luò)需要訓(xùn)練的參數(shù)在取值范圍內(nèi)求得有限個(gè)離散值[5]。PSO搜索尋優(yōu)實(shí)際上是粒子速度和位置的改變過程,所以FNN學(xué)習(xí)算法中的高斯隸屬函數(shù)均值、標(biāo)準(zhǔn)差和權(quán)重應(yīng)該與粒子的位置相對(duì)應(yīng)。

    第1步 初始化。蟻群初始化主要設(shè)定時(shí)間、迭代次數(shù)NC、最大迭代次數(shù)NCmax、螞蟻數(shù)目S、初始信息素ρ等,且讓全部螞蟻放置在蟻巢;粒子群初始化包括設(shè)置粒子數(shù)h、粒子初始位置、速度、慣性因子ωmax和ωmin、加速因子C1和C2、最大迭代次數(shù)Mp等。

    第2步 起動(dòng)全部螞蟻,第K(K=1,2,…,S)只螞蟻從集合IPi開始,按次序從IPi中選取元素,直至蟻群中所有螞蟻均搜尋到食物源。

    第3步 將每只螞蟻選取值作為FNN結(jié)構(gòu)參數(shù)訓(xùn)練樣本,然后按照式(15)求取其網(wǎng)絡(luò)輸出平均相對(duì)誤差,記錄目前FNN結(jié)構(gòu)參數(shù)最好解[5],且根據(jù)式(8)~式(11)對(duì)全部集合中各元素信息素更新調(diào)整。式(9)中eK定義為FNN期望輸出與實(shí)際值的平均相對(duì)誤差。定義公式:

    (15)

    式中:yi表示FNN的實(shí)際輸出值;Oi代表FNN期望輸出值;S為螞蟻樣本數(shù)量。因此,當(dāng)網(wǎng)絡(luò)訓(xùn)練輸出誤差eK減小時(shí),相應(yīng)信息素濃度就增大。式(10)和式(11)中的φ(m)=m/15,其中m為連續(xù)收斂次數(shù)。

    第4步 采用PSO改進(jìn)第3步中求得的最優(yōu)解。首先按式(15)計(jì)算粒子適應(yīng)度,并按式(12)和式(13)及時(shí)更新粒子的運(yùn)行速度與所處位置。假如當(dāng)前找到的粒子最好位置Pbest不能優(yōu)于粒子當(dāng)前位置Pnow,那么令Pbest=Pnow;假如當(dāng)前搜出的全局最好位置Gbest不能優(yōu)于粒子當(dāng)前位置Pnow,那么令Gbest=Pnow。如果沒有達(dá)到粒子群最大迭代次數(shù),繼續(xù)對(duì)粒子當(dāng)前位置與最優(yōu)位置Pbest及Gbest進(jìn)行比較執(zhí)行。

    第5步 循環(huán)次數(shù)增1。若目前最優(yōu)解得到改進(jìn),更新信息素濃度[5]。

    第6步 如果進(jìn)化滿足最大迭代次數(shù),便停止迭代過程,同時(shí)得出最優(yōu)解,否則將跳至第2步循環(huán)執(zhí)行。

    4 仿真分析

    利用MATLAB/Simulink對(duì)USM控制系統(tǒng)進(jìn)行仿真研究。應(yīng)用ACO與PSO結(jié)合(即ACO-PSO混合算法)訓(xùn)練優(yōu)化模糊神經(jīng)網(wǎng)絡(luò)。各參數(shù)選?。害?0.7,Q=20,NCmax=300,S=100;h=30,Mp=55,ωmax=1.6,ωmin=0.3,C1=C2=2。為對(duì)比不同算法訓(xùn)練模糊神經(jīng)網(wǎng)絡(luò)性能優(yōu)劣,本文采用同樣的訓(xùn)練樣本,先后通過ACO、PSO及ACO-PSO 3種算法對(duì)FNN加以學(xué)習(xí),訓(xùn)練結(jié)果如表1所示。由表1分析得到,采用ACO-PSO混合算法訓(xùn)練FNN,訓(xùn)練時(shí)間較短,訓(xùn)練精度較高,無論是在尋優(yōu)時(shí)間效率上還是尋優(yōu)精度方面,都取得較好的優(yōu)化性能。

    表1 三種算法的模糊神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果

    圖3、圖4分別是設(shè)定目標(biāo)轉(zhuǎn)速90 r/min且?guī)?.2 N·m負(fù)載時(shí),本文設(shè)計(jì)的經(jīng)ACO-PSO優(yōu)化FNN方法和按BP算法優(yōu)化的傳統(tǒng)FNN方法USM系統(tǒng)速度仿真特性曲線。在穩(wěn)態(tài)時(shí)圖4的速度波動(dòng)較圖3大且明顯。與傳統(tǒng)的BP算法優(yōu)化FNN控制相比,經(jīng)ACO-PSO優(yōu)化FNN控制器的控制系統(tǒng)的動(dòng)態(tài)響應(yīng)與靜態(tài)特性明顯優(yōu)于傳統(tǒng)FNN控制。

    圖3 ACO-PSO優(yōu)化FNN

    圖4 傳統(tǒng)BP訓(xùn)練FNN

    為驗(yàn)證經(jīng)ACO-PSO算法優(yōu)化的FNN控制器的抗干擾能力,t=0.15 s時(shí)刻又突加0.4 N·m負(fù)載,并同傳統(tǒng)BP算法訓(xùn)練的FNN控制進(jìn)行比較分析。從圖5、圖6對(duì)比可知,當(dāng)系統(tǒng)突加負(fù)載時(shí),由ACO-PSO算法優(yōu)化的FNN系統(tǒng)能使USM速度迅速穩(wěn)定在90 r/min,轉(zhuǎn)速響應(yīng)快,轉(zhuǎn)速波動(dòng)小,對(duì)外界負(fù)載擾動(dòng)能快速自動(dòng)調(diào)節(jié),抗干擾能力強(qiáng),穩(wěn)態(tài)后系統(tǒng)控制精度較高,而傳統(tǒng)BP算法訓(xùn)練的FNN,當(dāng)突加2倍負(fù)載時(shí),出現(xiàn)較大的轉(zhuǎn)速超調(diào)及脈動(dòng)。所以,綜合USM起動(dòng)、速度跟蹤及負(fù)載擾動(dòng)情況而言,經(jīng)ACO-PSO算法優(yōu)化的FNN控制器具有更好的動(dòng)態(tài)特性和魯棒性。

    圖5 ACO-PSO優(yōu)化FNN突加負(fù)載時(shí)速度響應(yīng)

    圖6 傳統(tǒng)BP訓(xùn)練FNN突加負(fù)載時(shí)速度響應(yīng)

    5 實(shí)驗(yàn)驗(yàn)證

    基于DSP的USM轉(zhuǎn)速控制系統(tǒng)結(jié)構(gòu)原理框圖及實(shí)驗(yàn)平臺(tái)如圖7、圖8所示[6]。

    圖7 USM速度控制結(jié)構(gòu)原理圖

    圖8 USM速度控制實(shí)驗(yàn)平臺(tái)截圖

    控制系統(tǒng)由PC機(jī)、DSP(數(shù)字信號(hào)處理器)、USM驅(qū)動(dòng)器、光電編碼器、示波器等功能設(shè)備組成。USM為直徑60 mm行波型電機(jī),DSP選用TMS320LF2407芯片。光電編碼器選用5 000線Model775A?,F(xiàn)場(chǎng)實(shí)驗(yàn)測(cè)試時(shí),通過光電編碼器將USM的旋轉(zhuǎn)角位移轉(zhuǎn)換成連續(xù)的電脈沖信號(hào),再通過F/V(頻率/電壓)變換和A/D(數(shù)/模)轉(zhuǎn)換后傳給DSP,DSP處理后轉(zhuǎn)換為速度值,作為輸入變量供給DSP進(jìn)行實(shí)時(shí)計(jì)算,同時(shí)通過DSP的串行口傳送至PC機(jī)進(jìn)行顯示和保存。數(shù)據(jù)采樣頻率選取為1 000 Hz。通過USM速度參考輸入量及實(shí)際轉(zhuǎn)速反饋量的比較運(yùn)算,在壓控振蕩器(采用CD4046)的控制端輸出控制電壓,從而改變USM驅(qū)動(dòng)電路頻率,有效驅(qū)動(dòng)USM。

    利用示波器觀察設(shè)計(jì)的控制系統(tǒng),分別跟蹤恒定目標(biāo)轉(zhuǎn)速80r/min和階躍目標(biāo)轉(zhuǎn)速80r/min→50 r/min→80 r/min→50 r/min時(shí)的實(shí)驗(yàn)結(jié)果,如圖9、圖10所示(USM負(fù)載力矩為0.5 N·m)。

    圖9 USM恒轉(zhuǎn)速跟蹤曲線

    圖10 USM階躍轉(zhuǎn)速跟蹤曲線

    圖9中USM跟隨目標(biāo)轉(zhuǎn)速的波動(dòng)很小,僅在±1.8 r/mim左右,近似一條直線,說明對(duì)負(fù)載的適應(yīng)性較強(qiáng)。圖10中當(dāng)USM目標(biāo)轉(zhuǎn)速設(shè)為80 r/min時(shí),實(shí)測(cè)的轉(zhuǎn)速波動(dòng)范圍為±1.7 r/mim;當(dāng)目標(biāo)轉(zhuǎn)速階躍為50 r/min時(shí),實(shí)測(cè)運(yùn)行的轉(zhuǎn)速波動(dòng)范圍為±2.1 r/mim。

    5 結(jié) 語

    本文用ACO和PSO相結(jié)合形成的ACO-PSO混合算法,優(yōu)化FNN控制器的高斯隸屬函數(shù)參數(shù)和神經(jīng)網(wǎng)絡(luò)權(quán)值,并應(yīng)用于USM的自適應(yīng)速度跟蹤系統(tǒng),電機(jī)轉(zhuǎn)速跟蹤性能強(qiáng),轉(zhuǎn)速脈動(dòng)明顯減小,控制精度也得到明顯改善。該系統(tǒng)呈現(xiàn)出較強(qiáng)的適應(yīng)性和魯棒性,取得較為滿意的優(yōu)化控制成效,為USM智能化控制提供一種新的技術(shù)方案,有一定的工程參考與應(yīng)用價(jià)值。

    [1] 黃帆,金龍,徐志科,等.基于卡爾曼濾波的超聲波電動(dòng)機(jī)測(cè)速算法研究[J].微特電機(jī),2015,43(11):1-4.

    [2] 喬維德.遺傳算法優(yōu)化的開關(guān)磁阻電動(dòng)機(jī)RFNN位置控制器設(shè)計(jì)[J].微特電機(jī),2016,44(2):75-77.

    [3] 喬維德.基于蟻群模糊神經(jīng)網(wǎng)絡(luò)的無刷直流電機(jī)仿真控制研究[J].電氣傳動(dòng)自動(dòng)化,2009,31(5):6-10.

    [4] 喬維德.改進(jìn)粒子群優(yōu)化神經(jīng)網(wǎng)絡(luò)的變壓器故障診斷[J].電工電氣,2015(12):24-27.

    [5] 于萬霞,杜太行,鄭宏興.基于蟻群粒子群的模糊神經(jīng)網(wǎng)絡(luò)交通流量預(yù)測(cè)[J].計(jì)算機(jī)工程與應(yīng)用,2007,43(25):197-199.

    [6] 喬維德.一種改進(jìn)的提升機(jī)同步電機(jī)直接轉(zhuǎn)矩控制[J].鹽城工學(xué)院學(xué)報(bào),2017,30(1):28-33.

    猜你喜歡
    螞蟻粒子神經(jīng)網(wǎng)絡(luò)
    神經(jīng)網(wǎng)絡(luò)抑制無線通信干擾探究
    電子制作(2019年19期)2019-11-23 08:42:00
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    我們會(huì)“隱身”讓螞蟻來保護(hù)自己
    螞蟻
    基于神經(jīng)網(wǎng)絡(luò)的拉矯機(jī)控制模型建立
    復(fù)數(shù)神經(jīng)網(wǎng)絡(luò)在基于WiFi的室內(nèi)LBS應(yīng)用
    基于支持向量機(jī)回歸和RBF神經(jīng)網(wǎng)絡(luò)的PID整定
    螞蟻找吃的等
    基于Matlab的α粒子的散射實(shí)驗(yàn)?zāi)M
    物理與工程(2014年4期)2014-02-27 11:23:08
    亚洲精品国产av蜜桃| 在线观看一区二区三区激情| 欧美日韩视频精品一区| 99精国产麻豆久久婷婷| 亚洲欧洲日产国产| 国产精品av久久久久免费| 性色av乱码一区二区三区2| 99国产综合亚洲精品| 1024视频免费在线观看| 天天操日日干夜夜撸| 黄片大片在线免费观看| 一区二区三区精品91| 亚洲精品乱久久久久久| 丰满迷人的少妇在线观看| 窝窝影院91人妻| netflix在线观看网站| 欧美国产精品va在线观看不卡| 99热全是精品| 亚洲av成人不卡在线观看播放网 | 超碰成人久久| 午夜免费鲁丝| 操美女的视频在线观看| 亚洲精品日韩在线中文字幕| 大香蕉久久网| 一二三四社区在线视频社区8| 精品视频人人做人人爽| 热re99久久精品国产66热6| 久久精品国产综合久久久| 免费一级毛片在线播放高清视频 | 亚洲国产精品成人久久小说| 久久久精品区二区三区| 手机成人av网站| 夜夜骑夜夜射夜夜干| 欧美 亚洲 国产 日韩一| 狠狠婷婷综合久久久久久88av| 免费在线观看视频国产中文字幕亚洲 | 日韩,欧美,国产一区二区三区| 亚洲九九香蕉| 一级片免费观看大全| 亚洲五月婷婷丁香| 亚洲精品中文字幕一二三四区 | 久久99一区二区三区| 汤姆久久久久久久影院中文字幕| 桃花免费在线播放| 青青草视频在线视频观看| 亚洲欧美日韩高清在线视频 | 国产av国产精品国产| 亚洲国产欧美在线一区| 国产精品 欧美亚洲| 一级毛片女人18水好多| 男男h啪啪无遮挡| 色婷婷久久久亚洲欧美| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情极品国产一区二区三区| 蜜桃在线观看..| 一个人免费在线观看的高清视频 | 51午夜福利影视在线观看| 永久免费av网站大全| 中国美女看黄片| 欧美成狂野欧美在线观看| 交换朋友夫妻互换小说| 超色免费av| 十分钟在线观看高清视频www| 国产日韩一区二区三区精品不卡| 国产野战对白在线观看| 亚洲第一欧美日韩一区二区三区 | 啦啦啦中文免费视频观看日本| 50天的宝宝边吃奶边哭怎么回事| 久久人人爽av亚洲精品天堂| 在线观看舔阴道视频| 亚洲va日本ⅴa欧美va伊人久久 | 国产深夜福利视频在线观看| 激情视频va一区二区三区| 国产一级毛片在线| 波多野结衣av一区二区av| 在线观看人妻少妇| 啦啦啦免费观看视频1| 久久精品熟女亚洲av麻豆精品| 亚洲专区中文字幕在线| 日日夜夜操网爽| 男女下面插进去视频免费观看| 高清视频免费观看一区二区| 亚洲国产看品久久| 国产一区二区三区综合在线观看| 天天添夜夜摸| 日韩大码丰满熟妇| 极品人妻少妇av视频| 极品人妻少妇av视频| 正在播放国产对白刺激| 黑人猛操日本美女一级片| 一区在线观看完整版| 免费观看av网站的网址| 免费在线观看完整版高清| 亚洲久久久国产精品| kizo精华| 女性生殖器流出的白浆| 国产av一区二区精品久久| 欧美另类一区| 欧美少妇被猛烈插入视频| 欧美大码av| av又黄又爽大尺度在线免费看| 后天国语完整版免费观看| 欧美人与性动交α欧美精品济南到| 性少妇av在线| av网站在线播放免费| 亚洲情色 制服丝袜| 欧美另类亚洲清纯唯美| 18禁国产床啪视频网站| 欧美激情 高清一区二区三区| 一二三四社区在线视频社区8| 亚洲精品中文字幕一二三四区 | 国产一区二区三区综合在线观看| 国产亚洲欧美在线一区二区| 国产av精品麻豆| 亚洲国产精品一区二区三区在线| 精品福利永久在线观看| 欧美成狂野欧美在线观看| 国产麻豆69| 97人妻天天添夜夜摸| 欧美国产精品一级二级三级| 国产精品 欧美亚洲| 久久久久久久久免费视频了| 两个人看的免费小视频| 精品一区二区三区av网在线观看 | 在线精品无人区一区二区三| 80岁老熟妇乱子伦牲交| 超碰97精品在线观看| 成人国产一区最新在线观看| 精品视频人人做人人爽| 亚洲中文av在线| 老鸭窝网址在线观看| 久久久久久亚洲精品国产蜜桃av| 黄色怎么调成土黄色| 亚洲专区国产一区二区| 欧美亚洲日本最大视频资源| 三上悠亚av全集在线观看| 欧美日韩福利视频一区二区| 黄色片一级片一级黄色片| 黑人猛操日本美女一级片| 中文字幕色久视频| 午夜91福利影院| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品美女久久av网站| 久久久久久久久久久久大奶| 精品熟女少妇八av免费久了| 美女扒开内裤让男人捅视频| netflix在线观看网站| 欧美 日韩 精品 国产| 亚洲 国产 在线| 国产一区有黄有色的免费视频| 亚洲第一青青草原| 俄罗斯特黄特色一大片| 日韩欧美免费精品| 亚洲av日韩精品久久久久久密| 亚洲精品国产精品久久久不卡| 在线观看人妻少妇| av视频免费观看在线观看| 亚洲国产毛片av蜜桃av| 天天躁夜夜躁狠狠躁躁| 考比视频在线观看| 国产亚洲欧美在线一区二区| 成人手机av| 国产日韩一区二区三区精品不卡| 国产在视频线精品| 亚洲国产欧美日韩在线播放| 久久毛片免费看一区二区三区| 精品国产一区二区三区久久久樱花| 男女下面插进去视频免费观看| 热99re8久久精品国产| 亚洲欧洲精品一区二区精品久久久| 精品国产国语对白av| 国产精品一二三区在线看| 久久av网站| 波多野结衣一区麻豆| av电影中文网址| av欧美777| 成年人黄色毛片网站| 国产成人av激情在线播放| a 毛片基地| tube8黄色片| 精品国产一区二区三区四区第35| 婷婷成人精品国产| 91国产中文字幕| 最黄视频免费看| svipshipincom国产片| 丁香六月欧美| 国产精品一区二区精品视频观看| 欧美日韩亚洲高清精品| 精品视频人人做人人爽| 久久国产精品大桥未久av| 女性生殖器流出的白浆| 美女大奶头黄色视频| 日韩制服骚丝袜av| 亚洲中文日韩欧美视频| 欧美日韩精品网址| 嫩草影视91久久| 女性生殖器流出的白浆| 久久九九热精品免费| 久久人人爽av亚洲精品天堂| 精品少妇一区二区三区视频日本电影| 久久人人97超碰香蕉20202| 黄片小视频在线播放| 久久久久国内视频| 亚洲第一欧美日韩一区二区三区 | 亚洲天堂av无毛| 少妇的丰满在线观看| 2018国产大陆天天弄谢| 91九色精品人成在线观看| 亚洲av成人不卡在线观看播放网 | 亚洲一区中文字幕在线| 日韩视频在线欧美| 欧美日韩成人在线一区二区| 亚洲专区字幕在线| 啦啦啦 在线观看视频| 日韩大片免费观看网站| 国产高清视频在线播放一区 | 婷婷色av中文字幕| 国产成人精品无人区| 青青草视频在线视频观看| 在线观看人妻少妇| 操出白浆在线播放| 亚洲中文字幕日韩| 永久免费av网站大全| 一区二区三区四区激情视频| 欧美在线黄色| 久久香蕉激情| 美女国产高潮福利片在线看| 亚洲九九香蕉| 久久人人爽av亚洲精品天堂| 欧美黑人精品巨大| 两人在一起打扑克的视频| 久久久久国产精品人妻一区二区| 俄罗斯特黄特色一大片| 91精品国产国语对白视频| 午夜福利乱码中文字幕| 99热国产这里只有精品6| 岛国在线观看网站| 大香蕉久久网| 亚洲精品国产色婷婷电影| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 老司机亚洲免费影院| 欧美激情极品国产一区二区三区| 久久影院123| 亚洲专区中文字幕在线| 黄色视频不卡| 操出白浆在线播放| 亚洲欧美日韩高清在线视频 | av天堂在线播放| www.999成人在线观看| 国产av又大| 国产人伦9x9x在线观看| 久久精品亚洲熟妇少妇任你| 日本91视频免费播放| 最近最新中文字幕大全免费视频| 亚洲国产中文字幕在线视频| 久久人人爽av亚洲精品天堂| 黄色 视频免费看| 久久狼人影院| 中国美女看黄片| 久久热在线av| 好男人电影高清在线观看| 亚洲国产精品一区二区三区在线| 中文欧美无线码| 看免费av毛片| 国产无遮挡羞羞视频在线观看| 国产不卡av网站在线观看| 18禁黄网站禁片午夜丰满| 丝袜美足系列| 午夜精品国产一区二区电影| 亚洲中文av在线| 日韩制服骚丝袜av| 一级黄色大片毛片| 久久亚洲国产成人精品v| 69av精品久久久久久 | √禁漫天堂资源中文www| 欧美激情久久久久久爽电影 | 97人妻天天添夜夜摸| av在线老鸭窝| 叶爱在线成人免费视频播放| 热re99久久精品国产66热6| 日本欧美视频一区| 人成视频在线观看免费观看| 婷婷丁香在线五月| 黄色毛片三级朝国网站| 亚洲精品第二区| 精品少妇黑人巨大在线播放| 满18在线观看网站| 免费不卡黄色视频| 国产在视频线精品| 久久精品国产亚洲av香蕉五月 | 日本91视频免费播放| 成年美女黄网站色视频大全免费| 中文字幕色久视频| 不卡av一区二区三区| 男女无遮挡免费网站观看| 免费观看av网站的网址| 可以免费在线观看a视频的电影网站| 一本大道久久a久久精品| 精品亚洲乱码少妇综合久久| 国产在线视频一区二区| 女警被强在线播放| 精品卡一卡二卡四卡免费| 国产高清视频在线播放一区 | 色播在线永久视频| 极品少妇高潮喷水抽搐| 叶爱在线成人免费视频播放| 最新在线观看一区二区三区| 欧美国产精品va在线观看不卡| 亚洲精品国产av蜜桃| 国产欧美亚洲国产| 日韩有码中文字幕| 美女福利国产在线| 日韩欧美国产一区二区入口| 欧美日韩亚洲综合一区二区三区_| 亚洲avbb在线观看| 亚洲av成人一区二区三| 欧美+亚洲+日韩+国产| 成在线人永久免费视频| 大香蕉久久成人网| 色综合欧美亚洲国产小说| 一级a爱视频在线免费观看| 三上悠亚av全集在线观看| 国产精品九九99| 夜夜夜夜夜久久久久| 男女高潮啪啪啪动态图| 亚洲精品在线美女| 性高湖久久久久久久久免费观看| 人妻久久中文字幕网| 日韩一区二区三区影片| 男女床上黄色一级片免费看| 精品一区二区三区四区五区乱码| 无限看片的www在线观看| 男女无遮挡免费网站观看| 后天国语完整版免费观看| 亚洲精品久久午夜乱码| 黑人巨大精品欧美一区二区蜜桃| 男女之事视频高清在线观看| 法律面前人人平等表现在哪些方面 | 91成人精品电影| 亚洲欧美精品自产自拍| 亚洲精品国产色婷婷电影| 国产欧美日韩一区二区三区在线| 美女高潮喷水抽搐中文字幕| 一级,二级,三级黄色视频| 美女扒开内裤让男人捅视频| 美女视频免费永久观看网站| 午夜福利在线免费观看网站| www.精华液| 亚洲精品国产区一区二| 欧美黑人精品巨大| 国产伦人伦偷精品视频| 欧美少妇被猛烈插入视频| 亚洲精品在线美女| 男女之事视频高清在线观看| 日韩制服骚丝袜av| 午夜激情av网站| 亚洲av日韩精品久久久久久密| 一区二区av电影网| 热re99久久国产66热| 搡老岳熟女国产| cao死你这个sao货| 国产人伦9x9x在线观看| 欧美日韩亚洲高清精品| 亚洲av欧美aⅴ国产| 黑人猛操日本美女一级片| 国产亚洲精品一区二区www | 女人精品久久久久毛片| 黄色视频在线播放观看不卡| 亚洲专区中文字幕在线| 亚洲精品粉嫩美女一区| 亚洲国产av新网站| 久久精品亚洲熟妇少妇任你| 两个人免费观看高清视频| 大香蕉久久成人网| 中文字幕av电影在线播放| 99久久人妻综合| 国产黄频视频在线观看| 9191精品国产免费久久| 久久天堂一区二区三区四区| 国产av又大| 欧美日韩国产mv在线观看视频| 波多野结衣av一区二区av| 后天国语完整版免费观看| 欧美激情久久久久久爽电影 | 精品人妻在线不人妻| 一区二区三区精品91| 91大片在线观看| 自线自在国产av| 亚洲激情五月婷婷啪啪| 欧美精品av麻豆av| 少妇粗大呻吟视频| 国产精品自产拍在线观看55亚洲 | 中文字幕人妻丝袜制服| 成人18禁高潮啪啪吃奶动态图| 美女午夜性视频免费| 欧美日韩视频精品一区| 国产成人免费无遮挡视频| 日韩电影二区| 午夜福利在线观看吧| 国产在视频线精品| 天天影视国产精品| 69精品国产乱码久久久| 精品卡一卡二卡四卡免费| 另类亚洲欧美激情| 自拍欧美九色日韩亚洲蝌蚪91| 黑人巨大精品欧美一区二区蜜桃| 久久精品国产亚洲av高清一级| xxxhd国产人妻xxx| 欧美在线黄色| 男男h啪啪无遮挡| 亚洲av成人不卡在线观看播放网 | 亚洲欧美一区二区三区久久| 老熟女久久久| 国产精品一区二区在线观看99| 人人妻人人添人人爽欧美一区卜| 国产又色又爽无遮挡免| 首页视频小说图片口味搜索| 亚洲精品国产区一区二| 中亚洲国语对白在线视频| 夜夜骑夜夜射夜夜干| 天天躁狠狠躁夜夜躁狠狠躁| 香蕉国产在线看| 亚洲自偷自拍图片 自拍| 国产精品二区激情视频| 亚洲国产欧美在线一区| 纯流量卡能插随身wifi吗| 免费高清在线观看视频在线观看| 国产成+人综合+亚洲专区| 99re6热这里在线精品视频| 国产激情久久老熟女| 韩国精品一区二区三区| 一区二区三区精品91| 9热在线视频观看99| 大陆偷拍与自拍| 精品国产乱子伦一区二区三区 | 欧美国产精品一级二级三级| av天堂在线播放| 国产熟女午夜一区二区三区| 97人妻天天添夜夜摸| 久久人妻福利社区极品人妻图片| 超碰成人久久| 国产精品国产av在线观看| 黑人欧美特级aaaaaa片| 国产欧美日韩精品亚洲av| 国产1区2区3区精品| 国产av精品麻豆| 欧美性长视频在线观看| 麻豆av在线久日| 在线观看www视频免费| 人人妻人人澡人人看| videos熟女内射| 2018国产大陆天天弄谢| 亚洲,欧美精品.| 青青草视频在线视频观看| www.av在线官网国产| 午夜福利影视在线免费观看| av网站免费在线观看视频| 在线观看www视频免费| 操出白浆在线播放| 搡老岳熟女国产| 一级毛片电影观看| 精品少妇黑人巨大在线播放| 91九色精品人成在线观看| 满18在线观看网站| 999久久久精品免费观看国产| 亚洲欧美成人综合另类久久久| 在线观看人妻少妇| 欧美另类一区| 成年人免费黄色播放视频| 黑人操中国人逼视频| 天天操日日干夜夜撸| 国产视频一区二区在线看| 美女主播在线视频| 天天影视国产精品| 制服人妻中文乱码| 国产一区二区三区在线臀色熟女 | 伊人亚洲综合成人网| 一区二区三区激情视频| 国产精品秋霞免费鲁丝片| tube8黄色片| 欧美另类一区| 女性生殖器流出的白浆| 女性被躁到高潮视频| 日韩大码丰满熟妇| 亚洲国产成人一精品久久久| 少妇的丰满在线观看| 久久天堂一区二区三区四区| 日日摸夜夜添夜夜添小说| 欧美xxⅹ黑人| 亚洲中文字幕日韩| 亚洲国产av新网站| 人妻 亚洲 视频| 淫妇啪啪啪对白视频 | 一级毛片女人18水好多| 久久人妻福利社区极品人妻图片| 麻豆乱淫一区二区| 国产成人av教育| tube8黄色片| 中国国产av一级| 亚洲av日韩精品久久久久久密| 欧美国产精品va在线观看不卡| 久久久久国内视频| 午夜视频精品福利| 大片电影免费在线观看免费| 日本精品一区二区三区蜜桃| 又黄又粗又硬又大视频| 亚洲一区二区三区欧美精品| 久久久久久久大尺度免费视频| 亚洲第一欧美日韩一区二区三区 | 精品久久久久久久毛片微露脸 | 日韩欧美免费精品| 黑人巨大精品欧美一区二区mp4| 久久香蕉激情| 青春草视频在线免费观看| 成人18禁高潮啪啪吃奶动态图| 视频区欧美日本亚洲| 久久久久网色| 一区在线观看完整版| 免费少妇av软件| 两个人看的免费小视频| 丝袜美腿诱惑在线| 韩国高清视频一区二区三区| 少妇人妻久久综合中文| 宅男免费午夜| 咕卡用的链子| 嫁个100分男人电影在线观看| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 欧美日本中文国产一区发布| av不卡在线播放| 欧美成狂野欧美在线观看| 91成人精品电影| 国产精品一区二区在线不卡| 国产男女超爽视频在线观看| 久热爱精品视频在线9| 精品一区二区三区四区五区乱码| 两个人免费观看高清视频| 国产成人a∨麻豆精品| 亚洲av欧美aⅴ国产| 久久久久久人人人人人| 熟女少妇亚洲综合色aaa.| 亚洲精品美女久久av网站| 色老头精品视频在线观看| 波多野结衣一区麻豆| 欧美av亚洲av综合av国产av| 50天的宝宝边吃奶边哭怎么回事| 免费观看av网站的网址| 亚洲精品在线美女| e午夜精品久久久久久久| 欧美激情 高清一区二区三区| 午夜91福利影院| 18在线观看网站| 97在线人人人人妻| av天堂久久9| 最近中文字幕2019免费版| 一级a爱视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美一区二区三区在线观看 | 最新在线观看一区二区三区| 欧美老熟妇乱子伦牲交| 一本久久精品| 久久性视频一级片| 国产一区有黄有色的免费视频| 一级毛片女人18水好多| 国产成人精品久久二区二区91| 午夜精品国产一区二区电影| www.精华液| 亚洲人成77777在线视频| 久久人人爽人人片av| 亚洲国产av新网站| 国产伦理片在线播放av一区| 一区二区三区乱码不卡18| 大香蕉久久成人网| 久久久久久久国产电影| 夫妻午夜视频| 成人免费观看视频高清| 一边摸一边做爽爽视频免费| 性少妇av在线| 亚洲第一青青草原| 国产av国产精品国产| 啪啪无遮挡十八禁网站| 日韩制服丝袜自拍偷拍| 欧美少妇被猛烈插入视频| 国产免费现黄频在线看| 国产视频一区二区在线看| 国产成人免费观看mmmm| 久久久久精品人妻al黑| 老熟妇乱子伦视频在线观看 | 欧美 日韩 精品 国产| 国产男女内射视频| 国产成人精品在线电影| 亚洲欧美日韩另类电影网站| 超碰成人久久| 久久香蕉激情| 少妇被粗大的猛进出69影院| 欧美大码av| 久久人妻福利社区极品人妻图片| 亚洲精品成人av观看孕妇| 久久久久视频综合| 在线天堂中文资源库| 丝袜脚勾引网站| av片东京热男人的天堂| 亚洲久久久国产精品| 欧美精品亚洲一区二区| 久久这里只有精品19| 国产精品二区激情视频| 美女国产高潮福利片在线看| 老司机午夜十八禁免费视频| 午夜福利,免费看| 9热在线视频观看99| 老司机午夜十八禁免费视频| 一二三四社区在线视频社区8| 老司机在亚洲福利影院| 国产一区二区在线观看av|