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

    基于Matlab的磁測(cè)量實(shí)驗(yàn)數(shù)據(jù)處理方法研究

    2016-07-07 02:25:03劉文峰張敏剛

    劉文峰,張敏剛,李 蒙

    (太原科技大學(xué)材料科學(xué)與工程學(xué)院,太原030024)

    ?

    基于Matlab的磁測(cè)量實(shí)驗(yàn)數(shù)據(jù)處理方法研究

    劉文峰,張敏剛,李 蒙

    (太原科技大學(xué)材料科學(xué)與工程學(xué)院,太原030024)

    摘 要:隨著科學(xué)技術(shù)的發(fā)展,材料分析方法日益豐富,材料測(cè)試設(shè)備不斷更新,儀器測(cè)量精度也越來(lái)越高。然而,與測(cè)試設(shè)備匹配的專業(yè)數(shù)據(jù)處理軟件的開發(fā)和應(yīng)用卻發(fā)展緩慢。基于量子設(shè)計(jì)公司的振動(dòng)樣品磁強(qiáng)計(jì)(VSM)所測(cè)NdFeB磁性材料的M-H數(shù)據(jù),利用MATLAB7.0自帶擬合工具對(duì)數(shù)據(jù)進(jìn)行了三角函數(shù)擬合,并通過擬合反函數(shù)對(duì)磁測(cè)量數(shù)據(jù)進(jìn)行了運(yùn)算,得到相應(yīng)的矯頑力Hc和剩磁Mr值,從而準(zhǔn)確地表示了材料的磁學(xué)性能。

    關(guān)鍵詞:MATLAB;磁測(cè)量數(shù)據(jù);三角函數(shù)擬合;反函數(shù)

    隨著科學(xué)技術(shù)的發(fā)展,材料分析方法日益豐富,材料測(cè)試設(shè)備不斷更新,儀器測(cè)量精度也越來(lái)越高。然而,與測(cè)試設(shè)備匹配的專業(yè)數(shù)據(jù)處理軟件的開發(fā)和應(yīng)用卻發(fā)展緩慢。MATLAB作為國(guó)際學(xué)術(shù)界公認(rèn)的一種準(zhǔn)確、可靠的科學(xué)計(jì)算標(biāo)準(zhǔn)軟件,具有編程效率高、用戶使用方便、擴(kuò)充能力強(qiáng)、矩陣和數(shù)組運(yùn)算高效方便等優(yōu)點(diǎn)。本文基于量子設(shè)計(jì)公司的振動(dòng)樣品磁強(qiáng)計(jì)(VSM)所測(cè)磁性材料的MH數(shù)據(jù),利用MATLAB7.0軟件進(jìn)行了三角函數(shù)擬合,并通過擬合反函數(shù)對(duì)磁測(cè)量數(shù)據(jù)進(jìn)行了運(yùn)算,得到相應(yīng)的矯頑力Hc和剩磁Mr值,從而準(zhǔn)確地表示了材料的磁學(xué)性能。

    1 數(shù)據(jù)擬合和數(shù)值計(jì)算思路

    1.1 數(shù)據(jù)擬合思路

    應(yīng)用MATLAB對(duì)數(shù)據(jù)擬合有很大的隨機(jī)性,選用什么函數(shù),函數(shù)怎么設(shè)定等都是要考慮的因素,通常采取下面的方法來(lái)擬合數(shù)據(jù)。

    1)用編程序的方法進(jìn)行擬合。

    2)利用MATLAB自帶擬合工具進(jìn)行擬合。

    3)擬合反函數(shù)。

    1.2 矯頑力Hc和最大磁能積(BH)max的計(jì)算思路

    從正面擬合函數(shù)和擬合出的反函數(shù)中選取相應(yīng)簡(jiǎn)單的表達(dá)式進(jìn)行計(jì)算,計(jì)算結(jié)果出來(lái)之后代回剛才未采用的函數(shù)中進(jìn)行驗(yàn)證計(jì)算。

    2 數(shù)據(jù)擬合結(jié)果及分析

    2.1 應(yīng)用線性多項(xiàng)式進(jìn)行擬合結(jié)果及分析

    表1為實(shí)驗(yàn)所測(cè)釹鐵硼材料的M-H(磁化強(qiáng)度-外加磁場(chǎng))數(shù)據(jù)。文中對(duì)所有測(cè)試數(shù)據(jù)進(jìn)行了擬合和運(yùn)算,但實(shí)驗(yàn)所測(cè)數(shù)據(jù)量太大,此處選擇顯示了部分有代表性的數(shù)據(jù)。

    基于實(shí)驗(yàn)所測(cè)釹鐵硼的M-H數(shù)據(jù),用MATLAB作圖得出釹鐵硼材料的磁滯回線,如圖1所示。

    計(jì)算材料的矯頑力和磁能積,需要先求出退磁曲線的具體方程。因此,選取退磁曲線對(duì)應(yīng)的原始數(shù)據(jù),分別進(jìn)行三次多項(xiàng)式擬合、六次多項(xiàng)式擬合、九次多項(xiàng)式擬合和十五次多項(xiàng)式擬合,擬合曲線如圖2所示。

    由圖2可以看出,采用三次多項(xiàng)式擬合和六次多項(xiàng)式擬合效果不是很好,故就九次多項(xiàng)式和十五次多項(xiàng)式擬合結(jié)果單獨(dú)進(jìn)行分析。

    表1 釹鐵硼的M-H數(shù)據(jù)Tab.1 M-H Date of NdFeB

    圖1 釹鐵硼材料的磁滯回線Fig.1 Magnetic hysteresis loop of NdFeB

    圖3(a)和(b)所示分別為采用九次多項(xiàng)式擬合和十五次多項(xiàng)式擬合的結(jié)果。

    圖2 采用多項(xiàng)式擬合結(jié)果Fig.2 Results of polynomial fitting

    由圖3可以看出,九次多項(xiàng)式擬合曲線中間與散點(diǎn)圖出現(xiàn)較大偏離,十五次多項(xiàng)式擬合雖然整體吻合較,但曲線端部擬合不是很好,出現(xiàn)了動(dòng)蕩,所以這兩種擬合方法所得函數(shù)均不利于精確計(jì)算,而用更高次擬合多項(xiàng)式進(jìn)行擬合后,曲線端部震蕩的程度和范圍會(huì)不斷擴(kuò)大,因此不能用線性多項(xiàng)式進(jìn)行精確擬合,則嘗試其他擬合方法,例如指數(shù)擬合,三角函數(shù)擬合等。

    選用其他擬合方法,由于選擇比較多,單一的編程序擬合方法會(huì)顯得比較盲目,工作量大,由于數(shù)據(jù)比較多,很容易出現(xiàn)運(yùn)算或程序的錯(cuò)誤,為克服以上缺點(diǎn),我們利用MATLAB自帶擬合工具cftool進(jìn)行多項(xiàng)式擬合。

    2.2 利用MATLAB自帶擬合工具cftool進(jìn)行多項(xiàng)式擬合結(jié)果及分析

    利用MATLAB自帶擬合工具cftool進(jìn)行九次多項(xiàng)式擬合結(jié)果,如圖4所示。

    圖3 多項(xiàng)式擬合結(jié)果Fig.3 Results of polynomial fitting

    圖4 利用MATLAB自帶擬合工具cftool進(jìn)行九次多項(xiàng)式擬合結(jié)果Fig.4 Result of nine times polynomial fitting by cftool

    由圖4可以看出,九次多項(xiàng)式擬合效果很顯然不是很好,因此采用MATLAB自帶擬合工具cftool進(jìn)行三角函數(shù)擬合。

    2.3 利用MATLAB自帶擬合工具cftool進(jìn)行三角函數(shù)擬合結(jié)果及分析

    利用MATLAB自帶擬合工具cftool進(jìn)行三角函數(shù)擬合,結(jié)果如圖5所示。

    由圖5可以看出,利用MATLAB自帶擬合工具cftool進(jìn)行三角函數(shù)擬合結(jié)果較好。于是,整理得到M與H的關(guān)系如式(1).

    M(H)= 1.784sin(0.000131H + 0.8901)+ 0.6026sin(0.0004586H -2.006)+28.33sin(1.875e -005H +0.1286)+2.323sin(0.0002417H +2.051)+ 0.4889sin(0.0007855H + 1.272)+ 0.4471sin (0.0005988H -0.1153)+0.2262sin(0.001001H + 3.037)+0.04257sin(0.001245H -0.614)(1)

    圖5 利用MATLAB自帶擬合工具cftool進(jìn)行三角函數(shù)擬合結(jié)果Fig.5 Result of trigonometric function fitting by cftool

    3 數(shù)值計(jì)算結(jié)果及分析

    3.1 利用擬合出的關(guān)系式求矯頑力Hc、剩磁Mr

    3.1.1 直接利用所求三角函數(shù)求解

    通過擬合可以得到式(1),利用此式可以求解矯頑力Hc(即當(dāng)M = 0時(shí),所對(duì)應(yīng)的H的值)。但是,若直接令M = 0,求解以下方程:

    0 =1.784sin(0.000131H +0.8901)+ 0.6026sin (0.0004586H - 2.006)+ 28.33sin(187500H + 0.1286)+2.323sin(0.0002417H +2.051)+0.4889sin (0.0007855H +1.272)+0.4471sin(0.0005988H -0.1153)+0.2262sin(0.001001H +3.037)+0.04257sin (0.001245H -0.614),則非常困難,因此采用擬合反函數(shù)法,導(dǎo)出H = f(M)的關(guān)系。

    3.1.2 擬合反函數(shù)法,導(dǎo)出H = f(M)的關(guān)系

    通過以上擬合發(fā)現(xiàn)擬合出的函數(shù)在指定區(qū)間是單調(diào)且是單值函數(shù),由反函數(shù)存在定理可知,在該區(qū)間內(nèi)一定存在反函數(shù)。因此,通過擬合反函數(shù)來(lái)求解矯頑力。

    具體步驟如下:

    1)對(duì)調(diào)H和M,即讓H為縱坐標(biāo),M為橫坐標(biāo),再做一次數(shù)據(jù)擬合,擬合出來(lái)的函數(shù)關(guān)系為H = f(M)的形式。

    2)令M = 0,可以得到H的值,此時(shí)的H即為矯頑力Hc.

    3)利用MATLAB自帶擬合工具進(jìn)行擬合;

    擬合結(jié)果如下:

    首先,出現(xiàn)數(shù)據(jù)散點(diǎn)如圖6所示。

    圖6 擬合反函數(shù)法所得結(jié)果Fig.6 Result of inverse function fitting

    圖7 選取合適函數(shù)擬合結(jié)果Fig.7 Result of selected function fitting

    擬合函數(shù)如式(2):

    f(x)= -0.000133x^9 -0.0005971x^8 +0.0338x7+ 0.1402x^6 -2.295x^5 -9.042x^4 +58.89x^3 +238.4x^2 + 289.6x -8889

    則H與M的關(guān)系如下:

    H(M)= -0.000133M9-0.0005971M8+0.0338M7+ 0.1402M6-2.295M5-9.042M4+58.89M3+238.4M2+ 289.6M -8889(2)3.1.3 利用反函數(shù)求解矯頑力Hc和剩磁Mr

    1)求解矯頑力Hc:

    令M =0,代入式(2),計(jì)算得到:

    H(0)= -8889 Oe

    因此所求的矯頑力Hc = -8889 Oe.

    2)求解剩磁Mr:

    令H = 0時(shí),代入式(2),利用MATLAB求解,結(jié)果如下:

    則根據(jù)經(jīng)驗(yàn)和實(shí)際我們?nèi)?shí)根:Mr =6.92209memu

    3.2 驗(yàn)證矯頑力Hc和剩磁Mr的正確性

    由于矯頑力Hc及剩磁Mr都是用擬合出來(lái)的反函數(shù)進(jìn)行求解,即式(2)求解出來(lái)的。在這里通過正擬合出來(lái)的三角函數(shù)式(1)驗(yàn)證所求的H與M,說明結(jié)果的正確性。

    將H = - 8889Oe代入式(1),求Mr的值,若Mr的值與0接近,說明矯頑力求解正確。

    利用MATLAB求解如下:

    通過驗(yàn)證發(fā)現(xiàn)Mr與0基本接近,說明矯頑力Hc的精確度比較高。

    同理當(dāng)Mr = 6.92209memu代入式(1),計(jì)算H的值,若H的值越接近0,說明剩磁Mr求解比較精確。

    利用MATLAB求解:

    在曲線擬合中,由于擬合函數(shù)本身存在誤差,故此誤差不算很大。

    4 結(jié)論

    利用MATLAB7.0自帶擬合工具cftool對(duì)磁性材料的M-H測(cè)試結(jié)果進(jìn)行三角函數(shù)擬合,可以得到較好的結(jié)果。在此基礎(chǔ)上,再通過擬合反函數(shù)進(jìn)行計(jì)算,可以得到磁性材料的矯頑力Hc和剩磁Mr的值,進(jìn)而可以準(zhǔn)確地表示材料的磁學(xué)性能。

    參考文獻(xiàn):

    [1] 楊云升.Matlab曲線擬合及其在試驗(yàn)數(shù)據(jù)處理中的應(yīng)用[J].電腦與信息技術(shù),2009,17(2):34-36.

    [2] 陳敏.應(yīng)用Matlab擬合傳感器特性曲線[J].南京師范大學(xué)學(xué)報(bào),2003,3(1):45-49.

    [3] 熊學(xué)輝.用改進(jìn)單純形法優(yōu)化由Matlab曲線擬合得到的表達(dá)式[J].黃岡師范學(xué)院學(xué)報(bào),2003,23(3):67-68.

    [4] 查艷萍,王剛.多目標(biāo)線性規(guī)劃MATLAB軟件求解[J].棗莊學(xué)院學(xué)報(bào),2010,27(5):25-27.

    [5] 徐春環(huán),包娟,高楊.利用MATLAB軟件處理X射線圖像[J].數(shù)理醫(yī)藥學(xué)雜志,2009,22(4):415-417.

    [6] 胡異丁,甘俊英.MATLAB軟件在數(shù)字圖像處理中的應(yīng)用[J].軟件導(dǎo)刊,2008,7(5):139-140.

    [7] 沈順成,何博,黎章杰.基于Matlab的C2,C3類氣體壓縮機(jī)功率求解[J].武漢理工大學(xué)學(xué)報(bào),2015,37(6):832-836.

    [8] 鄭君剛,陳彪,樊旭峰,等.大學(xué)物理實(shí)驗(yàn)數(shù)據(jù)處理的Matlab應(yīng)用[J].大學(xué)物理實(shí)驗(yàn),2015,28(2):116-117,124.

    Research on Magnetic Measurement Data Processing Method Based on Mtlab

    LIU Wen-feng,ZHANG Min-gang,LI Meng
    (School of Materials Science and Engineering,Taiyuan University of Science and Technology,Taiyuan 030024,China)

    Abstract:With the development of science and technology,the material analysis methods become increasingly rich,the material testing equipment constantly updates,and the testing accuracy gets higher and higher.However,the development and application of the professional data processing software that matches with the test equipment is developed slowly.In this work,based on the magnetic M-H data measured by the quantum design Vibrating Sample Magnetometer(VSM),the trigonometric function fitting and inverse function fitting were adopted to numerical calculation by MATLAB 7.0 software,thus the coercive force and remnant magnetization can be obtained,so the material magnetic properties were accurately analyzed.

    Key words:MATLAB,VSM data,trigonometric function fitting,inverse function

    中圖分類號(hào):TB31

    文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/ j.issn.1673 -2057.2016.03.012

    文章編號(hào):1673 -2057(2016)03 -0224 -05

    收稿日期:2016-04-01

    作者簡(jiǎn)介:劉文峰(1981 -),女,講師,主要研究方向?yàn)榇判圆牧系闹苽浜托阅苎芯俊?/p>

    极品少妇高潮喷水抽搐| 欧美黑人精品巨大| 国产精品av久久久久免费| 亚洲美女黄色视频免费看| 中文精品一卡2卡3卡4更新| 久久青草综合色| 少妇 在线观看| 精品视频人人做人人爽| 国产女主播在线喷水免费视频网站| 99热网站在线观看| 精品久久久精品久久久| 在线亚洲精品国产二区图片欧美| 欧美激情极品国产一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av电影在线观看一区二区三区| 美女国产高潮福利片在线看| 国产主播在线观看一区二区 | 王馨瑶露胸无遮挡在线观看| 嫩草影视91久久| www.999成人在线观看| 亚洲七黄色美女视频| 日本vs欧美在线观看视频| 高清不卡的av网站| 嫩草影视91久久| 亚洲中文av在线| 麻豆国产av国片精品| 久久国产精品人妻蜜桃| 999精品在线视频| 久久久欧美国产精品| 欧美精品一区二区大全| 欧美中文综合在线视频| 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| 中文字幕另类日韩欧美亚洲嫩草| 91麻豆av在线| 99国产精品一区二区三区| 91成人精品电影| 精品久久久久久久毛片微露脸 | 首页视频小说图片口味搜索 | 性色av乱码一区二区三区2| 久久久精品国产亚洲av高清涩受| 欧美乱码精品一区二区三区| 91成人精品电影| 成在线人永久免费视频| 成人国语在线视频| 老司机影院成人| videos熟女内射| 色精品久久人妻99蜜桃| 亚洲av电影在线观看一区二区三区| 少妇 在线观看| 午夜福利一区二区在线看| 国产福利在线免费观看视频| 久久久国产一区二区| 在现免费观看毛片| 两个人免费观看高清视频| 大型av网站在线播放| 久久久精品国产亚洲av高清涩受| 午夜福利影视在线免费观看| 精品视频人人做人人爽| 欧美日韩成人在线一区二区| 人人妻人人添人人爽欧美一区卜| 午夜精品国产一区二区电影| 欧美日韩综合久久久久久| 99热国产这里只有精品6| 啦啦啦在线免费观看视频4| 亚洲国产欧美日韩在线播放| 在线精品无人区一区二区三| 午夜av观看不卡| 人人妻人人添人人爽欧美一区卜| 乱人伦中国视频| 成人亚洲精品一区在线观看| 首页视频小说图片口味搜索 | 久久精品亚洲av国产电影网| 日日爽夜夜爽网站| 丁香六月天网| 精品免费久久久久久久清纯 | 欧美黑人欧美精品刺激| 国产成人av教育| 中文字幕高清在线视频| 新久久久久国产一级毛片| 精品国产一区二区三区久久久樱花| 操美女的视频在线观看| 丝袜脚勾引网站| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 亚洲欧美日韩高清在线视频 | 亚洲欧美精品综合一区二区三区| 午夜福利,免费看| 国产精品成人在线| 丰满少妇做爰视频| 午夜免费观看性视频| 麻豆av在线久日| 国产人伦9x9x在线观看| 亚洲精品国产av蜜桃| 黄色一级大片看看| 成人国语在线视频| 亚洲国产精品成人久久小说| 午夜福利,免费看| 欧美人与性动交α欧美精品济南到| 脱女人内裤的视频| 亚洲精品一区蜜桃| 精品福利永久在线观看| 青春草亚洲视频在线观看| 建设人人有责人人尽责人人享有的| 午夜福利,免费看| 麻豆国产av国片精品| 男男h啪啪无遮挡| 97在线人人人人妻| 男的添女的下面高潮视频| 久久久久精品国产欧美久久久 | 久久精品国产a三级三级三级| 国产精品一区二区在线观看99| 18禁裸乳无遮挡动漫免费视频| 国产精品熟女久久久久浪| 亚洲精品久久成人aⅴ小说| 亚洲精品自拍成人| 黄色毛片三级朝国网站| 新久久久久国产一级毛片| 男女免费视频国产| 欧美中文综合在线视频| 国产欧美日韩一区二区三区在线| 最近最新中文字幕大全免费视频 | 国产在线视频一区二区| 中文字幕色久视频| 亚洲欧美中文字幕日韩二区| 亚洲成av片中文字幕在线观看| 在线观看免费午夜福利视频| 免费日韩欧美在线观看| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久电影网| 2021少妇久久久久久久久久久| 美女中出高潮动态图| 日韩av在线免费看完整版不卡| √禁漫天堂资源中文www| 日韩一区二区三区影片| 91精品三级在线观看| 亚洲欧洲国产日韩| 日韩 亚洲 欧美在线| 免费在线观看完整版高清| 亚洲天堂av无毛| 国产免费一区二区三区四区乱码| 国产精品久久久人人做人人爽| 国产精品成人在线| av不卡在线播放| 午夜免费观看性视频| 香蕉国产在线看| 丝袜美足系列| 欧美在线黄色| 69精品国产乱码久久久| 一区二区三区乱码不卡18| 男女无遮挡免费网站观看| 亚洲一区二区三区欧美精品| 50天的宝宝边吃奶边哭怎么回事| 成人国产一区最新在线观看 | 少妇 在线观看| 在线天堂中文资源库| 国产成人91sexporn| 成年人午夜在线观看视频| 久久久欧美国产精品| 最黄视频免费看| 一级毛片我不卡| 美女中出高潮动态图| 超色免费av| 91精品伊人久久大香线蕉| 国产精品麻豆人妻色哟哟久久| 亚洲国产欧美一区二区综合| 亚洲,一卡二卡三卡| 91成人精品电影| 日韩一本色道免费dvd| 亚洲 欧美一区二区三区| 男女午夜视频在线观看| 十八禁人妻一区二区| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 国产精品偷伦视频观看了| 亚洲欧美日韩高清在线视频 | 欧美日韩亚洲高清精品| 交换朋友夫妻互换小说| 亚洲av片天天在线观看| 日韩精品免费视频一区二区三区| 18禁黄网站禁片午夜丰满| 在现免费观看毛片| 亚洲av日韩精品久久久久久密 | 99国产精品一区二区三区| a级片在线免费高清观看视频| 99热国产这里只有精品6| 日本欧美国产在线视频| 高清视频免费观看一区二区| 汤姆久久久久久久影院中文字幕| 狂野欧美激情性xxxx| 两个人看的免费小视频| 亚洲国产精品一区三区| 亚洲欧美一区二区三区国产| 激情视频va一区二区三区| 婷婷成人精品国产| 在线亚洲精品国产二区图片欧美| 青春草视频在线免费观看| 日韩人妻精品一区2区三区| 七月丁香在线播放| 中文字幕av电影在线播放| 女人被躁到高潮嗷嗷叫费观| 精品一区二区三卡| 久久人妻福利社区极品人妻图片 | www日本在线高清视频| 美女午夜性视频免费| 久久九九热精品免费| 91九色精品人成在线观看| 日韩电影二区| 精品一区二区三区四区五区乱码 | 叶爱在线成人免费视频播放| 午夜免费鲁丝| 日韩欧美一区视频在线观看| 国产成人免费观看mmmm| 丝袜在线中文字幕| 中国美女看黄片| 中文字幕高清在线视频| a级片在线免费高清观看视频| 两人在一起打扑克的视频| www日本在线高清视频| 日韩精品免费视频一区二区三区| 天天添夜夜摸| 国产一区二区在线观看av| 日日夜夜操网爽| 日韩一区二区三区影片| 午夜福利一区二区在线看| 日本a在线网址| 日本色播在线视频| 一区二区三区乱码不卡18| 免费在线观看完整版高清| 成人国产av品久久久| 9色porny在线观看| 久久狼人影院| 免费不卡黄色视频| 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 国精品久久久久久国模美| avwww免费| 亚洲国产av新网站| 亚洲国产欧美日韩在线播放| 亚洲人成电影观看| 人妻人人澡人人爽人人| 国产免费现黄频在线看| 青青草视频在线视频观看| 国产成人系列免费观看| 少妇人妻久久综合中文| 午夜精品国产一区二区电影| 亚洲色图 男人天堂 中文字幕| av不卡在线播放| 九草在线视频观看| 人人澡人人妻人| 亚洲国产看品久久| 亚洲免费av在线视频| 国产精品一二三区在线看| 国产成人精品久久久久久| 国产99久久九九免费精品| 日韩精品免费视频一区二区三区| 两人在一起打扑克的视频| 超碰97精品在线观看| 在线看a的网站| 男男h啪啪无遮挡| 尾随美女入室| 国产高清国产精品国产三级| av欧美777| 99久久精品国产亚洲精品| 亚洲精品在线美女| 国产精品国产三级专区第一集| 婷婷色综合www| 久久鲁丝午夜福利片| 国产成人啪精品午夜网站| 99久久人妻综合| 大香蕉久久成人网| 一边摸一边做爽爽视频免费| 男的添女的下面高潮视频| netflix在线观看网站| 一二三四社区在线视频社区8| 一边亲一边摸免费视频| 欧美在线黄色| av又黄又爽大尺度在线免费看| 国产一级毛片在线| 国产成人精品在线电影| 另类精品久久| 脱女人内裤的视频| 午夜福利视频在线观看免费| 久久九九热精品免费| 亚洲,一卡二卡三卡| 一区二区日韩欧美中文字幕| 免费在线观看黄色视频的| 亚洲欧美中文字幕日韩二区| 精品少妇黑人巨大在线播放| 亚洲欧洲国产日韩| 久久久久网色| 精品少妇黑人巨大在线播放| 欧美日韩一级在线毛片| 中文字幕精品免费在线观看视频| 国产精品国产av在线观看| 黄色片一级片一级黄色片| 国产精品久久久久久精品电影小说| 亚洲精品美女久久久久99蜜臀 | 99国产精品一区二区三区| 1024香蕉在线观看| e午夜精品久久久久久久| 一级黄色大片毛片| 欧美黑人欧美精品刺激| 欧美黄色淫秽网站| 韩国高清视频一区二区三区| 大陆偷拍与自拍| 久久久精品94久久精品| 久久久久久人人人人人| 亚洲七黄色美女视频| 久久人人爽av亚洲精品天堂| 精品亚洲乱码少妇综合久久| 美女大奶头黄色视频| 中文字幕av电影在线播放| 久久亚洲精品不卡| 亚洲人成电影免费在线| 中文乱码字字幕精品一区二区三区| 欧美黑人欧美精品刺激| xxx大片免费视频| 又粗又硬又长又爽又黄的视频| 国产日韩欧美视频二区| 欧美激情极品国产一区二区三区| 成人国语在线视频| 久久国产精品影院| 岛国毛片在线播放| av网站在线播放免费| 欧美亚洲日本最大视频资源| 人体艺术视频欧美日本| 国产成人啪精品午夜网站| 在线亚洲精品国产二区图片欧美| 永久免费av网站大全| 国产1区2区3区精品| 免费少妇av软件| 国产精品一国产av| 女性被躁到高潮视频| 亚洲,欧美精品.| 亚洲五月婷婷丁香| 亚洲精品自拍成人| 伦理电影免费视频| 欧美亚洲日本最大视频资源| 五月开心婷婷网| 一本色道久久久久久精品综合| 亚洲欧美中文字幕日韩二区| 一本一本久久a久久精品综合妖精| 国产精品秋霞免费鲁丝片| 日韩免费高清中文字幕av| 久久精品熟女亚洲av麻豆精品| 国产精品一国产av| 91九色精品人成在线观看| 少妇被粗大的猛进出69影院| 日本av免费视频播放| 无限看片的www在线观看| 国产一卡二卡三卡精品| 精品国产一区二区久久| 精品福利观看| 在线天堂中文资源库| 欧美日韩亚洲国产一区二区在线观看 | 久久久亚洲精品成人影院| 我要看黄色一级片免费的| 香蕉丝袜av| 91九色精品人成在线观看| 成人亚洲欧美一区二区av| 少妇猛男粗大的猛烈进出视频| www.av在线官网国产| 高潮久久久久久久久久久不卡| 两个人免费观看高清视频| √禁漫天堂资源中文www| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 日韩欧美一区视频在线观看| 国产高清国产精品国产三级| 国产精品二区激情视频| 久久 成人 亚洲| av有码第一页| 91九色精品人成在线观看| 国产精品av久久久久免费| 国语对白做爰xxxⅹ性视频网站| 欧美日韩亚洲综合一区二区三区_| 精品国产国语对白av| 90打野战视频偷拍视频| 又黄又粗又硬又大视频| 国产福利在线免费观看视频| 久久久久国产精品人妻一区二区| 免费人妻精品一区二区三区视频| 久久久精品区二区三区| 中文字幕人妻熟女乱码| 免费在线观看影片大全网站 | 国产黄色视频一区二区在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品国产区一区二| 国产欧美日韩综合在线一区二区| av线在线观看网站| 最新在线观看一区二区三区 | 国产国语露脸激情在线看| 91字幕亚洲| 日韩一卡2卡3卡4卡2021年| videos熟女内射| 少妇人妻 视频| 午夜免费成人在线视频| 十分钟在线观看高清视频www| 国产野战对白在线观看| 成人三级做爰电影| 亚洲视频免费观看视频| 啦啦啦中文免费视频观看日本| 啦啦啦 在线观看视频| 国产免费视频播放在线视频| 大陆偷拍与自拍| 免费在线观看视频国产中文字幕亚洲 | 国产精品 国内视频| 秋霞在线观看毛片| 99国产综合亚洲精品| 久久精品亚洲熟妇少妇任你| 亚洲激情五月婷婷啪啪| 丰满饥渴人妻一区二区三| 亚洲美女黄色视频免费看| 中文字幕亚洲精品专区| 亚洲av男天堂| 午夜福利乱码中文字幕| 中国国产av一级| 人妻一区二区av| 精品少妇久久久久久888优播| 久久鲁丝午夜福利片| 大片免费播放器 马上看| 国产有黄有色有爽视频| 777久久人妻少妇嫩草av网站| 国产精品国产三级专区第一集| 婷婷成人精品国产| 日韩一卡2卡3卡4卡2021年| 国产一区有黄有色的免费视频| 熟女少妇亚洲综合色aaa.| 蜜桃在线观看..| 黄色a级毛片大全视频| 19禁男女啪啪无遮挡网站| av国产久精品久网站免费入址| 丝袜美腿诱惑在线| 国产精品久久久久久精品电影小说| 18在线观看网站| 国产一级毛片在线| 两人在一起打扑克的视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品熟女亚洲av麻豆精品| 久久久久精品国产欧美久久久 | 久久人妻熟女aⅴ| 777久久人妻少妇嫩草av网站| 国产黄频视频在线观看| 自线自在国产av| 狂野欧美激情性xxxx| 久热这里只有精品99| 久久综合国产亚洲精品| 免费高清在线观看日韩| 欧美日韩精品网址| 成人黄色视频免费在线看| 亚洲av在线观看美女高潮| 成年女人毛片免费观看观看9 | 高清av免费在线| 久久亚洲国产成人精品v| 免费在线观看日本一区| 欧美激情高清一区二区三区| 国产精品一区二区精品视频观看| 香蕉国产在线看| 亚洲国产精品国产精品| 国产亚洲欧美精品永久| 久久av网站| 精品少妇内射三级| 国产免费福利视频在线观看| 秋霞在线观看毛片| 少妇人妻 视频| 欧美在线黄色| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 久久精品亚洲熟妇少妇任你| 日韩大码丰满熟妇| 狠狠精品人妻久久久久久综合| 国产成人欧美在线观看 | 色播在线永久视频| 国产精品久久久久久人妻精品电影 | 久久午夜综合久久蜜桃| 午夜福利免费观看在线| 国产老妇伦熟女老妇高清| 一区二区三区精品91| 成人国产一区最新在线观看 | 人人妻人人爽人人添夜夜欢视频| 大片免费播放器 马上看| 永久免费av网站大全| 另类精品久久| 最新在线观看一区二区三区 | 久久久精品国产亚洲av高清涩受| 天天躁狠狠躁夜夜躁狠狠躁| 久久ye,这里只有精品| 岛国毛片在线播放| 老司机亚洲免费影院| 久久精品亚洲av国产电影网| 日本五十路高清| 日本猛色少妇xxxxx猛交久久| 亚洲,欧美精品.| 国产免费一区二区三区四区乱码| 男女之事视频高清在线观看 | 国产熟女欧美一区二区| 欧美成人午夜精品| 极品人妻少妇av视频| 久久亚洲精品不卡| 丁香六月欧美| 欧美精品一区二区大全| 中文字幕制服av| 777米奇影视久久| 国产精品二区激情视频| 又黄又粗又硬又大视频| 青草久久国产| 丝瓜视频免费看黄片| 婷婷丁香在线五月| 捣出白浆h1v1| 亚洲精品日本国产第一区| 国产精品一国产av| 国产日韩欧美在线精品| 黄色怎么调成土黄色| 国产免费现黄频在线看| 亚洲精品一卡2卡三卡4卡5卡 | 国产亚洲欧美精品永久| 久久鲁丝午夜福利片| www日本在线高清视频| 国产精品一国产av| 巨乳人妻的诱惑在线观看| 国产一区二区 视频在线| 大陆偷拍与自拍| 欧美成人午夜精品| 最新在线观看一区二区三区 | e午夜精品久久久久久久| 男女午夜视频在线观看| 男女免费视频国产| 免费不卡黄色视频| 亚洲av成人不卡在线观看播放网 | 亚洲美女黄色视频免费看| 国产一级毛片在线| 精品熟女少妇八av免费久了| 中文字幕亚洲精品专区| 日韩伦理黄色片| 精品人妻一区二区三区麻豆| 大码成人一级视频| 国产精品麻豆人妻色哟哟久久| av片东京热男人的天堂| 亚洲 国产 在线| 精品久久久久久电影网| 新久久久久国产一级毛片| 各种免费的搞黄视频| 欧美精品高潮呻吟av久久| 成人三级做爰电影| 午夜福利乱码中文字幕| av一本久久久久| 欧美老熟妇乱子伦牲交| 久久热在线av| 在线天堂中文资源库| 亚洲国产欧美网| 成年美女黄网站色视频大全免费| 女性被躁到高潮视频| 日韩人妻精品一区2区三区| 亚洲欧美激情在线| 国产精品一国产av| 亚洲 国产 在线| 首页视频小说图片口味搜索 | 中文精品一卡2卡3卡4更新| 国产主播在线观看一区二区 | 国产不卡av网站在线观看| 成人18禁高潮啪啪吃奶动态图| 久久久久精品人妻al黑| 咕卡用的链子| 亚洲五月婷婷丁香| 天天影视国产精品| 久久 成人 亚洲| 亚洲一区二区三区欧美精品| 男女午夜视频在线观看| 午夜免费鲁丝| 亚洲黑人精品在线| 嫩草影视91久久| 成人18禁高潮啪啪吃奶动态图| 国产伦人伦偷精品视频| 在线亚洲精品国产二区图片欧美| 一级毛片电影观看| 一级毛片黄色毛片免费观看视频| 老汉色∧v一级毛片| www.精华液| 色网站视频免费| 欧美日韩黄片免| 亚洲国产中文字幕在线视频| 黑丝袜美女国产一区| 亚洲熟女毛片儿| 国产精品人妻久久久影院| 黄色视频不卡| 国产精品免费视频内射| 久久久久久亚洲精品国产蜜桃av| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 久久久精品国产亚洲av高清涩受| 国产伦理片在线播放av一区| 久久久久网色| 老司机午夜十八禁免费视频| 中文字幕人妻丝袜一区二区| 老司机在亚洲福利影院| e午夜精品久久久久久久| 国产欧美日韩精品亚洲av| 日韩 亚洲 欧美在线| 十分钟在线观看高清视频www| 国产一区二区三区综合在线观看| 水蜜桃什么品种好| 最新在线观看一区二区三区 | 欧美另类一区| 少妇裸体淫交视频免费看高清 | 国产精品久久久人人做人人爽| 夜夜骑夜夜射夜夜干| 日韩伦理黄色片| 国产伦理片在线播放av一区| 国产精品久久久久久精品电影小说| 日本vs欧美在线观看视频| 老汉色av国产亚洲站长工具| 亚洲欧美日韩另类电影网站|