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

    紅外光譜結(jié)合化學(xué)計量學(xué)檢驗(yàn)直液式走珠筆墨水

    2022-08-23 05:14:12劉新磊韓丹巖楊尚鵬
    紅外技術(shù) 2022年8期
    關(guān)鍵詞:判別函數(shù)均值紅外

    劉新磊,韓丹巖,肖 強(qiáng),楊尚鵬

    〈紅外應(yīng)用〉

    紅外光譜結(jié)合化學(xué)計量學(xué)檢驗(yàn)直液式走珠筆墨水

    劉新磊,韓丹巖,肖 強(qiáng),楊尚鵬

    (中國刑事警察學(xué)院,遼寧 沈陽 110035)

    檢驗(yàn)鑒定書寫墨水是法庭科學(xué)領(lǐng)域中一項(xiàng)重要的工作。利用傅里葉變換紅外光譜儀對市面上常見的20種直液式走珠筆墨水樣品進(jìn)行檢驗(yàn),采取化學(xué)計量學(xué)實(shí)現(xiàn)了對直液式走珠筆墨水的快速檢驗(yàn)。對傅里葉變換紅外光譜數(shù)據(jù)做標(biāo)準(zhǔn)化處理,采取自動基線矯正、峰面積歸一化、Savitzky-Golay 5點(diǎn)平滑3種方法對譜圖做預(yù)處理,通過誤差平方和指標(biāo)(sum of the squares errors, SSE)確定分類K的最優(yōu)取值,對樣品進(jìn)行K-均值聚類分析并解釋聚類結(jié)果。采取主成分分析方法驗(yàn)證K-均值聚類的結(jié)果,利用組均值均等檢驗(yàn)考察主成分變量對Fisher判別分析(Fisher discriminant analysis, FDA)模型的貢獻(xiàn)程度,構(gòu)建直液式走珠筆墨水的FDA判別模型。結(jié)果表明,K-均值聚類將全部墨水樣本聚類為3類,PCA-FDA(principal components analysis-fisher discriminant analysis)模型能夠?qū)Σ煌悇e直液式走珠筆墨水達(dá)到100%的預(yù)測分類,經(jīng)交叉驗(yàn)證后正確率達(dá)到100%。紅外光譜與PCA-FDA模型結(jié)合可用于直液式走珠筆墨水的快速、準(zhǔn)確檢驗(yàn)鑒定。

    直液式走珠筆墨水;傅里葉變換紅外光譜;化學(xué)計量學(xué);判別分析

    0 引言

    直液式走珠筆是一種借鑒了鋼筆結(jié)構(gòu)的走珠筆,與傳統(tǒng)的走珠筆相比,墨水流動性更好,其墨水的性質(zhì)與鋼筆相似,比中性筆的墨水更稀,目前成為了市場上主流書寫工具。因此建立直液式走珠筆的油墨種類鑒別方法對于司法鑒定工作具有重大意義。書寫墨水的種屬鑒定一直是法庭科學(xué)領(lǐng)域的研究重點(diǎn),不僅可以為文件真?zhèn)舞b定提供依據(jù)[1],還能縮小偵查范圍[2],成為認(rèn)定案件事實(shí)的證據(jù)。目前國內(nèi)法庭科學(xué)領(lǐng)域區(qū)分直液式走珠筆這種新型筆種墨水的文獻(xiàn)較少[3],而傳統(tǒng)檢驗(yàn)筆墨水成分的方法有光學(xué)檢驗(yàn)法、色譜法、質(zhì)譜法[4]、光譜法等方法[5]。其中質(zhì)譜法[6]、色譜法、溶解法均屬于有損檢驗(yàn),不符合當(dāng)今無損檢驗(yàn)的主流,破壞了檢材案發(fā)時的完整性,而光譜法具有快速、無損檢驗(yàn)的特點(diǎn),廣泛應(yīng)用于墨水檢驗(yàn)中[3]。

    國內(nèi)外學(xué)者針對墨水檢驗(yàn)展開了許多實(shí)驗(yàn)研究,證實(shí)了傅里葉變換紅外光譜法結(jié)合主成分分析可以檢驗(yàn)?zāi)甗7-9],也有學(xué)者嘗試將化學(xué)計量學(xué)[10-14]等技術(shù)應(yīng)用于墨水、紙張等物證的檢驗(yàn),然而這些檢驗(yàn)方法較為依賴樣本標(biāo)簽,大多根據(jù)人為預(yù)設(shè)標(biāo)簽構(gòu)建模型對樣品進(jìn)行大致分類,且缺少了評估分類變量對模型的貢獻(xiàn)程度,難以做到客觀、準(zhǔn)確地檢驗(yàn)?zāi)?。目前K-均值聚類與組均值均等檢驗(yàn)具有較好的穩(wěn)定性,已被用于藥物、細(xì)菌的檢驗(yàn),并且得到了較好的結(jié)果[15-16]。因此本文嘗試結(jié)合光譜數(shù)據(jù)與K-均值聚類預(yù)先得出分類標(biāo)簽,采集市面上常見的20個直液式走珠筆樣品,利用傅里葉變換紅外光譜儀對樣品進(jìn)行檢驗(yàn)并獲得紅外光譜譜圖[8],采取自動基線矯正、峰面積歸一化、Savitzky-Golay5點(diǎn)平滑3種方法對譜圖做預(yù)處理后,通過誤差平方和指標(biāo)(sum of the squares errors,SSE)確定最優(yōu)取值,對樣品進(jìn)行K-均值聚類得到樣品的初步聚類結(jié)果,將聚類結(jié)果作為分類標(biāo)簽,采用主成分分析[10]對原始數(shù)據(jù)降維,對樣品進(jìn)行分類,利用組均值均等檢驗(yàn)考察主成分降維后的變量對Fisher判別模型的貢獻(xiàn)程度,根據(jù)考察結(jié)果選取判別函數(shù)建立Fisher判別模型,從而為直液式走珠筆油墨的快速鑒定提供一種新的方式。

    1 理論基礎(chǔ)

    1.1 傅里葉變換紅外光譜

    傅里葉變換紅外光譜(Fourier transform infrared spectroscopy, FTIR)是一種通過數(shù)學(xué)處理將傅里葉變換,結(jié)合計算機(jī)與紅外光譜技術(shù)的分析方法。傅里葉紅外光譜法首先測量目標(biāo)干涉圖,通過對該圖進(jìn)行傅里葉變換獲得對應(yīng)的紅外光譜數(shù)據(jù),是一種無損、快速檢驗(yàn)樣品成分的方法[17]。

    1.2 組內(nèi)誤差平方和

    組內(nèi)誤差平方和代表數(shù)據(jù)誤差大小的偏差平方和。隨著聚類的類別值增加,樣品會被歸類到更加精細(xì)的簇中,SSE會隨之減少,期間SSE的下降幅度會驟減形成一個拐點(diǎn),然后隨著值的繼續(xù)增大而趨于平緩,該拐點(diǎn)就是聚類類別值的最優(yōu)取值。

    1.3 K-均值聚類

    K-均值聚類是一種劃分樣本為個集合的聚類算法,是無監(jiān)督學(xué)習(xí)類型的算法,通過迭代找到最佳聚類個數(shù)的一種劃分方案,使得用個聚類的均值來表示相應(yīng)各類樣本時所得到的總體誤差最小,所以K均值聚類是以最小誤差平方和為劃分集合依據(jù)的聚類算法。其原理是根據(jù)預(yù)先設(shè)定的值,將個樣本作為聚類的個聚類中心,計算所得其余樣品和每個聚類中心的歐幾里距離,將其余樣品分配給距離最近的聚類中心,通過多次迭代達(dá)到最優(yōu)結(jié)果。

    1.4 主成分分析

    主成分分析(principal components analysis, PCA)是一種降低光譜數(shù)據(jù)特征維度的方法,通過線性變換把原始數(shù)據(jù)從原始坐標(biāo)系統(tǒng)變換到一個新坐標(biāo)系統(tǒng)中,從而達(dá)到降低原始數(shù)據(jù)特征維度,同時獲得對原始數(shù)據(jù)信息方差貢獻(xiàn)最大的特征??赏ㄟ^主成分分析提取樣品光譜數(shù)據(jù)的最大主成分,進(jìn)一步構(gòu)建2維主成分得分圖,驗(yàn)證K-均值聚類結(jié)果的準(zhǔn)確性。

    1.5 組均值均等檢驗(yàn)

    組均值均等檢驗(yàn)的目的是評價分組變量對分組模型貢獻(xiàn)程度,將分組變量作為因子得出每個自變量的單因素ANOVA(analysis of variance)分析結(jié)果,根據(jù)分析結(jié)果得出對分組模型貢獻(xiàn)程度較高的變量,之后用該變量構(gòu)建Fisher判別模型。

    1.6 判別分析

    判別分析是在樣本分類確定的前提下,按照未知樣品特征所反映的信息判別其類別歸屬的一種統(tǒng)計分析法。按照相應(yīng)的規(guī)則,構(gòu)建一至多個判別函數(shù),根據(jù)已有樣品的大量數(shù)據(jù)信息計算待定系數(shù)從而確定判別函數(shù),計算判別指標(biāo),進(jìn)而推斷未知樣品的歸屬[18]。

    2 實(shí)驗(yàn)部分

    2.1 實(shí)驗(yàn)儀器及條件

    實(shí)驗(yàn)儀器為NICOLET5700傅里葉變換紅外光譜儀,分辨率4cm-1,掃描次數(shù)16次,掃描面積為100mm×100mm,掃描范圍4000~400cm-1。

    2.2 實(shí)驗(yàn)樣品

    從各地收集不同品牌、不同型號樣本共20個,均為黑色墨水,其中樣品1~樣品10為白雪品牌,樣品11~樣品13為晨光品牌,樣品14~樣品20為百樂品牌。

    2.3 實(shí)驗(yàn)方法

    用收集到的直液式走珠筆油墨樣品制作檢材,在上述實(shí)驗(yàn)條件下進(jìn)行檢測,為了保證數(shù)據(jù)的可復(fù)現(xiàn)性,每份檢測均檢測3次,以確保實(shí)驗(yàn)結(jié)果的可靠性[19]。光譜數(shù)據(jù)中存在大量冗余信息,會影響后續(xù)化學(xué)計量學(xué)方法的使用,因此采用The Unscrambler X 10.4(挪威CAMO公司)對光譜信號進(jìn)行自動基線矯正、峰面積歸一化、Savitzky-Golay5點(diǎn)平滑處理[20]。下面將通過化學(xué)計量學(xué)對樣品光譜數(shù)據(jù)進(jìn)行分析[21]。

    3 結(jié)果與分析

    3.1 K-均值聚類

    的取值很大程度上決定了K-均值聚類的準(zhǔn)確性,故引入手肘法(elbow method)確定的最優(yōu)取值,手肘法的核心指標(biāo)是SSE,如式(1)所示:

    式中:C是第個簇;是C中的樣本點(diǎn);mC的質(zhì)心(C中所有樣本的均值);SSE是所有樣本的聚類誤差,SSE值越小代表聚類效果越好。

    聚類數(shù)值的增大會使得樣本相應(yīng)簇的劃分更加精準(zhǔn),SSE會逐漸變小。因此從小到大逐漸增加的過程中,SSE會逐漸變小,期間會存在一個真實(shí)聚類數(shù)的臨界點(diǎn),當(dāng)值從臨界點(diǎn)左側(cè)增加到真實(shí)聚類數(shù)時,SSE會發(fā)生迅速地下降,然后隨著逐漸增大趨于平緩。SSE和關(guān)系圖會有一個拐點(diǎn),該拐點(diǎn)對應(yīng)的值就是數(shù)據(jù)的真實(shí)聚類數(shù)。從圖1可知,從1~3時下降得很快,當(dāng)取值≥3后,曲線下降逐漸趨于平緩,所以最佳聚類個數(shù)為3。

    圖1 SSE與簇的數(shù)量關(guān)系圖

    當(dāng)=3時,對樣品進(jìn)行K-均值聚類,聚類結(jié)果圖如圖2所示。

    由圖2可得,當(dāng)=3時,所有樣品在二維平面被顯著地分為了3類,且各類樣品在紅外譜圖中存在明顯的區(qū)別,表現(xiàn)為:第一類樣品(7#)中存在CS伸縮振動,CH面內(nèi)彎曲振動以及C=C伸縮振動,出峰位置分別為1172cm-1,1388cm-1以及1577cm-1;第二類樣品(13#)中存在972cm-1的CH的變形振動,以及1353cm-1和1470cm-1的CH的面內(nèi)彎曲振動;第三類樣品(18#)中存在1461cm-1的CH的面內(nèi)彎曲振動,結(jié)果如表1和圖3所示。為了驗(yàn)證K-均值聚類這一無監(jiān)督學(xué)習(xí)的準(zhǔn)確性,根據(jù)聚類結(jié)果劃分樣品的類別,采用主成分分析將各個類別的樣品投影到二維平面進(jìn)行驗(yàn)證。

    圖2 K-均值聚類結(jié)果

    表1 樣品聚類結(jié)果

    圖3 三種樣品紅外譜圖比較

    3.2 主成分分析

    主成分分析是一種無監(jiān)督學(xué)習(xí)的方法,通過將原始多維數(shù)據(jù)投影到新坐標(biāo)系統(tǒng)重新組合成一組新的線性無關(guān)數(shù)據(jù)來代替原始數(shù)據(jù),同時盡可能地反映原始數(shù)據(jù)特征信息[21]。PCA共提取了6個PC,圖4顯示了不同PC的方差貢獻(xiàn)率及累計貢獻(xiàn)率,其中PC1方差貢獻(xiàn)率最大,為80.466%,前6個PC的累計方差貢獻(xiàn)率為98.556%,包含了原始數(shù)據(jù)的絕大多數(shù)信息。以PC1、PC2和PC3繪制PCA的2維主成分得分圖,如圖5所示,樣本基于其在PC1、PC2、PC3上的得分被分為了3類,組內(nèi)距離較小,組間距離較大,驗(yàn)證了K-均值聚類的聚類結(jié)果??紤]到實(shí)際中需要預(yù)測未知樣本,本文采用Fisher判別分別建立3類樣品的預(yù)測模型。

    圖4 各成分貢獻(xiàn)率分布圖

    3.3 組均值均等檢驗(yàn)

    為了實(shí)現(xiàn)對未知樣品的預(yù)測,構(gòu)建Fisher判別模型前,使用組均值均等檢驗(yàn)考察分組變量對分組模型貢獻(xiàn)程度。組間均方與組內(nèi)均方的比值為F統(tǒng)計量,兩個自由度分別為自由度1(df1)和自由度2(df2),分別表示殘差平方和自由度和回歸平方和自由度,自由度1為分子,自由度2為分母,兩個自由度用來得到觀測顯著性(Sig值)。顯著性(Sig)和Wilks’lambda是評價分組變量的標(biāo)準(zhǔn)。Sig值表示區(qū)分不同組別的能力,如果Sig值較?。⊿ig<0.10)則表明組間差異較為顯著,如果Sig值較大(Sig>0.10)則表示組間差異不顯著。PC1的Sig值為0.000,表面PC1在各組間的差異不顯著。組內(nèi)平方和與總平方和的比值為Wilks’lambda,值的范圍在0~1之間,值越小表示組內(nèi)有很大差異,值接近1表示沒有組內(nèi)差異。根據(jù)K-均值聚類結(jié)果劃分樣品類別,用組均值均等性檢驗(yàn),結(jié)果如表2所示,PC3~PC6的Wilks’Lambda均大于0.97,表明這4個變量組內(nèi)差異很小,對模型影響不顯著,且sig值均大于0.10,表明這四者不能充分解釋各樣本的分類情況。PC1和PC2的Wilks’Lambda均接近0,表明這兩個變量對判別模型影響的顯著性極高,且Sig值均為0,可表明這兩個變量可以很好地解釋各樣本的分組信息。因此選取前兩個變量構(gòu)建Fisher判別模型。

    圖5 PCA得分圖

    表2 各組平均值的均等性檢驗(yàn)

    3.4 Fisher判別模型

    在構(gòu)建Fisher判別模型中使用前兩個典則判別函數(shù)1和2,其中1=-0.028PC1-24.94PC2-10.932,2=34.764PC1-39.654PC2-48.976。表3為Fisher判別函數(shù)摘要,判別函數(shù)1(1)的方差貢獻(xiàn)率為81.0%,判別函數(shù)2(2)的方差貢獻(xiàn)率為19.0%,兩者累計方差貢獻(xiàn)率達(dá)到了100%,表明這兩個函數(shù)能完全反映樣品光譜特征中的信息。1和2的相關(guān)性大于0.95,表明不同樣本類別與1和2的相關(guān)性很強(qiáng)。觀察到1和2在函數(shù)檢驗(yàn)(Function test)中的Wilks’Lambda數(shù)值分別為0.002和0.083,均接近0,表明這兩個判別函數(shù)對模型影響的顯著性很高,Sig值均為0,表明這兩個判別函數(shù)對變量的可解釋度很高,具有充分解釋樣本分類的能力,綜上選用1和2這兩個函數(shù)作為判別函數(shù),構(gòu)建Fisher判別模型進(jìn)行分類,用留一法(Leave-One-Out Cross Validation)作為驗(yàn)證模型的方法。選擇函數(shù)1(1)和函數(shù)2(2)繪制判別函數(shù)聯(lián)合分布圖,如圖6所示。由圖6和表4的結(jié)果可得,該判別模型對3類樣本實(shí)現(xiàn)了正確率100%的分類,同時經(jīng)過留一法驗(yàn)證后模型正確率仍為100%。如果需要區(qū)分未知樣品,只需將該樣品相應(yīng)的紅外光譜數(shù)據(jù)輸入模型,判別函數(shù)聯(lián)合分布圖中就會顯示出該未知樣品位置及與其相距最近的分組質(zhì)心,從而推斷出未知樣品屬于哪一類別。

    表3 Fisher判別函數(shù)摘要

    圖6 判別函數(shù)聯(lián)合分布圖

    表4 分類結(jié)果

    4 結(jié)論

    本文借助了傅里葉變換紅外光譜儀對20種直液式走珠筆墨水成分及含量進(jìn)行了快速無損檢驗(yàn),建立了一種鑒別直液式走珠筆墨水類別的方法。首先對樣品的傅里葉變換紅外光譜數(shù)據(jù)進(jìn)行自動基線矯正、峰面積歸一化、Savitzky-Golay5點(diǎn)平滑3種預(yù)處理方法。采取組內(nèi)誤差平方和指標(biāo)得出K-均值聚類的最優(yōu)取值,對樣品進(jìn)行K-均值聚類,實(shí)現(xiàn)了對直液式走珠筆墨水樣品初步分類。對各類樣品K-均值聚類結(jié)果進(jìn)行主成分分析,繪制主成分二維得分圖使各類樣品之間的關(guān)系更加明確,同時驗(yàn)證了K-均值聚類結(jié)果的準(zhǔn)確性。使用組均值均等檢驗(yàn)評價PCA降維后的分組變量對FDA模型的貢獻(xiàn)程度,得到兩個具有較高貢獻(xiàn)程度的判別函數(shù),用該判別函數(shù)構(gòu)建Fisher判別模型對3類樣本進(jìn)行分類識別,3類樣品均被正確地分類為實(shí)際對應(yīng)的3種品牌,分類正確率達(dá)到100%,交叉驗(yàn)證后預(yù)測正確率為100%。可通過該模型對案件涉及的未知直液式走珠筆墨水物證進(jìn)行快速分類,為司法檢驗(yàn)鑒定提供幫助。

    [1] 趙巍, 高樹輝. 基于超光譜技術(shù)檢驗(yàn)添改文件的實(shí)驗(yàn)研究[J]. 中國人民公安大學(xué)學(xué)報:自然科學(xué)版, 2019, 25(1):11-15.

    ZHAO Wei, GAO Shuhui. Experimental study on inspection of added and modified documents based on hyperspectral technology[J].: Science and Technology, 2019, 25(1): 11-15.

    [2] 林建成, 梁魯寧, 黃建同, 等. 顯微激光拉曼光譜技術(shù)檢驗(yàn)彩色激光打印機(jī)打印文件墨跡的實(shí)驗(yàn)研究[J]. 光散射學(xué)報, 2015, 27(4): 374-378(DOI:10.13883/j.issn1004-5929.201504012.).

    LIN Jiancheng, LIANG Luning, HUANG Jiantong, et al. Experimental study on ink marks of documents printed by color laser printer with micro Raman spectroscopy[J]., 2015, 27(4): 374-378.

    [3] 國中正, 陳維娜, 王曉賓, 等. 超光譜技術(shù)結(jié)合化學(xué)計量法區(qū)分直液式走珠筆墨跡種類[J]. 理化檢驗(yàn)(化學(xué)分冊), 2020, 56(8): 878-882.

    GUO Zhongzheng, CHEN Weina, WANG Xiaobin, et al. Hyperspectral technology combined with chemometrical method for distinction of ink type of straight liquid ballpoint pen[J].(Part B:Chemical Analysis), 2020, 56(8): 878-882.

    [4] 王舒超, 崔嵐, 宋輝, 等. 氣相色譜-質(zhì)譜法分析3種常見品牌噴墨打印機(jī)84種墨水制備的墨跡樣本中揮發(fā)性溶劑成分[J]. 理化檢驗(yàn)(化學(xué)分冊), 2021, 57(9): 781-787.

    WAN Shuchao, CUI Lan, SONG Hui, et al. Analysis of volatile solvent components in ink blot samples prepared by 84 inks from 3 common brands of ink-jet printers by GC-MS[J].(Part B:Chemical Analysis), 2021, 57(9): 781-787.

    [5] 何海洋, 趙鵬程, 何寧, 等. 基于氣相色譜技術(shù)研究黑色中性筆字跡中3種溶劑隨書寫時間的變化規(guī)律[J]. 理化檢驗(yàn)(化學(xué)分冊), 2021, 57(9): 769-775.

    HE Haiyang, ZHAO Pengcheng, HE Ning, et al. Study on the variation of three solvents in black gel pen handwriting with writing time based on gas chromatography technology[J].(Part B:Chemical Analysis) , 2021, 57(9): 769-775.

    [6] 趙歷, 馮偉科, 莫少芳, 等. 頂空-氣相色譜-質(zhì)譜法測定記號筆油墨中揮發(fā)性有機(jī)化合物[J]. 理化檢驗(yàn)(化學(xué)分冊), 2017, 53(12): 1375-1380.

    ZHAO Li, FENG Weike, MO Shaofang. Determination of volatile organic compounds in mark pen ink by HS-GC-MS[J].(Part B:Chemical Analysis), 2017, 53(12): 1375-1380.

    [7] 牛凡, 黃建同, 張云, 等. 基于拉曼光譜和多元統(tǒng)計學(xué)研究藍(lán)色中性筆油墨[J]. 光散射學(xué)報, 2017, 29(3): 234-238.

    NIU Fan, HUANG Jiantong, ZHANG Yun, et al. Study on blue gel pen inks based on Raman spectroscopy and multivariate statistical[J]., 2017, 29(3): 234-238.

    [8] 侯偉, 王繼芬. 基于紅外指紋光譜的快速鑒別黑色記號筆墨水[J]. 激光技術(shù), 2020, 44(4): 436-440.

    HOU Wei, WANG Jifen. Rapid identification of the black marker ink based on infrared fingerprint spectroscopy[J]., 2020, 44(4): 436-440.

    [9] LukáGál, Michal Oravec, Pavol Gemeiner. Principal component analysis for the forensic discrimination of black inkjet inks based on the Vis-NIR fibre optics reflection spectra[J]., 2015, 257: 285-292.

    [10] 姚付軍. 化學(xué)計量學(xué)方法對藍(lán)色圓珠筆和簽字筆油墨的分類研究[D]. 北京: 首都師范大學(xué), 2009: 4-13.

    YAO Fujun. Classification of Blue Ballpoint Pen Inks and Blue Gel Pen Inks Using Chemometrics Methods[D]. Beijing: Capital Normal University, 2009: 4-13.

    [11] Mohamad Asri Muhammad Naeim, Mat Desa Wan Nur Syuhaila, Ismail Dzulkiflee. Source determination of red gel pen inks using Raman spectroscopy and attenuated total reflectance Fourier transform infrared spectroscopy combined with Pearson's product moment correlation coefficients and principal component analysis[J]., 2018, 63(1): 285-291.

    [12] Senior Samir et al. Characterization and dating of blue ballpoint pen inks using principal component analysis of UV-Vis absorption spectra, IR spectroscopy, and HPTLC[J]., 2012, 57(4): 1087-93.

    [13] ZHANG Jin, JIANG Hong, DUAN Bin, et al. A rapid and nondestructive approach for forensic identification of cigarette inner liner papers using shift-excitation Raman difference spectroscopy and chemometrics[J]., 2021, 66(6): 2180-2189.

    [14] Michal Oravec, Anel Beganovi?, Luká? Gál, et al. Forensic classification of black inkjet prints using Fourier transform near-infrared spectroscopy and linear discriminant analysis[J]., 2019, 299: 128-134.

    [15] Maione Camila, de Oliveira Souza Vanessa Cristina, Togni Loraine Rezende, et al. Using cluster analysis and ICP-MS to identify groups of ecstasy tablets in Sao Paulo State, Brazil[J]., 2017, 62(6): 1479-1486.

    [16] Nurlaila Ika, Irawati Wahyu, Purwandari Kartika, et al. K-means clustering model to discriminate copper-resistant bacteria as bioremediation agents[J]., 2021, 179: 804-812.

    [17] 余靜, 張云, 龐松穎, 等. 紅外光譜技術(shù)在物證鑒定中的應(yīng)用[J]. 光譜學(xué)與光譜分析, 2016, 36(9): 2807-2811.

    YU Jing, ZHANG Yun, PANG Songying. Application of IR in the field of evidence identification[J]., 2016, 36(9): 2807-2811.

    [18] 孫家政, 姜紅, 段斌, 等. 差分拉曼光譜結(jié)合多元統(tǒng)計學(xué)對IC卡的分類研究[J]. 塑料工業(yè), 2021, 49(4): 102-106.

    SUN Jiazheng, JIANG Hong, DUAN Bin, et al. Classification of IC cards by differential Raman spectroscopy combined with multivariate statistics[J]., 2021, 49(4): 102-106.

    [19] 何欣龍, 王繼芬, 張倩, 等. 基于多分類模型的記號筆墨水紅外光譜分析[J].化學(xué)通報, 2019, 82(2): 169-174.

    HE Xinlong, WANG Jifen, ZHANG Qian, et al. Infrared spectroscopy analysis of marker ink based on multi-classification model[J]., 2019, 82(2): 169-174.

    [20] 何欣龍, 王繼芬, 于佳裔, 等. 基于判別分析的車用保險杠激光拉曼光譜鑒別研究[J]. 激光雜志, 2019, 40(10): 21-25.

    HE Xinlong, WANG Jifen, YU Jiayi, et al. Identification Raman spectroscopy of vehicle bumpers based on discrimination analysis[J]., 2019, 40(10): 21-25.

    [21] 趙昱萱, 曾樂洋子, 李開開. 紅外光譜法結(jié)合化學(xué)計量學(xué)對可擦筆油墨檢驗(yàn)的研究[J]. 光譜學(xué)與光譜分析, 2021, 41(8): 2420-2426.

    ZHAO Yuxuan, ZENG Leyangzi, LI Kaikai. Identification of different brands erasable pens by infrared spectroscopy combined with chemometrics methods[J]., 2021, 41(8): 2420-2426.

    Infrared Spectroscopy Combined with Chemometrics to Test Ink Type of Straight Liquid Ballpoint Pen

    LIU Xinlei,HAN Danyan,XIAO Qiang,YANG Shangpeng

    (,110035,)

    The inspection and identification of writing ink are important in the field of forensic science. A Fourier transform infrared spectrometer was used to test 20 types of direct liquid ball pen ink samples, and chemometrics was used to rapidly test the direct liquid ball pen ink. The Fourier transform infrared spectrum data were standardized, and the spectrum was pre-processed using three methods: automatic baseline correction, peak area normalization, and Savitzky-Golay 5-point smoothing. The optimal value of classification K was determined using the sum of squares error(SSE). The samples were analyzed using K-means clustering, and the clustering results were explained. The principal component analysis method was used to verify the results of K-means clustering. The group mean equality test was used to investigate the contribution of principal component variables to the Fisher discriminant analysis (FDA) model, and the FDA discriminant model of straight liquid ball pen ink was constructed. The results show that all the ink samples were clustered into three categories using K-means clustering. The principal components analysis–Fisher discriminant analysis (PCA–FDA) model achieved 100% prediction and classification of different categories of straight-liquid ball pen inks with an accuracy of 100% after cross-validation. Infrared spectroscopy combined withthe PCA–FDA model can be used for rapid and accurate inspection and identification of direct liquid ball pen inks.

    ink type of straight liquid ballpoint pen, Fourier transform infrared spectroscopy(FTIR), chemometrics, discriminant analysis

    O657.33

    A

    1001-8891(2022)08-0882-06

    2022-01-11;

    2022-02-11.

    劉新磊(1998-),男,碩士研究生,主要從事刑事科學(xué)技術(shù)方面的工作,E-mail:492270626@qq.com。

    韓丹巖(1965-),女,教授,主要從事刑事技術(shù)方面的工作,E-mail:13066585988@163.com。

    中國刑事警察學(xué)院研究生創(chuàng)新能力提升項(xiàng)目資助(2021YCYB33)。

    猜你喜歡
    判別函數(shù)均值紅外
    網(wǎng)紅外賣
    閃亮的中國紅外『芯』
    金橋(2021年4期)2021-05-21 08:19:20
    游樂設(shè)施事故與危險量化判別函數(shù)的構(gòu)建
    TS系列紅外傳感器在嵌入式控制系統(tǒng)中的應(yīng)用
    電子制作(2019年7期)2019-04-25 13:17:14
    探究上市公司財務(wù)預(yù)警的數(shù)學(xué)模型
    基于快速遞推模糊2-劃分熵圖割的紅外圖像分割
    均值不等式失效時的解決方法
    均值與方差在生活中的應(yīng)用
    關(guān)于均值有界變差函數(shù)的重要不等式
    對偶均值積分的Marcus-Lopes不等式
    国产女主播在线喷水免费视频网站 | 亚洲图色成人| 久久国内精品自在自线图片| 有码 亚洲区| 免费不卡的大黄色大毛片视频在线观看 | 欧美变态另类bdsm刘玥| 夫妻性生交免费视频一级片| avwww免费| 精品久久国产蜜桃| 男插女下体视频免费在线播放| 91精品国产九色| 在线国产一区二区在线| 又粗又爽又猛毛片免费看| 亚州av有码| 国产精品一及| 国产精品人妻久久久影院| 亚洲七黄色美女视频| 国产日本99.免费观看| 免费黄网站久久成人精品| 青春草国产在线视频 | 日韩欧美国产在线观看| 久久久久免费精品人妻一区二区| 夫妻性生交免费视频一级片| 久久精品91蜜桃| 又爽又黄a免费视频| 国产黄色视频一区二区在线观看 | 国产美女午夜福利| 美女 人体艺术 gogo| 午夜激情欧美在线| 久久精品国产清高在天天线| 亚洲精品成人久久久久久| 精品国产三级普通话版| 97超视频在线观看视频| 亚州av有码| 97在线视频观看| 欧美成人一区二区免费高清观看| 乱码一卡2卡4卡精品| 网址你懂的国产日韩在线| 日韩精品青青久久久久久| 国产黄a三级三级三级人| 性色avwww在线观看| 欧洲精品卡2卡3卡4卡5卡区| 婷婷精品国产亚洲av| 精华霜和精华液先用哪个| 国产视频首页在线观看| 非洲黑人性xxxx精品又粗又长| 国产成人精品一,二区 | 久久精品久久久久久久性| 99精品在免费线老司机午夜| 精品人妻偷拍中文字幕| 日韩强制内射视频| 国产一级毛片在线| 波野结衣二区三区在线| 亚州av有码| 日韩高清综合在线| 99精品在免费线老司机午夜| av天堂在线播放| 国产精品日韩av在线免费观看| 在线免费十八禁| 中文字幕久久专区| 欧美+亚洲+日韩+国产| 天堂√8在线中文| 亚洲av一区综合| 日韩 亚洲 欧美在线| 久久久精品大字幕| 天天躁夜夜躁狠狠久久av| 有码 亚洲区| 国产精品一区www在线观看| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线 | 69人妻影院| 少妇人妻精品综合一区二区 | 久久欧美精品欧美久久欧美| 2022亚洲国产成人精品| 18禁黄网站禁片免费观看直播| 久久精品国产亚洲av香蕉五月| 91aial.com中文字幕在线观看| 国产v大片淫在线免费观看| 99久国产av精品| av黄色大香蕉| 国产精品不卡视频一区二区| 成年版毛片免费区| av在线蜜桃| 国产精品福利在线免费观看| 丰满的人妻完整版| 亚洲内射少妇av| 成人午夜高清在线视频| 色播亚洲综合网| 插阴视频在线观看视频| 美女 人体艺术 gogo| 一卡2卡三卡四卡精品乱码亚洲| 久久久a久久爽久久v久久| 看片在线看免费视频| 人妻制服诱惑在线中文字幕| 内射极品少妇av片p| 久久久久久久久久成人| 在线a可以看的网站| 男人舔奶头视频| 日日撸夜夜添| 亚洲精品国产成人久久av| 中国美女看黄片| 久久精品国产亚洲网站| 日本三级黄在线观看| 小蜜桃在线观看免费完整版高清| 婷婷精品国产亚洲av| 2021天堂中文幕一二区在线观| 长腿黑丝高跟| 高清毛片免费观看视频网站| 少妇熟女欧美另类| 亚洲不卡免费看| 亚洲国产精品国产精品| 亚洲精品影视一区二区三区av| 国产精品嫩草影院av在线观看| 少妇人妻精品综合一区二区 | 亚洲国产日韩欧美精品在线观看| 免费看光身美女| 成人特级黄色片久久久久久久| 欧美精品一区二区大全| 日韩一本色道免费dvd| 少妇人妻一区二区三区视频| 高清日韩中文字幕在线| 天美传媒精品一区二区| 简卡轻食公司| 国产精品永久免费网站| 欧美日韩国产亚洲二区| ponron亚洲| avwww免费| 欧美性感艳星| 午夜爱爱视频在线播放| 又爽又黄a免费视频| 久久久久久久久久成人| 亚洲精品自拍成人| 国产国拍精品亚洲av在线观看| 一本久久精品| 直男gayav资源| 国产伦在线观看视频一区| 成人毛片a级毛片在线播放| 国产精品人妻久久久影院| 午夜爱爱视频在线播放| 搞女人的毛片| 黑人高潮一二区| 日韩 亚洲 欧美在线| 亚洲性久久影院| 男女下面进入的视频免费午夜| 一区二区三区四区激情视频 | 久久精品影院6| 久久精品国产亚洲网站| 国产精品久久久久久久电影| 欧美bdsm另类| 欧美xxxx黑人xx丫x性爽| 日韩大尺度精品在线看网址| 18禁黄网站禁片免费观看直播| 午夜老司机福利剧场| 免费观看精品视频网站| 午夜a级毛片| 男插女下体视频免费在线播放| 欧美变态另类bdsm刘玥| 欧美日韩精品成人综合77777| 一区二区三区免费毛片| 在线国产一区二区在线| av免费在线看不卡| 久久久国产成人精品二区| 免费av不卡在线播放| 亚洲av成人av| 欧美极品一区二区三区四区| 91午夜精品亚洲一区二区三区| 变态另类丝袜制服| 欧美色欧美亚洲另类二区| 国产精品国产高清国产av| 国产精品免费一区二区三区在线| 免费看av在线观看网站| 三级国产精品欧美在线观看| 亚洲国产精品成人综合色| 国产高清三级在线| 午夜福利在线在线| 美女内射精品一级片tv| 国产亚洲精品av在线| 午夜a级毛片| 亚洲av免费高清在线观看| 精品久久久噜噜| 亚洲激情五月婷婷啪啪| 国模一区二区三区四区视频| 久久久久久久久久成人| av又黄又爽大尺度在线免费看 | 直男gayav资源| 国产高潮美女av| 在线观看美女被高潮喷水网站| 亚洲天堂国产精品一区在线| 久久99蜜桃精品久久| 成人午夜高清在线视频| 精品人妻熟女av久视频| 国产精品一二三区在线看| 麻豆成人av视频| 99久久中文字幕三级久久日本| 精品无人区乱码1区二区| 免费无遮挡裸体视频| 国产av不卡久久| 亚洲精品456在线播放app| 国产亚洲精品久久久久久毛片| a级毛片免费高清观看在线播放| 午夜福利成人在线免费观看| 久久欧美精品欧美久久欧美| 日本欧美国产在线视频| 日本熟妇午夜| 国产av不卡久久| 精品日产1卡2卡| 国产成人精品久久久久久| 综合色丁香网| 日韩一区二区视频免费看| 久久精品国产自在天天线| 国产真实乱freesex| 看免费成人av毛片| 亚洲内射少妇av| 欧美极品一区二区三区四区| 2022亚洲国产成人精品| 欧美色视频一区免费| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品三级大全| 级片在线观看| 在线观看一区二区三区| 麻豆精品久久久久久蜜桃| 国产老妇女一区| 一级二级三级毛片免费看| 国产大屁股一区二区在线视频| 神马国产精品三级电影在线观看| 免费观看人在逋| 国产精品久久电影中文字幕| av天堂在线播放| 国产精华一区二区三区| 亚洲色图av天堂| 久久久久久久午夜电影| 天堂影院成人在线观看| 久久99精品国语久久久| 亚洲欧美中文字幕日韩二区| 亚洲电影在线观看av| 免费观看人在逋| 国产成人精品一,二区 | 成人鲁丝片一二三区免费| 欧美区成人在线视频| 18禁在线播放成人免费| 亚洲av中文字字幕乱码综合| 国产成年人精品一区二区| 国产v大片淫在线免费观看| 91久久精品电影网| 毛片女人毛片| 免费av不卡在线播放| 青春草亚洲视频在线观看| 蜜桃久久精品国产亚洲av| 夫妻性生交免费视频一级片| 日本av手机在线免费观看| 欧美最黄视频在线播放免费| 久久中文看片网| 成人鲁丝片一二三区免费| 一区二区三区高清视频在线| 午夜免费激情av| 少妇熟女aⅴ在线视频| 99在线人妻在线中文字幕| 99热这里只有精品一区| 91av网一区二区| 最近2019中文字幕mv第一页| 插逼视频在线观看| 男女视频在线观看网站免费| 日韩,欧美,国产一区二区三区 | 午夜激情福利司机影院| ponron亚洲| 亚洲va在线va天堂va国产| 欧美激情国产日韩精品一区| 青春草亚洲视频在线观看| 日韩大尺度精品在线看网址| 99热网站在线观看| 日本三级黄在线观看| 黄片无遮挡物在线观看| 中文字幕精品亚洲无线码一区| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| 欧美潮喷喷水| АⅤ资源中文在线天堂| 国产成人a区在线观看| 男女下面进入的视频免费午夜| 人妻系列 视频| 日本-黄色视频高清免费观看| 久久韩国三级中文字幕| 成人高潮视频无遮挡免费网站| 国产亚洲91精品色在线| 亚洲在线观看片| 热99在线观看视频| 中文字幕免费在线视频6| 自拍偷自拍亚洲精品老妇| 99在线视频只有这里精品首页| 亚洲经典国产精华液单| 免费搜索国产男女视频| 人妻制服诱惑在线中文字幕| 精品免费久久久久久久清纯| 国产伦在线观看视频一区| 国产精品人妻久久久影院| 99久国产av精品| 亚洲电影在线观看av| 免费观看人在逋| 赤兔流量卡办理| 黑人高潮一二区| 国产av一区在线观看免费| 亚洲精品成人久久久久久| 欧美色欧美亚洲另类二区| 男的添女的下面高潮视频| 午夜福利视频1000在线观看| 国产亚洲精品久久久久久毛片| 91久久精品国产一区二区成人| 成人欧美大片| 看非洲黑人一级黄片| 国产成人aa在线观看| 国产极品精品免费视频能看的| 亚洲精品影视一区二区三区av| 伊人久久精品亚洲午夜| 日本av手机在线免费观看| 97超碰精品成人国产| 国产精品三级大全| 看片在线看免费视频| 欧美精品国产亚洲| 成人美女网站在线观看视频| 99热只有精品国产| 又黄又爽又刺激的免费视频.| 久久亚洲精品不卡| 亚洲国产精品久久男人天堂| 久久精品国产亚洲av香蕉五月| 欧美成人精品欧美一级黄| 2021天堂中文幕一二区在线观| 午夜激情欧美在线| 99国产极品粉嫩在线观看| 最近2019中文字幕mv第一页| 亚洲精品亚洲一区二区| 久久人人爽人人片av| 久久久久久九九精品二区国产| 禁无遮挡网站| 成人亚洲欧美一区二区av| 女人十人毛片免费观看3o分钟| 内射极品少妇av片p| av在线天堂中文字幕| 在线播放国产精品三级| 色5月婷婷丁香| 久久久久久久久久久免费av| 国产精品福利在线免费观看| 日韩av在线大香蕉| 精品久久久久久久久久久久久| 亚洲精品日韩在线中文字幕 | 黄色欧美视频在线观看| 日本免费一区二区三区高清不卡| 噜噜噜噜噜久久久久久91| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 免费大片18禁| 丝袜喷水一区| 村上凉子中文字幕在线| 精品久久久久久久末码| 日韩一区二区视频免费看| 看非洲黑人一级黄片| 麻豆国产av国片精品| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品色激情综合| 免费看av在线观看网站| 边亲边吃奶的免费视频| 日日干狠狠操夜夜爽| 免费一级毛片在线播放高清视频| 变态另类丝袜制服| 久久久欧美国产精品| or卡值多少钱| 久久精品国产亚洲av涩爱 | 国产老妇女一区| 午夜老司机福利剧场| 亚洲成人av在线免费| 日韩欧美精品免费久久| 成年女人看的毛片在线观看| 麻豆国产av国片精品| 国产女主播在线喷水免费视频网站 | 丰满的人妻完整版| 久久久久久久亚洲中文字幕| 国产国拍精品亚洲av在线观看| 成人一区二区视频在线观看| 久久精品综合一区二区三区| 国产日韩欧美在线精品| 日本黄大片高清| 国产亚洲精品久久久久久毛片| 五月玫瑰六月丁香| 久99久视频精品免费| 亚洲成av人片在线播放无| 99久久中文字幕三级久久日本| 国产成人aa在线观看| 小说图片视频综合网站| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 久久精品久久久久久噜噜老黄 | 国产精品嫩草影院av在线观看| 欧美+日韩+精品| 久久久久久久久久黄片| 晚上一个人看的免费电影| 天天躁夜夜躁狠狠久久av| 能在线免费看毛片的网站| 亚洲在久久综合| 亚洲内射少妇av| 级片在线观看| 五月伊人婷婷丁香| 麻豆乱淫一区二区| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 精品国产三级普通话版| 国产高清三级在线| 国国产精品蜜臀av免费| 亚洲性久久影院| 69人妻影院| 国产白丝娇喘喷水9色精品| 久久午夜福利片| av福利片在线观看| 日韩亚洲欧美综合| 蜜臀久久99精品久久宅男| 女人被狂操c到高潮| 美女大奶头视频| 午夜精品国产一区二区电影 | 亚洲av不卡在线观看| 精品国产三级普通话版| 久久精品国产鲁丝片午夜精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品久久久久久亚洲av鲁大| 国产视频内射| 一区福利在线观看| 国产在线精品亚洲第一网站| 国产私拍福利视频在线观看| 美女高潮的动态| 麻豆久久精品国产亚洲av| 男插女下体视频免费在线播放| 在线a可以看的网站| 亚洲四区av| 在线观看av片永久免费下载| 成人二区视频| 精品一区二区免费观看| 成人一区二区视频在线观看| 我要搜黄色片| 国产精品.久久久| 亚洲欧美清纯卡通| 久久精品人妻少妇| 国产精品蜜桃在线观看 | 国产亚洲欧美98| 色噜噜av男人的天堂激情| 91精品国产九色| 草草在线视频免费看| 精品久久久久久久久久免费视频| 亚洲无线在线观看| 黄色一级大片看看| 男插女下体视频免费在线播放| 一进一出抽搐gif免费好疼| 久久人人爽人人片av| 97人妻精品一区二区三区麻豆| 国内久久婷婷六月综合欲色啪| 深夜精品福利| 日韩成人伦理影院| 嫩草影院新地址| 欧美日韩综合久久久久久| 久久久国产成人免费| 最好的美女福利视频网| 久久精品国产亚洲av天美| 乱码一卡2卡4卡精品| 成人毛片60女人毛片免费| 插阴视频在线观看视频| 白带黄色成豆腐渣| 狠狠狠狠99中文字幕| kizo精华| 久久久久久国产a免费观看| 国产综合懂色| 国产伦精品一区二区三区视频9| 国产精品蜜桃在线观看 | 一区二区三区四区激情视频 | 一夜夜www| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 欧美变态另类bdsm刘玥| 欧美性猛交黑人性爽| 亚洲丝袜综合中文字幕| 国产精品三级大全| 欧美激情国产日韩精品一区| 成人无遮挡网站| 亚洲一区二区三区色噜噜| 美女 人体艺术 gogo| 久久精品国产自在天天线| 国产久久久一区二区三区| 99久久中文字幕三级久久日本| 联通29元200g的流量卡| 最近中文字幕高清免费大全6| 欧美区成人在线视频| 最近2019中文字幕mv第一页| 26uuu在线亚洲综合色| 亚洲国产精品久久男人天堂| 男女边吃奶边做爰视频| 黄色一级大片看看| 国产老妇伦熟女老妇高清| 久久99蜜桃精品久久| 高清在线视频一区二区三区 | 亚洲aⅴ乱码一区二区在线播放| 精品熟女少妇av免费看| 亚洲一区高清亚洲精品| 亚洲一级一片aⅴ在线观看| 小蜜桃在线观看免费完整版高清| 久久韩国三级中文字幕| 国产爱豆传媒在线观看| 亚洲av电影不卡..在线观看| 99在线人妻在线中文字幕| 国产精品99久久久久久久久| 午夜福利视频1000在线观看| 村上凉子中文字幕在线| 国产高清三级在线| 青青草视频在线视频观看| 99视频精品全部免费 在线| 99热6这里只有精品| 婷婷六月久久综合丁香| 午夜精品一区二区三区免费看| 变态另类丝袜制服| 乱码一卡2卡4卡精品| 国产成年人精品一区二区| 国产欧美日韩精品一区二区| 日日摸夜夜添夜夜爱| av在线蜜桃| 蜜桃亚洲精品一区二区三区| 国产成年人精品一区二区| 麻豆成人av视频| 亚洲欧美清纯卡通| 亚洲自拍偷在线| 国产色爽女视频免费观看| 18禁在线播放成人免费| av.在线天堂| 免费看av在线观看网站| 亚洲自偷自拍三级| 国产成人freesex在线| 男人和女人高潮做爰伦理| 村上凉子中文字幕在线| 国产高清三级在线| 搡老妇女老女人老熟妇| 婷婷亚洲欧美| 一级毛片我不卡| 我的老师免费观看完整版| 中出人妻视频一区二区| 九色成人免费人妻av| 国产毛片a区久久久久| 国产片特级美女逼逼视频| 黄片wwwwww| 国产av麻豆久久久久久久| 毛片一级片免费看久久久久| 日韩精品青青久久久久久| 亚洲第一电影网av| 欧美性感艳星| 不卡一级毛片| 欧美另类亚洲清纯唯美| 丰满乱子伦码专区| 美女脱内裤让男人舔精品视频 | 寂寞人妻少妇视频99o| 亚洲成a人片在线一区二区| 日韩欧美三级三区| 天天一区二区日本电影三级| 午夜爱爱视频在线播放| 国产伦精品一区二区三区视频9| 人人妻人人看人人澡| 欧美成人一区二区免费高清观看| 2021天堂中文幕一二区在线观| 男女那种视频在线观看| 狠狠狠狠99中文字幕| 久久久久久九九精品二区国产| 亚洲国产精品久久男人天堂| 我要看日韩黄色一级片| 久久久久久久久久久免费av| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 麻豆精品久久久久久蜜桃| 黄片无遮挡物在线观看| 美女 人体艺术 gogo| 91aial.com中文字幕在线观看| 夜夜夜夜夜久久久久| 久久国产乱子免费精品| 亚洲美女搞黄在线观看| 亚洲精品色激情综合| 天堂av国产一区二区熟女人妻| 国产精品精品国产色婷婷| 国产精品国产高清国产av| 高清毛片免费观看视频网站| 日本免费一区二区三区高清不卡| 久久精品久久久久久久性| 成人国产麻豆网| 老女人水多毛片| 床上黄色一级片| 国产伦一二天堂av在线观看| 亚洲aⅴ乱码一区二区在线播放| 99视频精品全部免费 在线| 国产精品久久久久久精品电影小说 | 久久人人爽人人片av| 亚洲精品成人久久久久久| 亚洲天堂国产精品一区在线| 有码 亚洲区| 日韩国内少妇激情av| 午夜福利视频1000在线观看| 又粗又爽又猛毛片免费看| 如何舔出高潮| 久久久久久久久大av| 免费看光身美女| 欧美xxxx性猛交bbbb| 亚洲一级一片aⅴ在线观看| 久久久久国产网址| 成人性生交大片免费视频hd| 久久久久久伊人网av| 伊人久久精品亚洲午夜| 成人性生交大片免费视频hd| 久久久久久久久大av| 亚洲最大成人中文| 国产中年淑女户外野战色| 国产精品人妻久久久影院| 三级经典国产精品| 欧美丝袜亚洲另类| 免费不卡的大黄色大毛片视频在线观看 | 免费看日本二区| 欧美一区二区国产精品久久精品| 成人特级黄色片久久久久久久|