function ares=dcs(ab2,nlayer,res,h)
c=[0.003042 -0.001198 0.01284 0.0235 0.08688 0.2374 0.6194 1.1817 0.4248 -3.4507 2.7044 -1.1324 0.393 -0.1436 0.05812 -0.0252 0.01125 -0.004978 0.002072 -0.000318];%滤波系数C,共20个
rn=length(ab2); %极距的个数
d=log(10)/6; %取样间距
s=-2.1719; %位移系数
ps=zeros(rn,1); %视电阻率
for i=1:rn
for k=1:20
m=exp(k*d+s)/ab2(i);
T2=p(n); %电阻率转换函数
for j=n:-1:2
T1=p(j-1)*(p(j-1)*tanh(m*h(j-1))+T2)/(p(j-1)+T2*tanh(m*h(j-1)));%电阻率转换函数
T2=T1;
end
ps(i)=ps(i)+T1*c(k);
end
end
loglog(ab2,ps,'-o')
axis([1 1000 10 1000]) %设置坐标值范围,x[1,1000],y[1,1000]
grid on %打开网格
xlabel('AB/2 (m)');
ylabel('ρs (Ωm)');
title('对称四极电测深曲线');
return
评论1