MT_inversion 大地电磁最小二乘光滑约束反演

大地电磁最小二乘光滑约束反演，包含bostick反演对比

function [rho_aphase]=MT1D_FWD（rhoh）
%解析法，层状介质的大地电磁响应
%rho：各层电阻率
%h：层厚度
%rho_a：视电阻率
%phase：相位
mu=（4e-7）*pi;
T=logspace（-3440）;
k=zeros（size（rho2）size（T2））;
for N=1:size（rho2）
    k（N:）=sqrt（-i*2*pi*mu./（T.*rho（N）））;
end
m=size（rho2）;
z=-（i*mu*2*pi）./（T.*k（m:））;
for n=m-1:-1:1
    A=-（i*mu*2*pi）./（T.*k（n:））;
    B=exp（-2*k（n:）*h（n））;
    z=A.*（A.*（1-B）+z.*（1+B））./（A.*（1+B）+z.*（1-B））;
end
rho_a=（T./（mu*2*pi））.*（abs（z）.^2）;
phase=-atan（imag（z）./real（z））.*180/pi;
% subplot（211）;
% semilogx（Trho_a‘-*‘）
% xlabel（‘T（s）‘）
% ylabel（‘\rho_a（\Omega\cdotm）‘）
% subplot（212）;
% semilogx（Tphase‘-*‘）
% xlabel（‘T（s）‘）
% ylabel（‘\circ）‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         484  2013-10-13 11:31  MT_forward.m
     文件        1452  2013-10-13 11:11  nonlin_inv.m
     文件         709  2013-10-12 19:43  MT1D_FWD.m