大地电磁法MT一维反演Bostick方法Matlab程序

MT一维反演Bostick方法Matlab程序，分为两个Matlaab子程序，一个用于MT一维正演计算，另一个是Bostick反演程序。打开Bostick程序，点击运行，即可自动反演。修改程序中数据，即可完成MT一维反演问题。

%%%%%%  Bostick inversion   %%%%%%%
[bypcph]=compute（[20 200 10][500 2000]）;

T=logspace（-3440）;
u=（4e-7）*pi;
D=sqrt（（pc.*T）./（2*pi*u））;
pd=pc.*（180./（2*ph）-1）;
%loglog（pdD/1000‘r*‘）;

figure（1）

depth=[1 500 2000 100000];rho=[20 200 10 10];
stairs（depth/1000rho）;       %理论模型

hold on

stairs（D/1000pd）;             %Bostick反演模型
ylabel（‘电阻率/\Omega\cdotm‘）;
xlabel（‘深度/km‘）;
grid off


figure（2）
h=diff（D）;   %Bostick反演的层厚度
[b1y1pc1ph1]=compute（pdh）;    %Bostick反演结果的响应值

subplot（211）
semilogx（Tpc‘*‘）;  %理论响应---视电阻率
hold on
semilogx（Tpc1‘r‘）;  %Bostick反演模型的响应---视电阻率

subplot（212）
semilogx（Tph‘*‘）;  %理论响应---相位
hold on
semilogx（Tph1‘r‘）;  %Bostick反演模型的响应---相位

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        902  2012-05-24 18:22  compute.m

     文件        787  2012-05-24 18:22  Bostick.m

----------- ---------  ---------- -----  ----

                 1689                    2