clc
clear
format
%
%?一維大地電磁測(cè)深
%?2013-5-1
%-------------參數(shù)---------------------------------------------------------
nlayer=3;
res（1）?=?10;
res（2）?=?10000;
res（3）?=?100;

thick（1）?=?200;
thick（2）?=?200;
thick（3）?=?200;
%-------------常數(shù)u——真空中的磁導(dǎo)率（SI）-----------------------------------
u=4*pi*10^（-7）;?????
%-------------頻率freq及圓頻率w---------------------------------------------
f?=dlmread（‘freq2.dat‘）;
w=2*pi*f;
nf?=?length（f）;
%-------------計(jì)算視電阻率和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