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

    新型縮減矩陣構(gòu)造加快特征基函數(shù)法迭代求解*

    2019-09-21 05:50:24王仲根沐俊文林涵聶文艷
    物理學(xué)報 2019年17期
    關(guān)鍵詞:子域門限維數(shù)

    王仲根 沐俊文 林涵 聶文艷

    1) (安徽理工大學(xué)電氣與信息工程學(xué)院, 淮南 232001)

    2) (淮南師范學(xué)院機械與電氣工程學(xué)院, 淮南 232001)

    1 引 言

    矩量法[1]是求解頻域積分方程的一種有效數(shù)值方法, 被廣泛應(yīng)用于目標雷達散射截面計算、天線設(shè)計與分析、電磁環(huán)境預(yù)估、電磁兼容設(shè)計等領(lǐng)域.但矩量法需要對目標精細剖分, 隨著目標電尺寸的增大, 計算復(fù)雜度以及內(nèi)存需求都會急劇增大.為解決這個問題, 一些快速有效的矩量法被提出來, 如快速多極子法(fast multipole method,FMM)[2]、多層快速多極子法 (multilevel fast multipole method, MLFMM)[3,4]、預(yù)修正-快速傅立葉變換法(precorrected fast Fourier transform,P-FFT)[5]、自適應(yīng)積分法(adaptive integral method, AIM)[6]等, 這些方法可以降低矩陣向量積計算復(fù)雜度, 但不能減少未知數(shù)的數(shù)目.為降低未知數(shù)的數(shù)目, 有學(xué)者提出將宏基函數(shù)引入到矩量法中, 如子全域基函數(shù)法[7]、復(fù)合基函數(shù)法[8]、子域多層法[9]、特征模法[10,11]以及特征基函數(shù)法(characteristic basis function method, CBFM)[12-14],其中CBFM因考慮到各子域間的耦合作用而備受關(guān)注.為提高CBFM計算效率, 文獻[15]提出應(yīng)用物理光學(xué)法生成特征基函數(shù)(characteristic basis functions, CBFs), 但精確度不高; 文獻[16,17]分別應(yīng)用自適應(yīng)交叉近似-LU分解技術(shù)、自適應(yīng)交叉近似-奇異值分解來高效生成CBFs; 文獻[18,19]提出一種CBFs融合構(gòu)造方法, 提高了CBFM的計算精度; 文獻[20-23]應(yīng)用自適應(yīng)交叉近似(adaptive cross approximation, ACA)算法、快速偶極子法加快矩陣向量積運算, 提高縮減矩陣構(gòu)造效率; 文獻[24-26]將CBFM與MLFMM, AIM,P-FFT相結(jié)合, 通過迭代法求解縮減矩陣方程, 提高了CBFM分析電大目標電磁散射問題的能力,但是隨著目標電尺寸的增大, CBFs數(shù)目不斷增加,縮減矩陣維數(shù)會變得越來越大, 矩陣條件數(shù)變差,迭代求解縮減矩陣方程效率降低[27].

    本文提出一種新型縮減矩陣構(gòu)造方法, 應(yīng)用奇異值分解技術(shù)壓縮激勵源, 基于新激勵源求解出各子域的特征基函數(shù).運用伽略金方法構(gòu)造縮減矩陣時, 將新激勵源和特征基函數(shù)作為檢驗函數(shù)和基函數(shù), 得到一個對角子矩陣均為單位矩陣的縮減矩陣.新方法構(gòu)造的縮減矩陣與傳統(tǒng)方法構(gòu)造的縮減矩陣相比, 矩陣條件數(shù)得到了優(yōu)化, 迭代求解縮減矩陣方程的效率顯著提高, 并且該方法易于與MLFMM, AIM, P-FFT等算法相結(jié)合, 進一步提高了特征基函數(shù)法分析電大尺寸目標電磁散射問題的能力.

    2 特征基函數(shù)法

    CBFM首先將目標劃分為M個鄰接的子域,再將每個子域剖分成Ni個單元(i=1, 2,...,M).為獲得一組包含多角度電流信息的CBFs,CBFM采用不同入射方向和極化的激勵照射每個子域, 假設(shè)總的激勵數(shù)為Npws=2NθN?,Nθ,N?分別表示在θ,?方向上的激勵數(shù)目, 于是子域i上的主要特征基函數(shù)(primary characteristic basis functions, PCBFs)即可通過下式求得:

    式中,Ei表示擴展子域i的激勵矩陣, 維數(shù)為表示擴展子域i的自阻抗矩陣, 維數(shù)為為擴展子域i的電流系數(shù)矩陣,維數(shù)為為擴展子域i的未知數(shù)數(shù)目.通過直接求解(1)式, 得到擴展子域i的.由于采用多角度激勵源得到的必然含有冗余信息,故通過奇異值分解(singular value decomposition,SVD)壓縮矩陣去除冗余信息, 即

    式中,Ui和Vi均為酉矩陣, 維數(shù)分別為和為對角陣, 維數(shù)為通過設(shè)置合適的門限τ, 保留Ui中大于門限的前Ki個列向量并去除擴展部分作為子域i的最終CBFs.假設(shè)子域i經(jīng)過 SVD后得到Ki個 CBFs,則子域i的表面電流可由這Ki個CBFs線性組合表示:

    式中,ZR表示所有子域CBFs之間的相互作用,維數(shù)為是激勵向量;α為待求CBFs系數(shù).縮減矩陣方程構(gòu)造原理與矩量法構(gòu)造阻抗矩陣方程相似, 子域i上的第m個CBFs與子域j上第n個CBFs之間的相互作用表示為

    式中,Fi,m,Fj,n分 別表示子域i和j上的第m,n個 CBFs;fi,p(r) ,fj,q(r) 分別為子域i和j上的第p,q個 Rao-Wilton-Glisson (RWG)基 函 數(shù).Zij(p,q)=Zpq=〈fi,p(r),L(fj,q(r))〉,Zij表示子域i和j上所有RWG基函數(shù)之間的相互作用,Zpq是第p個RWG基函數(shù)和第q個RWG基函數(shù)之間的相互作用.Ji(p,m) 為聯(lián)系子域i中第m個CBFs和子域i中第p個RWG的線性標出系數(shù), 這些系數(shù)構(gòu)造矩陣Ji中的每一列對應(yīng)一個 CBFs;Jj(q,n)為聯(lián)系子域j中第n個 CBFs 和子域j中第q個RWG的線性標出系數(shù),Ni和Nj分別表示子域i和j所包含的RWG基函數(shù)的數(shù)目.因此子域i和j所有CBFs之間的相互作用可以表示為

    從(6)式和(7)式可以看出, 運用伽略金方法構(gòu)造縮減矩陣時, 使用Ji的共軛轉(zhuǎn)置同乘方程兩邊, 檢驗函數(shù)和基函數(shù)均采用CBFs(Ji).由(6)式和(7)式可得整個縮減矩陣方程的表達式為

    通過求解(8)式即可得到系數(shù)矩陣α, 通常求解(8)式可以選擇直接法求解, 但在分析電大復(fù)雜目標時, 縮減矩陣維數(shù)增大, 縮減矩陣方程需要通過迭代法求解.

    3 新型縮減矩陣構(gòu)造

    為提高縮減矩陣方程迭代求解效率, 本文提出一種新型縮減矩陣構(gòu)造方法, 首先應(yīng)用SVD對激勵矩陣進行壓縮:

    設(shè)定門限τ去除Ui中具有線性相關(guān)性的分量并將其表示為, 并將定義為激勵基函數(shù).假設(shè)每個子域經(jīng)過SVD后包含Li個激勵矢量, 將新的激勵源代入到(1)式, 求解出每個子域的CBFs:

    由于Li?Npws, 可以顯著減少方程求解次數(shù).通過求解(10)式, 每個子域可得到Li個將分別作為構(gòu)建縮減矩陣的檢驗函數(shù)和基函數(shù), 則縮減矩陣子矩陣可以表示為

    4 數(shù)值算例

    為驗證本文方法(novel characteristic basis function method, NCBFM)的有效性和精確性,分別對導(dǎo)體球、錐球帶縫體的雙站RCS以及杏仁體的單站RCS進行了計算.所有算例均在Intel(R) Core(TM) i5-6200U 2.30 GHz, 48 GB RAM 的PC 機上完成, 編譯器采用Visual studio 2013, BiCGStab迭代誤差為 0.001, 為了驗證NCBFM計算精度, 定義電流系數(shù)均方根誤差為

    算例1計算一個半徑為λ導(dǎo)體球的雙站RCS, 入射頻率為 300 MHz, 入射角度θ=0°,?=0°.應(yīng)用三角單元剖分球表面, 未知數(shù)為17278, 目標劃為8個子域.NCBFM和CBFM均為每個子域構(gòu)造800個激勵, 圖1(a)給出了2種方法在不同SVD門限下的電流誤差以及CBFs數(shù)目.從圖1(a)可以看出, 采用SVD壓縮激勵源, 電流誤差收斂速度更快; 另外, CBFM在分析電大復(fù)雜目標時, 激勵數(shù)目往往根據(jù)經(jīng)驗設(shè)定, 存在大量冗余計算, 而應(yīng)用SVD對激勵源進行壓縮, 只需設(shè)置合適的SVD門限, 就可以在保證精度的情況下減少冗余計算.根據(jù)電流誤差分析, CBFM和NCBFM的門限τ分別取0.005和0.008, 圖1(b)給出了左半球面4個子域SVD后奇異值的分布曲線.從圖1(b)可以看出, 2種方法在每個子域奇異值數(shù)目為80時即可達到門限設(shè)置要求, CBFM得到653個CBFs, 縮減矩陣維數(shù)為 6 53×653 , 矩陣條件數(shù)為5282, 采用BiCGStab迭代法求解縮減矩陣方程, 迭代26次即可收斂.NCBFM共得到649個CBFs, 縮減矩陣維數(shù)為 6 49×649 , 矩陣條件數(shù)為1785, 縮減矩陣方程迭代17次即可收斂,計算效率提高了34.6%.分別應(yīng)用CBFM和NCBFM計算了導(dǎo)體球HH極化雙站RCS, 計算結(jié)果如圖1(c)所示, 從圖1(c)可以看出,NCBFM與CBFM計算結(jié)果吻合較好, 計算精度較高.

    圖1 (a)不同SVD門限下2種方法的計算誤差及CBFs數(shù)目; (b)左半球面4個子域的奇異值分布曲線; (c)導(dǎo)體球雙站RCSFig.1.(a) Calculation error and numbers of CBFs under different SVD thresholds of two methods; (b) singular value distribution curve in four sub-domains of the left hemisphere; (c) bistatic RCS of PEC sphere.

    算例2計算一個錐球帶縫體的雙站RCS, 其幾何外形定義見文獻[28], 入射頻率為6 GHz, 入射角度θ=270°,?=0°.應(yīng)用三角單元剖分目標表面, 得到124685個未知數(shù), 目標被劃為48個子域.2種方法在每個子域上均設(shè)置1600個激勵, CBFM的SVD門限為0.001, 共得到7829個CBFs, 若采用直接法(LU分解)求解縮減矩陣方程需要耗時1239.6 s, 而采用迭代法求解縮減矩陣方程, 迭代86次即可收斂, 耗時104.5 s.NCBFM的SVD門限為0.002, 共得到7775個基函數(shù), 縮減矩陣方程迭代求解45次即可收斂, 耗時57.3 s, 計算效率提高了45.2%.2種方法計算的HH極化雙站RCS如圖2所示, 從圖2可以看出NCBFM計算的結(jié)果與CBFM和FEKO吻合較好.

    圖2 錐球帶縫體雙站RCSFig.2.Bistatic RCS of cone-sphere with gap.

    算例3計算一個252.3744 mm杏仁體的單站RCS, 入射頻率為20 GHz, 入射角為θ=90°,φ=0°- 1 80°.用三角單元對目標表面進行剖分,共得到153690個未知數(shù), 目標劃為52個子域.2種方法設(shè)置的激勵數(shù)均為1600, CBFM和NCBFM的SVD門限分別為0.001和0.002, 分別得到11410和11362個CBFs, 若采用直接法求解縮減矩陣方程, 需要耗時5387.3 s.應(yīng)用CBFM求解縮減矩陣方程, 迭代次數(shù)平均為128.2, 單次縮減矩陣方程求解平均耗時237.8 s; NCBFM迭代次數(shù)平均為63.5, 迭代次數(shù)減少了50.4%, 單次縮減矩陣方程求解平均耗時121.2 s.圖3給出了2種方法計算的HH極化單站RCS, 從圖3可以看出, NCBFM的計算結(jié)果與FEKO吻合較好, 具有較高的計算精度.

    圖3 杏仁體HH極化單站RCSFig.3.Monostatic RCS in HH polarization of NASA almond.

    表1給出了CBFM和NCBFM在阻抗矩陣填充、基函數(shù)構(gòu)造、縮減矩陣構(gòu)造以及縮減矩陣方程時的求解時間.從表1可以看出, NCBFM在基函數(shù)構(gòu)造方面計算效率有了小幅提高, 在縮減矩陣方程迭代求解方面計算效率都得到了顯著提高.

    5 結(jié) 論

    本文給出了一種新型縮減矩陣構(gòu)造方法, 該方法應(yīng)用奇異值分解技術(shù)壓縮激勵源, 并在新激勵源下求解出各子域的特征基函數(shù), 減少了冗余計算;在構(gòu)造縮減矩陣時, 選擇激勵基函數(shù)和特征基函數(shù)作為檢驗函數(shù)和基函數(shù), 將縮減矩陣的對角子矩陣優(yōu)化為單位矩陣, 提高了縮減矩陣方程的迭代求解效率.數(shù)值結(jié)果證明了本文方法在保證精度的前提下, 有效地提高了縮減矩陣方程的迭代求解效率.本文方法提高了縮減矩陣方程的迭代求解效率, 但在迭代過程中存在大量的矩陣向量積運算, 今后會進一步研究將多層快速多極子法、自適應(yīng)積分法、預(yù)修正-快速傅里葉變換法、快速偶極子法等算法引入到該方法中, 加快矩陣向量積運算, 以提高特征基函數(shù)法分析電大目標電磁散射特性的效率.

    表1 計算時間比較Table 1.Comparison of computation time.

    猜你喜歡
    子域門限維數(shù)
    β-變換中一致丟番圖逼近問題的維數(shù)理論
    基于鏡像選擇序優(yōu)化的MART算法
    基于規(guī)則的HEV邏輯門限控制策略
    地方債對經(jīng)濟增長的門限效應(yīng)及地區(qū)差異研究
    中國西部(2021年4期)2021-11-04 08:57:32
    基于子域解析元素法的煤礦疏降水量預(yù)測研究
    煤炭工程(2021年7期)2021-07-27 09:34:20
    一類齊次Moran集的上盒維數(shù)
    隨機失效門限下指數(shù)退化軌道模型的分析與應(yīng)用
    一種基于壓縮感知的三維導(dǎo)體目標電磁散射問題的快速求解方法
    關(guān)于齊次Moran集的packing維數(shù)結(jié)果
    涉及相變問題Julia集的Hausdorff維數(shù)
    精品国产超薄肉色丝袜足j| 欧美黄色淫秽网站| 满18在线观看网站| 国产亚洲欧美在线一区二区| 久久精品影院6| 1024香蕉在线观看| 久久久国产欧美日韩av| 国产精品98久久久久久宅男小说| 色婷婷久久久亚洲欧美| netflix在线观看网站| 一区二区三区国产精品乱码| 两个人看的免费小视频| 成熟少妇高潮喷水视频| 夜夜夜夜夜久久久久| 日韩高清综合在线| 国内毛片毛片毛片毛片毛片| 一个人免费在线观看的高清视频| 99热只有精品国产| 国产91精品成人一区二区三区| 宅男免费午夜| 亚洲精品久久成人aⅴ小说| www日本黄色视频网| 啦啦啦免费观看视频1| www日本黄色视频网| 男女做爰动态图高潮gif福利片| avwww免费| 色综合亚洲欧美另类图片| 久久精品国产亚洲av香蕉五月| 男女下面进入的视频免费午夜 | 精品不卡国产一区二区三区| 老司机在亚洲福利影院| 成年版毛片免费区| 白带黄色成豆腐渣| 波多野结衣巨乳人妻| 国产成+人综合+亚洲专区| 激情在线观看视频在线高清| 国产熟女xx| 成人国产综合亚洲| 免费电影在线观看免费观看| 国产又爽黄色视频| 日本 av在线| 午夜福利一区二区在线看| 18禁国产床啪视频网站| 亚洲成人国产一区在线观看| 国产精品九九99| 日韩欧美免费精品| 国产一区二区三区在线臀色熟女| 日韩欧美一区视频在线观看| 国产av一区二区精品久久| 亚洲国产欧洲综合997久久, | 国产成人av激情在线播放| 人妻丰满熟妇av一区二区三区| 波多野结衣巨乳人妻| 老司机深夜福利视频在线观看| 又黄又爽又免费观看的视频| 国产亚洲av高清不卡| 99久久久亚洲精品蜜臀av| 欧美又色又爽又黄视频| 亚洲成人精品中文字幕电影| 在线观看一区二区三区| 99热只有精品国产| 日韩三级视频一区二区三区| 变态另类丝袜制服| a级毛片a级免费在线| 中文字幕另类日韩欧美亚洲嫩草| 51午夜福利影视在线观看| 男女视频在线观看网站免费 | 亚洲一区中文字幕在线| 国产精品电影一区二区三区| 亚洲午夜理论影院| 欧美色视频一区免费| 最近最新免费中文字幕在线| 亚洲精华国产精华精| 国产精品久久久久久人妻精品电影| 日本a在线网址| 免费人成视频x8x8入口观看| 久久久久九九精品影院| 热99re8久久精品国产| 国产区一区二久久| 国产v大片淫在线免费观看| 在线观看午夜福利视频| 成年人黄色毛片网站| 亚洲一卡2卡3卡4卡5卡精品中文| 日本一区二区免费在线视频| 这个男人来自地球电影免费观看| 男人舔女人下体高潮全视频| 欧美亚洲日本最大视频资源| 日韩免费av在线播放| 美女午夜性视频免费| 亚洲国产毛片av蜜桃av| 久久国产精品男人的天堂亚洲| 亚洲男人天堂网一区| 欧美激情极品国产一区二区三区| 男女床上黄色一级片免费看| 亚洲男人天堂网一区| 村上凉子中文字幕在线| 国产主播在线观看一区二区| 国产男靠女视频免费网站| 成人三级黄色视频| 天天躁夜夜躁狠狠躁躁| 国产亚洲精品第一综合不卡| 波多野结衣av一区二区av| 自线自在国产av| 人人妻人人澡欧美一区二区| 999精品在线视频| 亚洲成人免费电影在线观看| 欧美激情久久久久久爽电影| 精品高清国产在线一区| 中文在线观看免费www的网站 | 久久天堂一区二区三区四区| 日韩精品青青久久久久久| 又紧又爽又黄一区二区| 精品一区二区三区av网在线观看| www.精华液| 欧美成人免费av一区二区三区| 麻豆av在线久日| 国产一区二区在线av高清观看| 亚洲激情在线av| 国产av在哪里看| 国产午夜福利久久久久久| 亚洲av片天天在线观看| 一级毛片女人18水好多| 中文在线观看免费www的网站 | 一二三四社区在线视频社区8| 久久久久久免费高清国产稀缺| 天天一区二区日本电影三级| 一级毛片精品| 欧美日韩亚洲国产一区二区在线观看| 制服丝袜大香蕉在线| 亚洲熟妇中文字幕五十中出| 免费在线观看完整版高清| 成年女人毛片免费观看观看9| 国产精华一区二区三区| 母亲3免费完整高清在线观看| 久久久久久久精品吃奶| 91成人精品电影| 日韩精品青青久久久久久| 一区二区三区高清视频在线| 久久精品91蜜桃| 国产精品国产高清国产av| 一区二区三区国产精品乱码| 精品久久久久久久久久免费视频| 国产日本99.免费观看| 免费在线观看黄色视频的| 久久久久九九精品影院| 18美女黄网站色大片免费观看| 免费人成视频x8x8入口观看| 亚洲aⅴ乱码一区二区在线播放 | 大型黄色视频在线免费观看| 非洲黑人性xxxx精品又粗又长| 黑丝袜美女国产一区| 男女床上黄色一级片免费看| 国产精品久久久久久精品电影 | 精品久久久久久久毛片微露脸| 19禁男女啪啪无遮挡网站| 久久中文字幕一级| 成在线人永久免费视频| 国产亚洲欧美在线一区二区| 亚洲国产精品sss在线观看| 麻豆久久精品国产亚洲av| 亚洲国产欧美日韩在线播放| 老熟妇乱子伦视频在线观看| 亚洲人成网站高清观看| 精品久久久久久久久久久久久 | 侵犯人妻中文字幕一二三四区| 欧美激情久久久久久爽电影| 国产精品乱码一区二三区的特点| 久久久久免费精品人妻一区二区 | 日本成人三级电影网站| 久久 成人 亚洲| 国产精品久久电影中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成伊人成综合网2020| 叶爱在线成人免费视频播放| 久久久精品欧美日韩精品| 亚洲av美国av| 亚洲精品国产区一区二| 午夜影院日韩av| 亚洲真实伦在线观看| 亚洲av五月六月丁香网| 精品日产1卡2卡| 在线av久久热| 亚洲第一欧美日韩一区二区三区| 别揉我奶头~嗯~啊~动态视频| 欧美成人性av电影在线观看| 久久久水蜜桃国产精品网| 黄色视频不卡| 亚洲一码二码三码区别大吗| 亚洲人成电影免费在线| 一区二区三区精品91| 精品熟女少妇八av免费久了| 久久久国产成人免费| 无人区码免费观看不卡| 久9热在线精品视频| 后天国语完整版免费观看| 国内精品久久久久久久电影| 嫁个100分男人电影在线观看| 最近最新中文字幕大全免费视频| 黄频高清免费视频| 满18在线观看网站| 99国产精品一区二区蜜桃av| 亚洲精华国产精华精| 麻豆av在线久日| 母亲3免费完整高清在线观看| 1024视频免费在线观看| 亚洲精品美女久久av网站| 免费高清视频大片| 久久久久久大精品| 日韩三级视频一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 成人免费观看视频高清| 一个人免费在线观看的高清视频| 亚洲av电影在线进入| 国产av不卡久久| 老汉色av国产亚洲站长工具| 久久久久久人人人人人| 精品一区二区三区视频在线观看免费| 777久久人妻少妇嫩草av网站| 精品第一国产精品| 午夜久久久在线观看| 制服人妻中文乱码| 999久久久精品免费观看国产| 国产成人欧美| 99久久久亚洲精品蜜臀av| 国产男靠女视频免费网站| 美女免费视频网站| 欧美黑人精品巨大| 久久国产乱子伦精品免费另类| 91成人精品电影| 国产av不卡久久| 十八禁网站免费在线| 国产一区二区激情短视频| 日本免费a在线| 99国产精品一区二区蜜桃av| av超薄肉色丝袜交足视频| 制服诱惑二区| 在线观看免费午夜福利视频| 成年免费大片在线观看| 亚洲av电影不卡..在线观看| 亚洲天堂国产精品一区在线| 黄色毛片三级朝国网站| 国产亚洲精品一区二区www| 女人爽到高潮嗷嗷叫在线视频| 99国产精品一区二区三区| 亚洲色图av天堂| 999精品在线视频| 亚洲欧美日韩无卡精品| 亚洲成av片中文字幕在线观看| 88av欧美| 久久久精品欧美日韩精品| 18美女黄网站色大片免费观看| av福利片在线| 熟女少妇亚洲综合色aaa.| 亚洲第一电影网av| 看片在线看免费视频| 亚洲va日本ⅴa欧美va伊人久久| 色婷婷久久久亚洲欧美| 国产乱人伦免费视频| 国产精品免费视频内射| 久久精品91蜜桃| 久久亚洲精品不卡| 成人午夜高清在线视频 | www日本在线高清视频| 婷婷精品国产亚洲av在线| 在线天堂中文资源库| 国产精品,欧美在线| 亚洲精品粉嫩美女一区| 又大又爽又粗| 国产片内射在线| 久久精品国产综合久久久| 国产精品99久久99久久久不卡| 男女做爰动态图高潮gif福利片| 免费搜索国产男女视频| 国产97色在线日韩免费| 久久久久久久久中文| 一进一出抽搐gif免费好疼| 久久国产乱子伦精品免费另类| 琪琪午夜伦伦电影理论片6080| 国产精品美女特级片免费视频播放器 | 大型黄色视频在线免费观看| 此物有八面人人有两片| 国产成年人精品一区二区| 性色av乱码一区二区三区2| 真人一进一出gif抽搐免费| 黄色丝袜av网址大全| 国产成+人综合+亚洲专区| 婷婷精品国产亚洲av| 12—13女人毛片做爰片一| 亚洲精品久久成人aⅴ小说| 巨乳人妻的诱惑在线观看| 国产亚洲欧美在线一区二区| 精品熟女少妇八av免费久了| 美女扒开内裤让男人捅视频| 欧美乱色亚洲激情| 亚洲真实伦在线观看| 国产三级黄色录像| a在线观看视频网站| 欧美黄色淫秽网站| 国内精品久久久久精免费| 人妻丰满熟妇av一区二区三区| 欧美一区二区精品小视频在线| 免费电影在线观看免费观看| 美女高潮喷水抽搐中文字幕| 亚洲五月天丁香| 91麻豆av在线| 国产99白浆流出| 在线观看一区二区三区| 极品教师在线免费播放| 一本久久中文字幕| 一进一出抽搐动态| 精品久久久久久久人妻蜜臀av| 美国免费a级毛片| 此物有八面人人有两片| 国产一区二区激情短视频| 欧美日韩一级在线毛片| 首页视频小说图片口味搜索| 午夜久久久在线观看| 在线国产一区二区在线| 日韩中文字幕欧美一区二区| 12—13女人毛片做爰片一| 99精品欧美一区二区三区四区| 国产欧美日韩精品亚洲av| 国产麻豆成人av免费视频| 麻豆成人午夜福利视频| 欧美另类亚洲清纯唯美| 欧美乱码精品一区二区三区| 国产av在哪里看| 久久精品国产亚洲av香蕉五月| 久久久久久亚洲精品国产蜜桃av| 国产又黄又爽又无遮挡在线| 女性被躁到高潮视频| 999久久久精品免费观看国产| 久久婷婷人人爽人人干人人爱| 亚洲国产精品成人综合色| a在线观看视频网站| 丰满人妻熟妇乱又伦精品不卡| 欧美精品亚洲一区二区| 久久婷婷人人爽人人干人人爱| 一个人观看的视频www高清免费观看 | 男女午夜视频在线观看| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 九色国产91popny在线| 久久久久久人人人人人| 国产精品二区激情视频| 麻豆av在线久日| 亚洲男人的天堂狠狠| 亚洲最大成人中文| 少妇粗大呻吟视频| 久久欧美精品欧美久久欧美| 国产免费av片在线观看野外av| 欧美日韩亚洲综合一区二区三区_| 国产成人av教育| 免费在线观看亚洲国产| 人妻丰满熟妇av一区二区三区| 熟女电影av网| 人人妻,人人澡人人爽秒播| 露出奶头的视频| 大香蕉久久成人网| 久久伊人香网站| 国产亚洲精品综合一区在线观看 | 国内精品久久久久精免费| 麻豆成人av在线观看| 亚洲人成网站在线播放欧美日韩| 国产野战对白在线观看| 久久久久久大精品| 国产成人一区二区三区免费视频网站| 国产又黄又爽又无遮挡在线| 在线免费观看的www视频| 中文字幕另类日韩欧美亚洲嫩草| 久久性视频一级片| 淫妇啪啪啪对白视频| 中文字幕人妻熟女乱码| 国产亚洲精品av在线| 波多野结衣高清无吗| 婷婷六月久久综合丁香| 99久久无色码亚洲精品果冻| 亚洲aⅴ乱码一区二区在线播放 | 女同久久另类99精品国产91| 黄色 视频免费看| 啦啦啦免费观看视频1| 亚洲国产欧美一区二区综合| 国产成人一区二区三区免费视频网站| 亚洲va日本ⅴa欧美va伊人久久| 久热这里只有精品99| 国产又爽黄色视频| 日日夜夜操网爽| 国产精品乱码一区二三区的特点| 午夜福利高清视频| 国产又色又爽无遮挡免费看| 久久精品国产亚洲av高清一级| aaaaa片日本免费| 亚洲av成人一区二区三| 日日爽夜夜爽网站| 一进一出抽搐gif免费好疼| 久久久国产欧美日韩av| av欧美777| 午夜福利在线在线| 亚洲第一青青草原| 亚洲五月天丁香| 日韩欧美国产在线观看| 在线观看www视频免费| av欧美777| 一本久久中文字幕| 一级毛片精品| 色在线成人网| 国产精品 国内视频| 亚洲免费av在线视频| 免费在线观看日本一区| 丝袜在线中文字幕| 亚洲精品美女久久av网站| 日本一区二区免费在线视频| 狂野欧美激情性xxxx| 韩国av一区二区三区四区| 此物有八面人人有两片| 久久中文字幕一级| av福利片在线| 欧美最黄视频在线播放免费| 国产亚洲av高清不卡| 成人国产一区最新在线观看| 91成年电影在线观看| 非洲黑人性xxxx精品又粗又长| 免费观看精品视频网站| 99riav亚洲国产免费| 国产三级黄色录像| av福利片在线| www.精华液| 成人特级黄色片久久久久久久| 国产av不卡久久| x7x7x7水蜜桃| 亚洲 欧美一区二区三区| 色老头精品视频在线观看| 一本综合久久免费| avwww免费| 国产成人精品无人区| 99久久精品国产亚洲精品| 国产成人欧美在线观看| 久久精品亚洲精品国产色婷小说| 亚洲国产毛片av蜜桃av| 免费女性裸体啪啪无遮挡网站| 中出人妻视频一区二区| 国产精品1区2区在线观看.| 草草在线视频免费看| 久久天躁狠狠躁夜夜2o2o| 精品国产乱子伦一区二区三区| 欧美又色又爽又黄视频| 亚洲七黄色美女视频| 中文字幕人妻熟女乱码| 91老司机精品| 国产成人欧美| 美女午夜性视频免费| 99久久无色码亚洲精品果冻| 日韩欧美在线二视频| 久久久久免费精品人妻一区二区 | 免费观看精品视频网站| 亚洲熟女毛片儿| 国产av在哪里看| 国产黄色小视频在线观看| 亚洲天堂国产精品一区在线| 一级黄色大片毛片| 最近最新中文字幕大全免费视频| 黑人操中国人逼视频| 国产一卡二卡三卡精品| 老司机福利观看| 国产又爽黄色视频| 精品日产1卡2卡| 国产真人三级小视频在线观看| 久久久国产成人精品二区| 99国产极品粉嫩在线观看| 久久久国产成人免费| 亚洲三区欧美一区| 午夜福利视频1000在线观看| 国语自产精品视频在线第100页| 国产成人精品久久二区二区91| 久久久久免费精品人妻一区二区 | 亚洲av日韩精品久久久久久密| 午夜免费鲁丝| 国产精品免费一区二区三区在线| 99热这里只有精品一区 | 51午夜福利影视在线观看| 91在线观看av| 91成年电影在线观看| 窝窝影院91人妻| 国产欧美日韩精品亚洲av| 亚洲男人的天堂狠狠| 精华霜和精华液先用哪个| 禁无遮挡网站| 亚洲人成网站在线播放欧美日韩| 天天添夜夜摸| 国产精品自产拍在线观看55亚洲| 人人澡人人妻人| 国产一区在线观看成人免费| 99在线人妻在线中文字幕| av片东京热男人的天堂| 国产精品国产高清国产av| 99精品欧美一区二区三区四区| 色哟哟哟哟哟哟| 亚洲成人国产一区在线观看| 最近在线观看免费完整版| 久9热在线精品视频| 中文资源天堂在线| 欧美黄色淫秽网站| 国产精品自产拍在线观看55亚洲| 一级作爱视频免费观看| 午夜免费激情av| 可以免费在线观看a视频的电影网站| 成人18禁在线播放| 人妻久久中文字幕网| 精品欧美国产一区二区三| 成年免费大片在线观看| 哪里可以看免费的av片| 长腿黑丝高跟| 在线观看免费日韩欧美大片| 午夜激情福利司机影院| 免费看a级黄色片| 高清毛片免费观看视频网站| 国产视频内射| 国产激情久久老熟女| 亚洲天堂国产精品一区在线| 免费无遮挡裸体视频| 中出人妻视频一区二区| 欧美日韩瑟瑟在线播放| 国产精品影院久久| 久久久久久国产a免费观看| 免费在线观看黄色视频的| 婷婷亚洲欧美| 美女高潮到喷水免费观看| 最新在线观看一区二区三区| 欧美性长视频在线观看| 久久精品国产亚洲av香蕉五月| 亚洲成a人片在线一区二区| 熟女电影av网| 免费观看人在逋| 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 好男人在线观看高清免费视频 | 国产男靠女视频免费网站| 国产精品久久久久久人妻精品电影| 日韩欧美国产在线观看| 日韩精品青青久久久久久| 久久久久久久午夜电影| 伦理电影免费视频| 黄网站色视频无遮挡免费观看| 国产精品爽爽va在线观看网站 | 大香蕉久久成人网| 日日摸夜夜添夜夜添小说| 久久亚洲精品不卡| 啦啦啦免费观看视频1| 国产欧美日韩一区二区精品| 757午夜福利合集在线观看| 精品国产美女av久久久久小说| 99热6这里只有精品| 国产精品国产高清国产av| 国产成年人精品一区二区| 免费在线观看亚洲国产| 亚洲一区中文字幕在线| 嫩草影院精品99| 中文字幕最新亚洲高清| 精华霜和精华液先用哪个| 婷婷丁香在线五月| 亚洲av五月六月丁香网| 大型黄色视频在线免费观看| 天堂√8在线中文| 熟女少妇亚洲综合色aaa.| 色综合婷婷激情| 美女免费视频网站| 免费在线观看成人毛片| 精品国产亚洲在线| 亚洲精品中文字幕一二三四区| xxx96com| 人人妻人人澡欧美一区二区| 嫁个100分男人电影在线观看| 国产精品国产高清国产av| 俄罗斯特黄特色一大片| 一级毛片高清免费大全| 我的亚洲天堂| 欧美又色又爽又黄视频| 男人操女人黄网站| 国产成+人综合+亚洲专区| 天堂√8在线中文| 中文字幕最新亚洲高清| 在线观看66精品国产| 日本a在线网址| 国产1区2区3区精品| 国产精品免费一区二区三区在线| 国产成人啪精品午夜网站| 欧美av亚洲av综合av国产av| 操出白浆在线播放| 精品国产国语对白av| 日韩三级视频一区二区三区| svipshipincom国产片| 精品国产美女av久久久久小说| 色av中文字幕| 村上凉子中文字幕在线| 人人妻人人澡人人看| 黄色 视频免费看| 国产不卡一卡二| 亚洲九九香蕉| 日韩高清综合在线| 国产v大片淫在线免费观看| 美女免费视频网站| 美女高潮到喷水免费观看| 夜夜躁狠狠躁天天躁| 老司机靠b影院| 麻豆av在线久日| 国产精品国产高清国产av| 国产伦人伦偷精品视频| 成人18禁在线播放| 国产三级在线视频| 免费在线观看亚洲国产| 波多野结衣高清无吗| 桃红色精品国产亚洲av| 久热爱精品视频在线9| 又大又爽又粗| 手机成人av网站|