資源簡介
應(yīng)用地球物理勘探中的一維電阻率反演程序,基于matlab編寫的,
代碼片段和文件信息
function?[rhos]?=?s1fwd(nlayerresthickr)
%
%?c系數(shù),來自程志平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
????%
????%?求電阻率轉(zhuǎn)換函數(shù)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??電阻率\正演數(shù)據(jù).dat
?????目錄??????????0??2013-05-10?00:35??電阻率
-----------?---------??----------?-----??----
????????????????43600????????????????????12
- 上一篇:圖書借閱管理信息系統(tǒng)
- 下一篇:qt圖像處理
評論
共有 條評論