莊 為李臻峰,2宋飛虎,2李 靜,2汪迪松朱冠宇
(1. 江南大學(xué)機(jī)械工程學(xué)院,江蘇 無錫 214122;2. 江蘇省食品先進(jìn)制造裝備技術(shù)重點(diǎn)實(shí)驗(yàn)室,江蘇 無錫 214122)
中國(guó)西瓜產(chǎn)量一直穩(wěn)居世界第一,在世界總產(chǎn)量中始終占據(jù)60%以上[1]。西瓜傳統(tǒng)的鑒別方法主要是看、聽,效率極低。而振動(dòng)聲學(xué)檢測(cè)方法簡(jiǎn)單,抗干擾能力強(qiáng),檢測(cè)裝置成本低廉,易實(shí)現(xiàn)小型化和智能化,被廣泛用于瓜果品質(zhì)檢測(cè)[2]。
國(guó)內(nèi)外利用振動(dòng)特性對(duì)西瓜的研究備受關(guān)注, Nourain等[3]研究發(fā)現(xiàn)其楊氏模量與固有頻率有著較好的相關(guān)關(guān)系;Abbasadeh等[4-5]利用激光多普勒測(cè)振儀非接觸測(cè)量西瓜的振動(dòng)特性,利用頻域信號(hào)對(duì)西瓜的成熟度進(jìn)行預(yù)測(cè);Mao等[6]利用聲學(xué)裝置研究硬度指標(biāo)與mf2的關(guān)系;Taniwaki等[7-8]利用激光多普勒測(cè)振儀測(cè)量西瓜的振動(dòng)頻譜,計(jì)算彈性參數(shù)EI=f22m2/3,研究通過EI值預(yù)測(cè)2種甜瓜達(dá)到最佳成熟的時(shí)間范圍;方漢良等[9]開發(fā)出一套通過對(duì)敲打西瓜產(chǎn)生的振動(dòng)頻率及聲音頻率進(jìn)行測(cè)量、分析判斷西瓜是否成熟的測(cè)試儀;浦宏杰等[10]基于振動(dòng)特性用子帶頻譜質(zhì)心描述西瓜內(nèi)部品質(zhì)中的糖度及瓤色;陸勇等[11]檢測(cè)西瓜的振動(dòng)頻率隨貯藏時(shí)間的變化規(guī)律;高宗梅等[12]利用激光多普勒測(cè)振系統(tǒng)研究了加速度振幅、掃頻速率、檢測(cè)點(diǎn)位置3個(gè)因素對(duì)西瓜振動(dòng)頻譜響應(yīng)特性的影響;魏巍[13]分析了西瓜品質(zhì)指標(biāo)與振動(dòng)參數(shù)的相關(guān)性。國(guó)內(nèi)外的研究多數(shù)提取振動(dòng)信號(hào)的某些主要頻率值,研究這些頻率值與品質(zhì)的關(guān)系,但是沒有充分考慮頻率間的相互影響,并且這些研究是針對(duì)某特定品種的西瓜而言的,沒有充分考慮西瓜品種的不同對(duì)建模的影響以及對(duì)檢測(cè)效果的影響。
目前,國(guó)內(nèi)外應(yīng)用聲學(xué)振動(dòng)特性鑒別西瓜品種鮮有報(bào)道,主要是研究某特定品種瓜果的品質(zhì)與振動(dòng)參數(shù)的關(guān)系[6]、用振動(dòng)參數(shù)描述品質(zhì)的變化趨勢(shì)[11]以及研究不同品種瓜果的某品質(zhì)能否用同一類型的振動(dòng)參數(shù)來描述[14]。建模時(shí),主要是利用振動(dòng)信號(hào)的主要共振峰頻率與西瓜品質(zhì)的關(guān)系。但是,不同品種西瓜的內(nèi)部品質(zhì)預(yù)測(cè)模型有所差異,判斷其品質(zhì)時(shí)會(huì)因品種的不同而造成誤判,因此,實(shí)現(xiàn)對(duì)未知西瓜的品種鑒別到品質(zhì)分級(jí)的過程倍受關(guān)注,對(duì)檢測(cè)精度的提升度具有一定的意義。目前,利用聲學(xué)振動(dòng)特性研究西瓜從品種鑒別到品質(zhì)分級(jí)過程鮮有報(bào)道。本研究擬以市面上常見的皇冠、早佳8424、京欣和硒砂4個(gè)西瓜品種為樣本,基于振動(dòng)特性分析不同品種西瓜振動(dòng)頻譜的3個(gè)主要共振峰的關(guān)系,利用神經(jīng)網(wǎng)絡(luò)建立西瓜品種鑒別模型。同時(shí),研究4個(gè)品種西瓜的振動(dòng)峰頻率和糖度的關(guān)系。
西瓜:早佳8424、皇冠、京欣和硒砂,適熟瓜,樣本數(shù)分別為44,20,20,20個(gè),分別對(duì)每種西瓜編號(hào)。早佳8424西瓜樣本采摘于無錫北塘區(qū)西瓜大棚內(nèi),其它品種樣本均于無錫濱湖某農(nóng)貿(mào)市場(chǎng)購買。采購時(shí),盡量控制西瓜樣本瓜形、大小、色澤均勻一致。所有測(cè)試均在采摘后1 d內(nèi)完成,樣本保存環(huán)境為室溫25 ℃,相對(duì)濕度60%~70%。
數(shù)字手持折射儀:PAL-1型,日本ATAGO公司;
信號(hào)激勵(lì)力錘:086C01型,美國(guó)PCB公司;
信號(hào)接收加速度傳感器:352C68型,美國(guó)PCB公司;
信號(hào)調(diào)節(jié)儀:482C05型,美國(guó)PCB公司;
NI數(shù)據(jù)采集卡:USB-6259型,美國(guó)NI公司;
西瓜的敲擊振動(dòng)檢測(cè)系統(tǒng)見圖1。
圖1 西瓜敲擊振動(dòng)檢測(cè)系統(tǒng)圖
1.3.1 頻率響應(yīng)信號(hào)的采集 本研究先確定試驗(yàn)因素,包括沖擊錘頭材料、測(cè)試過程中加權(quán)平均次數(shù)和敲擊位置。試驗(yàn)時(shí),將西瓜沿著瓜臍和瓜蒂方向水平放置在工作臺(tái)上,力錘在相應(yīng)部位垂直于西瓜表面敲擊,激勵(lì)信號(hào)和振動(dòng)信號(hào)見圖2,圖3為L(zhǎng)abview程序界面上顯示的頻率響應(yīng)函數(shù),計(jì)算機(jī)利用Labview編寫的程序處理數(shù)據(jù)并保存,以便后續(xù)處理。
圖2 力錘沖擊力信號(hào)和加速度沖擊響應(yīng)信號(hào)
Figure 2 Force hammer impact force signal and acceleration impulse response signal
圖3 西瓜的頻率響應(yīng)函數(shù)FRF
1.3.2 糖度檢測(cè) 本研究通過破壞性試驗(yàn)測(cè)量西瓜的糖度。試驗(yàn)時(shí),將西瓜沿著瓜蒂和瓜臍對(duì)半切開,將同一西瓜樣本的瓜臍部位、瓜蒂部位和瓜中心部位的瓜瓤取出榨汁,利用折射儀測(cè)量樣本的糖度,每個(gè)西瓜樣本取3次進(jìn)行測(cè)量,取平均值。折射儀使用前需要清洗清零,每次測(cè)量需用清水洗凈擦干后再進(jìn)行下次測(cè)量。
本研究選取頻率響應(yīng)函數(shù)的共振峰頻率作為振動(dòng)特征值,如圖3所示,選取3個(gè)主要共振峰頻率f1、f2和f3。
2.2.1 不同材質(zhì)錘頭對(duì)頻率響應(yīng)函數(shù)的影響 本研究以早佳8424西瓜樣本為例,將加速度傳感器固定在瓜臍位置,力錘敲擊瓜蒂位置,沖擊錘頭選用鋁、塑料、硬橡膠和軟橡膠4種材質(zhì),采用單次敲擊。圖4表明,使用鋁質(zhì)錘頭的頻率響應(yīng)函數(shù)在250 Hz之前的曲線更加平滑,出現(xiàn)了4個(gè)較平滑明顯的共振峰,故試驗(yàn)采用鋁質(zhì)錘頭。
圖4 不同材質(zhì)錘頭的頻率響應(yīng)函數(shù)
采用鋁質(zhì)錘頭敲擊西瓜表面,在同一位置敲擊3次。如圖5所示,3次頻率響應(yīng)函數(shù)的重合度非常高,表明系統(tǒng)的單次敲擊穩(wěn)定性較高,能夠采集到穩(wěn)定數(shù)據(jù)滿足本研究的要求。
圖5 鋁質(zhì)錘頭3次敲擊的頻率響應(yīng)函數(shù)
Figure 5 The frequency response function of three strokes of aluminum hammer head
2.2.2 加權(quán)平均次數(shù)對(duì)頻率響應(yīng)函數(shù)的影響 以早佳8424西瓜樣本為例,將加速度傳感器固定在瓜臍位置,力錘敲擊位置為瓜蒂位置,采用單次敲擊。分別采用軟橡膠、硬橡膠、塑料和鋁質(zhì)4種錘頭敲擊西瓜表面,得到1次敲擊和10次敲擊加權(quán)平均后的頻率響應(yīng)函數(shù)。如圖6~9所示,不同錘頭經(jīng)過10次平均后,頻響曲線相比單次平滑??梢姸啻纹骄螅盘?hào)中的噪聲降低,信號(hào)平滑度更高。本試驗(yàn)采用的加權(quán)平均次數(shù)為10次。
2.2.3 敲擊位置對(duì)頻率響應(yīng)函數(shù)的影響 以早佳8424西瓜樣本為例,采用鋁質(zhì)錘頭進(jìn)行10次均勻、游動(dòng)敲擊,如圖10所示,將加速度傳感器固定于瓜臍處,沖擊力錘在圖中點(diǎn)1~8處逐一敲擊,得到各點(diǎn)處的頻響見圖11。圖11中可見,A處共振峰的幅值隨各點(diǎn)大小不一,但是頻率基本相同。B處共振峰有部分明顯平滑,其頻率也基本相同,B處后面的共振峰不明顯,但共振峰處頻率也有相同的趨勢(shì)。由此可見,敲擊位置對(duì)頻響函數(shù)的幅值有影響,但是對(duì)頻率的影響基本可以忽略。本試驗(yàn)將加速度傳感器固定在瓜臍位置,力錘敲擊位置為瓜蒂位置。
圖6 軟橡膠錘頭單次和10次平均的頻率響應(yīng)函數(shù)
圖7 硬橡膠錘頭單次和10次平均的頻率響應(yīng)函數(shù)
圖8 塑料錘頭單次和10次平均的頻率響應(yīng)函數(shù)
圖12為不同品種西瓜樣本的3個(gè)共振峰頻率關(guān)系圖,第1共振峰頻率為x軸,第2共振峰頻率為y軸,第3共振峰頻率為z軸,可以看出,早佳8424、京欣和皇冠樣本的3個(gè)共振峰頻率點(diǎn)的線性關(guān)系明顯,并且明顯分布在不同區(qū)域。但是硒砂樣本的共振峰頻率點(diǎn)分布沒有其它3種西瓜明顯,但是其頻率點(diǎn)分布最低,并且有部分硒砂樣本的頻率點(diǎn)分布于京欣樣本的頻率點(diǎn)分布區(qū)域。
圖9 鋁質(zhì)錘頭單次和10次平均的頻率響應(yīng)函數(shù)
圖10 敲擊-測(cè)點(diǎn)組合分布圖
圖11 點(diǎn)1~8處逐一敲擊的頻率響應(yīng)函數(shù)
根據(jù)這一分布特點(diǎn),利用PNN神經(jīng)網(wǎng)絡(luò)模型對(duì)該4種西瓜進(jìn)行分類。分類變量分別為1(早佳8424)、2(京欣)、3(硒砂)、4(皇冠)。選取3個(gè)共振峰頻率f1、f2、f3作為PNN神經(jīng)網(wǎng)絡(luò)的輸入,建立網(wǎng)絡(luò)模型。每種西瓜樣本取20個(gè),共80個(gè)樣本,其中每種西瓜取10個(gè)樣本用于建立網(wǎng)絡(luò)模型,余下的樣本用于預(yù)測(cè)。
圖12 不同品種西瓜的第1共振峰頻率、第2共振峰頻率和
Figure 12 The diagram of the first formant frequency, the second formant frequency and third formant frequency resonance peak frequency of different varieties of watermelon
該模型是通過Matlab R2015a建立的,首先將采集到的3個(gè)峰值頻率輸入,利用40個(gè)樣本作為訓(xùn)練樣本,得到西瓜品種識(shí)別的PNN網(wǎng)絡(luò)模型,再用40個(gè)樣本作為測(cè)試樣本,進(jìn)行網(wǎng)絡(luò)性能檢驗(yàn),判斷該網(wǎng)絡(luò)能否用來對(duì)未知樣本的類別進(jìn)行預(yù)測(cè)。
由表1可知,京欣和皇冠樣本的真實(shí)類別和預(yù)測(cè)類別完全一致,早佳8424樣本中有1個(gè)被鑒別為皇冠,硒砂樣本中有2個(gè)樣本被鑒別為京欣,預(yù)測(cè)集的準(zhǔn)確率為92.5%,模型效果可以接受,說明采用聲學(xué)振動(dòng)技術(shù)可以解決不同品種西瓜的分類問題。
表1 PNN神經(jīng)網(wǎng)絡(luò)模型對(duì)未知樣本的判別結(jié)果
王書茂等[15]研究發(fā)現(xiàn)西瓜的第1階固有頻率與含糖量有較好的關(guān)系,還發(fā)現(xiàn)含糖量不同,瓜瓤的彈性模量變化較大。瓜瓤彈性模量的變化關(guān)系到西瓜的振動(dòng)特性,瓜瓤的彈性模量變小,質(zhì)量一定的西瓜固有頻率必然降低[15]。危艷君[16]、浦宏杰[10]等利用聲學(xué)特征來建立糖度檢測(cè)模型。聲學(xué)檢測(cè)技術(shù)應(yīng)用于西瓜糖度檢測(cè)[17]等領(lǐng)域,西瓜品質(zhì)的變化主要表現(xiàn)為含糖量等參數(shù)的變化[18],含糖量是評(píng)價(jià)西瓜品質(zhì)的最核心指標(biāo)[19]。因此,研究糖度和聲學(xué)特征之間的相關(guān)關(guān)系,可以實(shí)現(xiàn)對(duì)不同含糖量西瓜的無損檢測(cè)。
由于4種西瓜糖度預(yù)測(cè)模型建立的方法一致,故僅以早佳8424瓜為例進(jìn)行建模,校正集34個(gè)。如圖13~15所示,分別為早佳8424樣本的3個(gè)主要共振峰頻率f1、f2、f3與糖度的關(guān)系圖,表2為各共振峰頻率與糖度建立線性模型,模型的確定系數(shù)R2均偏低。
在自變量較多時(shí),其中有的自變量可能對(duì)因變量的影響不是很大,而且3個(gè)共振峰頻率之間可能不完全相互獨(dú)立的。在這種情況下可用逐步多元線性回歸分析,對(duì)3個(gè)共振峰頻率進(jìn)行篩選,以此建立的多元回歸模型預(yù)測(cè)效果會(huì)更好。
圖13 早佳8424樣本的第1共振峰頻率f1與糖度的關(guān)系圖
Figure 13 The relationship of the first peak frequencyf1and the sugar degree of “Zaojia 8424”samples
圖14 早佳8424樣本的第2共振峰頻率f2與糖度的關(guān)系圖
Figure 14 The relationship of the second peak frequencyf2and the sugar degree of “Zaojia 8424”samples
圖15 早佳8424樣本的第3共振峰頻率f3與糖度的關(guān)系圖
Figure 15 The relationship of the third peak frequencyf3and the sugar degree of “Zaojia 8424”samples
表2早佳8424樣本峰值頻率與糖度之間的線性模型
Table 2 The linear model of the peak frequency and the sugar degree of “Zaojia 8424”samples
峰值頻率f擬合方程R2f1f1=-5.712 16x+233.264 80.732 50f2f2=-6.821 2x+317.0790.713 71f3f3=-9.114x+401.30.566 55
統(tǒng)計(jì)數(shù)據(jù)計(jì)算使用軟件SPSS 21.0,在每一步,檢驗(yàn)統(tǒng)計(jì)量F的概率小于設(shè)置值的自變量被進(jìn)入,檢驗(yàn)統(tǒng)計(jì)量F的概率大于設(shè)置值的自變量被移除,當(dāng)沒有更多的變量被納入或移除的時(shí)候逐步多元線性回歸就結(jié)束了。如表3所示,為早佳8424樣本糖度預(yù)測(cè)模型,T表示樣本糖度,模型的確定系數(shù)依次升高,標(biāo)準(zhǔn)估計(jì)的誤差依次降低,因此,選用模型:
T=-0.079×f1-0.022×f2-0.032×f3+38.528。
(1)
其它品種西瓜的糖度預(yù)測(cè)模型建立的過程與早佳8424樣本模型的一致,并且模型的確定系數(shù)均在0.86 以上,預(yù)測(cè)模型見表4。陳春雷等[20]研究發(fā)現(xiàn)不同品種西瓜的皮厚、中糖、邊糖、外觀和風(fēng)味等參數(shù)特征有所差異。并且不同品種西瓜果實(shí)的瓤色不同[21]。王書茂等[15]研究發(fā)現(xiàn)含糖量和瓜瓤不同,振動(dòng)特性不同。因此各品種西瓜的振動(dòng)頻率和固有頻率會(huì)受到影響,糖度預(yù)測(cè)模型有所差異。當(dāng)有其它品種西瓜時(shí),利用上述方法建模,增加相應(yīng)的模型,檢測(cè)時(shí)便可以根據(jù)該品種的模型先鑒別后檢測(cè)。
表3 早佳8424樣本糖度預(yù)測(cè)模型
表4 西瓜糖度預(yù)測(cè)模型
利用分類模型鑒別出的10個(gè)早佳8424瓜作為預(yù)測(cè)集,其中9個(gè)為早佳8424瓜,1個(gè)為皇冠瓜。圖16為早佳8424樣本的預(yù)測(cè)糖度與實(shí)際糖度關(guān)系圖,其中實(shí)際糖度為10%、預(yù)測(cè)糖度為12.279%的坐標(biāo)點(diǎn)為誤判的皇冠樣本預(yù)測(cè)糖度與實(shí)際糖度關(guān)系點(diǎn),剔除該點(diǎn)后,對(duì)余下的樣本建立線性回歸模型,模型確定系數(shù)為0.859 96,標(biāo)準(zhǔn)估計(jì)的誤差為0.631 81??梢?,該糖度預(yù)測(cè)模型可以用來對(duì)西瓜糖度進(jìn)行預(yù)測(cè)。表5為早佳8424預(yù)測(cè)集樣本的實(shí)際糖度值與預(yù)測(cè)糖度值的比較及誤差研究,2號(hào)樣本為誤判的皇冠樣本,剔除皇冠樣本后,9個(gè)預(yù)測(cè)集樣本平均誤差為6.1%,預(yù)測(cè)集樣本的標(biāo)準(zhǔn)差(RMSEP)為1.08,沒有剔除皇冠樣本的10個(gè)預(yù)測(cè)集樣本平均誤差為7.56%,預(yù)測(cè)集樣本的標(biāo)準(zhǔn)差(RMSEP)為1.60,并且對(duì)于誤判的皇冠樣本,其相對(duì)誤差達(dá)到20.72%。校正集樣本的標(biāo)準(zhǔn)差(RMSEC)為0.986。對(duì)于其中9個(gè)品種鑒別準(zhǔn)確的瓜,其糖度測(cè)量誤差小于7%;而對(duì)于品種鑒別錯(cuò)誤的瓜,由于模型選用了錯(cuò)誤預(yù)測(cè)模型,測(cè)量誤差較大,可見,對(duì)未知市售西瓜先鑒別后測(cè)量糖度是非常重要的。
本研究自行搭建的振動(dòng)檢測(cè)系統(tǒng)完成基于振動(dòng)特性的西瓜品種鑒別以及西瓜糖度預(yù)測(cè)研究。通過試驗(yàn)分析皇冠、早佳8424、京欣和硒砂4個(gè)西瓜品種頻率響應(yīng)函數(shù)的3個(gè)主要共振峰頻率,發(fā)現(xiàn)4種西瓜樣本的頻率點(diǎn)分布區(qū)域明顯不同,利用PNN神經(jīng)網(wǎng)絡(luò)可有效地對(duì)不同品種西瓜進(jìn)行鑒別。
圖16 早佳8424樣本的實(shí)際糖度與預(yù)測(cè)糖度關(guān)系圖
Figure 16 The diagram of the actual sugar degree and predicted sugar degree of “Zaojia 8424”samples
表5早佳8424樣本實(shí)際糖度值與預(yù)測(cè)糖度值的比較及誤差
Table 5 The comparison of actual sugar degree and predicted sugar degree of “Zaojia 8424”samples
樣本實(shí)際測(cè)量值模型測(cè)量值相對(duì)誤差/%剔除皇冠瓜樣本后的平均誤差/%112.5011.03911.69210.0012.27920.72312.8012.2474.3249.209.7345.80511.2010.0282.60611.409.22012.04711.3010.8274.19810.509.21312.2699.509.3271.821010.009.9830.176.1
同時(shí),在此基礎(chǔ)上研究4種西瓜的糖度預(yù)測(cè)模型,以本地早佳8424樣本為例,分別建立3個(gè)主要共振峰頻率與糖度的線性模型,模型的確定系數(shù)R2分別為0.732 50,0.713 70,0.566 55,相關(guān)性均較低。這種情況下采用逐步多元線性回歸分析,對(duì)每種西瓜的3個(gè)共振峰頻率進(jìn)行篩選,得到相應(yīng)糖度預(yù)測(cè)模型的確定系數(shù)均大于0.86。經(jīng)驗(yàn)證,該模型可以用來對(duì)西瓜糖度進(jìn)行預(yù)測(cè)。對(duì)預(yù)測(cè)集樣本先鑒別品種再預(yù)測(cè)糖度,對(duì)品種鑒別準(zhǔn)確的瓜,其糖度測(cè)量誤差較小,而對(duì)品種鑒別錯(cuò)誤的瓜,測(cè)量誤差較大,由此可知,先鑒別西瓜的品種,進(jìn)而選取相應(yīng)模型測(cè)量糖度是非常重要的。
[1] 王琛, 張琳, 趙姜, 等. 中國(guó)西瓜市場(chǎng)形勢(shì)分析與展望[J]. 農(nóng)業(yè)展望, 2013, 9(4): 27-30.
[2] 呂飛玲. 基于聲學(xué)特性的西瓜品質(zhì)無損檢測(cè)方法研究[D]. 杭州: 浙江大學(xué), 2003: 10.
[3] NOURAIN Jamal, YING Yi-bin, WANG Jian-ping, et al. Firmness evaluation of melon using its vibration characteristic and finite element analysis[J]. Journal of Zhejiang University Science, 2005(6): 483-490.
[4] ABBASZADEH R, RAJABIPOUR A, DELSHAD M, et al. Evaluation of watermelons texture using their vibration responses[J]. Biosystems Engineering, 2013, 115(1): 102-105.
[5] ABBASZADEH R, RAJABIPOUR A, AHMADI H, et al. Prediction of watermelon quality based on vibration spectrum[J]. Postharvest Biology and Technology, 2013, 86: 291-293.
[6] MAO Jian-hua, YU Yong-hua, RAO Xiu-qin, et al. Firmness prediction and modeling by optimizing acoustic device for watermelons[J]. Journal of Food Engineering, 2016, 168: 1-6.
[7] TANIWAKI M, TOHRO M, SAKURAI N. Measurement of ripening speed and determination of the optimum ripeness of melons by a nondestructive acoustic vibration method[J]. Postharvest Biology and Technology, 2010, 56(1): 101-103.
[8] TANIWAKI M, TAKAHASHI M, SAKURAI N. Determina-tion of optimum ripeness for edibility of postharvest melons using nondestructive vibration[J]. Food Research International, 2009, 42(1): 137-141.
[9] 方漢良, 任志英, 鄭志強(qiáng), 等. 西瓜成熟度測(cè)試儀設(shè)計(jì)研究[J]. 赤峰學(xué)院學(xué)報(bào): 自然科學(xué)版, 2015, 31(12): 53-54.
[10] 浦宏杰, 汪迪松, 李臻峰. 基于振動(dòng)頻譜子帶質(zhì)心的西瓜內(nèi)部品質(zhì)檢測(cè)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2017, 48(1): 281-288.
[11] 陸勇, 李臻峰, 浦宏杰, 等. 基于聲振法的西瓜貯藏時(shí)間檢測(cè)[J]. 浙江農(nóng)業(yè)學(xué)報(bào), 2016, 28(4): 682-687.
[12] 高宗梅, 張文, 任夢(mèng)佳, 等. 西瓜品質(zhì)LDV測(cè)振響應(yīng)影響因素分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2015, 46(S1): 134-140.
[13] 魏巍. 西瓜振動(dòng)特性及振動(dòng)對(duì)其品質(zhì)影響研究[D]. 呼和浩特市: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2015: 24-32.
[14] 屠康, 馬龍, 潘秀娟. 敲擊振動(dòng)無損檢測(cè)3種梨果實(shí)品質(zhì)參數(shù)的研究[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報(bào), 2005(1): 50-53.
[15] 王書茂, 焦群英, 籍俊杰. 西瓜成熟度無損檢驗(yàn)的沖擊振動(dòng)方法[J]. 農(nóng)業(yè)工程學(xué)報(bào), 1999, 15(3): 241-245.
[16] 危艷君, 饒秀勤, 漆兵. 基于聲學(xué)特性的西瓜糖度檢測(cè)系統(tǒng)[J]. 農(nóng)業(yè)工程報(bào), 2012, 28(3): 283-287.
[17] 汪迪松, 浦宏杰, 李臻峰, 等. 基于聲振法的西瓜內(nèi)部糖度檢測(cè)研究[J]. 上海農(nóng)業(yè)報(bào), 2017, 33(2): 125-130.
[18] 趙洪衛(wèi), 韓東海, 宋曙輝, 等. 小型西瓜果實(shí)成熟度表征因子篩選[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2012, 28(17): 281-286.
[19] 任毅. 西瓜果實(shí)含糖量QTL定位及糖轉(zhuǎn)運(yùn)蛋白基因功能初析[D]. 北京: 中國(guó)農(nóng)業(yè)科學(xué)院, 2016: 6.
[20] 孫菊華, 陳春雷, 陳建才. 甜瓜品比試驗(yàn)簡(jiǎn)報(bào)[J]. 上海農(nóng)業(yè)科技, 2011(1): 50-51.
[21] 朱迎春, 劉君璞, 鄧云, 等. 河南地區(qū)中、小果型西瓜品種比較試驗(yàn)[J]. 中國(guó)瓜菜, 2016, 29(12): 51-54.