论文-利用超声检测西瓜的成熟度

分类号： O426.9本科生毕业论文（设计）题目：西瓜成熟度的超声检测法作者单位物理学与信息技术学院作者姓名张俊贤专业班级应用电子技术教育（1）班指导教师（职称）郭建中（教授）论文（设计）完成时间二〇一四年五月西瓜成熟度的超声检测法张俊贤陕西师范大学物理学与信息技术学院，陕西西安 710062摘要: 超声波法无损检测技术以其操作简单、快速、适应性强、检测精度较高等特点而在农产品品质无损检测中得到了广泛应用。

为达无损检测西瓜成熟度的目的，本研究在自主搭建的超声学检测台上采集透过西瓜的超声学信号，提取超声从脉冲发射点传到接收点过程中声透过率参数（δ）。

从中选取相关系数前30位的声透过率值，借助软件中多元线性回归函数建立多元线性回归方程。

研究不同发射点和信号接收点组合所建模型的预测值和实测值。

试验结果表明，发射点和接收点位于赤道部位能获得较好的测定模型。

利用Matlab软件中的仿真技术，对超声衰减与西瓜硬度进行仿真。

成熟度指数与坚硬度有较为接近的线性关系。

因此,用西瓜成熟度指数来检验西瓜的成熟程度是一种比较有效的无损检测方法。

关键词：西瓜；超声衰减；Matlab仿真；坚硬度1前言随着现代食品工业的发展以及为了满足市场竞争的需求，人们越来越多地把主要精力投向提高产品品质。

因此，建立有效的农产品无损检测方法对促进食品产业健康蓬勃发展以及提高国民生活质量发挥着重要作用。

在农产品加工中，水果的成熟度以及随后的质量评估都非常重要。

要求有一种可靠、快速、非破坏性的检测技术，用来测量水果中与成熟和品质有关的一些物理性质。

西瓜作为世界上的重要水果之一，为人们所喜爱，给人们带来可观的经济效益。

在判断西瓜成熟度方面,人们积累了丰富的经验,比如“计算生长期”,“观察外部形态”,“触、拍、压、闻”等无损检验法,为消费者提供成熟的西瓜起到了关键作用。

但这些方法依赖于个人经验等主观因素,尤其消费者难以掌握。

然而按西瓜成熟度分级销售,对销售者和消费者都是有益的。

因此,科学地、客观地对西瓜进行无损检验变得尤其重要,国内外学者通过检测西瓜瓜瓤硬度、颜色、糖度来判断西瓜成熟度的方法已经开始了研究。

为了达到检测西瓜成熟度且又不破坏西瓜的目的，国内外众多学者利用西瓜的声学特性来实现西瓜内部品质的研究，常用的声学参数是声振荡频率。

1975年，Clark[1]研究表明声波在西瓜内部的传递和反射时间随着瓜瓤成熟、变红而延长，从而实现西瓜成熟度的检测。

1968年Abbot等人[2]研究利用苹果声学特性检测苹果硬度，建立了弹性模量与苹果固有基频关系方程，后来的研究对此方程做了改进，现已成为声学参数提取的重要等式；有研究表明，众多声学参数与西瓜质量、瓜瓤糖度、弹性模量、瓜瓤颜色及空心的相关系数并不高，而且不同品种的西瓜相关性也不一致，所以还有待进一步研究[3]。

西瓜越成熟，质地越软，所以可用西瓜的声振动频率来计算瓜瓤弹性模量，与西瓜的成熟度相关性比较好。

同时由声谱图的峰值可以判定西瓜空心存在[4]。

Farabee等人[5]在频谱图中提取含50%能量波段的中心频率参数以及85～160Hz内的能量值，研究这2个声学参数与西瓜成熟度之间的关系，通过线性回归分析，发现两者与西瓜含糖量相关系数分别是0.67和0.56。

JSugiyama[6]于1998年开发便携式仪器测定甜瓜硬度；2009年肖珂等人[7]根据声谱图中峰值出现的频率范围，可以判定西瓜的成熟度。

其中频率范围大于189Hz 的是未成熟瓜，而在160～189Hz内是成熟瓜，在133～160Hz的是过熟瓜。

BDiezma-Iglesias[8-9]利用声学检测西瓜内部空心，并在后续研究中建立西瓜的有限元模型，与理论模型比较发现此方法更适用于长型西瓜检测。

该项技术在国内发展仍处于初步理论研究阶段[10-11]。

如何找到更好的超声学参数以及适用于不同品种的检测模型一直是该领域研究的重点，为此本实验利用超声信号在西瓜中传递容易衰减的特性，提取超声衰减特征参数值（δ），与糖度做相关性分析，实现西瓜糖度检测，判断西瓜成熟度。

2 超声波法检测技术原理超声波在被检测物料中传播时，物料的声学特性和内部组织的变化会对超声波的传播产生一定的影响，通过对超声波受影响程度和状况的探测了解物料性能和结构变化的技术称为超声波法检测技术。

超声波无损检测技术是一门涉及多种学科的综合技术，利用低能量超声波对处理对象进行无损伤探测（NDT ），其特点是频率较高（介于0.1～0.2MHz ）、能量较低（往往小于100mW ）和大多采用脉冲式操作。

因为能量低，所以声波通过时不会根本改变介质的物理或化学性质。

超声波通过介质时大致表现为三种形式：压缩波、表面波和切变波。

在NDT 应用中，较之其它两种形式，压缩波是最重要的考虑因素。

这种压缩波在介质中的传递是通过介质的压缩和膨胀进行的，但这种介质质点在声波作用下以原始位置为原点所发生的振荡仍服从Hooke 定律，换言之，介质的结构在声波传递过程中未发生任何根本性的破坏。

超声波检测技术最为常用的两个测量参数是通过介质的声速和振幅衰减。

超声波检测技术分为连续式和脉冲式。

前者操作较为复杂，仪器和技术要求较高，测量精度也相应较高，主要使用在一些专门的研究场所。

工业上比较实用的是脉冲式技术，其优点是操作简单，快速，易于自动化，测量精度对农产品品质检测也是足够的。

3 西瓜超声检测系统超声检测西瓜成熟度的系统见图1图1 超声检测系统简图接收 探头发射 探头计算机高频脉冲发射器样本宽频 接收器 数字存储 示波器绘图仪此装置主要包括信号发生装置（高频脉冲发射装置，超声频率为25kHz ）、信号接收装置（宽频接收器）、信号处理装和数据分析处理系统（个人电脑，数据采集控制系统）。

软件部分采用C 语言编程设计数据采集窗口，输入编号和信号发射点，与高频脉冲发射装置连接，从而完成数据采集。

采用MATLAB 完成数据处理。

4 实验设计4.1西瓜选样样本选自产于陕西西安地区的西瓜，共30个。

选取成熟度分别为60%，80%和100%的西瓜样本。

每个西瓜编号后，称质量，测量西瓜的质量，半径。

西瓜梗部和蒂部2个点分别标号1，2，赤道上等距标定4个点3，4，5，6，其中3号点为靠地点，5号点向阳。

西瓜标号示意图见下图。

图2 西瓜标号示意图图2 西瓜示意图4.2超声信号采集将发射探头按标号对应贴在西瓜上，记录试验日期、样品编号、发射点、接收点、试验次数、压电式加速度传感器编号等信息。

声波信号由计算机记录在硬盘上。

按编号顺序依次单独发射，在每一个点发射5次，其余编号的每个点单独接收一次。

4.3超声特征参数提取超声波在介质中传递过程复杂，有阻尼振动情况下，振幅的衰减很快，可用于表征声波在介质中传播情况。

本研究发射—接收点组合ij 间的透过率计算方法：①采用傅里叶变换将时域信号转换成频域信号，得到信号幅频谱；②对发射点和接收点幅频谱上每个频率（下称频率点）对应的幅值采用公式（1）计算透过率：1325赤道 4（6）梗部蒂部δij=Aj/Ai（1）（1）式中：δij为点i发射，点j接收时点组合ij间的声透过率；Ai为接收点i的幅值Aj为发射点j的幅值表1 为超声检测西瓜的超声透过率表1 超声检测西瓜的超声透过率ijδij ijδij12 0.6822 41 0.702513 0.6847 42 0.728914 0.7176 43 0.891015 0.8435 45 0.711816 0.7534 46 0.758921 0.7417 51 0.737823 0.8356 52 0.772324 0.8630 53 0.828325 0.7986 54 0.801526 0.7963 56 0.789531 0.7044 61 0.792032 0.7750 62 0.788434 0.8019 63 0.818535 0.8640 64 0.808836 0.8090 65 0.7998由表1可见，由1号点发射脉冲时，赤道部位（3，4，5，6）接收信号，获得的超声透过率建立模型，预测值与真实值之间相关系数比2号点大；由2号点发射脉冲时，根据数据可以得出一个结论，赤道部位接收所获模型相关性比1号点大；由赤道部位发射脉冲时，1号和2号接收点相关系数相对赤道部位接收信号要低。

由此可见，信号接收点最好选择在赤道部位。

同为赤道部位接收信号，选择较好的信号发射点。

比较发现，由赤道部位发射脉冲，获得的相关系数都较高。

其中由4号点发射脉冲时，各组合获得的相关系数普遍都高。

所以最后选定脉冲信号发射点和信号接收点都在赤道部位能达到更好的坚硬度预测效果。

分析出现这种状况的原因：①在1号发射或接收脉冲信号时，由于梗部像橄榄球两端，与小球接触面积太小，产生信号及感应信号随机性大，获得检测效果稍差；②在2号及赤道部位，信号接收和感应面积较大，能充分反映西瓜的状况。

并且瓜瓤在泛瓤期间，是横向增长的，赤道部位是整个西瓜瓜瓤主体，能较好反映西瓜的成熟状况。

4.4西瓜硬度检测超声波的衰减公式：l e A A X α-=0（2）式中，A x 为检测西瓜是的振幅，A 0为不检测西瓜，空载时的振幅，e 为自然对数，α为衰减系数，不同的介质中衰减系数是不同的，L 为测距（西瓜赤道直径）。

而衰减系数α由3个部分组成：42cf bf af ++=α（3）a 代表介质弹性摩擦吸收系数，与频率的一次方成正比。

b 代表介质粘滞性与热传导的吸收系数，与频率的二次方成正比。

c 代表介质内散射体的瑞利散射吸收系数，与频率的四次方成正比。

由于衰减系数与实际的声强衰减成正比，因此只需计算衰减系数即可反映超声波的衰减变化情况。

在实际数据处理中，只需用（2）式计算即可。

（2）式经化简可得：α=8.686/L式中L 为西瓜赤道直径。

然后利用Matlab 对超声在西瓜的传播过程中的衰减特性进行仿真。

根据前人研究的陈国发现，由于超声波的衰减与西瓜的坚硬度有着直接的联系，当超声衰减增加的同时西瓜的硬度会减小。

西瓜坚硬度FE(N)与超声衰减A （db/mm ）之间有较好的相关性[12]。

其关系为：25.1334.73167.67315.18123-+-=A A A FE (4)根据（4）式超声衰减与西瓜坚硬度的关系，利用Matlab 进行仿真。

仿真图见下图。

图2 西瓜坚硬度随超声波衰减变化的情况 表2 超声衰减与西瓜坚硬度的数值5 数据分析首先，利用Matlab 中相关性函数分析1024个衰减值与坚硬度间的相关性，。

然后按照需要，比如选取相关系数前3位的衰减值，借助软件中多元线性回归函数建立3个衰减值与坚硬度间的三元线性回归方程。

利用matlab 对（2）式进行仿真，得到下图。