MATLAB与C＋＋混合编程实现航磁异常提取的小波分析方法研究

Research on Wavelet Analysis Method of Aeromagnetic Anomaly Extraction based on MATLAB and C++ Mixed Programming

阐述了小波多尺度分析原理和Mallat算法。针对以往MATLAB文件执行速度较慢的问题,提出了一种新的通过MEX文件来实现在MATLAB环境中调用C++源码程序的混合编程技术,并对球体磁异常叠加理论模型进行正演计算和对比研究,最后利用该方法对某测区的航磁异常进行分解,并展示两个阶次典型的航磁异常细节特征图。结果表明:用小波分解得到的航磁异常细节特征图,同样起到了功率谱分解的作用,不但可以从异常细节图上清晰地圈闭出深大断裂带的赋存位置,而且在速度上较以往的单一计算方法也有了明显的提高,对今后类似地区数据处理具有较重要的参考价值。

We explain a multi-scale analysis of wavelet theory and Mallat Algorithm.In view of the previous slow pace of implementation for MATLAB files,We propose a new mixed programming technique to call C＋＋ source process using MEX file in MATLAB environment.We make a calculation and comparison of magnetic anomalies superposition of the sphere theoretical model,and using the method to decompose the aeromagnetic anomalies of a certain measured area and its surrounding areas,we get two orders of the typical characteristics of aeromagnetic anomaly map.The results show that detail of aeromagnetic anomaly map obtained through wavelet decomposition also plays a role in the power spectrum analysis.The deep faults can be located on the anomalous map,hence,this map provide a reference for the similar regional data processing in future.