基于自适应傅里叶分解的船舶轴频电磁场信号提取.pdf
《基于自适应傅里叶分解的船舶轴频电磁场信号提取.pdf》由会员分享,可在线阅读,更多相关《基于自适应傅里叶分解的船舶轴频电磁场信号提取.pdf(7页珍藏版)》请在咨信网上搜索。
1、基于自适应傅里叶分解的船舶轴频电磁场信号提取徐震寰1,吴永飞1,裴建新2,3(1.太原理工大学计算机科学与技术学院(大数据学院),山西太原,030024;2.中国海洋大学海洋地球科学学院,山东青岛,266100;3.中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛,266100)摘要:船舶轴频电磁场是水面船舶及水下航行器非常关键的特征,然而由于静态电磁场的存在,大大降低了轴频电磁场的信噪比。为实现低信噪比情况下微弱轴频电磁场信号的有效检测,文中提出一种基于自适应傅里叶分解的信号提取方法,可以将复杂非平稳信号自低频到高频自适应地分解为一系列具有瞬时频率的单分量之和的形式。通过分别处理仿真
2、和实测数据,结果表明,该算法能够克服短时傅里叶变换、小波变换及经验模态分解等方法的缺点,可快速有效地提取到轴频电磁场信号,进而为后续船舶及水下航行器的定位追踪提供参考。关键词:船舶;轴频电磁场;自适应傅里叶分解;信号提取中图分类号:U674.7;TJ630.34文献标识码:A文章编号:2096-3920(2023)04-0593-07DOI:10.11993/j.issn.2096-3920.2023-0049ExtractionofShaft-RateElectromagneticFieldofShipsBasedonAdaptiveFourierDecompositionXUZhenhua
3、n1,WUYongfei1,PEIJianxin2,3(1.CollegeofComputerScienceandTechnology(CollegeofDataScience),TaiyuanUniversityofTechnology,Taiyuan030024,China;2.CollegeofMarineGeosciences,OceanUniversityofChina,Qingdao266100,China;3.KeyLaboratoryofSubmarineGeosciencesandProspectingTechniquesofMinistryofEducation,Ocean
4、UniversityofChina,Qingdao266100,China)Abstract:The shaft-rate electromagnetic field is a crucial feature of ships and undersea vehicles.However,due to thepresenceofstaticelectromagneticfields,thesignal-to-noiseratiooftheshaft-rateelectromagneticfieldsisgreatlyreduced.Toeffectivelydetectweakshaft-rat
5、eelectromagneticfieldsignalsunderlowsignal-to-noiseratioconditions,thispaperproposedasignal extraction method based on adaptive Fourier decomposition(AFD).It could adaptively decompose complex non-stationarysignalsfromlowfrequencytohighfrequencyintoaseriesofsinglecomponentswithinstantaneousfrequency
6、.Throughsimulationandmeasureddataprocessing,theresultsshowthatthisalgorithmcanovercometheshortcomingsofmethodssuchasshort-timeFouriertransform,wavelettransform,andempiricalmodedecomposition.Itcanquicklyandeffectivelyextractshaft-rateelectromagneticfieldsignals,therebyprovidingareferenceforpositionin
7、gandtrackingothershipsandunderseavehicles.Keywords:ship;shift-rateelectromagneticfield;adaptiveFourierdecomposition;signalextraction收稿日期:2023-05-15;修回日期:2023-06-16.基金项目:山西省基础研究计划(自由探索类青年科学研究项目)项目资助(202103021223050).作者简介:徐震寰(1989-),男,博士,讲师,主要研究方向为海洋电磁信号处理.第31卷第4期水下无人系统学报Vol.31 No.42023年8月JOURNALOFUNM
8、ANNEDUNDERSEASYSTEMSAug.2023引用格式 徐震寰,吴永飞,裴建新.基于自适应傅里叶分解的船舶轴频电磁场信号提取 J.水下无人系统学报,2023,31(4):593-599.水下无人系统学报sxwrxtxb.xml-5930引言水面船舶及水下航行器的钢质船壳会因不同金属材料间产生的电化学反应不断被腐蚀。为保护船壳,现代船舶普遍采用牺牲阳极的阴极保护系统和外加电流的阴极保护系统来产生保护电流进行防腐。不同金属材料间产生的腐蚀电流与阴极保护系统的保护电流均会通过海水流向螺旋桨,然后通过轴承返回船壳形成回路。螺旋桨轴承船体回路中的电阻抗会随着螺旋桨轴承规律性的旋转而发生周期性变
9、化,其在船舶周围产生以轴承转动速率为基频的极低频电磁场信号,即轴频电磁场信号1-3。轴频电磁场信号的基频一般位于 17Hz 频段内,它的信号谱线特征明显,可作为有效的非声探测手段对船舶等进行探测、识别与定位4-5。由于舰船产生的轴频电磁场信号幅值低,易受环境噪声及船舶静态电磁场影响,如何有效提取轴频电磁场以提高数据信噪比,是国内外学者长期致力的研究课题。目前在提取轴频电磁场信号的时频特征时,主要采用短时傅里叶变换6-7、小波分析8-12和以经验模态分解(empiricalmodedecom-position,EMD)为核心的希尔伯特-黄变换(Hilbert-Huangtransform,HHT
10、)等13-14电磁信号重构方法。然而,在基于短时傅里叶变换的时频分析方法中,窗函数的类型和大小都会影响信号的时频分析效果;在小波域时频分析方法中,首先面临的是选择何种基函数的问题,其次是所选基函数的适用性问题,导致该方法在信号时频分析的自适应性方面表现欠佳;HHT 虽然能精确刻画信号能量随时间和频率的分布,但由于 EMD 方法是根据人为经验提出的算法,缺乏系统严谨的数学理论框架作为基础,算法本身也仍存在模态混叠、端点效应等制约其工程应用效果的瓶颈问题。近年来,一种新的信号重构方法自适应傅里叶分解(adaptiveFourierdecomposition,AFD)算法被提出,它有着完备的数学证明
11、过程15,可以将一个复杂非平稳信号自适应地分解为若干个具有瞬时频率的傅里叶本征模态函数之和,更有利于信号时频特征的提取16,该算法已在饮用水污染预警17、心脏信号压缩与疾病分类18、无线信道处理19、机械故障诊断20等领域取得了不错的应用效果。因此,文中将 AFD 算法用于电磁场信号的重构,不仅可以拓展该算法的应用领域,而且也可为轴频电磁场信号的提取提供一种新的思路。1AFD 算法1.1基本理论f(t)f(z)z=eitH2(D)D=z C:|z|1CAFD 算法是在函数的复 Hardy 空间中建立的,已知实数信号,可以将其转化为解析信号,其中,接着在空间中进行分解,其中,表示复平面中以原点为
12、中心的开单位圆,表示复平面15。f+H2(D)f+=f1f+对于只含有正频率的信号,令,则可以表示为f+=f1=f1,ea1ea1+R1(z)(1)ea1f1,ea1ea1R1(z)a1式中,表示内积运算,可通过柯西公式求解。上式表示的是:一个信号可以通过一个有理正交系统来分解,第 1 次分解把信号分解为与方向一致的和 与方 向 垂 直 的 残 差。是Hardy 空间单位圆盘字典中的原子,一般取值为 0。R1(z)对于残差可以进一步写成如下形式R1(z)=f2(z)za11a1z(2)其中f2(z)=(f1(z)f1,ea1ea1(z)1a1zza1(3)z=a1f1(z)f1,ea1ea1(
13、z)=0f2H2(D)当时,有。事实上,仍属于 Hardy 空间,因此可以对其进行与式(1)类似的分解过程,即f2(z)=f2,ea2ea2+f3(z)za21a2z(4)联立式(1)、(2)和(4)得f+=f1,ea1ea1+(f2,ea2ea2+f3(z)za21a2z)za11a1z=f1,ea1B1+f2,ea2B2+f3(z)2k=1zak1akz=2k=1fk,eakBk+f3(z)2k=1zak1akz(5)重复上述过程,经过 n 次分解之后,可得2023年8月水下无人系统学报第31卷594JournalofUnmannedUnderseaSystemssxwrxtxb.xml-
14、f+=nk=1fk,eakBk(z)+Rn(z)=nk=1fk,eakBk(z)+fn+1(z)nk=1zak1akz(6)其中残差可通过递归公式求得fn+1(z)=(fk(z)fk,eakea1k(z)1akzzak(7)ak式中,可以通过能量最大选择准则获得ak=argmaxA2fk(a):a D(8)A2fk(ak)=?fk,eak?2=(1|ak|2)|fk(ak)|2(9)limnRn(z)=0由于是可以证明的16,因此在L2规范意义(L2-normsense)下可得如下表达式f+=k=1fk,eakBk(z)(10)Bk(z)式中,为由 k 阶加权 Blaschke 乘积所构成的有
15、理正交基函数Bk(z)=eak(z)k1q=1zaq1aqz(11)akak D k=1,2,eak(z)AFD 算法中起关键作用的就是的选择,是复平面单位圆盘 D 内的复数。为单位模化的 Szeg核,其表达式为eak(z)=1|ak|21akz(12)a1=a2=am=0Bk(eit)k=1如果,则为下述序列1,eit,ei(m1)t,1|am+1|2eimt1am+1eit,1|am+r|2eimt1am+reitm+r1j=m+1eitaj1ajeit,(13)ak=0eak(z)=1而当所有的时,式(10)就是傅里叶级数正频率展开的形式。1.2分解流程f(t)f1,f2,fn基于 AF
16、D 对实数信号进行分解,可得到一系列单一分量,分解过程如图 1 所示。G(t)G1(t)=G(t)a1=0首先对原始信号去均值,然后将其转换至复Hardy 空间得到。令,那么分解得到的第 1 个分量为f1=G1,ea1B1(z)(14)同理,分解得到的第 2 个分量为f2=G2,ea2B2(z)(15)fn重复以上步骤,即可得到第 n 次分解的单分量成分。在每一步分解的过程中,AFD 算法都是从给定的信号中尽可能抓取能量大的部分,然后通过 n 层抓取来实现信号重构。最后通过式(16)恢复出原信号fr=2Ren1fic0(16)fic0=12w20f(eit)dtf(t)式中:为复 Hardy
17、空间中的单分量解析信号;Re 表示实部;是对进行傅里叶分解的第 1 项系数。Eerr 通常情况下,可以采用相对能量误差 Eerr来表征算法在迭代过程中的有效性,若迭代过程中,则认为满足分解次数。Eerr=?f frf?2(17)2AFD 算法性能分析2.1仿真信号分析为说明 AFD 算法重构信号的有效性及特点,首先设计了 2 种简单的振动信号:一是主频为 50Hz。EerrEerr图1AFD 信号分解过程Fig.1SignaldecompositionbyAFD2023年8月徐震寰,等:基于自适应傅里叶分解的船舶轴频电磁场信号提取第4期水下无人系统学报sxwrxtxb.xml-595T=sin
18、(t)+cos(t)=2f的调幅信号(见图 2(a),该信号的幅值随时间变化逐渐增强;二是由 3 个主频分别为 0.3、1和 3Hz 的简谐信号叠加而成的信号(见图 2(b),其表达式 为:,其 中为 角 频率。然后利用提出的 AFD 算法,分别对这 2 种信号进行了分解重构,黑色实线为仿真信号的波形,红色虚线为 AFD 重构的结果。由图 2 可知:AFD算法经过一定的分解次数后,可以很好地恢复原始信号;随着信号复杂程度的增加,若要恢复原始信号,所需的分解层数将大大增加。如恢复单频调幅信号需要分解 46 次,而恢复三频简谐信号则需要分解 99 次,这是因为,对于能量相对较弱的高频成分,AFD
19、需要更多的迭代次数进行逼近。2.2实测数据分析实测轴频电磁场信号的干扰主要来自于能量较大的低频静态电磁场,而 AFD 的分解过程是由低频向高频依次分解的,这意味着仅需将前几次分解的低频成分舍去,即可提取得到轴频电磁场信号。文中采用于 2015 年 6 月在南黄海实测的船舶电磁场数据进行分析,测试区域水深 37m,海底为泥沙底质。测试船只长 68m,宽 15.6m,排水量2650t,螺旋桨转动频率为 220r/min(约为 3.67Hz,对应轴频电磁场的基频),航行速度为 6kn。a1=0ea10 2a1=0a2a2ea2图 3 给出了基于 AFD 算法提取轴频电场(Ey 分量)的过程:图 3(
20、a)表示在分解初始状态,当时(红点),分别在复平面和实平面范围内的取值(蓝色实线),红色实线表示复平面中的单位圆;图 3(b)中的黑色线为实测 Ey 分量的水平电场信号,红色虚线为初始重构结果,由于初始,所以重构曲线是 1 条幅值为 0 的直线;图 3(c)为根据能量最大选择准则(式(8)和式(9)获得的;图 3(d)给出了根据该计算所得的;由图 3(e)可知,当仅分解 1 次时,AFD 就已经较好地逼近了有静态电场引起的信号阶跃;图 3(f)与图 3(c)类似,a3an是用于下一次分解的;由于电场信号的抖动程度相对磁场信号剧烈,文中对其进行了 10 次分解,图 3(g)中的红色星号表示这 1
21、0 次分解所使用在复平面单位圆上的位置;图 3(h)中的红色虚线表示基于 AFD 重构的低频信号,主要包含能量较大的静态电场部分;此时,分解残差中主要包含轴频电场信号(见图 3(i)。同样,文中对实测的磁场数据(Bx 分量)进行了处理,由于磁场信号相对于电场信号更为平稳,仅需分解 3 次即可提取到轴频磁场信号。00.20.40.60.811 ea1 实部 ea2 实部 ea1 绝对值ea2 绝对值0虚部 实部虚部实部1.001 100 2相位(02)相位(02)02464 60.51.01.52.0(a)a1=00200 400 600实测电场 EyAFD 重构结果实测电场 EyAFD 重构结
22、果实测电场 EyAFD 重构结果800 1 0003020100102030Ey/(V/m)(b)分解 1 次重构结果(e)分解 2 次重构结果时间/s01214060110123456(d)a2=0.921 53+0.125 25i0200 400 600时间/s800 1 0003020100102030Ey/(V/m)1.0 0.500.51.0an 实部1.00.500.51.0an 虚部(g)an 在复平面单位圆上的位置0200 400时间/s600 800 1 0003020100102030Ey/(V/m)轴频电场/(V/m)(h)分解 10 次重构结果0200 400 600
23、800 1 0003020100102030(i)R10(z)时间/s0642a3 相位a3 绝对值0 00.20.40.60.80.5目标函数1.01.52.0105(f)根据能量最大选择准则获得的 a30642a2 相位a2 绝对值0 00.20.40.60.80.5目标函数1.02.01.52.53.0105(c)根据能量最大选择准则获得的 a200.05调幅信号AFD 重构结果0.10时间/s0.150.201.51.00.5振幅00.51.01.5(a)分解 46 次重构结果(b)分解 99 次重构结果时间/s108642AFD 重构结果包含 3 个主频成分的叠加信号050振幅5图2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 自适应 傅里叶 分解 船舶 电磁场 信号 提取
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。