15 桂林理工大学硕士学位论文r口lL图4 . 3模拟信号及其各层分解信号重构频域图塞毫召k●黛 桂林理工大学硕士学位论文第四章高频大地电磁法( E H _ 4 )数据信息提取及处理方法4 . 1高频大地电磁法( E H - 4 )数据读取通过阅读E H _ 4 说明书及查阅其它参考文献, 了解了E H _ 4 时间序列Y 文件的存储方式以及存储格式, 利用M a tla b 软件编写正确读取时间序列文件数据的程序, 并能够画出各信号道的时间序列图形( H y 、 E x 、 H x 、 E y )和傅里叶变换频谱图( H y r 、 E x r 、H xr、 E yr、 Iqyi、 E xi、 H xi、 聊)。 由本文第二章中2. 4节的内容对E 聃数据存储格式的研究可知, E H - 4 时间序列数据都保存在Y 文件当中, 该数据文件主要由两个部分组成13tl: 一是头文件信息; 二是时间序列数据。 在实际数据采集过程中, 仪器每次叠加采集都分成3段进行, 每段采集有4 0 9 6 个数据, 于是有4 0 9 6 × 3共1228 8 个数据, 总共4个信号道则一共有1228 8 x 4 个数据, 再加上头文件中采集信息的7 个数据, 这样便构成了一次叠加的完整数据块。编写的数据读取程序首先将头文件信息数据单独读取出来, 然后将时间序列数据按照Y 文件存储格式分别读取出}Iy 、 E x 、 H x 、 E y 四道数据。 头文件中的第一个数是决定了该时间数据块的模式, 当D ecim a tio nex p o n en r- = o 时, 表示该时间数据块为高频采集模块: D ecim a tio nex p o n en t- - 4 时, 表示该时间数据块为低频采集模块。 在该读取程序当中设置一个k l气1228 8 * 4 + 7 )× ( k - 1)和m o d e( k )= A ( k l+ 1)语句( k 为叠加次数), 这样便可以让程序根据头文件信息自动判断采集模式。 如图4 . 1所示, 为Im a g em 读取低频数据及频谱图:}㈧图4. Ilm a g em 软件读取低频数据及其频谱圈17 如图4 . 2所示, 为M a tla b 程序读取低频数据及频谱图:E xH xE y竹日№印H 11日 H i印匾强匝旺15 01∞雹——i————砧UU555●j●j‘j444U3j¨333"2jU2221j1j1j11l¨_ j¨-■ L【坠●L●■It■U5●j●U3拈2lj1¨-图4. 2 M a tla b程序读取低频数据及其凝谱图如图4 . 3所示, 为Im a g cm 读取高频数据及频谱图:P O R n}—‘・- t _}, ●}●●1 07 ●pl'r—P1 5t‘P‘●PS ●:_ l,qP柚|JLl1‘,}1一z ●‘1●}1■i}■’J }、 1。 1。