MT一维正演程序

clc
clear
format
%
% 一维大地电磁测深
% 2013-5-1
%-------------参数---------------------------------------------------------
nlayer=3;
res（1） = 10;
res（2） = 10000;
res（3） = 100;

thick（1） = 200;
thick（2） = 200;
thick（3） = 200;
%-------------常数u——真空中的磁导率（SI）-----------------------------------
u=4*pi*10^（-7）;     
%-------------频率freq及圆频率w---------------------------------------------
f =dlmread（‘freq2.dat‘）;
w=2*pi*f;
nf = length（f）;
%-------------计算视电阻率和Phi---------------------------------------------
for n=1:nf
    Z = w（n）*u/（sqrt（1i*w（n）*u/res（nlayer）））; 
    if nlayer==1
        Z = w（n）*u/（sqrt（1i*w（n）*u/res（nlayer）））; 
    else
        for ii = nlayer-1:-1:1
        k（ii）=sqrt（1i*u*w（n）*（1/res（ii）））;       %%%%
        zz（ii）=w（n）*u/k（ii）;                     %%%% 
        a=1-exp（2*1i*k（ii）*thick（ii））;           %%%%     
        b=1+exp（2*1i*k（ii）*thick（ii））;           %%%%  
        Z=zz（ii）*（zz（ii）*a+Z*b）/（zz（ii）*b+Z*a）;
        end 
    end 
        rhos（n1）=（abs（Z）*abs（Z））/（w（n）*u）;
        %rhos（n1）=0.2/f（n）*（abs（Z）*abs（Z））;
        Phi（n1）=abs（angle（Z））/pi*180;
end

delta=503*s