CN111916089B - 基于声信号特征分析的冰雹检测方法和装置 - Google Patents
基于声信号特征分析的冰雹检测方法和装置 Download PDFInfo
- Publication number
- CN111916089B CN111916089B CN202010729978.7A CN202010729978A CN111916089B CN 111916089 B CN111916089 B CN 111916089B CN 202010729978 A CN202010729978 A CN 202010729978A CN 111916089 B CN111916089 B CN 111916089B
- Authority
- CN
- China
- Prior art keywords
- signal
- identified
- hail
- acoustic
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 41
- 238000004458 analytical method Methods 0.000 title claims abstract description 40
- 238000012549 training Methods 0.000 claims abstract description 41
- 239000013598 vector Substances 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000006870 function Effects 0.000 claims description 15
- 230000005236 sound signal Effects 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 12
- 230000001681 protective effect Effects 0.000 claims description 5
- 230000000087 stabilizing effect Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 11
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 230000009467 reduction Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L17/00—Speaker identification or verification techniques
- G10L17/26—Recognition of special voice characteristics, e.g. for use in lie detectors; Recognition of animal voices
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L15/00—Speech recognition
- G10L15/06—Creation of reference templates; Training of speech recognition systems, e.g. adaptation to the characteristics of the speaker's voice
- G10L15/063—Training
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L15/00—Speech recognition
- G10L15/08—Speech classification or search
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L15/00—Speech recognition
- G10L15/08—Speech classification or search
- G10L15/14—Speech classification or search using statistical models, e.g. Hidden Markov Models [HMMs]
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Multimedia (AREA)
- Data Mining & Analysis (AREA)
- Acoustics & Sound (AREA)
- Human Computer Interaction (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Signal Processing (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于声信号特征分析的冰雹检测方法,确定多个训练样本,对各个训练样本进行特征分析,得到训练样本的特征向量,对各组特征向量进行聚类运算,得到第一聚类中心、第二聚类中心,以及隶属度函数,根据待识别信号的特征参数确定待识别向量,根据所述隶属度函数和待识别向量确定待识别信号的聚类中心,得到待识别聚类中心,在所述第一聚类中心和第二聚类中心中查找距所述待识别聚类中心近的聚类中心,将查找得到的聚类中心所表征的信号类型确定为待识别信号的信号类型,以确定待识别信号为冰雹信号或者雨声信号,实现对待识别信号的高效检测,相应的检测过程具有更高的准确性。
Description
技术领域
本发明涉及电气象检测技术领域,尤其涉及一种基于声信号特征分析的冰雹检测方法和装置。
背景技术
冰雹是一种常见的自然现象,一旦发生会产生极大的破坏力。长期以来,由于冰雹具有突发性强、破坏性大以及发生范围小的特点,所以冰雹的特征和冰雹降落时的天气状况一直都是气象研究者们的重要研究内容。在青藏高原、云南高原以及一些山地更是冰雹频发地区,对当地农作物及人民的生产生活均有不同程度的影响,故而对冰雹监测技术的发展尤为重视。
目前对于冰雹的监测主要有两种方式,一种是人工观测,人工观测通常是通过冰雹落下时砸在冰雹板上的印迹得出冰雹的粒径和密度,但冰雹砸在冰雹板上的印迹往往会重叠在一起,所以辨别冰雹的重复落点次数以及冰雹的粒径对于人工观测人员来说难度很大,而且人工观测的方法不仅耗时耗力,并且得出的冰雹粒径、历时、密度也会存在很大的误差,导致无法准确得出冰雹的量级。另一种则是通过天气雷达监测冰雹云,根据雷达回波的统计特征识别冰雹等强对流天气。虽然利用天气雷达探测冰雹云的方式在一定程度上提升了检测效率,然而仍然不能够十分准确地判断冰雹的量级等信息,可见传统方案通常存在检测结果片面,准确度低的问题。
发明内容
针对以上问题,本发明提出一种基于声信号特征分析的冰雹检测方法和装置,以获得冰雹的相关信息与天气雷达的回波图像进行统计应对,从而来实现全天候智能监测、准确率更高,更是节省人力、物力、财力等方面的消耗。
为实现本发明的目的,提供一种基于声信号特征分析的冰雹检测方法,包括如下步骤:
S10,将标注类别的多个冰雹信号和多个雨声信号确定为各个训练样本,对各个训练样本进行特征分析,得到各个训练样本分别对应的一组特征向量;所述特征向量包括相应训练样本的多个特征参数;
S20,对各组特征向量进行聚类运算,得到表征冰雹信号的第一聚类中心,表征雨声信号的第二聚类中心,以及隶属度函数;
S30,根据待识别信号的特征参数确定待识别向量,根据所述隶属度函数和待识别向量确定待识别信号的聚类中心,得到待识别聚类中心;
S40,在所述第一聚类中心和第二聚类中心中查找距所述待识别聚类中心近的聚类中心,将查找得到的聚类中心所表征的信号类型确定为待识别信号的信号类型;所述信号类型包括冰雹信号和雨声信号。
在一个实施例中,上述基于声信号特征分析的冰雹检测方法,还包括:
若所述待识别信号的信号类型为冰雹信号,提取所述待识别信号的能量大小,根据所述待识别信号的能量大小确定所述待识别信号的直径。
具体地,根据所述待识别信号的能量大小确定所述待识别信号的直径包括:
若所述待识别信号的能量大于1且小于或者等于300,则所述待识别信号的直径小于或者等于1.5cm,若所述待识别信号的能量大于300且小于540,则所述待识别信号的直径大于1.5cm且小于3cm,若所述待识别信号的能量大于或等于540,则所述待识别信号的直径大于或等于3cm。
在一个实施例中,所述待识别聚类中心包括两个聚类中心。
在一个实施例中,上述特征参数包括均方根、波形因子、峭度因子、最大带宽能量比和/或平均振幅值。
一种基于声信号特征分析的冰雹检测装置,包括拾音板、声波传感器、数采装置和计算机;
所述计算机通过所述数采装置连接声波传感器,所述声波传感器放置于拾音板的正下方,垂直地面放置且不与上方的拾音板相接触,以采集冰雹和雨水落在拾音板上产生的声信号;当声波传感器接收到声信号时,数采装置开始工作,数采装置将声波传感器采集到的声信号传输给计算机;所述计算机以当前接收的声信号为待识别信号执行所述的基于声信号特征分析的冰雹检测方法。
在一个实施例中,上述拾音板为方形,拾音板的四角由四个长度为10cm的弹簧支撑悬空。
在一个实施例中,上述数采装置的最外侧是壳体,是整个数采装置的保护壳,所述壳体内部分别布置采集卡、存储器、控制器和蓄电池,控制器通过导线分别与存储器和采集卡相连,控制器与蓄电池通过稳压电路进行连接,为整个数采装置供电。
上述基于声信号特征分析的冰雹检测方法,将标注类别的多个冰雹信号和多个雨声信号确定为各个训练样本,对各个训练样本进行特征分析,得到各个训练样本分别对应的一组特征向量,对各组特征向量进行聚类运算,得到表征冰雹信号的第一聚类中心,表征雨声信号的第二聚类中心,以及隶属度函数,根据待识别信号的特征参数确定待识别向量,根据所述隶属度函数和待识别向量确定待识别信号的聚类中心,得到待识别聚类中心,在所述第一聚类中心和第二聚类中心中查找距所述待识别聚类中心近的聚类中心,将查找得到的聚类中心所表征的信号类型确定为待识别信号的信号类型,以确定待识别信号为冰雹信号或者雨声信号,实现对待识别信号的高效检测,相应的检测过程具有更高的准确性。
附图说明
图1是一个实施例的基于声信号特征分析的冰雹检测方法流程图;
图2是一个实施例的基于声信号特征分析的冰雹检测装置结构示意图;
图3是一个实施例的数采装置结构示意图;
图4是一个实施例的基于声信号特征分析的冰雹检测装置工作过程示意图;
图5为一个实施例中降雹和雨声信号的时域波形图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
在本文中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是,本文所描述的实施例可以与其它实施例相结合。
参考图1所示,图1为一个实施例的基于声信号特征分析的冰雹检测方法流程图,包括如下步骤:
S10,将标注类别的多个冰雹信号和多个雨声信号确定为各个训练样本,对各个训练样本进行特征分析,得到各个训练样本分别对应的一组特征向量;所述特征向量包括相应训练样本的多个特征参数。
在一个实施例中,所述特征参数包括均方根、波形因子、峭度因子、最大带宽能量比和/或平均振幅值。
S20,对各组特征向量进行聚类运算,得到表征冰雹信号的第一聚类中心,表征雨声信号的第二聚类中心,以及隶属度函数。
上述聚类运算采用的聚类算法可以为基于马氏距离的模糊聚类算法(fuzzy c-means algorithm based on Mahalanobis distance,M-FCM)。
具体地,采用基于马氏距离的模糊聚类算法对冰雹进行分类识别,可以预先采集200个降雹和雨声信号,采用统计分析的方法对其进行时域和频域的特征分析,提取出均方根、波形因子、峭度因子、最大带宽能量比和平均振幅值等特征,具体计算见公式(1)至公式(5),这里平均振幅值为声信号在频域上频率范围[0Hz,15000Hz]的振幅平均值,再利用分类算法对降雹和雨声信号的特征向量进行训练学习。
在一个示例中,特征参数的确定公式包括:
其中,X={x1,x2,…,xN}为降雹或雨声信号样本,N为数据采样点的个数。xn(m)为对声信号进行分帧加窗处理得到的第n帧信号,M为帧长。dk为声信号最大能量衰减到-3dB时的带宽,表示求X的平均值后再求绝对值,()max表示求最大值。
S30,根据待识别信号的特征参数确定待识别向量,根据所述隶属度函数和待识别向量确定待识别信号的聚类中心,得到待识别聚类中心。
具体地,待识别聚类中心的确定过程可以包括:
式中,t的取值为1或者2,当t的取值为2时,待识别聚类中心包括两个聚类中心,C1t表示第t个聚类中心,xj表示第j个待识别样本,n为待识别样本的数目,b为加权参数,U(xj)表示隶属度函数,n表示待识别信号的信号个数。
在一个实施例中,所述待识别聚类中心包括两个聚类中心,如表征冰雹信号的聚类中心,和表征雨声信号的聚类中心,此时,待识别信号包括两类信号,相应的天气为雨加冰雹。
S40,在所述第一聚类中心和第二聚类中心中查找距所述待识别聚类中心近的聚类中心,将查找得到的聚类中心所表征的信号类型确定为待识别信号的信号类型;所述信号类型包括冰雹信号和雨声信号。
上述基于声信号特征分析的冰雹检测方法,将标注类别的多个冰雹信号和多个雨声信号确定为各个训练样本,对各个训练样本进行特征分析,得到各个训练样本分别对应的一组特征向量,对各组特征向量进行聚类运算,得到表征冰雹信号的第一聚类中心,表征雨声信号的第二聚类中心,以及隶属度函数,根据待识别信号的特征参数确定待识别向量,根据所述隶属度函数和待识别向量确定待识别信号的聚类中心,得到待识别聚类中心,在所述第一聚类中心和第二聚类中心中查找距所述待识别聚类中心近的聚类中心,将查找得到的聚类中心所表征的信号类型确定为待识别信号的信号类型,以确定待识别信号为冰雹信号或者雨声信号,实现对待识别信号的高效检测,相应的检测过程具有更高的准确性。
在一个实施例中,上述基于声信号特征分析的冰雹检测方法,还包括:
若所述待识别信号的信号类型为冰雹信号,提取所述待识别信号的能量大小,根据所述待识别信号的能量大小确定所述待识别信号的直径。
具体地,根据所述待识别信号的能量大小确定所述待识别信号的直径包括:
若所述待识别信号的能量大于1且小于或者等于300,则所述待识别信号的直径小于或者等于1.5cm,若所述待识别信号的能量大于300且小于540,则所述待识别信号的直径大于1.5cm且小于3cm,若所述待识别信号的能量大于或等于540,则所述待识别信号的直径大于或等于3cm。
具体地,所采集到的降雹信号分别包含直径d≥3cm、1.5cm<d<3cm、d≤1.5cm三种大中小冰雹声信号,分别对大中小冰雹做频谱分析,得到其信号主要集中在[0Hz,15000Hz]频率范围内,再计算[0Hz,15000Hz]频率上信号的平均能量,并进行统计分析得到冰雹直径d≤1.5cm时,其能量在0到300范围内;冰雹直径为1.5cm<d<3cm时,其能量在300到540范围内;冰雹直径d≥3cm时,其能量大于540,为后续确定冰雹量级做准备。
本实施例能够识别冰雹信号的直径,可以提高相应检测方案的完整性。
在一个实施例中,以待识别信号包括两类信号为例对上述基于声信号特征分析的冰雹检测方法进行说明,此时上述基于声信号特征分析的冰雹检测方法主要分为两步,一是训练,这里所述训练样本的类别已知,先提取选定训练样本的特征参数组成特征向量,再利用M-FCM算法对提取的特征值进行迭代运算,计算得到降雹和雨声信号的聚类中心C01(第一聚类中心)、C02(第二聚类中心)和隶属度函数U(xi),如公式(6)所示,式中C0t为训练样本聚类中心,xi表示第i个训练样本,m为类别数目,b为加权参数。
二是识别,依照训练过程的步骤,对待识别样本提取特征,并根据训练过程中得到的隶属度函数U,由公式(7)计算得到待识别样本的聚类中心C11、C12,公式(7)中,C1t为待识别样本的聚类中心,xj表示第j个待识别样本,n为待识别样本的数目。同时分别计算C11、C12到C01、C02的距离,分别将待识别样本的聚类中心与训练样本的聚类中心进行比较以及C01、C02到C11、C12的距离进行比较,待识别样本的聚类中心坐标最接近训练样本的聚类中心以及到训练样本聚类中心距离最短的一类,即该未知样本所属类别与训练样本所属类别一致,从而完成对降雹信号的识别。依据直径不同的冰雹产生的能量大小确定冰雹的量级,后再提取出降雹信号的能量大小判断冰雹的量级大小。经实验表明,降雹识别准确率平均可达93.333%,当降雹能量高于300时,识别准确率高达100%,但降雹能量低于300时,识别准确率为88.889%,相对较低。
在一个实施例中,提供一种基于声信号特征分析的冰雹检测装置,包括拾音板、声波传感器、数采装置和计算机;
所述计算机通过所述数采装置连接声波传感器,所述声波传感器放置于拾音板的正下方,垂直地面放置且不与上方的拾音板相接触,以采集冰雹和雨水落在拾音板上产生的声信号;当声波传感器接收到声信号时,数采装置开始工作,数采装置将声波传感器采集到的声信号传输给计算机;所述计算机以当前接收的声信号为待识别信号执行上述任一实施例所述的基于声信号特征分析的冰雹检测方法。
上述基于声信号特征分析的冰雹检测装置,首次提出从“声学”的角度对冰雹进行监测,能够判断降雨的过程中是否伴随冰雹降落,若有则可检测到冰雹的量级,同时将数据实时传输给后台监测终端。本装置整体数据传输为无线传输,综合成本更低,性能稳定性更高并且传输速度更快。与现有冰雹监测技术相比,本实施例所提出的设计结构更加简单,测量数据误差更小,准确性更高。本实施例对人工除雹工作以及冰雹灾情的损失估算提供数据依据,具有一定的价值。
在一个实施例中,所述拾音板为方形,拾音板的四角由四个长度为10cm的弹簧支撑悬空。
在一个实施例中,所述数采装置的最外侧是壳体,是整个数采装置的保护壳,所述壳体内部分别布置采集卡、存储器、控制器和蓄电池,控制器通过导线分别与存储器和采集卡相连,控制器与蓄电池通过稳压电路进行连接,为整个数采装置供电。
具体地,数采装置的最外侧是壳体,是整个数采装置的保护壳,而在壳体内部分别布置有采集卡、存储器、控制器和蓄电池四部分。控制器通过导线分别与存储器和采集卡相连,控制器与蓄电池通过稳压电路进行连接,为整个数采装置供电。采集卡通过数据连接线外连一个声波传感器,并将声波传感器放置于拾音板的正下方,垂直地面放置且不与上方的拾音板相接触,以便于采集冰雹和雨水落在拾音板上产生的声信号。整个数采装置采用蓄电池供电,当声波传感器接收到声信号时,数采装置开始工作。采集卡将声波传感器采集到的声信号直接送给控制器进行处理,同时放入存储器保存,控制器对实时数据处理完成后一并将采集的数据和处理结果通过存储器传输给后台计算机。
在一个实施例中,上述基于声信号特征分析的冰雹检测装置可以参考图2所示,其中计算机11执行上述基于声信号特征分析的冰雹检测方法,能够判断降雨的过程中是否有冰雹降落,以及若存在冰雹降落时,判断出冰雹的量级大小。本装置主要包含拾音板1、声波传感器3、数采装置、计算机11。拾音板1为方形,拾音板1的四角由四个长度约为10cm的弹簧2支撑悬空,悬空的目的是为了采集到的降雹和雨声信号特征更加的明显,有利于对降雹声信号的识别。数采装置的最外侧是壳体9,是整个数采装置的保护壳,而在壳体内部分别布置有采集卡5、存储器8、控制器7和蓄电池6。采集卡5通过数据连接线4外连一个声波传感器3,并将声波传感器3放置于拾音板1的正下方,垂直地面放置且不与上方的拾音板1相接触,以便于采集冰雹和雨水落在拾音板1上产生的声信号。
如图3所示,整个数采装置采用蓄电池6供电,采集卡将声波传感器3采集到的声信号直接送给控制器7进行处理,同时放入存储器8保存,控制器8对实时数据处理完成后一并将采集的数据和处理结果传输给后台计算机11。
如图4所示,本实施例的分类识别算法主要分为两步,一是训练,先提取选定训练样本的特征参数组成特征向量,再利用M-FCM算法对提取的特征值进行迭代运算,计算得到降雹和雨声信号的聚类中心C01、C02和隶属度函数U;二是识别,依照训练过程的步骤,对待识别样本提取特征,并根据训练过程中得到的隶属度函数U,计算得到待识别样本的聚类中心C11、C12,并分别计算C11、C12到C01、C02的距离,将计算得到的聚类中心与第二步中通过训练样本计算得到的聚类中心以及C01、C02到C11、C12的距离进行比较,聚类中心最接近以及距离最短的一类就是未知样本所属的类别,从而完成对降雹信号的识别结果,后续再提取出降雹信号的能量大小判断冰雹的量级大小。经实验表明,降雹识别准确率平均可达93.333%,当降雹能量高于300时,识别准确率高达100%,但降雹能量低于300时,识别准确率为88.889%,相对较低。
如图5所示为本实施例所设计的声信号采集装置在实验过程中所真实采集的降雹和雨声信号的时域波形图,图5中,(a)表示集的降雹声信号,(b)表示集的降雨声信号。
表1特征值统计表
如表1所示,预先采集了200个降雹和雨声信号,采用统计分析的方法对其进行时域和频域的特征分析,提取出均方根、波形因子、峭度因子、最大带宽能量比、平均振幅值五个特征组成特征向量,这里分别选取了降雹和雨声信号各5个样本提取了上述五个特征。所得到的试验结果,分别选取了96组数据作为训练样本和55组数据作为测试样本,已知训练样本1是雨声信号和训练样本2是降雹信号,利用M-FCM算法得到测试样本1是雨声信号、测试样本2是降雹信号。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
需要说明的是,本申请实施例所涉及的术语“第一\第二\第三”仅仅是区别类似的对象,不代表针对对象的特定排序,可以理解地,“第一\第二\第三”在允许的情况下可以互换特定的顺序或先后次序。应该理解“第一\第二\第三”区分的对象在适当情况下可以互换,以使这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
本申请实施例的术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或模块的过程、方法、装置、产品或设备没有限定于已列出的步骤或模块,而是可选地还包括没有列出的步骤或模块,或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或模块。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (8)
1.一种基于声信号特征分析的冰雹检测方法,其特征在于,包括如下步骤:
S10,将标注类别的多个冰雹信号和多个雨声信号确定为各个训练样本,对各个训练样本进行特征分析,得到各个训练样本分别对应的一组特征向量;所述特征向量包括相应训练样本的多个特征参数;
S20,对各组特征向量进行聚类运算,得到表征冰雹信号的第一聚类中心,表征雨声信号的第二聚类中心,以及隶属度函数;
S30,根据待识别信号的特征参数确定待识别向量,根据所述隶属度函数和待识别向量确定待识别信号的聚类中心,得到待识别聚类中心;
S40,在所述第一聚类中心和第二聚类中心中查找距所述待识别聚类中心近的聚类中心,将查找得到的聚类中心所表征的信号类型确定为待识别信号的信号类型;所述信号类型包括冰雹信号和雨声信号。
2.根据权利要求1所述的基于声信号特征分析的冰雹检测方法,其特征在于,还包括:
若所述待识别信号的信号类型为冰雹信号,提取所述待识别信号的能量大小,根据所述待识别信号的能量大小确定所述待识别信号的直径。
3.根据权利要求2所述的基于声信号特征分析的冰雹检测方法,其特征在于,根据所述待识别信号的能量大小确定所述待识别信号的直径包括:
若所述待识别信号的能量大于1且小于或者等于300,则所述待识别信号的直径小于或者等于1.5cm,若所述待识别信号的能量大于300且小于540,则所述待识别信号的直径大于1.5cm且小于3cm,若所述待识别信号的能量大于或等于540,则所述待识别信号的直径大于或等于3cm。
4.根据权利要求1至3任一项所述的基于声信号特征分析的冰雹检测方法,其特征在于,所述待识别聚类中心包括两个聚类中心。
5.根据权利要求1至3任一项所述的基于声信号特征分析的冰雹检测方法,其特征在于,所述特征参数包括均方根、波形因子、峭度因子、最大带宽能量比和/或平均振幅值。
6.一种基于声信号特征分析的冰雹检测装置,其特征在于,包括拾音板、声波传感器、数采装置和计算机;
所述计算机通过所述数采装置连接声波传感器,所述声波传感器放置于拾音板的正下方,垂直地面放置且不与上方的拾音板相接触,以采集冰雹和雨水落在拾音板上产生的声信号;当声波传感器接收到声信号时,数采装置开始工作,数采装置将声波传感器采集到的声信号传输给计算机;所述计算机以当前接收的声信号为待识别信号执行权利要求1至5任一项所述的基于声信号特征分析的冰雹检测方法。
7.根据权利要求6所述的基于声信号特征分析的冰雹检测装置,其特征在于,所述拾音板为方形,拾音板的四角由四个长度为10cm的弹簧支撑悬空。
8.根据权利要求6所述的基于声信号特征分析的冰雹检测装置,其特征在于,所述数采装置的最外侧是壳体,是整个数采装置的保护壳,所述壳体内部分别布置采集卡、存储器、控制器和蓄电池,控制器通过导线分别与存储器和采集卡相连,控制器与蓄电池通过稳压电路进行连接,为整个数采装置供电。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010729978.7A CN111916089B (zh) | 2020-07-27 | 2020-07-27 | 基于声信号特征分析的冰雹检测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010729978.7A CN111916089B (zh) | 2020-07-27 | 2020-07-27 | 基于声信号特征分析的冰雹检测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111916089A CN111916089A (zh) | 2020-11-10 |
CN111916089B true CN111916089B (zh) | 2022-11-04 |
Family
ID=73281437
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010729978.7A Active CN111916089B (zh) | 2020-07-27 | 2020-07-27 | 基于声信号特征分析的冰雹检测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111916089B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113095442B (zh) * | 2021-06-04 | 2021-09-10 | 成都信息工程大学 | 基于半监督学习在多维度雷达数据下的冰雹识别方法 |
DE102022201680A1 (de) | 2022-02-17 | 2023-08-17 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung eingetragener Verein | System mit einem akustischen sensor und verfahren zur echtzeitfähigen erfassung meteorologischer daten |
CN117074258A (zh) * | 2023-10-17 | 2023-11-17 | 北京中科技达科技有限公司 | 基于声学原理的冰雹测量装置及方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5528224A (en) * | 1994-05-03 | 1996-06-18 | Scientific Technology, Inc. | Optical and acoustic weather identification system |
KR100755677B1 (ko) * | 2005-11-02 | 2007-09-05 | 삼성전자주식회사 | 주제 영역 검출을 이용한 대화체 음성 인식 장치 및 방법 |
CN104616020A (zh) * | 2014-09-22 | 2015-05-13 | 国家电网公司 | 一种基于聚类分析的输电线路覆冰检测方法和系统 |
US10475447B2 (en) * | 2016-01-25 | 2019-11-12 | Ford Global Technologies, Llc | Acoustic and domain based speech recognition for vehicles |
CN110161506B (zh) * | 2019-07-01 | 2023-03-31 | 江苏省气象科学研究所 | 一种基于多源气象观测资料的分类型冰雹落区识别方法 |
-
2020
- 2020-07-27 CN CN202010729978.7A patent/CN111916089B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111916089A (zh) | 2020-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111916089B (zh) | 基于声信号特征分析的冰雹检测方法和装置 | |
CN102959354B (zh) | 用于利用LiDAR数据来分析树冠层的方法和装置 | |
CN102445640B (zh) | 基于向量机和人工鱼群优化的gis设备智能识别方法 | |
CN108564005A (zh) | 一种基于卷积神经网络的人体跌倒辨识方法 | |
US20150185270A1 (en) | Method for recognizing transformer partial discharge pattern based on singular value decomposition algorithm | |
CN111048114A (zh) | 一种设备异常声音检测的设备及方法 | |
CN105841961A (zh) | 一种基于Morlet小波变换和卷积神经网络的轴承故障诊断方法 | |
CN108648764B (zh) | 基于雨水敲击声识别的雨量测量系统及其测量方法 | |
CN112735473B (zh) | 基于声音识别无人机的方法及系统 | |
CN112331231B (zh) | 基于音频技术的肉鸡采食量检测系统 | |
CN109164450B (zh) | 一种基于多普勒雷达数据的下击暴流预测方法 | |
CN108694953A (zh) | 一种基于Mel子带参数化特征的鸟鸣自动识别方法 | |
CN112727710B (zh) | 一种基于音频信号的风场落雷密度统计方法及系统 | |
CN112068120A (zh) | 基于二维傅里叶变换的微多普勒时频平面单兵小分队识别方法 | |
CN117572837B (zh) | 一种智慧电厂ai主动运维方法及系统 | |
CN104749249B (zh) | 一种基于超声波技术的检测种子净度的方法 | |
CN108064745A (zh) | 动物嚎叫监测系统以及基于机器学习的嚎叫状态识别方法 | |
CN106979998B (zh) | 苹果新鲜度的仿生嗅觉快速无损检测装置及检测方法 | |
CN113361383A (zh) | 一种基于遗传算法特征参数优化的瓷绝缘子损伤神经网络识别方法 | |
CN112380903A (zh) | 一种基于WiFi-CSI信号增强的人体活动识别方法 | |
CN114783462A (zh) | 一种基于cs-music矿井提升机故障源定位分析方法 | |
CN108924753A (zh) | 用于场景的实时检测的方法和装置 | |
CN114118168A (zh) | 多站联合的电磁脉冲事件识别方法、系统和设备 | |
CN114121025A (zh) | 一种面向变电站设备的声纹故障智能检测方法及装置 | |
CN112347917A (zh) | 一种燃气轮机故障诊断方法、系统、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |