CN103870686A - 一种基于信息融合的超声相控阵方法 - Google Patents
一种基于信息融合的超声相控阵方法 Download PDFInfo
- Publication number
- CN103870686A CN103870686A CN201410082361.5A CN201410082361A CN103870686A CN 103870686 A CN103870686 A CN 103870686A CN 201410082361 A CN201410082361 A CN 201410082361A CN 103870686 A CN103870686 A CN 103870686A
- Authority
- CN
- China
- Prior art keywords
- array
- function
- obtaining
- array element
- fusion
- 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.)
- Granted
Links
- 230000004927 fusion Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000000523 sample Substances 0.000 claims abstract description 27
- 239000011159 matrix material Substances 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 12
- 239000000203 mixture Substances 0.000 claims description 9
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000003491 array Methods 0.000 claims description 2
- 230000001174 ascending effect Effects 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 238000012545 processing Methods 0.000 abstract description 4
- 239000006185 dispersion Substances 0.000 abstract 1
- 230000008030 elimination Effects 0.000 abstract 1
- 238000003379 elimination reaction Methods 0.000 abstract 1
- 238000001514 detection method Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于信息融合的超声相控阵方法,属于超声信号处理技术,所述的算法包括分位图法和贝叶斯估计两部分,其中,分位图法的建立包括以下步骤:取出超声相控阵探头各阵元依次单独采集信号的中位数、上四分位数和下四分位数,得到四分位离散度,求出淘汰点,得到有效数据的融合集;贝叶斯估计的建立包括以下步骤:得出各有效数据的特征函数,求出决策函数,得到风险函数,提取决策函数的贝叶斯风险,得出参数的最有估计值,应用本发明将提高各超声相控阵阵元所采集信息的利用率,能够得到对现实环境更为准确、可靠的描述。
Description
技术领域
本发明属于超声信号处理领域,特别涉及一种基于信息融合的超声相控阵方法。
背景技术
超声相控阵技术于上世纪八十年代被提出,初期主要应用于医学超声成像诊断。近年来国内外超声相控阵技术发展很快,在医学诊断和工业检测等方面的研究很受欢迎。另外,在相控阵的系统设计、系统的模拟以及实际应用等方面也有所进展,相控阵技术开始走向数字化。随着超声相控阵技术越来越受人们的重视,也出现了以下难点:超声相控阵各阵元得到的回波信息具有不同的信号特征,单个信息不能够全面描述被测体;超声波在复合材料中幅值衰减大,所以超声阵列得到的回波幅度很小,后续信息不容易处理;当多个阵元同时发射时,各阵元的发射波、回波之间会发生干涉作用,当各波相位不同时,信号强度会减弱。
根据上述技术难点,中国发明专利申请号200810140346.6采用D一S证据理论融合算法对超声相控阵多个阵元得到的信息进行数据融合,对数字信号处理器DSP数据融合模块数据的特征描述更加准准确。然而,证据理论具有很强的主观性,不区分不确定数据信息和不知道的信息,往往会把不知道的信息认为成是等可能得信息。中国发明专利申请号201210334509.0采用反向传播神经网络(BP)对相控阵各阵元采集的信息进行处理分析,高效完成凸优化处理。由于在网络训练期间,当权值改变时,会出现围绕最佳解的振荡情况,因此BP算法不能在训练期间实现全局误差最小,只可保证实现局部误差最小。中国发明专利申请号97109099.8采用傅里叶变换对超声信号进行频谱分析;中国发明专利申请号200410011403.2通过提升小波包变换对超声信号进行滤波,然后利用希尔伯特变换和正交解调方法对滤波信号进行联合时频分析提取信号在不同频率段的信息特征。但以上这些特征提取方法都是基于希尔伯特变换和傅里叶变换,所以特征提取速度慢,而又由于相控阵超声信号是一种非稳态的时变信号,不满足傅里叶变换的条件,不能准确反映超声信号的特征。我们要严格控制各相邻发射阵元之间具有恒定的相位差,通过硬件实现成本很高,比较困难,我们可以用软件来实现,本发采用信息融合技术来解决上述难点。
信息融合(Information Fusion,或称为数据融合Data Fusion)是电子战、机器人、系统导航、柔性制造、故障诊断以及数字图像处理等领域的一个重要的课题。多传感器信息融合是一种处理多个或多类传感器系统的新方法,又可被称为多元合成、多源关联、混合传感器或多传感器融合等,但是有一种更广泛的说法,即多传感器信息融合,简称为信息融合。它减少了单一传感器的信息盲区,提高了多源信息处理结果的质量,有利于检测人员对检测情况做出更好的判断和决策。目前,尚没有把信息融合技术用于研究超声相控阵的相关文献和实际应用系统。将两者集合,可以更加灵活而有效地控制声束,所以,选择最优信息融合算法运用于实际超声相控阵系统具有广阔的应用和发展前景。因此,本发明采用信息融合技术领域中的其他方法来解决超声相控阵检测中信号处理方面存在的技术难点的可能性。
发明内容
为了解决上述问题,本发明的目的是提供一种基于信息融合的超声相控阵方法,该方法能够提高各超声相控阵阵元所采集信息的利用率,能够得到对现实环境更为准确、可靠的描述。
本发明是通过下述技术方案加以实现的,一种基于信息融合的超声相控阵方法,其特征在于,包括以下步骤:
(一)超声相控阵探头的各个阵元依次单独采集信号,每个阵元采集到的所有信号组成一个数组,得到多个数组;
(二)分位图法的建立
1)各阵元采集得到的信号组成的各数组元素都由小到大排列,求出每个数组的中位数M、上四分位数F1、下四分位数F2;
2)求出信号保留区间(P1,P2),其中dF=F2-F1,P1=F1-a/2*dF,P2=F2+a/2*dF,其中a为经验值常数;各数组元素在区间内的保留,在区间外的剔除,得到比原来数组元素个数少的新数组,数组的大小就是最佳融合数;
(三)贝叶斯估计的建立
1)经过分位图法可以得到新的数组,求出每个新数组的m个元素彼此之间的置信距离测度 (i,j=1,2,...,n),组成置信矩阵D,;其中,,是在第i个样本值取的条件下的概率密度,是第i个样本值的方差;令,经过数学推导得,,其中,是标准正态分布的概率,b表示变量t的上限;
3)由关系矩阵R,得到每个数组参加信息融合的最佳融合数据集;
4)依新数组各元素的pdf(概率密度)曲线作为各特征函数,记为p(x);
7)对风险函数求平均,得到决策函数的贝叶斯风险B,其中,
本发明的优点在于:
1、经过分位图法能在很大程度上排除离异值的干扰,具有很强的抗干扰性。
2、通过分位图法,减少了参与计算的测量数据的数据量,从而降低了运算量,加快了运算速度。
3、若,则认为第i个样本值与第j个样本值相差比较大,相容性差。若,则认为第i个样本值与第j个样本值相差较小,相容性好,第i个样本值是支持第j个样本值的。若,则认为第i个样本值与第j个样本值是相互支持的。如果某个样本值与另一个样本值相容性差,应把后者删除。这样,降低了计算空间的维数,提高运算效率。
应用本发明将提高各超声相控阵阵元所采集信息的利用率,能够得到对现实环境更为准确、可靠的描述。
附图说明
图1:实现本发明的总体流程图;
图2:实现本发明步骤一分位图法的流程图;
图3:实现本发明步骤二贝叶斯估计的流程图;
图4:实施例1超声相控阵1#阵元单独发射的波形图;
图5:实施例1超声相控阵2#阵元单独发射的波形图;
图6:实施例1超声相控阵3#阵元单独发射的波形图;
图7:实施例1超声相控阵4#阵元单独发射的波形图;
图8:实施例1超声相控阵5#阵元单独发射的波形图;
图9:实施例1超声相控阵6#阵元单独发射的波形图;
图10:实施例1超声相控阵7#阵元单独发射的波形图;
图11:实施例1超声相控阵8#阵元单独发射的波形图;
图12:实施例1本发明通过软件实现的波形图。
具体实施方式
为了更清楚的理解本发明,结合附图和实施例详细描述本发明:
实施例1 :
如图 1 至图12所示: 实施例1的对象是不锈钢块。超声相控阵各阵元1#—8#单独发射的波形图如图4—11,可以看出八个图形并不完全相同,即各阵元采集的信号特征有差别,为了综合利用各信号提供的有用信息特征,摒除错误信息特征,采用本发明提供的一种基于信息融合的超声相控阵方法。实施例1的具体算法出图1—3的流程图。
实施例1的步骤如下:
(1)把1#—8#阵元的文本数据导入MATALB,组成8个11088维1列的数组A1、A2、A3、A4、A5、A6、A7、A8。
(2)分位图法的建立
a,把各数组转置并组成8维1088列的矩阵X。
b,分别求出矩阵X各列的中位数、上四分位数、下四分位数,组成1维1088列矩阵M、F1、F2。
c,求出有效区间(P1,P2)。
d,判断各数据是否在有效区间内,在有效区间内的数据保留,不在有效区间内的数据剔除。
(3)贝叶斯估计的建立
把矩阵X每一列得到的有效数据利用贝叶斯估计进行融合,每一列最后都得到一个最有估计,各列最有估计组成1维1088列的数组u,使用MATLAB进行仿真,得到图12。
由于8个阵元同时发射时,各阵元之间的发射波、回波之间存在干涉作用,当相位不同时,信号强度会减弱,因此得到的结果并不理想。而通过本发明一种超声相控阵多阵元的多源信息融合的算法,得到的仿真图形如图12。由图12,能看出通过分位图法后再由贝叶斯估计信息融合算法得到的图形包含的信号特征更加明显,充分利用了各阵元信息的有效性,使有用信号得到加强,没用的信号得到减弱。
实施例2
如表1所示:如下为超声相控阵八个阵元对同一点采集得到的数据特征。
表1 超声相控阵八个阵元对同一点采集得到的数据特征
实施例2步骤如下:
(1)分位图法的建立
(2)贝叶斯估计的建立
应用分位图法后,
融合结果为50.763。
通过以上的仿真结果我们可以得到:
1,利用分位图法之前通过R可以看出最佳融合数为4,融合集为:{ }。而这与通过分位图法得到的传感器组是完全一致的。
2,实施例2是基于贝叶斯估计的数据融合方法,通过引入数据探测技术中的分位图法对其加以改进。首先利用分位图法排除尽可能多的离异值,得到一致阵元组,然后从数学方法的角度,以置信距离测度作为数据融合的融合度,利用置信距离矩阵、关系矩阵得到多传感器的最佳融合数、融合集,以贝叶斯估计理论为基础得到最优融合数据,最后通过仿真证明了本发明的有效性。
根据上述说明,结合本领域技术可实现本发明的方案。
Claims (1)
1.一种基于信息融合的超声相控阵方法,其特征在于,包括以下步骤:
(一)超声相控阵探头的各个阵元依次单独采集信号,每个阵元采集到的所有信号组成一个数组,得到多个数组;
(二)分位图法的建立
1)各阵元采集得到的信号组成的各数组元素都由小到大排列,求出每个数组的中位数M、上四分位数F1、下四分位数F2;
2)求出信号保留区间(P1,P2),其中dF=F2-F1,P1=F1-a/2*dF,P2=F2+a/2*dF,其中a为经验值常数;各数组元素在区间内的保留,在区间外的剔除,得到比原来数组元素个数少的新数组,数组的大小就是最佳融合数;
(三)贝叶斯估计的建立
1)经过分位图法可以得到新的数组,求出每个新数组的m个元素彼此之间的置信距离测度 (i,j=1,2,...,n),组成置信矩阵D,;其中,,是在第i个样本值取的条件下的概率密度,是第i个样本值的方差;令,经过数学推导得,,其中,是标准正态分布的概率,b表示变量t的上限;
3)由关系矩阵R,得到每个数组参加信息融合的最佳融合数据集;
4)依新数组各元素的pdf(概率密度)曲线作为各特征函数,记为p(x);
7)对风险函数求平均,得到决策函数的贝叶斯风险B,其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410082361.5A CN103870686B (zh) | 2014-03-07 | 2014-03-07 | 一种基于信息融合的超声相控阵方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410082361.5A CN103870686B (zh) | 2014-03-07 | 2014-03-07 | 一种基于信息融合的超声相控阵方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103870686A true CN103870686A (zh) | 2014-06-18 |
CN103870686B CN103870686B (zh) | 2015-08-12 |
Family
ID=50909211
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410082361.5A Active CN103870686B (zh) | 2014-03-07 | 2014-03-07 | 一种基于信息融合的超声相控阵方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103870686B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107179102A (zh) * | 2017-06-14 | 2017-09-19 | 国网湖北省电力公司孝感供电公司 | 一种多传感器数据融合的输电线路走廊烟火识别系统及方法 |
CN107274006A (zh) * | 2017-05-23 | 2017-10-20 | 中国地质大学(武汉) | 一种基于贝叶斯推断的多源气象风速融合方法 |
CN108444556A (zh) * | 2018-04-16 | 2018-08-24 | 西南石油大学 | 一种高精度电磁流量计的设计方法 |
CN110346004A (zh) * | 2019-08-16 | 2019-10-18 | 杭州山科智能科技股份有限公司 | 一种双声道超声时差法的流量测量数据融合方法 |
CN115399790A (zh) * | 2022-10-31 | 2022-11-29 | 深圳联友科技有限公司 | 基于贝叶斯数据融合脑电ers_erd的分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101210859A (zh) * | 2007-12-24 | 2008-07-02 | 哈尔滨工业大学 | 高超声速进气道起动/不起动模式的融合分类测定方法 |
-
2014
- 2014-03-07 CN CN201410082361.5A patent/CN103870686B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101210859A (zh) * | 2007-12-24 | 2008-07-02 | 哈尔滨工业大学 | 高超声速进气道起动/不起动模式的融合分类测定方法 |
Non-Patent Citations (1)
Title |
---|
徐兵 等: "《基于贝叶斯估计的超声红外复合测距系统》", 《解放军理工大学学报(自然科学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107274006A (zh) * | 2017-05-23 | 2017-10-20 | 中国地质大学(武汉) | 一种基于贝叶斯推断的多源气象风速融合方法 |
CN107179102A (zh) * | 2017-06-14 | 2017-09-19 | 国网湖北省电力公司孝感供电公司 | 一种多传感器数据融合的输电线路走廊烟火识别系统及方法 |
CN108444556A (zh) * | 2018-04-16 | 2018-08-24 | 西南石油大学 | 一种高精度电磁流量计的设计方法 |
CN110346004A (zh) * | 2019-08-16 | 2019-10-18 | 杭州山科智能科技股份有限公司 | 一种双声道超声时差法的流量测量数据融合方法 |
CN115399790A (zh) * | 2022-10-31 | 2022-11-29 | 深圳联友科技有限公司 | 基于贝叶斯数据融合脑电ers_erd的分析方法 |
CN115399790B (zh) * | 2022-10-31 | 2023-02-03 | 深圳联友科技有限公司 | 基于贝叶斯数据融合脑电ers_erd的分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103870686B (zh) | 2015-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103870686A (zh) | 一种基于信息融合的超声相控阵方法 | |
CN104698837B (zh) | 一种时变线性结构工作模态参数识别方法、装置及应用 | |
CN109649432B (zh) | 基于导波技术的云端平台钢轨完整性监测系统及方法 | |
CN111880157B (zh) | 一种雷达图像中目标的检测方法及系统 | |
CN104331583B (zh) | 一种基于实测海杂波数据的多重分形建模方法 | |
CN107015205B (zh) | 一种分布式mimo雷达检测的虚假目标消除方法 | |
CN101893698B (zh) | 噪声源测试分析方法及其装置 | |
CN106680825A (zh) | 一种声学阵列成像系统与方法 | |
CN104634872A (zh) | 高铁钢轨伤损在线监测方法 | |
CN107202989A (zh) | 一种适用于被动拖曳线列阵声呐的复杂弱目标检测和跟踪方法 | |
CN107942308A (zh) | 复合高斯环境下基于伽玛纹理的自适应Rao检测方法 | |
Zhao et al. | A frequency-domain adaptive matched filter for active sonar detection | |
CN109471091A (zh) | 中重频pd雷达同时解模糊与检测跟踪方法 | |
CN112668527A (zh) | 一种超声导波半监督成像检测方法 | |
CN110702792A (zh) | 一种基于深度学习的合金组织超声检测分类方法 | |
Wan et al. | Parameter estimation of multi frequency hopping signals based on space-time-frequency distribution | |
CN111310719A (zh) | 一种未知辐射源个体识别及检测的方法 | |
Wang et al. | Co-channel interference suppression for LTE passive radar based on spatial feature cognition | |
CN113608193A (zh) | 一种基于UNet的雷达多目标距离和速度估计方法 | |
Kabir et al. | CSI-DeepNet: A Lightweight Deep Convolutional Neural Network Based Hand Gesture Recognition System Using Wi-Fi CSI Signal | |
Zhang et al. | An ultrasonic reverse time migration imaging method based on higher-order singular value decomposition | |
CN112881812A (zh) | 基于机器学习编码的全闪实时定位方法及装置 | |
Huang et al. | LPI waveform recognition using adaptive feature construction and convolutional neural networks | |
Yun et al. | Radar signal recognition algorithms based on neural network and grey relation theory | |
CN116010798A (zh) | 一种基于YOLOv5的复杂典型通信信号快速识别分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
OL01 | Intention to license declared | ||
OL01 | Intention to license declared |