２０２１年４月上
第５０卷㊀第７期
施㊀工㊀技㊀术
ＣＯＮＳＴＲＵＣＴＩＯＮＴＥＣＨＮＯＬＯＧＹ
７７
㊀ＤＯＩ:１０ ７６７２/ｓｇｊｓ２０２１０７００７７
∗
王晓曙１ꎬ赵湖潮１ꎬ张延杰１ꎬ许利东１ꎬ丁红岗１ꎬ任孟德２
(１.云南建投基础工程有限责任公司ꎬ云南㊀昆明㊀６５００００ꎻ２.西安建筑科技大学土木工程学院ꎬ陕西㊀西安㊀７１００５５)[摘要]采用探地雷达进行强风化地层隧道超前地质预报时ꎬ受现场探测环境㊁强风化地层的不均匀性以及解译技术制约ꎬ造成富水区雷达反射波的解译困难ꎮ将探地雷达瞬时属性和频 谱属性应用到富水强风化地层研究中ꎬ介绍了瞬时属性和频谱属性的原理ꎬ并阐述了探地雷达属性在富水强风化地层预报中的可行性ꎮ结合工程实例分析研究表明：富水地层的瞬时振幅增强ꎬ瞬时相位轴稀疏ꎬ瞬时频率降低ꎮ根据１００ＭＨｚ雷达天线实测结果ꎬ良好地层的主频值在７０ＭＨｚ左右ꎬ富水地层的主频值在３０ＭＨｚ左右ꎮ
国际组织的作用[关键词]隧道ꎻ地质预报ꎻ探地雷达ꎻ强风化地层ꎻ三瞬属性ꎻ频谱属性[中图分类号]ＴＵ７５３
[文献标识码]Ａ
[文章编号]１００２￣８４９８(２０２１)０７￣００７７￣０４
ＡｐｐｌｉｃａｔｉｏｎＡｎａｌｙｓｉｓｏｆｔｈｅＧｒｏｕｎｄＰｅｎｅｔｒａｔｉｎｇＲａｄａｒＡｔｔｒｉｂｕｔｅｓｉｎ
Ｗａｔｅｒ￣ｒｉｃｈａｎｄＳｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄＳｔｒａｔｕｍ
ＷＡＮＧＸｉａｏｓｈｕ１ꎬＺＨＡＯＨｕｃｈａｏ１ꎬＺＨＡＮＧＹａｎｊｉｅ１ꎬＸＵＬｉｄｏｎｇ１ꎬ
ＤＩＮＧＨｏｎｇｇａｎｇ１ꎬＲＥＮＭｅｎｇｄｅ２
(１.ＹｕｎｎａｎＣｏｎｓｔｒｕｃｔｉｏｎａｎｄＩｎｖｅｓｔｍｅｎｔＦｏｕｎｄａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇＣｏ.ꎬＬｔｄ.ꎬＫｕｎｍｉｎｇꎬＹｕｎｎａｎ㊀６５００００ꎬＣｈｉｎａꎻ
２.ＳｃｈｏｏｌｏｆＣｉｖｉｌＥｎｇｉｎｅｅｒｉｎｇꎬＸｉ ａｎＵｎｉｖｅｒｓｉｔｙｏｆＡｒｃｈｉｔｅｃｔｕｒｅａｎｄＴｅｃｈｎｏｌｏｇｙꎬＸｉ ａｎꎬＳｈａａｎｘｉ㊀７１００５５ꎬＣｈｉｎａ)Ａｂｓｔｒａｃｔ:Ｉｔｉｓｄｉｆｆｉｃｕｌｔｔｏｉｎｔｅｒｐｒｅｔｔｈｅｒａｄａｒｒｅｆｌｅｃｔｉｏｎｗａｖｅｓｉｎｗａｔｅｒ￣ｒｉｃｈｓｔｒａｔｕｍｗｈｅｎｕｓｉｎｇｇｒｏｕｎｄｐｅｎｅｔｒａｔｉｎｇｒａｄａｒ(ＧＰＲ)ｔｏｆｏｒｅｃａｓｔｔｈｅｓｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄｓｔｒａｔｕｍｔｕｎｎｅｌｄｕｅｔｏｎｕｍｅｒｏｕｓｆａｃｔｏｒｓꎬｓｕｃｈａｓｔｈｅｆｉｅｌｄｄｅｔｅｃｔｉｏｎｅｎｖｉｒｏｎｍｅｎｔꎬｔｈｅｉｎｈｏｍｏｇｅｎｅｉｔｙｏｆｓｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄｓｔｒａｔｕｍａｎｄｔｈｅｉｎｔｅｒｐｒｅｔａｔｉｏｎｔｅｃｈｎｏｌｏｇｙａｎｄｓｏｏｎ.ＴｈｅｉｎｓｔａｎｔａｎｅｏｕｓａｔｔｒｉｂｕｔｅｓａｎｄｓｐｅｃｔｒｕｍａｔｔｒｉｂｕｔｅｓｏｆＧＰＲａｒｅａｐｐｌｉｅｄｔｏｔｈｅｓｔｕｄｙｏｆｔｈｅｗａｔｅｒ￣ｒｉｃｈａｎｄｓｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄｓｔｒａｔｕｍ.ＴｈｅｐｒｉｎｃｉｐｌｅｓｏｆＧＰＲａｔｔｒｉｂｕｔｅｓａｒｅｉｎｔｒｏｄｕｃｅｄꎬａｎｄｔｈｅｆｅａｓｉｂｉｌｉｔｙｏｆｔｈｅｗａｔｅｒ￣ｒｉｃｈａｎｄｓｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄｓｔｒａｔｕｍｐｒｅｄｉｃｔｉｏｎｗｉｔｈＧＰＲｂａｓｅｄｏｎｔｈ
ｅａｔｔｒｉｂｕｔｅｓｉｓｅｘｐｏｕｎｄｅｄ.Ｆｉｎａｌｌｙꎬｔｈｅｅｘａｍｐｌｅｓａｒｅａｎａｌｙｚｅｄｉｎｔｈｉｓｐａｐｅｒａｎｄｉｔｓｈｏｗｓｔｈａｔｔｈｅｗａｔｅｒ￣ｒｉｃｈａｎｄｓｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄｓｔｒａｔｕｍｅｎｈａｎｃｅｔｈｅｉｎｓｔａｎｔａｎｅｏｕｓａｍｐｌｉｔｕｄｅｏｆｒｅｆｌｅｃｔｅｄｓｉｇｎａｌꎬｔｈｅｉｎｓｔａｎｔａｎｅｏｕｓ
ｐｈａｓｅａｘｉｓｉｓｓｐａｒｓｅꎬａｎｄｒｅｄｕｃｅｔｈｅｉｎｓｔａｎｔａｎｅｏｕｓｆｒｅｑｕｅｎｃｙ.Ｂａｓｅｄｏｎｔｈｅｒｅｓｕｌｔｓｍｅａｓｕｒｅｄｂｙａ１００ＭＨｚｒａｄａｒａｎｔｅｎｎａꎬｔｈｅｄｏｍｉｎａｎｔｆｒｅｑｕｅｎｃｙｖａｌｕｅｏｆｔｈｅｇｏｏｄｓｔｒａｔｕｍｉｓａｐｐｒｏｘｉｍａｔｅｌｙ７０ＭＨｚꎬａｎｄｔｈｅｄｏｍｉｎａｎｔｆｒｅｑｕｅｎｃｙｖａｌｕｅｏｆｔｈｅｗａｔｅｒ￣ｒｉｃｈｓｔｒａｔｕｍｉｓａｂｏｕｔ３０ＭＨｚ.Ｋｅｙｗｏｒｄｓ:ｔｕｎｎｅｌｓꎻｇｅｏｌｏｇｉｃａｌｐｒｅｄｉｃｔｉｏｎꎻｇｒｏｕｎｄｐｅｎｅｔｒａｔｉｎｇｒａｄａｒꎻｓｔｒｏｎｇ￣ｗｅａｔｈｅｒｅｄｓｔｒａｔｕｍꎻｔｈｒｅｅｉｎｓｔａｎｔａｎｅｏｕｓａｔｔｒｉｂｕｔｅｓꎻｓｐｅｃｔｒｕｍａｔｔｒｉｂｕｔｅ∗国家十三五重点研发计划(２０１７ＹＦＣ０８０４６０１)ꎻ国家自然科学基金地区项目(５１７６４０２０)
[作者简介]王晓曙ꎬ高级工程师ꎬ董事长ꎬ总经理
[通信作者]张延杰ꎬ工程师ꎬ博士ꎬＥ￣ｍａｉｌ:ｙａｎｊｉｅ＿ｔｍ＠１６３.ｃｏｍ[收稿日期]２０２０￣１２￣２１
０㊀引言
在隧道施工中ꎬ常常会遇到富水地质ꎬ强风化砂岩遇水软化ꎬ黏聚力大幅度降低ꎬ导致隧道大变形㊁坍塌事故时有发生[１￣２]ꎮ提高强风化地层隧道掌子面前方地质含水情况预报的准确性ꎬ是确保当
前隧道安全施工亟待解决的问题ꎮ
与其他地球物理探测方法相比ꎬ探地雷达具有高分辨率㊁结果直观㊁扫描迅速等优势ꎬ对含水体非常敏感ꎬ非常适合含水地质体的探测[３￣５]ꎮ但受隧道强风化地质㊁岩体破碎和探测环境的干扰ꎬ探测人员仅依靠探地雷达剖面图中的强反射轴来对地质体含水情况进行推断和解释ꎬ很容易造成误判和漏判ꎮ为从雷达信号中挖掘更多的信息来刻画和解释探测目标ꎬ将探地雷达的三瞬属性和频谱属性应用到地质解译中ꎮ赵文轲等[６]对探地雷达属性
７８㊀施工技术第５０卷
分析技术进行较全面的总结㊁评价和分析ꎬ认为探地雷达属性技术拥有良好的应用前景ꎮ刘东坤等[７]对混合型充填溶洞㊁空溶洞㊁富水破碎带３种不良地质体的雷达反射波进行频谱分析ꎬ得出了３种不
同的主频值ꎬ为推测不良地质体提供依据ꎮ刘宗辉等[８]利用雷达属性提取技术ꎬ对岩溶不良地质体多角度定性分析属性特征ꎬ在定性基础上实现属性特征定量表示ꎮ李尧等[９]基于钻孔雷达对典型不良地质体进行正演模拟并进行复信号分析ꎬ利用复信号技术可提高钻孔雷达解译的精度和准确性ꎮ刘斌等[１０]将复信号分析技术引入地质雷达预报岩溶裂隙水的研究与应用中ꎬ多角度㊁多标准地对岩溶裂隙水进行预报ꎮ刘成禹等[１１]基于探地雷达单道信号处理方法确定岩溶的位置㊁大小和形状ꎬ对岩溶不良地质体多角度定性分析属性特征ꎬ在定性基础上实现属性特征定量表示ꎮ上述研究主要对岩溶地层的雷达属性分析ꎬ都取得了良好的效果ꎬ但目前尚不存在针对富水强风化地层探地雷达属性分析工作ꎮ
本文从风口山隧道中选取富水区段和良好区段作为对比分析案例ꎬ首先对雷达原始数据进行常规处理ꎬ然后利用三瞬属性和频谱属性分别在时间域和频率域上对探地雷达反射波瞬时振幅㊁瞬时相位㊁瞬时频率和能量㊁频率进行分析检测ꎬ在前者定性分析基础上ꎬ根据主频值进行量化分析ꎬ从而达到对探地雷达数据深入分析ꎬ提高对富水地质体的预报准确度ꎮ
１㊀属性分析原理
１ １㊀三瞬属性
㊀㊀三瞬属性指瞬时振幅㊁瞬时相位与瞬时频率３种属性ꎮ瞬时振幅反映了反射信号的强度随着传播深度㊁被探测介质的电磁性质等因素的改变而变化ꎬ与该时刻地质雷达信号总能量的平方根成正比ꎮ当
地质体有明显变化时ꎬ特别遇到富水段ꎬ地质的介电常数差就会很大ꎬ瞬时振幅明显增强ꎮ瞬时相位反映了地质雷达剖面图上同相轴的连续性ꎬ无论反射波的能量强弱ꎬ其相位都能显示出来ꎮ当电磁波在各向同性均匀介质中传播时ꎬ其相位是连续的ꎻ当电磁波在不均匀介质中传播时ꎬ其相位在剖面图中不连续ꎮ瞬时频率是相位的时间变化率ꎬ反映了组成地层的物性变化ꎮ当电磁波通过不同介质时ꎬ电磁波反射信号的频率将会发生变化ꎮ由于水会吸收掉电磁波中的高频部分ꎬ当地质体富水时瞬时频率会明显偏低ꎮ
１ ２㊀频谱属性
㊀㊀雷达反射波频谱分析ꎬ就是把信号从时域转换到频域中去分析ꎮ设雷达信号时域函数为ｆ(ｔ)ꎬ经
过傅里叶变换后的频率域函数为Ｆ(ｗ)ꎬ频率域函
数Ｆ(ｗ)为时域函数ｆ(ｔ)的频谱密度函数ꎬ简称为
频谱[１２]ꎮ进行频谱分析后ꎬ可以得到雷达波频率与能量分布的关系曲线ꎮ能量与振幅相关ꎬ设能量谱
为Ｓ(ｗ)ꎬ根据能量定理可以得到能量谱和振幅谱
之间的关系ꎬ即:
维纳滤波
Ｓ(ｗ)＝｜Ｆ(ｗ)｜２(１)㊀㊀因地质的完整性㊁破碎性情况可以间接通过雷达反射波频率的高低呈现出来ꎬ所以可以通过对雷达反射波频谱分析得到该地层的频率值ꎬ将其值用作判断地质体的定量指标ꎮ
在电磁场的作用下ꎬ水会导致围岩介质的电导
率显著增大ꎬ电场电荷在高电导率的介质中更容易
做功而产生大量的电荷能量消耗ꎬ雷达更容易接受
低频反射信号ꎬ而高频振荡电磁波则更易被吸收ꎮ
因此ꎬ通过建立频率与能量的关系曲线进行频谱分
析ꎬ寻求地质体含水情况的主频值差异ꎬ进而定量
分析地质体含水情况ꎮ
２㊀工程概况与探测方法㊀
２ １㊀工程概况
㊀㊀风口山隧道是一座双向四车道双连拱隧道ꎬ全长５１５ｍꎮ隧道区属低中山构造剥蚀地貌ꎬ地层以强风化砂岩为主ꎬ节理裂隙发育ꎬ岩体破碎ꎬ多呈碎石状碎裂结构ꎻ地下水类型主要为基岩裂隙水ꎬ主要受大气降雨补给ꎬ水量受区内降雨及季节性影响较大ꎮ强风化砂岩遇水软化ꎬ黏聚力大幅度降低ꎬ隧道大变形㊁坍塌风险增加ꎮ因此ꎬ对其富水区进行超前地质预报识别是一项重要内容ꎮ
２ ２㊀预报方式
㊀㊀本次探测的雷达设备选用瑞典ＭＡＬＡ/Ｘ３Ｍ型探地雷达ꎬ雷达脉冲重复频率１０~２００ｋＨｚ(标准１００ｋＨｚ)ꎮ采用发射频率为１００ＭＨｚ的屏蔽天线ꎬ样点数设为６６０ꎬ时窗为５００~７００ｎｓꎬ采样频率１０２４ＭＨｚꎬ叠加次数１２８次ꎬ道间距为０ １ｍꎬ以点测方式采集数据ꎮ测线布置如图１所示ꎮ为减少外界环境杂波的干扰ꎬ探测前应将掌子面附近的开挖台车㊁机械设备撤离ꎬ并暂停掌子面附近下台阶钢拱立架㊁锚杆钻孔㊁焊接等施工作业ꎮ
采用ＭＡＴＧＰＲ软件对雷达反射波数据进行分
析预处理ꎬ预处理步骤主要包括背景去除㊁带通滤
波㊁反褶积等ꎬ通过预处理来去除反射波中的杂波㊁
增强有效波ꎬ提高信噪比ꎮ
３㊀探地雷达波形图像分析
㊀㊀采用ＭＡＴＧＰＲ软件对现场探测雷达反射波数据预处理后ꎬ得出富水地层与良好地层的探地雷达
2012全国百强县２０２１Ｎｏ.７王晓曙等：探地雷达属性在富水强风化地层中的应用分析７９
㊀
图１㊀探地雷达测线布置
反射波波形堆积图ꎬ如图２所示ꎮ由富水地层与良好地层的现场实际开挖掌子面情况可知ꎬ富水地层的隧道上台阶左侧壁有股状水流出ꎬ导致掌子面局部滑塌ꎻ良好地层的隧道掌子面未见有水ꎬ呈干燥状态ꎮ由于该隧道地层为强风化砂砾岩ꎬ节理裂隙发育ꎬ岩体较破碎ꎬ地层的不均匀性导致富水区段与良好区段的波形图均存在局部波形较杂乱ꎬ同相轴错断ꎬ振幅变化大ꎬ在地质解释时往往都解释为岩体破碎ꎬ却很难推测前方围岩含水情况ꎮ可见ꎬ针对强风化地层ꎬ仅仅依靠反射波形很难准确解译前方地质情况ꎬ需进一步挖掘雷达反射波的深层信息来提高超前地质预报的准确性
ꎮ
图２㊀探地雷达反射波波形堆积图
４㊀探地雷达属性分析４ １㊀三瞬属性４ １ １㊀瞬时振幅
㊀㊀瞬时振幅只呈现雷达反射波的振幅信息ꎬ可排除相位信息的干扰ꎬ使雷达反射波信号的强弱更容易被区分ꎮ富水地层和良好地层的瞬时振幅分别如图３所示ꎮ富水地层在２５０~３５０ｎｓ位置呈现深红ꎬ说明电磁波遇到富水地层时产生强烈反射信号ꎬ即能量大ꎬ振幅强ꎻ相比良好地层瞬时振幅ꎬ却没有出现红区域ꎬ即反射信号一般ꎮ因此ꎬ通过瞬时振幅图的颜判断反射信号的强弱ꎬ定性推测地层的含水情况
ꎮ
北京亚运会开幕式图３㊀探地雷达瞬时振幅
４ １ ２㊀瞬时相位
㊀㊀瞬时相位只呈现雷达反射波的相位信息ꎬ可排除振幅信息的干扰ꎬ微弱的反射波也能清晰呈现出来ꎮ富水地层和良好地层的瞬时相位分别如图４所示ꎮ由图４可知ꎬ由于岩体破碎ꎬ富水地层和良好地层的同相轴均错断不连续ꎻ富水地层相位轴较稀疏ꎬ良好地层相位轴较密集ꎬ说明富水地层频率较低ꎬ良好地层频率较高
ꎮ
图４㊀探地雷达瞬时相位
４ １ ３㊀瞬时频率
㊀㊀富水地层和良好地层的瞬时频率分别如图５所示ꎮ由图５可知ꎬ富水地层反射波频率较低ꎬ良好地层反射波频率相对较高ꎬ这是因为水会吸收掉反射波中的高频成分ꎬ雷达更容易接收到低频成分ꎮ因此ꎬ通过分析瞬时频率图的颜可以反映反射波频率高低ꎬ间接定性判断地层含水情况ꎮ
４ ２㊀频谱属性
硫铁矿制硫酸㊀㊀对富水区段的雷达反射波其中２道单道波做傅里叶变换进行频谱分析ꎬ在频谱分析图中建立振幅和频率的分布关系曲线ꎬ如图６所示ꎮ两幅图均呈现高频成分的振幅较低ꎬ低频成分的振幅高达１７０左右ꎬ这说明电磁波在含水量较大的围岩介质中传播时ꎬ雷达反射波中大量高频部分被吸收ꎬ且该部
８０㊀施工技术第５０
卷
图５㊀探地雷达瞬时频率
分高频反射波能量耗费较高ꎬ因此所对应的能量较小ꎮ通过对频谱分析后得到的主频值是围岩富水的一种综合反映ꎬ从图６中可以看出该富水区段的主频值大约为３０ＭＨｚꎬ因此可以把该主频值作为推断围岩富水的量化指标判据
ꎮ
图６㊀富水地层频谱分析
作为对比ꎬ对良好区段的雷达反射波中２道单道波做傅里叶变换进行频谱分析ꎬ在频谱分析图中建立振幅和频率的分布关系曲线ꎬ如图７所示ꎮ两幅图均呈现出两个峰值ꎬ峰值振幅５０左右ꎬ主频值
７０ＭＨｚ左右ꎻ从而７０ＭＨｚ左右的主频值作为一个量化指标可对推测围岩良好
ꎮ
图７㊀良好地层频谱分析
５㊀结语㊀㊀本文利用探地雷达三瞬属性和频谱属性ꎬ对富
㊀㊀㊀
水区段和良好区段的雷达反射波信号进行对比分析ꎬ得出以下主要结论ꎮ
１)强风化地层岩体破碎ꎬ地质的不均匀性导致
富水区段与良好区段的波形图均存在局部波形较杂乱㊁同相轴错断㊁振幅变化大ꎬ根据波形图很难推测前方围岩含水情况ꎮ
２)通过对富水地层与良好地层的三瞬属性对比分析ꎬ富水地层的瞬时振幅明显增强ꎬ瞬时相位轴较稀疏ꎬ瞬时频率明显降低ꎮ
３)通过对富水地层与良好地层的频谱对比分
析ꎬ得到电磁波在富水地层传播过程中高频的吸收程度比在良好地层高ꎬ良好地层的主频值在７０ＭＨｚ左右ꎬ富水地层的主频值在３０ＭＨｚ左右ꎮ
鞍山信托
４)将探地雷达属性应用到富水强风化地层中
能提高预报的准确度ꎬ但探地雷达易受探测环境干扰ꎮ为确保施工安全ꎬ建议结合地质勘查资料㊁掌子面地质素描及外界降雨量综合分析ꎬ必要时可同时采用其他探测方法进行综合超前地质预报ꎮ
参考文献:
[１]㊀邹翀ꎬ雷胜友ꎬ岳喜军ꎬ等.富水全㊁强风化砂岩强度特性试验
及本构关系探讨[Ｊ].中国工程科学ꎬ２０１１ꎬ１３(１):７４￣８０.[２]㊀张卫中ꎬ陈从新ꎬ余明远.风化砂岩的力学特性及本构关系研究[Ｊ].岩土力学ꎬ２００９ꎬ３０(Ｓ１):３３￣３６.[３]㊀刘志生.地铁施工区段地下空洞探测及病害处理研究[Ｊ].施工技术ꎬ２０１９ꎬ４８(１３):１０４￣１０７.
[４]㊀刘澜波ꎬ钱荣毅.探地雷达：浅表地球物理科学技术中的重要
工具[Ｊ].地球物理学报ꎬ２０１５ꎬ５８(８):２６０６￣２６１７.
[５]㊀李术才ꎬ刘斌ꎬ孙怀凤ꎬ等.隧道施工超前地质预报研究现状
及发展趋势[Ｊ].岩石力学与工程学报ꎬ２０１４ꎬ３３(６):
１０９０￣１１１３.
[６]㊀赵文轲ꎬ陈国顺ꎬ田钢ꎬ等.探地雷达属性技术进展[Ｊ].地球物
理学进展ꎬ２０１２ꎬ２７(３):１２６２￣１２６７.
[７]㊀刘东坤ꎬ巨能攀ꎬ霍宇翔.地质雷达在不同介质填充下的频谱
差异分析[Ｊ].现代隧道技术ꎬ２０１３ꎬ５０(５):２３￣２８.
[８]㊀刘宗辉ꎬ刘毛毛ꎬ周东ꎬ等.基于探地雷达属性分析的典型岩
溶不良地质识别方法[Ｊ].岩土力学ꎬ２０１９(８):１￣９.
[９]㊀李尧ꎬ李术才ꎬ刘斌ꎬ等.钻孔雷达探测地下不良地质体的正
演模拟及其复信号分析[Ｊ].岩土力学ꎬ２０１７ꎬ３８(１):３００￣３０８.
[１０]㊀刘斌ꎬ李术才ꎬ李树忱ꎬ等.复信号分析技术在地质雷达预报
岩溶裂隙水中的应用研究[Ｊ].岩土力学ꎬ２００９ꎬ３０(７):
２１９１￣２１９６.
[１１]㊀刘成禹ꎬ余世为.基于探地雷达单道信号处理的岩溶分析方
法[Ｊ].岩土力学ꎬ２０１６ꎬ３７(１２):３６１８￣３６２６.[１２]㊀李录明.信号与系统分析[Ｍ].成都：电子科技大学出版
社ꎬ２００２.
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议