地球物理电法反演

应用地球物理勘探中的一维电阻率反演程序，基于matlab编写的，

function [rhos] = s1fwd（nlayerresthickr）

%
% c系数，来自程志平1997，p
%
c=[0.003042-0.0011980.012840.02350.086880.23740.61941.18170.4248-3.4507...
    2.7044-1.13240.393-0.14360.05812-0.025210.01125-0.0049780.002072-0.000318];

sum = 0;
for k=1:20
    %
    % 求电阻率转换函数t（k）
    %
    T = res（nlayer）;
    m= exp（k*log（10）/6-2.1719）/r;
    for i = nlayer-1:-1:1
        tt1 = 1 - exp（-2*m*thick（i））;
        tt2 = 1 + exp（-2*m*thick（i））;
        T   = res（i）*（res（i）*tt1+T*tt2）/（res（i）*tt2+T*tt1）;
    end
    t（k） = T;
    
    sum = sum + t（k）*c（k）;
end
rhos=sum;
return

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

     文件        402  2010-05-25 09:18  电阻率\31.dat

     文件        188  2010-05-16 14:33  电阻率\ab.dat

     文件        625  2013-05-09 18:43  电阻率\core.m

     文件        837  2013-05-09 23:42  电阻率\ew.dat

     文件       1360  2013-05-02 19:18  电阻率\pqfile.txt

     文件        177  2013-05-05 21:41  电阻率\re.txt

     文件        536  2013-04-29 15:24  电阻率\tmz.m

     文件      15551  2013-05-10 13:06  电阻率\tmzf.fig

     文件      22209  2013-05-10 13:09  电阻率\tmzf.m

     文件        828  2013-05-09 23:22  电阻率\we.dat

     文件        887  2013-05-05 23:21  电阻率\正演数据.dat

     目录          0  2013-05-10 00:35  电阻率

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

                43600                    12