异常处理转换的基础是
基于富立叶变换与反变换，将磁异常频谱 与频域处理转换因子相乘后再反变换即可 实现，计算简便，速度快。
四、实验原理
1、频率域向上延拓的计算公式：
T e T 2 1
延拓后的场：
h u 2 v 2
( h >0)
2 i ( ux vy ) T2 T2 e dudv 反付立叶变换
三、实验要求
球体磁异常上延计算、分量转换、导数计算、化磁极可任选 其一，具体要求如下： 1、上延计算：对ΔT、Za、Ha各分量向上延拓5m、10m，画 出对应结果图； 2、分量计算：进行ZaHa，HaZa分量转换，画出对应结 果图； 3、导数计算：利用Za异常计算垂向一、二阶导数，画出对应 结果图； 4、化磁极：利用ΔT数据进行化极处理，画出对应结果图； 5、观察向上延拓、分量转换、导数计算、化磁极前后的异常 特征，分析这些处理转换的作用；
一、实验目的
1、加深对磁性体磁异常在频率域处理转换 原理与作用的认识
2、用Matlab语言编程实现球体磁异常(包括
Za、Ha、Δt)的向上延拓、分量转换、导数
计算、化磁极处理，培养数据处理的实际动
手能力。
二、实验内容
利用两个大小与埋深不同的球体产生的磁异
常（ ΔT 、 Za 、 Ha ），利用频率域处理转换 因子实现向上延拓、分量转换、导数计算、 化磁极处理，异常数据利用实验一中球体的 正演程序计算而来。
五、结果分析
观察向上延拓、分量转换、导数计算、化
磁极前后的异常特征，分析这些处理转换 的作用；
六、实验报告
报告要有封面、要装订，同时要电子版
报告内容要包括实验目的、实验内容、
实验原理、计算程序代码、实验结果、 结果分析或小结
数据准备
球体模型示意图
30 20 10
(nT)
0 -10 -20 -30 5 0 10 15 25 35 45 55 65 75 85 95 105 115 125 135 145 155
(m)
20 30
数据准备
设在地下有两个磁性体，一个浅而小
（半径r=5m，埋深R=25m），一个深而 大（半径r=50m，埋深R=100m ），有效 磁化强度 Ms=0.2A/m, is=90°情况下， 正演其磁异常（ΔT、Za、Ha）。
四、实验原理
2、频率域磁异常分量换算的计算公式：
(u,0) Z (u,0) i Za Ha ： H a a
3、频率域磁异常n阶垂向导数的计算公式
n T 2 2 2 1 ( u v ) T 1 n z
n
四、实验原理
4、频率域化磁极的计算公式：
Z a u 2 v2 i (u l v m) n u 2 v 2 u 2 v2 i (u l0 v m0 ) n0 u 2 v 2 T