频率域位场处理与转换1225-.docx

第PAGE \* MERGEFORMAT16页 重磁实验报告 实验名称：频率域位场处理与转换 学生姓名：刘国强 专业名称：勘查技术与工程 学生学号：201226020126 指导老师：巨鹏 王万银  目录  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc438832376 一．实验内容与要求  PAGEREF _Toc438832376 \h 2  HYPERLINK \l _Toc438832377 （1） 实验内容：  PAGEREF _Toc438832377 \h 2  HYPERLINK \l _Toc438832378 （2） 要求：  PAGEREF _Toc438832378 \h 2  HYPERLINK \l _Toc438832379 二．基本原理  PAGEREF _Toc438832379 \h 2  HYPERLINK \l _Toc438832380 （1）空间域位场延拓：  PAGEREF _Toc438832380 \h 2  HYPERLINK \l _Toc438832381 （2）频率域位场延拓：  PAGEREF _Toc438832381 \h 3  HYPERLINK \l _Toc438832382 三.输入输出数据格式设计  PAGEREF _Toc438832382 \h 3  HYPERLINK \l _Toc438832383 1.输入数据：  PAGEREF _Toc438832383 \h 3  HYPERLINK \l _Toc438832384 2.重要变量名  PAGEREF _Toc438832384 \h 3  HYPERLINK \l _Toc438832385 四．总体设计  PAGEREF _Toc438832385 \h 4  HYPERLINK \l _Toc438832386 五．测试结果  PAGEREF _Toc438832386 \h 5  HYPERLINK \l _Toc438832387 六．结论与建议  PAGEREF _Toc438832387 \h 6  一．实验内容与要求 实验内容： 对一个网格化数据进行向上延拓，得到延拓后的结果，并计算延拓后的均方误差。 要求： 进行软件设计、代码编写、软件测试、编写实验报告。 二．基本原理 分别介绍空间域和频率域位场延拓的计算公式。 （1）空间域位场延拓： 延拓分为向上延拓和向下延拓。 向上延拓是只将实测平面上的数据换到平面之上的平面上的计算。对于二度体（令z向下为正）： U(x,z)=-zπ-∞+∞Uξ,0ξ-x2+z2dξ （z0） 其中U（ξ，0）为剖面上各点的实测值。Z为延拓的高度，即转换后的平面和观测平面的距离，由于z坐标向下为正，因此z0。空间域的延拓实际是积分的计算。 向下延拓是由实测磁场向磁源方向延拓。 （2）频率域位场延拓： 设场源位于z=H平面以下（H0），则磁场为在z=H平面以上对x、y、z连续的函数，若z=0观测平面的磁场T（x，y，0）为已知，则由外部的狄利克莱问题： T（x，y，z）=-z2π-∞+∞-∞+∞T（ξ，η，0）[(x-ξ)2+(y-η)2+z2]3/2dξdη 对其进行变量x，y进行傅里叶变换成 STu，v，z=-∞∞-∞∞T（x，y，z）e-2πi(ux+vy)dxdy 因此： STu，v，z= STu，v，0e2π(u2+v2)1/2z 在频率域中，延拓变成了对观测数据的傅里叶变换乘以延拓因子。 三.输入输出数据格式设计 1.输入数据： （1）建立文件：wave.par （2）观测面高度之差：height=3.3m （3）低高度网格化数据：A20_mag.grd （4）特征值：eigval=1.701411*1038 （6）输出文件名：expound_a20_2D.grd 2.重要变量名 （1）输入文件：input_filename (2) 输出文件：out_filename (3) 点数：mpoint (4)线数：nline (5)异常值：filed 四．总体设计 结束 调用子程序，输出文件 调用子程序，计算延拓因子 调用子程序，对扩边数据进行快速傅里叶变换 调用子程序，进行扩边 调用程序，读入 A20_mag.grd场值 调用程序，计算扩边点位 调用子程序，读入 A20_mag.grd数据 开始 调用子程序，读入wave.par文