Pig state audio recognition based on underdetermined blind source separation and deep learning
-
摘要:目的
为解决群养环境下生猪音频难以分离与识别的问题，提出基于欠定盲源分离与ECA-EfficientNetV2的生猪状态音频识别方法。
方法以仿真群养环境下4类生猪音频信号作为观测信号，将信号稀疏表示后，通过层次聚类估计出信号混合矩阵，并利用lp范数重构算法求解lp范数最小值以完成生猪音频信号重构。将重构信号转化为声谱图，分为进食声、咆哮声、哼叫声和发情声4类，利用ECA-EfficientNetV2网络模型识别音频，获取生猪状态。
结果混合矩阵估计的归一化均方误差最低为3.266×10−4，分离重构的音频信噪比在3.254~4.267 dB之间。声谱图经ECA-EfficientNetV2识别检测，准确率高达98.35%；与经典卷积神经网络ResNet50和VGG16对比，准确率分别提升2.88和1.81个百分点；与原EfficientNetV2相比，准确率降低0.52个百分点，但模型参数量减少33.56%，浮点运算量(FLOPs)降低1.86 G，推理时间减少9.40 ms。
结论基于盲源分离及改进EfficientNetV2的方法，轻量且高效地实现了分离与识别群养生猪音频信号。
Abstract:ObjectiveIn order to solve the problem of difficult separation and recognition of pig audio under group rearing environment, we propose a method of pig state audio recognition based on underdetermined blind source separation and ECA-EfficientNetV2.
MethodFour types of pig audio signals were simulated as observation signals in group rearing environment. After the signals were sparsely represented, the signal mixing matrix was estimated by hierarchical clustering, and the lp-paradigm reconstruction algorithm was used to solve for the minimum of lp-paradigm to complete the reconstruction of pig audio signals. The reconstructed signals were transformed into acoustic spectrograms, which were divided into four categories, namely, eating sound, roar sound, hum sound and estrous sound. The audio was recognized using the ECA-EfficientNetV2 network model to obtain the state of the pigs.
ResultThe normalized mean square error of the hybrid matrix estimation was as low as 3.266×10−4, and the signal-to-noise ratios of the separated reconstructed audio ranged from 3.254 to 4.267 dB. The acoustic spectrogram was recognized and detected by ECA-EfficientNetV2 with an accuracy of up to 98.35%, and the accuracy improved by 2.88 and 1.81 percentage points compared with the classical convolutional neural networks ResNet50 and VGG16, respectively. Compared with the original EfficientNetV2, the accuracy decreased by 0.52 percentage points, but the amount of the model parameters reduced by 33.56%, the floating-point operations (FLOPs) reduced by 1.86 G, and inference time reduced by 9.40 ms.
ConclusionThe method based on blind source separation and improvement of EfficientNetV2 lightly and efficiently realizes separating and recognizing audio signals of group-raised pigs.
-
保鲜环境温湿度调控是延长果蔬保鲜周期的有效方式之一[1-2]。然而，温度超调量过大可能会导致冷害[3]，环境湿度过高会加速微生物生长[4]，保鲜环境参数波动震荡也会影响果蔬保鲜效果[5]。因而，提高保鲜环境控制系统的控制性能对于保障果蔬储运品质有重要作用。我们前期开发了基于双限值的果蔬保鲜控制系统，该系统能实时调控保鲜环境参数，但存在环境参数超调量大、波动严重、抗干扰能力不强等问题。因此，有必要开展果蔬保鲜环境参数调控策略优化试验，提高系统控制性能。
国内外学者对环境参数的调控进行了一系列研究，王广海等[6]以双限值作为系统控制算法，控制系统温度、相对湿度、氧气和二氧化碳等环境参数。虽然双限值控制逻辑简单，但存在超调量大、环境参数频繁波动等问题[7]。Barros 等[8]设计了基于比例−积分−微分（Proportional-integral-derivative，PID）控制器的温度加热系统。传统PID控制结构易于实现，在工业中最为常用[9]，然而传统PID需要人为调整控制器参数，缺少自适应能力，鲁棒性较低[10]，难以实现环境参数的高质量控制。赵鑫鑫[11]设计了车厢温度模糊规则，应用于冷藏车箱温度控制。相较于传统PID，超调量有所降低，但仍存在自适应能力弱、抗干扰能力较差等问题。近年来，人工智能算法发展迅速，许多学者对自适应控制技术进行了研究。Aftab等[12]提出采用基于Lyapunov函数的人工神经网络对PID控制器的比例、积分和微分项进行在线整定，相较于传统PID取得较好的跟踪效果，但系统的每个控制参数都需使用一个独立的神经网路，计算量大且复杂。Salcedo等[13]设计了一种结合Smith预估器的状态反馈控制器，其控制性能和精度高于传统PID控制器，但需要建立控制对象动态行为的数学模型，该模型缺乏自适应能力，控制环境或对象发生变化时，控制精度将大幅下降。Silveira等[14]设计了LPPT控制算法用于制冷系统，通过2个具有自适应能力的非零开关量，实现降低系统能耗的同时具备较强的抗干扰能力，但与传统PID相比该系统温度有明显波动，稳态误差较大，无法保持恒定。
目前，针对果蔬保鲜环境参数智能调控方面研究较少，大多数冷藏、冷冻等控制系统仍采用传统双限值和PID控制技术[15]。因此，本文将根据传统果蔬保鲜环境控制系统超调量大、控制精度低、波动严重、鲁棒性差等现状，结合神经网络和PID控制，设计一种基于BP神经网络-PID（Back-propagation neural network-PID，BPNN-PID）的控制策略，使PID控制器具备自学习、自适应能力，确保果蔬保鲜环境控制系统在运行过程中处于优化状态，从而提高系统控制性能，为果蔬保鲜环境参数调控的实际应用提供一定参考。
1. 果蔬保鲜系统的基本组成
如图1所示，果蔬保鲜环境调控试验平台[6]由制冷系统、加湿系统和控制系统3个部分组成。制冷系统主要由变频器(型号为SC-650，由河南尚川电子科技有限公司生产)、变频压缩机(型号为DTH356LDPC9FQ，由上海日立电器有限公司生产)、直流蒸发风机(额定电压24 V、额定电流8 A)、冷凝风机(额定电压220 V，额定电流0.42 A)和制冷管路组成。果蔬保鲜环境控制系统通过变频器调节变频器压缩机工作频率改变压缩机的转速，再利用直流蒸发风机带动气流，穿过蒸发盘管形成冷空气，经开孔隔板(开孔率16.11%)[16]均匀流入保鲜室，降低货物温度后从回风道返回，实现果蔬保鲜环境温度控制。
![]()
图 1 果蔬保鲜运输车箱体结构示意图1：冷凝器；2：冷凝风机；3：蒸发风机；4：回风道；5：保鲜室；6： 传感器盒；7：货物；8：支撑架；9：气流轨道；10：排水管；11：三通接头；12：积水槽；13：温度传感器；14：超声波雾化器；15：开孔隔板；16：加湿风机；17：蒸发器；18：补水箱；19：排水阀；20：变频器；21：变频压缩机；22：制冷管路；23：压力室；24：控制器；25：触摸屏；26：电子计算机Figure 1. Schematic diagram of the box structure of the transport vehicle for fruit and vegetable fresh-keeping1: Condenser; 2: Condensing fan; 3: Evaporating fan; 4: Return air duct; 5: Keep-freshing room; 6: Sensor box; 7: Cargo; 8: Support frame; 9: Air flow track; 10: Drain pipe; 11: Tee connector; 12: Standing water tank; 13: Temperature sensor; 14: Ultrasonic atomizer; 15: Open partition; 16: Humidifying fan; 17: Evaporator; 18: Water replenishment tank; 19: Drain valve; 20: Inverter; 21: Inverter compressor; 22: Refrigeration line; 23: Pressure chamber; 24: Controller; 25: Touch screen; 26: Electronic computer加湿系统主要由水槽、超声波雾化器(额定电压24 V、额定电流1 A)和加湿风机(额定电压24 V、额定电流0.55 A)组成。超声波雾化器置于液面下2~3 cm处，将液态水雾化形成微小雾粒，通过加湿风机吹进保鲜区，提高保鲜区环境相对湿度。
控制系统总体结构如图2所示，该系统以Cortex-M3架构的STM32F103C8T6微处理器为核心，配合数据采集模块、执行机构、触摸屏、PC电子计算机、电路保护装置等对果蔬保鲜环境进行控制。
![]()
图 2 果蔬保鲜环境控制系统结构框架图Figure 2. Structural framework diagram of the environmental control system for fruit and vegetable fresh-keeping2. 控制系统硬件与软件设计
本研究将设计一款基于STM32F103C8T6的果蔬保鲜环境调控系统，该系统微处理器工作频率最高为72 MHz, 内置高达128 K字节的Flash存储器和20 K字节的SRAM, 具备足够的通用I/O端口[17]，满足系统开发需求。
2.1 硬件电路
2.1.1 电源电路
电源电路如图3所示，考虑到本系统的运行环境复杂，信息通讯频繁且运行时间长等一系列问题，电路设计需降低干扰源、阻断耦合以及提高敏感设备的阈值。系统器件需要的电压等级分别为24.0、5.0和3.3 V，其中24.0 V由外部开关电源提供，隔离5.0 V电源电路采用DC-DC电源隔离模块B0505S-1WR2，为系统提供信号电源，3.3 V的电源采用AMS1117-3.3电压转换芯片。
![]()
图 3 果蔬保鲜环境控制系统的电源电路Figure 3. Power supply circuit of the environmental control system for fruit and vegetable fresh-keeping2.1.2 CAN总线驱动电路
控制器局域网总线(Controller area network，CAN)驱动电路如图4所示，ISO1050是一款将隔离通道和CAN收发器集成在一个封装内的隔离型CAN总线收发器。与隔离式电源一起使用，可防止数据总线或者其他电路上的噪音电流进入本地接地而产生的干扰或损坏敏感电路。为抵消电信号的反射，CAN总线输出两端增加了1个120 Ω的终端电阻。
2.1.3 RS485接口电路
RS485的通信电路如图5所示，MCU根据Modbus-RTU通讯协议，通过RS485收发器与从机进行异步通讯。SP3485芯片是3.3 V 低功耗半双工的收发器，将其RO及DI引脚分别与USART的RX和TX引脚相连，将RE和DE引脚直接用普通IO口来控制数据传输方向，采用轮询发送和中断接收的数据传输方式。A与B之间接120 Ω电阻避免信号发射问题。
2.1.4 中间继电器驱动电路
中间继电器驱动电路如图6所示，当NPN型三极管导通时，继电器吸合，并联在继电器两端的发光二极管被点亮，表明继电器正在工作。与线圈并联的续流二极管可以吸收线圈断电时产生的感应电动势，防止晶体管被击穿。继电器输出端串联1个保险丝，当电路出现温度异常时可以迅速切断电路。
2.1.5 数据采集模块
数据采集模块由温度传感器(型号：TH600NXC，供电：10~30 V，精度：±0.3 ℃，范围：−40~80 ℃，通讯方式：RS485)、相对湿度传感器(型号：TH600NXT，供电：10~30 V，精度：±3%，范围：0%~100%，通讯方式：RS485)组成。各传感器单独作为1个节点接入RS485总线，微处理器利用Modbus-RTU通讯协议定时采集数据，数据经过微处理器解析后得到实时的保鲜室环境参数。
2.2 软件设计
果蔬保鲜控制系统软件设计主要由硬件驱动程序、数据采集及处理模块、控制算法程序、历史数据存储模块和触摸屏驱动程序5个部分组成。
硬件驱动程序使控制系统及各执行器正常运行；数据采集及处理模块对箱内传感器信号进行采集和处理；控制算法程序将数学模型转换为机器语言，实现对果蔬保鲜环境参数的自动控制；历史数据存储模块使处理器将采集并处理后的环境数据以txt文件格式保存于SD卡中，方便用户后期对果蔬保鲜环境数据进行分析；触摸屏驱动程序使用户通过人机交互页面对果蔬保鲜系统下达控制指令、设置控制参数，并实时显示箱内环境参数。配合BPNN-PID控制算法得到基于BPNN-PID的果蔬保鲜环境控制策略，系统运行程序流程如图7所示。
![]()
图 7 基于BPNN-PID的果蔬保鲜环境控制策略Figure 7. Environmental control strategy of fruit and vegetable fresh-keeping based on BPNN-PID3. 传统PID和BPNN-PID控制算法原理及实现
3.1 常规PID原理
常规PID控制系统，由PID控制器和被控对象组成[18]。PID控制器是一种线性控制器，它根据给定值r(k)与被控对象输出值y(k)构成控制偏差e(k)，通过对系统当前状态的监测和反馈，实现对系统行为的控制。
$$ e\left( k \right) = r\left( k \right) - y\left( k \right) ，$$ (1) 式中：k为采样序号；r(k) 为系统期望值；y(k) 为系统实际值；e(k) 为系统偏差。
PID控制算法的核心是3个参数：比例、积分和微分系数，它们分别控制着系统的响应速度、稳定性和抗干扰性能。其中，比例环节反应了控制系统的偏差信号e(k)，偏差一旦产生，控制器立即反应，以减小偏差。积分环节主要用于消除静态误差，提高系统的误差度。微分环节主要反映偏差信号的变化趋势，并在偏差信号值过大之前，在系统中引用一个有效的早期修正信号，从而减小系统震荡。常规PID控制又分为位置式PID与增量式PID，两者的表达式分别为[19]：
$$ u(k) = {K_{\rm{P}}}e(k) + {K_{\rm{I}}}\displaystyle\sum\limits_{k = 0}^k {e(k) + {K_{\rm{D}}}[e(k) - e(k - 1)]} ，$$ (2) $$ \Delta u(k) ={K_{\rm{P}}}[e(k) - e(k - 1)] + {K_{\rm{I}}}e(k)+ $$ $$ {K_{\rm{D}}}[e(k) - 2e(k - 1) + e(k - 2)] ，$$ (3) 式中：u(k)为第k次采样时刻的计算机输出值；KP 为比例系数；KI为积分系数；KD为微分系数。
由位置式PID表达式可知，控制器每次输出都与过去的状态有关，导致计算机运算量过大。如MCU出现故障，输出会大幅度变化，这种情况在实际生产中是不允许的。而增量式PID只与过去2次的状态有关，大大增加了系统的容错率，因此本系统选择增量式PID作为控制器的基础。
3.2 BPNN-PID原理
3.2.1 BPNN-PID控制器设计
常规的果蔬保鲜环境温度PID控制器无法根据被控对象等因素自适应调整控制参数，在PID控制器的基础上增加BP神经网络，构成一个具有自适应能力的果蔬保鲜环境温度BPNN-PID控制器，如图8所示。
![]()
图 8 果蔬保鲜环境温度BPNN-PID控制器C：外部偏置常量；r(k) ：系统期望值；y(k) ：系统实际值；e(k) ：系统偏差；u(k)：计算机输出值；KP：比例系数；KI ：积分系数；KD：微分系数Figure 8. Ambient temperature BPNN-PID controller for fruit and vegetable fresh-keepingC: External bias constant; r(k): System expected value; y(k): Actual value of the system; e(k): System deviation; u(k):Computer output; KP: Proportional coefficient; KI: Integration coefficient; KD: Differential coefficient3.2.2 BPNN-PID控制器自适应过程
本文设计的果蔬保鲜环境控制系统的BP神经网络输入层由目标温度、实际温度、温度误差和外部偏置常量共4个神经元组成[20]，由于控制参数的取值范围为0~1，并根据前期试验进行调试，确定当外部偏置常量C=1时控制效果最好。输出层由比例、积分、微分系数组成，共3个神经元，再经过调试后确定隐藏层一共5个神经元，最终建立的BP神经网络结构如图9所示。
![]()
图 9 BP神经网络结构x1：目标温度；x2：实际温度；x3：温度误差；x4：外部偏置常量；l：输入层神经元序号；m隐含层神经元序号；n：输出层神经元序号；KP ：比例系数；KI：积分系数；KD：微分系数Figure 9. BP neural network structurex1: Target temperature; x2: Actual temperature; x3: Temperature error; x4: External bias constant; l: Input layer neuron number; m: Hidden layer neuron number; n: Output layer neuron number; KP: Proportional coefficient; KI: Integration coefficient; KD: Differential coefficient1)误差正向传播：根据图9所示的神经网络结构图可得，输入层的输出
$O_l^{(1)} $ 为$$ O_{l}^{(1)}=x_{l}, \;l=1,2,3,4 \; ,$$ (4) 式中：l为输入层神经元序号；xl为输入层第l个输入。
隐含层的输入和输出分别为
$$ \left\{\begin{array}{l}{\rm{n e t}}_m^{(2)}(k)=\displaystyle\sum_{l=0}^4 w_{l m}^{(2)} O_l^{(1)} \\ O_m^{(2)}(k)=f\left[ { {\rm{net}} }_m^{(2)}(k)\right]\end{array}, \;m=1,2,3,4,5\; ,\right. $$ (5) 式中：m为隐含层神经元序号；
${\rm{net}}_m^{(2)} $ 为隐含层第m个神经元输入；${w_{lm}^{(2)}}$ 为隐含层权值系数；$O_m^{(2)} $ 为隐含层第m个神经元输出。输出层的输入和输出分别为
$$ \left\{\begin{array}{l}{\rm{n e t}}_n^{(3)}(k)=\displaystyle\sum_{m=0}^5 w_{m n}^{(3)} O_m^{(2)} \\ O_n^{(3)}(k)=g\left[{\rm{n e t}}_n^{(3)}(k)\right]\end{array},\;n=1,2,3\right. ,$$ (6) 式中：n为输出层神经元序号；
${\rm{net}}_n^{(3)} $ 为输出层第n个神经元输入；${w_{mn}^{(3)}}$ 为输出层权值系数；$O_n^{(3)} $ 为输出层第n个神经元输出。隐含层的激活函数
$ f(x) $ 采用正负对称的Sigmoid函数，输出层激活函数$ g(x) $ 采用非负的Sigmoid函数。$$ f(x)=\dfrac{{\text{e}}^{x}-{\text{e}}^{-x}}{{\text{e}}^{x}+{\text{e}}^{-x}}\text{，}g(x)=\dfrac{1}{1+{\text{e}}^{-x}} 。$$ (7) 2)误差反向传播：本系统设计所用的神经网络的输出性能指标函数为
$$ E(k)=\dfrac{1}{2}\left[r(k)-y(k)\right]^{2}。 $$ (8) 引入学习效率η和惯性系数α，根据式(4)~(8)通过梯度下降法，得到输出层和隐藏层的权重更新和误差项，经调试取η=0.1，α=0.25，系统控制效果最佳。
$$ \left\{\begin{array}{l}\Delta w_{m n}^{(3)}(k)=0.1 \Delta w_{m n}^{(3)}(k-1)+0.25 \delta_n^{(3)} O_n^{(2)}(k) \\ \delta_n^{(3)}=\operatorname{sgn}\left[\dfrac{\partial \operatorname{yout}(k)}{\partial u(k)}\right] e(k) \dfrac{\partial u(k)}{\partial O_l^{(3)}(k)} g^{\prime}\left[{\rm{n e t}}_n^{(3)}(k)\right]\end{array}\right., $$ (9) 式中：
$ \delta_{n}^{(3)} $ 为输出层神经元误差项。$$ \left\{\begin{array}{l}\Delta w_{l m}^{(2)}(k)=0.1 \Delta w_{l m}^{(2)}(k-1)+0.25 \delta_m^{(2)} O_l^{(1)}(k) \\ \delta_m^{(2)}=f^{\prime}\left[{\rm{n e t}}_m^{(2)}(k)\right] \displaystyle\sum_{n=1}^3 \delta_n^{(3)} w_{m n}^{(3)}(k)\end{array}\right., $$ (10) 式中：
$ \delta_{m}^{(2)} $ 为隐藏层神经元误差项。由于
$\dfrac{\partial {\rm{y o u t}}(k)}{\partial u(k)}$ 无法直接计算，用$\operatorname{sgn}\left[\dfrac{\partial {\rm{yout}}(k)}{\partial u(k)}\right]$ 代替[21]。以上为果蔬保鲜环境BPNN训练1次的全部过程。4. 试验与结果分析
研究发现，果蔬保鲜环境温度与相对湿度存在较强的耦合关系[22]，且环境温度对厚表皮类果蔬品质的影响更为明显[23]。因此本研究在温度控制方面采用BPNN-PID控制算法，实现对温度的精准控制；湿度方面搭配双限值控制算法，组成BPNN-PID控制策略，可满足果蔬保鲜的基本环境要求。
本研究通过自主搭建的果蔬保鲜系统，采用常规PID和BPNN-PID 2种不同控制策略进行果蔬保鲜试验，分析果蔬保鲜系统在不同控制策略下环境参数的超调量、稳定时间和稳态误差。根据试凑法[24]确定常规PID控制器控制参数KP=0.8、KI=0.75、KD=0.5，果蔬保鲜环境BPNN-PID控制系统的隐含层、输出层的初始权重均为−0.5~0.5的随机值，根据前期试验，学习率η=0.1、惯性系数α=0.25时控制性能最优。
试验材料采用赣南脐橙，总质量40 kg，购于水果批发市场，果实新鲜，无病虫害，无机械损伤。脐橙采用塑料筐装放，塑料筐呈网状结构，尺寸规格(长×宽×高)为600 mm×425 mm×350 mm。脐橙以每框20 kg装箱，共2框放置保鲜室，外界环境温度(26±1) ℃，保鲜室初始温度(26±1) ℃，据相关文献[25]，脐橙适宜的保鲜参数为：温度为6 ℃，相对湿度为80%~90%。因此设置果蔬保鲜环境目标温度为6 ℃、目标湿度为80%~90%。为分析果蔬保鲜系统在不同控制策略下的抗干扰能力，在实载试验中，待果蔬保鲜系统运行稳定(240 min)后，对系统施加干扰(模拟实际取货时打开箱门5 min)，每5 min记录1次数据，运行时间为6 h。
试验完成后计算系统的超调量、稳态误差和稳定时间。超调量计算公式如下：
$$ \sigma = {T_{{\rm{max}}}}(t) - {T_{{\rm{out}}}}(\infty ) ，$$ (11) 稳态误差计算公式如下：
$$ {e_{{\rm{ss}}}} = {T_{{\rm{out}}}}(t) - {T_{{\rm{out}}}}(\infty ) ，$$ (12) 稳定时间评判公式如下：
$$ \dfrac{{|{T_{{\rm{out}}}}(t) - {T_{{\rm{out}}}}(\infty )|}}{{|{T_{{\rm{out}}}}(0 + ) - {T_{{\rm{out}}}}(\infty )|}} \leqslant \Delta ，$$ (13) 式中：σ为超调量；ess为稳态误差；Tmax(t)为环境参数最大值；Tout(t)为环境参数瞬态值；Tout(0+)为环境参数初始值；Tout(∞)为环境参数稳态值；Δ为稳定时间评判系数，取2%[26]。
4.1 空载试验与结果分析
4.1.1 空载常规PID控制策略试验
空载状态下，果蔬保鲜环境系统采用常规PID控制策略的控制响应曲线如图10所示，环境温度的超调量为3.5 ℃、稳定时间为100 min、稳态误差为±0.2 ℃；相对湿度的超调量为9.5%、稳定时间为100 min。
![]()
图 10 基于常规PID控制策略的空载保鲜环境参数变化曲线Figure 10. Variation curve of no-load fresh-keeping environment parameters based on conventional PID control strategy从该试验结果可以看出，以STM32F103C8T6微处理器为核心的控制系统利用RS485通讯方式可实时采集箱内环境参数；利用CAN总线可实时调节变频压缩机频率从而控制系统制冷量；通过继电器驱动电路控制环境相对湿度。系统按照设定的控制流程执行，实现了果蔬保鲜环境参数的自动控制。因此，本文设计的控制系统软硬件满足果蔬保鲜要求。
4.1.2 空载BPNN-PID控制策略试验
空载状态下，果蔬保鲜环境系统采用BPNN-PID控制策略的控制响应曲线如图11、12所示，果蔬保鲜环境系统在运行过程中通过自学习、自适应能力动态调整控制参数，最终收敛于KP=0.588，KI=0.666，KD=0.617，环境温度的超调量为2.5 ℃、稳定时间为70 min、稳态误差为±0.1 ℃；相对湿度的超调量为4.3%、稳定时间为70 min。
![]()
图 11 基于BPNN-PID控制策略的空载保鲜环境参数变化曲线Figure 11. Variation curve of no-load fresh-keeping environment parameters based on BPNN-PID control strategy空载试验中，与常规PID控制策略相比，基于BPNN-PID控制策略的果蔬保鲜系统环境温度超调量减小了1 ℃、稳态误差减小了0.1 ℃、稳定时间缩短了30 min；环境相对湿度超调量减小了5.2%、稳定时间缩短了30 min。由此可知，采用BPNN-PID控制策略的果蔬保鲜系统控制性能均有明显提升。
4.2 实载试验与结果分析
4.2.1 实载常规PID控制策略试验
由图13可知，实载状态下，果蔬保鲜环境温度控制采用常规PID控制算法，相对湿度控制采用双限值控制算法时，温度的超调量为3.8 ℃、稳定时间为105 min、稳态误差为±0.5 ℃；相对湿度的超调量为5%、稳定时间为80 min。
![]()
图 13 基于常规PID控制策略的实载保鲜环境参数变化曲线Figure 13. Variation curves of real-load fresh-keeping environment parameters based on conventional PID control strategy与双限值控制策略相比，温、湿度控制精度有所提高，温、湿度波动减小，但仍存在超调量较大、稳定时间慢等问题。针对干扰，箱内环境温、湿度经过105 min的上下波动才重新恢复稳定。因此常规PID控制策略缺少动态调节能力，存在抗干扰能力不强、适应性较弱等缺陷，系统控制性能仍有提升空间。
4.2.2 实载BPNN-PID控制策略试验
由图14、15可知，实载状态下，果蔬保鲜环境温度控制采用BPNN-PID控制算法，相对湿度控制采用双限值控制算法时，果蔬保鲜环境系统控制参数最终收敛于KP=0.717，KI=0.682，KD=0.656，温度超调量为1.7 ℃，稳定时间为80 min，稳态误差为±0.2 ℃，相对湿度超调量为2.8%、稳定时间为55 min。
![]()
图 14 基于BPNN-PID控制策略的实载保鲜环境参数变化曲线Figure 14. Variation curves of real-load fresh-keeping environment parameters based on BPNN-PID control strategy实载试验中，与常规PID控制策略相比，基于BPNN-PID控制策略的果蔬保鲜系统环境温度超调量减小了2.1 ℃、稳态误差减小了0.3 ℃、稳定时间缩短了25 min；环境相对湿度超调量减小了2.2%、稳定时间缩短了25 min。BPNN-PID控制策略在超调量、稳定时间、稳态误差和控制精度等方面都有了较大的提升。针对干扰，BPNN-PID控制策略有效地抑制了保鲜室内环境参数的剧烈波动，经过80 min的动态调节后，最终收敛于目标参数。因此BPNN-PID控制策略的抗干扰能力更强、鲁棒性优异、自适应性好，具有良好的动态调节能力，其控制性能明显优于常规PID控制策略的，能够更有效地完成果蔬保鲜环境的控制工作。
5. 结 论
本研究根据果蔬保鲜运输的温、湿度控制要求，以STM32为核心处理器，设计了果蔬保鲜环境控制系统硬件和软件，并通过试验对比常规PID和BPNN-PID 2种控制策略的环境参数调控效果，研究结果对于果蔬保鲜环境参数调控有重要意义。经研究获得以下结论：
1)所开发的控制系统能准确采集保鲜室各环境参数，并根据控制策略准确控制相应的执行机构，实现保鲜环境参数的实时控制。
2)基于BPNN-PID控制策略的果蔬保鲜环境控制系统具有响应速度快、超调量小、控制精度高、抗干扰能力强等优点，系统控制性能明显提升。
3)保鲜环境温、湿度调控过程中，温度和相对湿度具有较强的耦合关系，其中温度占主导地位，相对湿度的波动趋势与温度基本一致，随着系统温度控制性能的提高，相对湿度的控制性能有了明显的改善。
值得注意的是，加湿器性能、果蔬生理特性、堆码方式等因素也可能会对研究结果产生一定的影响，这些因素将在后续研究中进行探讨。
-
![]()
图 2 基于多窗谱估计的改进谱减法算法的流程图
Figure 2. Flow chart of improved spectral subtraction algorithm based on multi-window spectral estimation
![]()
图 4 不同状态的生猪音频声谱图示例
图中每个坐标点的颜色代表此时间点此频率上的能量大小，颜色越深，能量越高
Figure 4. Examples of audio spectrograms of pigs in different states
The color of each coordinate point in the graph represents the energy magnitude at this frequency at this time point，and the darker the color, the higher the energy
![]()
图 5 同音频下宽带声谱图与窄带声谱图
图中每个坐标点的颜色代表此时间点此频率上的能量大小，颜色越深，能量越高
Figure 5. Wideband and narrowband spectrograms under the same audio
The color of each coordinate point in the graph represents the energy magnitude at this frequency at this time point，and the darker the color, the higher the energy
![]()
图 6 盲源分离处理模型
s(t)为源信号矢量，F( · )为混合系统，v(t)为噪声矢量，x(t)为观测信号，H( · )为解混系统， $ \hat {\boldsymbol{S}}{\text{(}}t{\text{)}} $为解混信号，即重构后的信号
Figure 6. Blind source separation processing model
s(t) is the source signal vector, F( · ) is the hybrid system, v(t) is the noise vector, x(t) is the observed signal, H( · ) is the demixed system, $ \hat {\boldsymbol{S}}{\text{(}}t{\text{)}} $ is the demixed signal, that is, the reconstructed signal
![]()
图 7 欠定盲源分离流程图
x(t)为观测信号，X(t)为稀疏表示后的稀疏向量，$ \hat {\boldsymbol{A}} $为混合矩阵估计值，$ \hat {\boldsymbol{S}}{\text{(}}t{\text{)}} $为分离出的信号，s(t)为最终重构出的信号
Figure 7. Flow chart of underdetermined blind source separation
x(t) is the observed signal, X(t) is the sparse vector after sparse representation, $ \hat {\boldsymbol{A}} $ is the mixing matrix estimate, $ \hat {\boldsymbol{S}}{\text{(}}t{\text{)}} $ is the separated signal, and s(t) is the final reconstructed signal
![]()
图 8 ECA-EffficientNetV2网络模型结构图
Figure 8. Structure of the ECA-EffficientNetV2 network model
![]()
图 10 能熵比法端点检测波形图
红色实线和虚线分别为当前有声段起始点和终止点
Figure 10. Waveforms of energy entropy ratio endpoint detection
The solid and dashed red lines are the start point and end point of the current audible segment respectively
![]()
图 11 低信噪比下的能熵比端点检测波形图
红色实线和虚线分别为当前有声段起始点和终止点
Figure 11. Waveforms of energy entropy ratio endpoint detection at low signal-to-noise ratio
The solid and dashed red lines are the start point and end point of the current audible segment respectively
![]()
图 14 单源点检测前后的观测信号散点图
Figure 14. Scatter plots of observation signals before and after single-source point detection
![]()
图 17 训练集与验证集的准确率和损失值对比曲线
Figure 17. Comparison curves of accuracy and loss value between training set and verification set
表 1 ECA-EfficientNetV2架构表
Table 1 ECA-EfficientNetV2 architecture table
阶段
Stage操作
Operator层数
Layers步长
Stride0 Conv3×3 1 2 1 Fused-MBConv1, k3×3 1 1 2 Fused-MBConv2, k3×3 4 2 3 Fused-MBConv2, k3×3 4 2 4 MBConv4, k3×3, ECA 5 2 5 MBConv6, k3×3, ECA 7 1 6 MBConv6, k3×3, ECA 12 2 7 Conv1×1, Pooling, FC 1 
下载: 导出CSV
表 2 谱减法与改进谱减法降噪效果对比
Table 2 Comparison of noise reduction effects of spectral subtraction and improved spectral subtraction
降噪前信噪比/dB
SNR before noise reduction谱减法降噪后信噪比/dB
SNR after spectral subtraction改进谱减法降噪后信噪比/dB
SNR after improved spectral subtraction for noise reduction0 2.31 4.67 5 5.97 7.42 15 15.21 15.58
下载: 导出CSV
表 3 不同模型在生猪声谱图上的预测性能比较
Table 3 Comparison of prediction performance of different models on pig spectrograms
模型
Model准确率/%
Accuracy精确率/%
Precision召回率/%
RecallF1分数/%
F1-score参数量/M
Parameter quantityFLOPs/
G推理时间/ms
Inference timeResNet50 95.47 95.35 95.56 95.45 25.56 8.21 11.16 VGG16 96.54 96.73 96.58 96.54 138.30 30.97 10.11 EfficientNetV2-S 98.87 98.69 98.86 98.85 21.46 5.73 34.34 ECA-EfficientNetV2 98.35 98.33 98.37 98.35 14.26 3.87 24.94
下载: 导出CSV
-
[1] 杨亮, 王辉, 陈睿鹏, 等. 智能养猪工厂的研究进展与展望[J]. 华南农业大学学报, 2023, 44(1): 13-23. doi: 10.7671/j.issn.1001-411X.202209050 [2] GHANI S H, KHAN W. Extraction of UAV sound from a mixture of different sounds[J]. Acoustics Australia, 2020, 48(3): 363-373. doi: 10.1007/s40857-020-00197-z
[3] HE P, QI M, LI W, et al. A general nonstationary and time-varying mixed signal blind source separation method based on online Gaussian process[J]. International Journal of Pattern Recognition and Artificial Intelligence, 2020, 34(11): 2058015. doi: 10.1142/S021800142058015X
[4] ADAM A M, FAROUK R M, EL-DESOUKY B S. Generalized gamma distribution for biomedical signals denoising[J]. Signal, Image and Video Processing, 2023, 17(3): 695-704. doi: 10.1007/s11760-022-02277-1
[5] 张振华. 基于音频分析技术的猪异常检测[D]. 太原: 太原理工大学, 2017. [6] 沈明霞, 王梦雨, 刘龙申, 等. 基于深度神经网络的猪咳嗽声识别方法[J]. 农业机械学报, 2022, 53(5): 257-266. doi: 10.6041/j.issn.1000-1298.2022.05.026 [7] JI N, SHEN W Z, YIN Y L, et al. Investigation of acoustic and visual features for pig cough classification[J]. Biosystems Engineering, 2022, 219: 281-293.
[8] 彭硕, 陶亮, 查文文, 等. 基于稀疏分量分析的生猪音频欠定盲源分离研究[J]. 畜牧兽医学报, 2023, 54(7): 2794-2809. [9] HU Y, LOIZOU P C. Speech enhancement based on wavelet thresholding the multitaper spectrum[J]. IEEE Transactions on Speech and Audio Processing, 2004, 12(1): 59-67. doi: 10.1109/TSA.2003.819949
[10] 杨稷, 沈明霞, 刘龙申, 等. 基于音频技术的肉鸡采食量检测方法研究[J]. 华南农业大学学报, 2018, 39(5): 118-124. doi: 10.7671/j.issn.1001-411X.2018.05.017 [11] 赵欢, 王纲金, 赵丽霞. 一种新的对数能量谱熵语音端点检测方法[J]. 湖南大学学报(自然科学版), 2010, 37(7): 72-77. [12] 朱国俊, 唐振博, 冯建军, 等. 启动方式对混流泵噪声特性的影响[J]. 农业工程学报, 2023, 39(13): 34-42. doi: 10.11975/j.issn.1002-6819.202304091 [13] 杜晓冬, 滕光辉, 刘慕霖, 等. 基于轻量级卷积神经网络的种鸡发声识别方法[J]. 农业机械学报, 2022, 53(10): 271-276. doi: 10.6041/j.issn.1000-1298.2022.10.029 [14] CAI W, LI M. A unified deep speaker embedding framework for mixed-bandwidth speech data [EB/OL]. arXiv: 2012.00486[2023-12-01]. https://ui.adsabs.harvard.edu/abs/2020arXiv201200486C.
[15] 丁楠. 基于特征学习的语音情感识别研究[D]. 南京: 南京邮电大学, 2023. [16] CANDAN C. An automated window selection procedure For DFT based detection schemes to reduce windowing SNR loss[EB/OL]. arXiv: 1710.10200[2023-12-01]. https://ui.adsabs.harvard.edu/abs/2017arXiv171010200C.
[17] 张逸, 陈书畅, 刘必杰, 等. 基于盲源分离的工业谐波源负荷分类识别方法[J]. 中国电机工程学报, 2024, 44(10): 3850-3862. [18] 邹亮, 张鹏, 陈勋. 基于三阶统计量的欠定盲源分离方法[J]. 电子与信息学报, 2022, 44(11): 3960-3966. doi: 10.11999/JEIT210844 [19] 凌康杰, 岳学军, 刘永鑫, 等. 基于移动互联的农产品二维码溯源系统设计[J]. 华南农业大学学报, 2017, 38(3): 118-124. doi: 10.7671/j.issn.1001-411X.2017.03.019 [20] 王晶, 李炜, 洪心睿, 等. 基于改进密度聚类算法的语音信号欠定盲分离[J]. 信息与控制, 2023, 52(6): 784-796. [21] KEMIHA M, KACHA A. Single-channel blind source separation using adaptive mode separation-based wavelet transform and density-based clustering with sparse reconstruction[J]. Circuits Systems and Signal Processing, 2023, 42(9): 5338-5357. doi: 10.1007/s00034-023-02350-1
[22] WANG M, CAI X X, ZHU K F. Underdetermined mixed matrix estimation of single source point detection based on noise threshold eigenvalue decomposition[C]//Proceedings of the 8th International Conference on Communications, Signal Processing, and Systems. Singapore: Springer, 2020: 704-711.
[23] LI Y B, NIE W, YE F, et al. A complex mixing matrix estimation algorithm in under-determined blind source separation problems[J]. Signal, Image and Video Processing, 2017, 11(2): 301-308. doi: 10.1007/s11760-016-0937-y
[24] 冯建英, 石岩, 王博, 等. 基于聚类分析的数据挖掘技术及其农业应用研究进展[J]. 农业机械学报, 2022, 53(S1): 201-212. doi: 10.6041/j.issn.1000-1298.2022.S1.021 [25] XIE Y, XIE K, XIE S L. Underdetermined blind separation of source using l p-norm diversity measures[J]. Neurocomputing, 2020, 411: 259-267.
[26] XU P F, JIA Y J, JIANG M X. Blind audio source separation based on a new system model and the Savitzky-Golay filter[J]. Journal of Electrical Engineering, 2021, 72(3): 208-212. doi: 10.2478/jee-2021-0029
[27] FISLI S, DJENDI M. Hybrid PSO-NLMS (HPSO-NLMS) algorithm for blind speech quality enhancement in time domain[J]. Applied Acoustics, 2021, 177: 107936.
[28] TAN M, LE Q V. EfficientNetV2: Smaller models and faster training[EB/OL]. arXiv: 2104.00298[2023-12-01]. https://ui.adsabs.harvard.edu/abs/2021arXiv210400298T.
[29] WANG X K, JIA X, ZHANG M Y, et al. Object detection in 3D point cloud based on ECA mechanism[J]. Journal of Circuits, Systems and Computers, 2023, 32(5): 2350080.
-
期刊类型引用(5)
1. 高旻. 视频分析技术在港口交通流量监控中的应用. 电视技术. 2025(01): 223-225 . 
百度学术
2. 陈钇果. 基于FPGA架构的人工智能加速器设计与实现. 软件. 2025(01): 59-61 . 百度学术
3. 沃英达. 基于深度学习算法的激光散斑音频信号提取与增强. 电声技术. 2025(02): 52-54+77 . 百度学术
4. 蓝壮青. 基于欠定盲源分离的双路音频信号噪声自适应分离. 现代电子技术. 2024(24): 68-72 . 百度学术
5. 吴源源. 基于机器学习的音频信号处理系统设计. 电声技术. 2024(11): 182-184 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 372
- HTML全文浏览量: 50
- PDF下载量: 22
- 被引次数: 5
