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

    基于傅里葉變換紅外光譜的金銀花產(chǎn)地鑒別

    2018-10-30 10:14皇甫義靜楊春艷毛貽龍劉美楊洋
    湖北農(nóng)業(yè)科學(xué) 2018年14期
    關(guān)鍵詞:微分金銀花二階

    皇甫義靜 楊春艷 毛貽龍 劉美 楊洋

    摘要:基于傅里葉變換紅外光譜技術(shù),結(jié)合光譜檢索的方法進(jìn)行金銀花(Lonicera japonica Thunb.)產(chǎn)地鑒別研究。測(cè)試了4個(gè)省、6個(gè)種植區(qū)共130份金銀花樣本的紅外光譜,利用Omnic 8.2軟件建立了由各類樣品的平均紅外光譜、一階導(dǎo)數(shù)平均譜和二階導(dǎo)數(shù)平均譜組成的光譜庫(kù)。比較了基于不同算法、不同波數(shù)范圍和不同類型光譜數(shù)據(jù)建立的模型檢索鑒別結(jié)果表明,基于1 500~1 200 cm-1波數(shù)范圍的二階導(dǎo)數(shù)光譜數(shù)據(jù)的絕對(duì)微分差算法檢索的鑒別效果相對(duì)較好。有119個(gè)樣品的光譜檢索匹配得分最高值對(duì)映的光譜為自身所在類別,占樣品總數(shù)的91.5%。紅外光譜結(jié)合光譜檢索的方法是金銀花產(chǎn)地鑒別的一種簡(jiǎn)便、快捷方法。

    關(guān)鍵詞:金銀花(Lonicera japonica Thunb.);傅里葉變換紅外光譜;光譜檢索;產(chǎn)地;鑒別

    中圖分類號(hào):R284 文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):0439-8114(2018)14-0105-04

    DOI:10.14088/j.cnki.issn0439-8114.2018.14.025 開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Abstract: Fourier transform infrared Spectroscopy (FTIR) was used to identify the origin of honeysuckle. The infrared spectra of 130 samples from 6 growing areas in 4 provinces were collected, and the first-derivate spectra and second-derivate spectra were taken by the software Omnic 8.2,and 3 spectra libraries were constructed. The correlation search, square differential search, absolute differential search and expert search of the spectra, the first-derivate spectra and second-derivate spectra for all samples are performed with corresponding library in different spectrum range, severally. The recognition effect were compared, and found that the model based on absolute differential reach of second-derivate spectra in the specified range of 1 500~1 200 cm-1 get a better identification effect, with the correct rate of 91.5%. The results showed that FTIR combined with the spectral retrieval method can be an effective and easy way for distinguishing the origin of honeysuckle herb.

    Key words: Lonicera japonica Thunb.; fourier transform infrared spectroscopy; spectra retrieval; honeysuckle; origin; identification

    金銀花(Lonicera japonica Thunb.)為忍冬科忍冬屬植物忍冬的干燥花蕾,屬常用藥材,既能宣散風(fēng)熱,又具有清解血毒、殺菌消炎等功效[1],醫(yī)學(xué)上廣泛應(yīng)用于急性熱病及外科感染性疾病,如上呼吸道感染、急性扁桃體炎、咽炎、癤瘡?fù)茨[、乳腺炎等,是眾多藥品,如金銀花顆粒、銀黃片/膠囊、銀翹解毒片/顆粒等的主要成分[2]。金銀花主要化學(xué)成分為揮發(fā)油、環(huán)烯醚萜類、黃酮類、三萜及三萜皂苷類、有機(jī)酸類及無機(jī)元素等化合物[3]。忍冬在全國(guó)各地均有分布和種植,而各地氣候、土壤等生態(tài)、地理?xiàng)l件的差異,導(dǎo)致不同產(chǎn)區(qū)金銀花質(zhì)量有明顯差異[4-7],傳統(tǒng)上以河南新密、封丘,山東平邑、費(fèi)縣為道地產(chǎn)區(qū)。在銷售過程中,有的商家以次充好,將非道地產(chǎn)區(qū)的金銀花作為道地藥材出售,損害消費(fèi)者利益,因此對(duì)金銀花產(chǎn)地的鑒別具有重要意義。紅外光譜具有快速、簡(jiǎn)便,且能全面反映樣品化學(xué)成分信息的優(yōu)點(diǎn),符合中藥整體使用的原則[8],在中藥的鑒別和質(zhì)量分析方面得到廣泛應(yīng)用[9-12]。本研究采用傅里葉變換紅外光譜結(jié)合光譜檢索對(duì)7個(gè)產(chǎn)地的金銀花進(jìn)行研究,以期為藥材市場(chǎng)金銀花質(zhì)量鑒別與控制提供依據(jù)。

    1 材料與方法

    1.1 材料與試劑

    試驗(yàn)樣品于2016年7月采自山東省平邑縣和費(fèi)縣,河南省原陽(yáng)縣和封丘縣,云南省瀘水縣,山西省汾西縣共6個(gè)產(chǎn)區(qū),均為栽培品種。采集時(shí),同一產(chǎn)區(qū)盡量選擇花蕾大小均勻的樣本,樣品詳細(xì)信息見表1。

    試劑為KBr,分析純,購(gòu)于天津市風(fēng)船化學(xué)試劑科技有限公司。

    1.2 主要儀器

    試驗(yàn)所用紅外光譜儀為Perkin-Elmer Frontier型傅里葉變換紅外光譜儀,裝配DTGS檢測(cè)儀,光譜掃描范圍為4 000~400 cm-1,分辨率為4 cm-1,累計(jì)掃描次數(shù)16次。

    1.3 樣品紅外光譜的采集

    光譜采集中,室溫保持在22 ℃,空氣濕度保持在40%~45%。測(cè)試時(shí),將每份樣品的花蕾剪下,在45 ℃恒溫下烘干至恒重,將樣品在瑪瑙研缽中磨為細(xì)粉后按樣品與溴化鉀1∶70的質(zhì)量加入溴化鉀攪磨均勻,壓制成透明薄片,放入光譜儀中采集光譜。

    1.4 光譜預(yù)處理和試驗(yàn)方法

    采集所得光譜均已扣除背景,用Nicolet Omnic 8.2 軟件進(jìn)行吸光度、自動(dòng)基線校正、自動(dòng)平滑和縱坐標(biāo)歸一化等預(yù)處理,計(jì)算了每份樣品紅外光譜的一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜以及每類樣品前10份樣本的平均紅外光譜、一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜。在Omnic 8.2軟件中分別建立了由每類樣品前10份樣本的平均紅外光譜、一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜組成的Lib1、Lib2和Lib3光譜數(shù)據(jù)庫(kù)。將所有樣品的紅外光譜、一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜分別在相應(yīng)光譜庫(kù)中進(jìn)行檢索,軟件自動(dòng)得出匹配結(jié)果。

    2 結(jié)果與分析

    2.1 不同產(chǎn)地金銀花樣品的紅外光譜特征分析

    圖1是不同產(chǎn)地金銀花樣品的平均紅外光譜,表2是不同產(chǎn)地金銀花的紅外光譜主要吸收帶。不同產(chǎn)地金銀花的紅外光譜非常相似,且特征吸收峰的峰位相近。在4 000~700 cm-1波數(shù)范圍內(nèi)共呈現(xiàn)了18個(gè)吸收峰,其中3 360 cm-1附近強(qiáng)而寬的吸收峰來自O(shè)-H伸縮振動(dòng);2 925和2 855 cm-1附近吸收峰來自CH2的反對(duì)稱伸縮振動(dòng)和對(duì)稱伸縮振動(dòng)吸收;1 730 cm-1附近吸收峰來自C=O伸縮振動(dòng),1 640 cm-1附近強(qiáng)的吸收峰主要來自烯醚雙鍵伸縮振動(dòng)[13];1 540 cm-1附近吸收峰來自C=C伸縮振動(dòng)吸收;1 440 cm-1附近吸收峰歸屬C-H彎曲振動(dòng);1 400和1 380 cm-1附近吸收峰歸屬中C-H剪式彎曲振動(dòng),1 250 cm-1附近吸收峰來自C-OH伸縮振動(dòng)[14];1 150、1 070、1 050 cm-1附近階梯增強(qiáng)的吸收峰歸屬多糖中醇、酯基、醚基的C-O伸縮振動(dòng)[15];930 cm-1附近吸收峰來自苯環(huán)上C-H面外彎曲振動(dòng)[16]。

    雖然不同產(chǎn)地金銀花樣品的紅外光譜非常相似,但吸收峰的位置、相對(duì)峰強(qiáng)和形狀仍稍有差異,具體表現(xiàn)在:①只有山東省費(fèi)縣樣品的光譜在1 730 cm-1附近C=O伸縮振動(dòng)吸收峰處顯示了弱的吸收峰,且在1 070 cm-1附近未呈現(xiàn)明顯吸收峰;②在1 640 cm-1附近烯醚雙鍵伸縮振動(dòng)吸收峰處,雖然6個(gè)產(chǎn)地樣品都呈現(xiàn)了吸收峰,但相對(duì)峰強(qiáng)和峰位差異明顯。封丘縣和原陽(yáng)縣產(chǎn)的金銀花光譜中呈現(xiàn)了極強(qiáng)的吸收峰,且峰位向短波方向移動(dòng),平邑縣、費(fèi)縣和汾西縣樣品的光譜中呈現(xiàn)了強(qiáng)吸收峰,云南產(chǎn)金銀花的紅外光譜中呈現(xiàn)了中等強(qiáng)度的吸收峰,且峰位向長(zhǎng)波方向移動(dòng);③在1 540 cm-1附近C=C伸縮振動(dòng)吸收峰處和1 440 cm-1附近C-H彎曲振動(dòng)吸收峰處,只有云南省瀘水縣樣品的光譜未顯示吸收峰;④在1 380和705 cm-1附近的吸收峰處,只有封丘縣和原陽(yáng)縣的樣品未顯示吸收峰,且在1 150 cm-1附近C-O伸縮振動(dòng)吸收峰處,只有原陽(yáng)縣的樣品未顯示吸收峰;⑤在930 cm-1附近吸收峰處,只有平邑縣和費(fèi)縣樣品顯示了弱吸收峰。由以上分析,可將不同產(chǎn)地的樣品初略區(qū)分開來。

    2.2 光譜檢索對(duì)金銀花產(chǎn)地的鑒別

    在Omnic 8.2軟件的光譜檢索模塊中,設(shè)置有絕對(duì)差、相關(guān)性、平方差、平方微分差、絕對(duì)微分差和專家檢索等6種算法。本研究選擇相關(guān)性、平方微分差、絕對(duì)微分差和專家檢索算法進(jìn)行光譜檢索鑒別比較[17]。由于不同產(chǎn)地金銀花樣品的紅外光譜非常相似,為了放大樣品光譜的差異,通過Omnic 8.2軟件計(jì)算了各樣品光譜的一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜;同時(shí)為了更客觀地進(jìn)行鑒別,分別取各產(chǎn)地前10份樣品紅外光譜的平均譜、一階導(dǎo)數(shù)光譜的平均譜和二階導(dǎo)數(shù)光譜的平均譜依次組成光譜庫(kù)Lib1、Lib2和Lib3。視全部130份樣品為未知產(chǎn)地的樣品,將其光譜分別在相應(yīng)光譜庫(kù)中進(jìn)行檢索匹配。比較了基于4 000~400 cm-1、1 800~1 200 cm-1和1 500~1 200 cm-1波數(shù)范圍的檢索鑒別分類效果。

    2.2.1 不同算法和不同類型光譜的檢索鑒別分析 表3為未知產(chǎn)地樣品的紅外光譜與相應(yīng)光譜庫(kù)進(jìn)行不同算法的檢索鑒別分析結(jié)果。將全部130份金銀花樣品的紅外光譜分別在光譜庫(kù)Lib1中進(jìn)行全譜范圍的相關(guān)性、平方微分差、絕對(duì)微分差和專家檢索。結(jié)果表明,大部分未知分類樣品的紅外光譜在光譜庫(kù)Lib1中檢索匹配得分最高值對(duì)應(yīng)的光譜均為自身所在類別,4種算法的檢索正確率均達(dá)85%以上,其中絕對(duì)微分差算法的檢索正確率最高,達(dá)88.5%,說明4種算法都可以鑒別不同產(chǎn)地的金銀花,其中基于絕對(duì)微分差算法建立的模型更有利于不同產(chǎn)地金銀花樣品的鑒別。

    將全部樣品的紅外光譜、一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜分別在Lib1、Lib2和Lib3光譜庫(kù)中進(jìn)行全譜范圍的絕對(duì)微分差算法的光譜檢索。結(jié)果表明,紅外光譜和導(dǎo)數(shù)光譜與相應(yīng)光譜庫(kù)的檢索匹配正確率均較高,未知產(chǎn)地樣品的紅外光譜和一階導(dǎo)數(shù)光譜分別與光譜庫(kù)Lib1和Lib2的檢索,匹配正確的樣品數(shù)均有114個(gè),正確率為87.7%;未知產(chǎn)地樣品的二階導(dǎo)數(shù)光譜與光譜庫(kù)Lib3的檢索,匹配正確的樣品數(shù)有115個(gè),正確率為88.5%,這與導(dǎo)數(shù)光譜具有放大光譜差異的特征相吻合。因此3種類型光譜數(shù)據(jù)均可用于鑒別不同產(chǎn)地的金銀花,其中基于二階導(dǎo)數(shù)光譜的檢索鑒別效果相對(duì)更好。

    2.2.2 不同范圍光譜數(shù)據(jù)的檢索鑒別分析 將全部未知產(chǎn)地樣品的二階導(dǎo)數(shù)光譜與光譜庫(kù)Lib3分別在4 000~400 cm-1、1 800~1 200 cm-1和1 500~1 200 cm-1波數(shù)范圍內(nèi)進(jìn)行絕對(duì)微分差算法檢索,結(jié)果見表4。表4中僅列出鑒別分類不正確的樣品及匹配得分高的前3個(gè)類別和得分,3個(gè)波數(shù)范圍均鑒別正確的108份樣品的檢索信息未列出。

    在基于4 000~400 cm-1波數(shù)范圍光譜數(shù)據(jù)的檢索中,有115個(gè)樣品的光譜檢索匹配得分最高值對(duì)映的光譜為自身所在類別,占樣品總數(shù)的88.5%;在1 800~1 200 cm-1波數(shù)范圍的檢索中,有113個(gè)樣品的光譜檢索匹配得分最高值對(duì)映的光譜為自身所在類別產(chǎn)地,占樣品總數(shù)的86.9%;在1 500~1 200 cm-1波數(shù)范圍的檢索中,有119個(gè)樣品的光譜檢索匹配得分最高值對(duì)映的光譜為自身所在類別產(chǎn)地,占樣品總數(shù)的91.5%。說明基于3個(gè)波數(shù)范圍的光譜數(shù)據(jù)均可以鑒別不同產(chǎn)地的金銀花,其中1 500~1 200 cm-1波數(shù)范圍的光譜數(shù)據(jù)鑒別效果更好。

    比較分析未知產(chǎn)地樣品的二階導(dǎo)數(shù)光譜與光譜庫(kù)Lib3進(jìn)行絕對(duì)微分差算法檢索的結(jié)果,主要有以下特點(diǎn):①在基于4 000~400 cm-1、1 800~1 200 cm-1和1 500~1 200 cm-1范圍的光譜數(shù)據(jù)的檢索中,瀘水縣和原陽(yáng)縣金銀花樣品產(chǎn)地的匹配正確率均為100%。②基于4 000~400 cm-1范圍光譜數(shù)據(jù)的檢索中,鑒別分類錯(cuò)誤的樣品C8、C9和C19,在1 800~1 200 cm-1和1 500~1 200 cm-1范圍的檢索中均歸類正確;基于4 000~400 cm-1和1 800~1 200 cm-1范圍光譜數(shù)據(jù)的檢索中,鑒別分類錯(cuò)誤的樣品A16、C11、D23、E1和E9在1 500~1 200 cm-1范圍的檢索中均歸類正確。③基于1 500~1 200 cm-1范圍光譜數(shù)據(jù)的檢索匹配得分均高于基于4 000~400 cm-1和1 800~1 200 cm-1范圍的匹配得分,其中基于4 000~400 cm-1范圍的匹配得分最低。④匹配分類錯(cuò)誤的情況主要集中在封丘縣、平邑縣、費(fèi)縣和原陽(yáng)縣,如在基于1 500~1 200 cm-1范圍的光譜數(shù)據(jù)的檢索中,鑒別錯(cuò)誤的A21、A26、A28、A29、A34、A38、C18、C26、D13、D14和D17共11份樣品均采集自封丘縣、平邑縣和費(fèi)縣,鑒別結(jié)果中除樣品C26的產(chǎn)地被鑒別為汾西縣外,其余10份均被鑒別為平邑縣、封丘縣或原陽(yáng)縣。

    3 結(jié)論與討論

    基于傅里葉變換紅外光譜技術(shù),利用光譜檢索法對(duì)6個(gè)種植區(qū)的130份金銀花樣品進(jìn)行鑒別研究,利用Omnic 8.2軟件對(duì)所有測(cè)試到的光譜進(jìn)行自動(dòng)基線校正等預(yù)處理,計(jì)算了各樣品紅外光譜的一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜,建立了分別由各產(chǎn)地樣品的紅外光譜、一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜的平均光譜組成的光譜庫(kù);光譜檢索中選擇相關(guān)性、平方微分差、絕對(duì)微分差和專家檢索4種算法,以及4 000~400 cm-1、1 800~1 200 cm-1和1 500~1 200 cm-1范圍的光譜數(shù)據(jù)計(jì)算所有樣品的匹配得分。結(jié)果顯示,各產(chǎn)地金銀花樣品的紅外光譜整體非常相似,通過吸收峰比對(duì),不同產(chǎn)地的樣品在幾個(gè)特殊吸收峰處具有差異,如只有山東省費(fèi)縣樣品的光譜在1 730 cm-1附近顯示了弱的吸收峰,且在1 070 cm-1附近未顯示明顯吸收峰;相關(guān)性、平方微分差、絕對(duì)微分差和專家檢索法均能鑒別樣品的產(chǎn)地,絕對(duì)微分差檢索算法對(duì)樣品的鑒別效果較好;比較了基于紅外光譜、一階導(dǎo)數(shù)光譜和二階導(dǎo)數(shù)光譜以及3個(gè)不同波數(shù)范圍的鑒別效果,結(jié)果表明基于1 500~1 200 cm-1范圍的二階導(dǎo)數(shù)光譜數(shù)據(jù)的絕對(duì)微分差法對(duì)金銀花樣品產(chǎn)地的鑒別效果最優(yōu),鑒別正確率達(dá)91.5%。

    本研究中瀘水縣樣品的檢索正確率達(dá)100%,且在其他產(chǎn)地所有樣品的檢索中,瀘水類別均未在匹配得分前五中出現(xiàn),這與瀘水縣的地理位置和氣候等生態(tài)環(huán)境與其他5個(gè)產(chǎn)地有較大差異的特征相一致。匹配分類錯(cuò)誤的情況主要集中在封丘縣、平邑縣、費(fèi)縣和原陽(yáng)縣,這與這4個(gè)產(chǎn)地有著相似的地理背景、降雨和氣候條件等特征相一致。雖然封丘縣和原陽(yáng)縣屬于河南省,平邑縣和費(fèi)縣屬于山東省,但它們都是金銀花的道地產(chǎn)區(qū),屬于同一生態(tài)區(qū)劃,所以金銀花的品質(zhì)特征差異甚小。若將封丘縣、平邑縣、費(fèi)縣和原陽(yáng)縣四個(gè)產(chǎn)地的樣品作為道地樣品歸為一類,則基于1 500~1 200 cm-1范圍的光譜數(shù)據(jù)檢索中,除C26外其余樣品的鑒別結(jié)果均與其實(shí)際分類相吻合,正確率將達(dá)99.2%。

    參考文獻(xiàn):

    [1] 國(guó)家藥典委員會(huì).中華人民共和國(guó)藥典(一部)[M].北京:中國(guó)藥典科技出版社,2015.221.

    [2] 王天志,李永梅.金銀花的研究進(jìn)展[J].華西藥學(xué)雜志,2000, 15(4):292-298.

    [3] 張小娜,童 杰,周衍晶,等.忍冬屬藥材藥效成分及藥理作用研究進(jìn)展[J].中國(guó)藥理學(xué)通報(bào),2014,30(8):1049-1054.

    [4] 王玲娜,張永清.不同種植區(qū)域忍冬植株形態(tài)特征比較[J].四川農(nóng)業(yè)大學(xué)學(xué)報(bào),2014,32(4):393-397.

    [5] 張重義,李 萍,李會(huì)軍,等.道地與非道地產(chǎn)區(qū)金銀花質(zhì)量的比較[J].中國(guó)中藥雜志,2007,32(9):786-788.

    [6] 王玲娜,劉紅燕,李 佳,等.不同種植點(diǎn)忍冬生長(zhǎng)形狀的變異及其與地理氣候因子的相關(guān)性分析[J].植物資源與環(huán)境學(xué)報(bào),2015,24(1):36-41.

    [7] 辛立紅,呂玲霞,管仁煒,等.金銀花質(zhì)量的印象因素及控制措施[J].山東農(nóng)業(yè)科學(xué),2017,49(3):154-157.

    [8] 孫素琴.中藥紅外光譜分析的關(guān)鍵方法和技術(shù)[C].全國(guó)第20屆光譜儀器與分析監(jiān)測(cè)學(xué)術(shù)研討會(huì),2013.48-52.

    [9] 包華音.壁虎藥材紅外光譜指紋圖譜研究[J].光散射學(xué)報(bào),2015, 27(3):276-280.

    [10] 黃冬蘭,陳小康,徐永群.溪黃草四種基源物質(zhì)的二維相關(guān)紅外光譜鑒別研究[J].光散射學(xué)報(bào),2015,27(1):64-68.

    [11] 申云霞,趙艷麗,張 霽,等.林藥復(fù)合種植滇龍膽紅外光譜鑒別研究[J].光譜學(xué)與光譜分析,2016,36(8):2468-2473.

    [12] GUO Y Z,BEI R,WANG J J,et al. Analysis of Chuanxiong Rhizoma and its active components by Fourier transform infrared spectroscopy combined with two-dimensional correlation infrared spectroscopy[J].Spectrochimica Acta Part A:Molecular and Biomolecular Spectroscopy,2016,153:550-559.

    [13] 張 芳,張永清,于 曉,等.忍冬不同種質(zhì)花蕾紅外光譜分析[J].中國(guó)實(shí)驗(yàn)方劑學(xué)雜志,2012,18(9):83-86.

    [14] 趙花榮,溫樹敏,王曉燕,等.金銀花道地品與非道地品的FTIR光譜研究[J].光譜學(xué)與光譜分析,2005,25(5):705-707.

    [15] 邱 璐,李曉勇,劉 鵬,等.九種薔薇科植物葉片的傅里葉紅外光譜與親緣關(guān)系分析[J].光譜學(xué)與光譜分析,2014,34(2):344-349.

    [16] 楊春艷,劉 飛,杜蓮英,等.基于紅外光譜的蘋果也病害的快速診斷[J].光散射學(xué)報(bào),2017,29(2):181-186.

    [17] 楊春艷,劉 飛,王元忠,等.紅外光譜結(jié)合判別分析對(duì)瑪咖產(chǎn)地的鑒別[J].江蘇農(nóng)業(yè)科學(xué),2017,45(5):170-173.

    猜你喜歡
    微分金銀花二階
    金銀花“香溢”致富路
    擬微分算子在Hp(ω)上的有界性
    一類二階迭代泛函微分方程的周期解
    金銀花
    上下解反向的脈沖微分包含解的存在性
    一類二階中立隨機(jī)偏微分方程的吸引集和擬不變集
    二階線性微分方程的解法
    一類二階中立隨機(jī)偏微分方程的吸引集和擬不變集
    借助微分探求連續(xù)函數(shù)的極值點(diǎn)
    金銀花又開
    欧美变态另类bdsm刘玥| 日本wwww免费看| 国产男女内射视频| 丁香六月欧美| 亚洲欧美一区二区三区国产| 国产真人三级小视频在线观看| 亚洲久久久国产精品| 纵有疾风起免费观看全集完整版| 婷婷丁香在线五月| 午夜福利影视在线免费观看| 极品少妇高潮喷水抽搐| 日本黄色日本黄色录像| 日本色播在线视频| 丰满迷人的少妇在线观看| www日本在线高清视频| 亚洲精品中文字幕在线视频| 欧美黑人欧美精品刺激| 国产一区二区三区av在线| 一边亲一边摸免费视频| 热99国产精品久久久久久7| 日本五十路高清| 久久鲁丝午夜福利片| 又大又爽又粗| 日本猛色少妇xxxxx猛交久久| 一本大道久久a久久精品| 丝瓜视频免费看黄片| 欧美黄色片欧美黄色片| 母亲3免费完整高清在线观看| 桃花免费在线播放| 亚洲精品av麻豆狂野| 天堂俺去俺来也www色官网| 亚洲第一av免费看| 大话2 男鬼变身卡| 亚洲精品乱久久久久久| 欧美 亚洲 国产 日韩一| 亚洲国产精品国产精品| 老司机影院毛片| 日韩伦理黄色片| 青春草亚洲视频在线观看| 国产精品一区二区精品视频观看| 热99久久久久精品小说推荐| 欧美亚洲 丝袜 人妻 在线| 男女高潮啪啪啪动态图| 亚洲三区欧美一区| 考比视频在线观看| 国产成人影院久久av| 午夜久久久在线观看| 国产精品熟女久久久久浪| 国产成人系列免费观看| cao死你这个sao货| 伦理电影免费视频| 妹子高潮喷水视频| 黄色一级大片看看| 日韩伦理黄色片| xxxhd国产人妻xxx| www.自偷自拍.com| 精品国产一区二区三区久久久樱花| 一区二区av电影网| 精品高清国产在线一区| 国产欧美日韩一区二区三区在线| 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 99久久综合免费| 欧美 日韩 精品 国产| 欧美黑人欧美精品刺激| 又大又黄又爽视频免费| 国产成人精品久久二区二区免费| 精品亚洲乱码少妇综合久久| 女性生殖器流出的白浆| 日本vs欧美在线观看视频| 亚洲免费av在线视频| 国产99久久九九免费精品| 大片免费播放器 马上看| 激情视频va一区二区三区| 亚洲成人国产一区在线观看 | 亚洲欧美清纯卡通| 青春草视频在线免费观看| 久久人妻熟女aⅴ| 一级毛片 在线播放| 精品少妇内射三级| 少妇精品久久久久久久| 丝袜美腿诱惑在线| 丝袜人妻中文字幕| 国产精品成人在线| 国产极品粉嫩免费观看在线| 精品熟女少妇八av免费久了| 久久精品久久久久久噜噜老黄| 老司机午夜十八禁免费视频| av片东京热男人的天堂| 国产成人免费观看mmmm| 2018国产大陆天天弄谢| 夫妻午夜视频| 黄色a级毛片大全视频| 国产免费一区二区三区四区乱码| 老司机深夜福利视频在线观看 | 国产精品免费视频内射| 免费在线观看影片大全网站 | www.999成人在线观看| 亚洲一码二码三码区别大吗| 丰满人妻熟妇乱又伦精品不卡| 黄色视频在线播放观看不卡| 考比视频在线观看| 一区二区三区精品91| 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 高清视频免费观看一区二区| 99热国产这里只有精品6| 日本av手机在线免费观看| 飞空精品影院首页| 亚洲国产精品999| 男女边摸边吃奶| 国产高清国产精品国产三级| 亚洲av成人不卡在线观看播放网 | 黄色怎么调成土黄色| 久久久精品区二区三区| 国产免费现黄频在线看| 亚洲国产欧美在线一区| a 毛片基地| 国产午夜精品一二区理论片| 婷婷色麻豆天堂久久| 久久人人爽av亚洲精品天堂| 免费在线观看影片大全网站 | 亚洲伊人久久精品综合| 欧美性长视频在线观看| 欧美av亚洲av综合av国产av| 电影成人av| 欧美大码av| 水蜜桃什么品种好| 国产精品九九99| 国产一区二区在线观看av| 99久久综合免费| av在线老鸭窝| 亚洲国产欧美一区二区综合| 啦啦啦中文免费视频观看日本| 韩国精品一区二区三区| 满18在线观看网站| 黄色a级毛片大全视频| 久久久久久久国产电影| 最近最新中文字幕大全免费视频 | 久久久久精品人妻al黑| 七月丁香在线播放| 亚洲欧美日韩另类电影网站| 人人妻人人澡人人爽人人夜夜| 亚洲av男天堂| 99热全是精品| 成人亚洲精品一区在线观看| 国产亚洲欧美精品永久| 国产精品久久久久久精品古装| 精品久久久久久久毛片微露脸 | 久久久久精品国产欧美久久久 | videos熟女内射| 成年美女黄网站色视频大全免费| 中国国产av一级| 国产成人欧美| 一级毛片黄色毛片免费观看视频| 成人影院久久| 成年av动漫网址| 欧美日韩综合久久久久久| 国产精品亚洲av一区麻豆| 久久精品亚洲熟妇少妇任你| 亚洲久久久国产精品| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 国产精品香港三级国产av潘金莲 | 十八禁高潮呻吟视频| 永久免费av网站大全| 久久久久精品国产欧美久久久 | 乱人伦中国视频| 亚洲欧美精品自产自拍| 校园人妻丝袜中文字幕| 亚洲伊人色综图| 高清黄色对白视频在线免费看| 国产精品 欧美亚洲| 久热爱精品视频在线9| 丁香六月天网| avwww免费| 人体艺术视频欧美日本| 高清黄色对白视频在线免费看| 最黄视频免费看| 亚洲国产欧美网| 国产又爽黄色视频| 中文欧美无线码| 少妇猛男粗大的猛烈进出视频| 精品人妻一区二区三区麻豆| 亚洲精品自拍成人| 亚洲av片天天在线观看| 国产又色又爽无遮挡免| 免费不卡黄色视频| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 国产日韩一区二区三区精品不卡| 91成人精品电影| 一级,二级,三级黄色视频| 欧美日韩视频高清一区二区三区二| 韩国精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 免费在线观看日本一区| 久久久久久久精品精品| 天天躁狠狠躁夜夜躁狠狠躁| 久久99精品国语久久久| 免费日韩欧美在线观看| 久久鲁丝午夜福利片| 久久久久视频综合| 999精品在线视频| 丝袜美腿诱惑在线| 人人妻人人添人人爽欧美一区卜| 丝瓜视频免费看黄片| 久久精品国产亚洲av高清一级| 久久ye,这里只有精品| 韩国精品一区二区三区| 一区二区三区精品91| 大片免费播放器 马上看| 少妇精品久久久久久久| 国产熟女午夜一区二区三区| a级毛片黄视频| 国产精品国产av在线观看| 欧美日韩福利视频一区二区| 黄色视频在线播放观看不卡| 免费高清在线观看日韩| 久久午夜综合久久蜜桃| av在线老鸭窝| 99香蕉大伊视频| av不卡在线播放| 欧美老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆| 精品少妇一区二区三区视频日本电影| 欧美日韩亚洲综合一区二区三区_| 亚洲欧洲日产国产| 老熟女久久久| 在线观看免费日韩欧美大片| 一本久久精品| 欧美激情极品国产一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 中文字幕制服av| 51午夜福利影视在线观看| 人人妻人人添人人爽欧美一区卜| av欧美777| 国产熟女午夜一区二区三区| 亚洲欧美一区二区三区国产| 最近中文字幕2019免费版| 777米奇影视久久| 亚洲精品国产区一区二| 一级片免费观看大全| 欧美久久黑人一区二区| 男人添女人高潮全过程视频| 亚洲欧洲国产日韩| 日韩中文字幕欧美一区二区 | 97精品久久久久久久久久精品| 国产欧美亚洲国产| 丰满人妻熟妇乱又伦精品不卡| 日韩制服丝袜自拍偷拍| 久久九九热精品免费| 日韩免费高清中文字幕av| 亚洲精品一二三| 一本色道久久久久久精品综合| 亚洲五月色婷婷综合| 国产成人精品在线电影| 黄色a级毛片大全视频| 免费黄频网站在线观看国产| 婷婷色综合大香蕉| 久久人人爽人人片av| 久久人妻福利社区极品人妻图片 | 久久女婷五月综合色啪小说| 9色porny在线观看| 欧美黑人精品巨大| 中文字幕色久视频| 好男人电影高清在线观看| 国产精品 国内视频| 亚洲成人手机| 久久av网站| 久久久久视频综合| 悠悠久久av| 欧美精品人与动牲交sv欧美| 欧美乱码精品一区二区三区| www.av在线官网国产| svipshipincom国产片| 女人高潮潮喷娇喘18禁视频| 亚洲美女黄色视频免费看| 亚洲精品日本国产第一区| 欧美激情极品国产一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 另类亚洲欧美激情| 午夜激情久久久久久久| 精品亚洲成a人片在线观看| 成人黄色视频免费在线看| 9热在线视频观看99| 成人亚洲欧美一区二区av| 国产99久久九九免费精品| av视频免费观看在线观看| 国产精品亚洲av一区麻豆| 纯流量卡能插随身wifi吗| 国产精品 国内视频| av电影中文网址| 十八禁网站网址无遮挡| 久久国产精品人妻蜜桃| 两个人免费观看高清视频| 在线av久久热| 午夜免费成人在线视频| 国产极品粉嫩免费观看在线| 悠悠久久av| 丝袜人妻中文字幕| 狂野欧美激情性bbbbbb| 亚洲一码二码三码区别大吗| 又黄又粗又硬又大视频| 亚洲视频免费观看视频| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 母亲3免费完整高清在线观看| 国产高清国产精品国产三级| 一区二区日韩欧美中文字幕| 成年女人毛片免费观看观看9 | 亚洲激情五月婷婷啪啪| av电影中文网址| 亚洲精品美女久久av网站| √禁漫天堂资源中文www| av在线老鸭窝| 两人在一起打扑克的视频| 蜜桃国产av成人99| 中文欧美无线码| 欧美日韩精品网址| 久久99一区二区三区| 国产成人精品久久二区二区免费| 老汉色av国产亚洲站长工具| 久久久久久久久免费视频了| 免费在线观看完整版高清| 亚洲国产中文字幕在线视频| 国产免费又黄又爽又色| 人人妻人人澡人人看| 天堂俺去俺来也www色官网| 欧美精品人与动牲交sv欧美| 国产精品久久久久成人av| 看十八女毛片水多多多| xxxhd国产人妻xxx| 欧美黑人欧美精品刺激| 制服人妻中文乱码| 久久天堂一区二区三区四区| 两性夫妻黄色片| 久久国产亚洲av麻豆专区| 国产一级毛片在线| 亚洲av日韩精品久久久久久密 | 亚洲,一卡二卡三卡| 亚洲图色成人| 精品免费久久久久久久清纯 | 亚洲国产欧美网| 国产免费又黄又爽又色| 日韩av不卡免费在线播放| 亚洲国产精品一区三区| 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 一级片'在线观看视频| 久久精品久久久久久久性| 国产在线视频一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久久精品古装| 日韩人妻精品一区2区三区| 午夜福利影视在线免费观看| 天天添夜夜摸| 男女之事视频高清在线观看 | 亚洲成av片中文字幕在线观看| 男女边吃奶边做爰视频| 午夜福利乱码中文字幕| 中文字幕色久视频| 欧美av亚洲av综合av国产av| 18禁裸乳无遮挡动漫免费视频| 老熟女久久久| 两人在一起打扑克的视频| 亚洲国产精品国产精品| 99热网站在线观看| 国产亚洲精品第一综合不卡| 日韩大片免费观看网站| 丝瓜视频免费看黄片| 精品熟女少妇八av免费久了| 亚洲 国产 在线| 女警被强在线播放| 亚洲av在线观看美女高潮| 少妇人妻 视频| 国产无遮挡羞羞视频在线观看| 这个男人来自地球电影免费观看| 午夜两性在线视频| 欧美成狂野欧美在线观看| 一级黄片播放器| 亚洲精品自拍成人| 韩国高清视频一区二区三区| 欧美成狂野欧美在线观看| 妹子高潮喷水视频| 日本av免费视频播放| 久久久久久免费高清国产稀缺| 国产黄频视频在线观看| 亚洲av综合色区一区| 国产精品成人在线| 侵犯人妻中文字幕一二三四区| 天天躁夜夜躁狠狠久久av| 黄色视频在线播放观看不卡| 亚洲精品日韩在线中文字幕| 激情视频va一区二区三区| 91国产中文字幕| 亚洲熟女毛片儿| 免费人妻精品一区二区三区视频| 热re99久久精品国产66热6| 国产黄色免费在线视频| 亚洲欧美激情在线| 亚洲av国产av综合av卡| 黄网站色视频无遮挡免费观看| 色精品久久人妻99蜜桃| 久久人人爽人人片av| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美一区二区三区国产| 免费看十八禁软件| 精品少妇一区二区三区视频日本电影| 最新的欧美精品一区二区| 亚洲欧美日韩另类电影网站| 欧美亚洲 丝袜 人妻 在线| 国产精品一二三区在线看| 天天躁夜夜躁狠狠久久av| 久久久久国产精品人妻一区二区| 制服诱惑二区| 人人妻,人人澡人人爽秒播 | 亚洲精品一区蜜桃| 中文精品一卡2卡3卡4更新| 日韩欧美一区视频在线观看| 欧美久久黑人一区二区| 在线天堂中文资源库| 亚洲国产欧美网| 黄频高清免费视频| 亚洲精品日本国产第一区| 男女边摸边吃奶| 热re99久久精品国产66热6| av福利片在线| 国产亚洲av片在线观看秒播厂| 久久久久国产一级毛片高清牌| 午夜免费观看性视频| 中文字幕最新亚洲高清| 亚洲国产日韩一区二区| 欧美日韩国产mv在线观看视频| 男女之事视频高清在线观看 | 欧美xxⅹ黑人| 欧美成人午夜精品| 日本一区二区免费在线视频| 丝袜喷水一区| 国产精品 欧美亚洲| 亚洲av美国av| 国产成人精品久久二区二区免费| 男女免费视频国产| 少妇 在线观看| 亚洲美女黄色视频免费看| 亚洲 国产 在线| 亚洲精品乱久久久久久| 另类亚洲欧美激情| 久久 成人 亚洲| 香蕉丝袜av| 9热在线视频观看99| 大香蕉久久成人网| www.自偷自拍.com| 熟女少妇亚洲综合色aaa.| a级片在线免费高清观看视频| 久久99一区二区三区| 日日夜夜操网爽| 亚洲精品美女久久久久99蜜臀 | 男女高潮啪啪啪动态图| 日韩,欧美,国产一区二区三区| 亚洲av片天天在线观看| 一区二区日韩欧美中文字幕| 爱豆传媒免费全集在线观看| 免费日韩欧美在线观看| 午夜福利免费观看在线| 性少妇av在线| 免费看av在线观看网站| 中文字幕人妻丝袜制服| 丰满少妇做爰视频| 国产av一区二区精品久久| 欧美精品人与动牲交sv欧美| 欧美成狂野欧美在线观看| 日本wwww免费看| 欧美日韩视频高清一区二区三区二| 成人国产一区最新在线观看 | 国产成人影院久久av| 男人爽女人下面视频在线观看| 另类亚洲欧美激情| 女警被强在线播放| 一级黄片播放器| 男女床上黄色一级片免费看| 亚洲人成电影免费在线| 在现免费观看毛片| 国产有黄有色有爽视频| 精品一区在线观看国产| 亚洲精品成人av观看孕妇| 国产成人精品久久二区二区91| 亚洲欧洲日产国产| 啦啦啦啦在线视频资源| 亚洲激情五月婷婷啪啪| 黄色怎么调成土黄色| 亚洲精品国产av蜜桃| 老司机影院毛片| 中文字幕av电影在线播放| 人人妻,人人澡人人爽秒播 | 亚洲精品日韩在线中文字幕| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产看品久久| 免费人妻精品一区二区三区视频| 中文字幕高清在线视频| 狠狠精品人妻久久久久久综合| 成人手机av| 99九九在线精品视频| 日韩欧美一区视频在线观看| 中文欧美无线码| 久久精品国产综合久久久| 国产xxxxx性猛交| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| 国产亚洲av片在线观看秒播厂| 免费女性裸体啪啪无遮挡网站| 亚洲av电影在线观看一区二区三区| 欧美精品亚洲一区二区| 亚洲国产精品一区二区三区在线| 国产成人av激情在线播放| 女人高潮潮喷娇喘18禁视频| 麻豆国产av国片精品| 欧美日韩国产mv在线观看视频| 青春草视频在线免费观看| 制服诱惑二区| 超碰97精品在线观看| 天天添夜夜摸| 色综合欧美亚洲国产小说| 波多野结衣av一区二区av| 午夜福利免费观看在线| 亚洲精品日韩在线中文字幕| 欧美国产精品va在线观看不卡| 久久精品国产亚洲av涩爱| 99久久综合免费| 国产精品一区二区在线观看99| 亚洲欧美中文字幕日韩二区| 欧美精品人与动牲交sv欧美| 亚洲 国产 在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人欧美| 久久精品久久久久久噜噜老黄| 亚洲少妇的诱惑av| 两人在一起打扑克的视频| 国产91精品成人一区二区三区 | 99国产精品免费福利视频| 中文字幕最新亚洲高清| 大陆偷拍与自拍| 国产成人a∨麻豆精品| 亚洲av欧美aⅴ国产| 一级a爱视频在线免费观看| 男女下面插进去视频免费观看| 国产高清视频在线播放一区 | 少妇人妻久久综合中文| 午夜免费观看性视频| 国产免费视频播放在线视频| 女性生殖器流出的白浆| 久久国产精品影院| 亚洲免费av在线视频| 国产野战对白在线观看| 精品一区二区三区av网在线观看 | 亚洲 欧美一区二区三区| 欧美日韩综合久久久久久| 国产精品一二三区在线看| 高清黄色对白视频在线免费看| 中文字幕亚洲精品专区| 免费久久久久久久精品成人欧美视频| 午夜免费观看性视频| 午夜日韩欧美国产| 国产黄色免费在线视频| 国产无遮挡羞羞视频在线观看| 婷婷丁香在线五月| av有码第一页| 一本综合久久免费| 欧美精品一区二区免费开放| 在线观看一区二区三区激情| 免费观看a级毛片全部| 国产成人系列免费观看| 国产精品一区二区免费欧美 | 只有这里有精品99| 建设人人有责人人尽责人人享有的| 少妇的丰满在线观看| 高清视频免费观看一区二区| 久久 成人 亚洲| 久久精品成人免费网站| 少妇精品久久久久久久| 国产亚洲欧美精品永久| 人人妻人人澡人人爽人人夜夜| 水蜜桃什么品种好| 麻豆av在线久日| 亚洲少妇的诱惑av| 9色porny在线观看| 久久鲁丝午夜福利片| 一级a爱视频在线免费观看| 一级毛片 在线播放| 国产精品香港三级国产av潘金莲 | 亚洲国产精品国产精品| 男女边吃奶边做爰视频| 最近中文字幕2019免费版| 少妇 在线观看| 最黄视频免费看| 日本欧美视频一区| 久久热在线av| videosex国产| 热re99久久精品国产66热6| a级毛片在线看网站| 日本91视频免费播放| 各种免费的搞黄视频| 国产不卡av网站在线观看| 婷婷色综合www| 国产精品 国内视频| 午夜福利视频在线观看免费| 亚洲av男天堂| 男人添女人高潮全过程视频| 欧美人与性动交α欧美精品济南到| 老熟女久久久| 69精品国产乱码久久久| 亚洲 国产 在线| 国产精品一区二区在线不卡| www.999成人在线观看| a 毛片基地|