function [Di]=banjin(Rm)
Rm=1.222225;% 泥漿電阻率值
A=importdata('xing10-013.txt');%根據(jù)需要修改輸入數(shù)據(jù)文件名
hh=1;%深度數(shù)據(jù)在數(shù)據(jù)文件xing10-013.txt（即矩陣A）中的列數(shù)
h=A(:,hh);
ss=3;%淺視電阻率在數(shù)據(jù)文件（即矩陣A）中的列數(shù)
Rlls=A(:,ss);
dd=2;%深視電阻率在數(shù)據(jù)文件xing10-013.txt（即矩陣A）中的列數(shù)
Rlld=A(:,dd);
oo=4;%沖洗帶電阻率在數(shù)據(jù)文件xing10-013.txt（即矩陣A）中的列數(shù)
Rxo=A(:,oo);
DD=7;%井徑在數(shù)據(jù)文件xing10-013.txt（即矩陣A）中的列數(shù)
Dh=A(:,DD);
Rxo=Rxo';
Rlls=Rlls';
Rlld=Rlld';
Sd=load('jiaozhengD.txt');%讀入深視電阻率校正曲線函數(shù)的系數(shù)值
Ss=load('jiaozhengS.txt');%讀入深視電阻率校正曲線函數(shù)的系數(shù)值
S=load('banjinXX.txt');%讀入侵入半徑校正曲線函數(shù)的系數(shù)值
x=log10(Rlld/Rm);
x0=log10(Rlls/Rm);
for t=1:9434%在數(shù)據(jù)文件xing10-013.xls中數(shù)據(jù)總行數(shù)
k=[6 8 10 12 14];
m=k-Dh(t);
if find(m==0)~=zeros(1,0)%當(dāng)井徑為[6 8 10 12 14]中的某個值時
    n=find(m==0);
    p1=Sd(:,n);
    p2=Ss(:,n);
    Rlldcc(t)=polyval(p1,x(t))*Rlld(t);%計算得到深視電阻率校正后井徑8in下的深視電阻率
    Rllscc(t)=polyval(p2,x(t))*Rlls(t);%計算得到淺視電阻率校正后井徑8in下的淺視電阻率 
else%當(dāng)井徑不為[6 8 10 12 14]中的任何值
  
    r1=find(m<0);
    [n11 n0]=size(r1);
   switch n0
       case{0}%當(dāng)井徑小于6in時
           n1=n0+1;
           p11=Sd(:,n1);
    p12=Sd(:,n1+1);
    y11=polyval(p11,x(t));
    y12=polyval(p12,x(t));
    y1=[y11 y12];
    x1=[k(n1) k(n1+1)];
    Rlldc(t)=interp1(x1,y1,Dh(t),'linear','extrap');
    Rlldcc(t)=Rlldc(t)*Rlld(t); %利用插值方法計算得到深視電阻率校正后井徑8in下的深視電阻率
   
   p21=Ss(:,n1);
    p22=Ss(:,n1+1);
    y21=polyval(p21,x0(t));
    y22=polyval(p22,x0(t));
    y2=[y21 y22];
    Rllsc(t)=interp1(x1,y2,Dh(t),'linear','extrap');
    Rllscc(t)=Rllsc(t)*Rlls(t);%利用插值方法計算得到淺視電阻率校正后井徑8in下的淺視電阻率
       case{5}%當(dāng)井徑大于14in時
           n1=n0-1;
          p11=Sd(:,n1);
    p12=Sd(:,n1+1);
    y11=polyval(p11,x(t));
    y12=polyval(p12,x(t));
    y1=[y11 y12];
    x1=[k(n1) k(n1+1)];
    Rlldc(t)=interp1(x1,y1,Dh(t),'linear','extrap');
    Rlldcc(t)=Rlldc(t)*Rlld(t);%利用插值方法計算得到深視電阻率校正后井徑8in下的深視電阻率
    
   p21=Ss(:,n1);
    p22=Ss(:,n1+1);
    y21=polyval(p21,x0(t));
    y22=polyval(p22,x0(t));
    y2=[y21 y22];
    Rllsc(t)=interp1(x1,y2,Dh(t),'linear','extrap');
    Rllscc(t)=Rllsc(t)*Rlls(t);%利用插值方法計算得到淺視電阻率校正后井徑8in下的淺視電阻  
       otherwise%當(dāng)井徑在[6 8 10 12 14]之間但不等于6或8或10或12或14
           n1=n0;
            p11=Sd(:,n1);
    p12=Sd(:,n1+1);
    y11=polyval(p11,x(t));
    y12=polyval(p12,x(t));
    y1=[y11 y12];
    x1=[k(n1) k(n1+1)];
    Rlldc(t)=interp1(x1,y1,Dh(t));
    Rlldcc(t)=Rlldc(t)*Rlld(t);%利用插值方法計算得到深視電阻率校正后井徑8in下的深視電阻率
    
   p21=Ss(:,n1);
    p22=Ss(:,n1+1);
    y21=polyval(p21,x0(t));
    y22=polyval(p22,x0(t));
    y2=[y21 y22];
    Rllsc(t)=interp1(x1,y2,Dh(t));
    Rllscc(t)=Rllsc(t)*Rlls(t);%利用插值方法計算得到淺視電阻率校正后井徑8in下的淺視電阻
   end
   
end
end
Xj=log10(Rlldcc./Rllscc);%轉(zhuǎn)化為侵入半徑校正曲線中的橫坐標(biāo)值
Yj=log10(Rlldcc./Rxo);%轉(zhuǎn)化為侵入半徑校正曲線中的縱坐標(biāo)值
b=length(Xj);
for j=1:b
for i=1:8
    p=S(:,i)';
    Y(j,i)=polyval(p,Xj(j));
    p=zeros(1,3);
end
end
X=[20 30 40 50 60 80 100 120];
for i=1:b
    Di(i)=interp1(Y(i,:),X,Yj(i),'linear','extrap');%插值計算對應(yīng)深度下的泥漿侵入半徑

end

Di=Di*25.4/1000;%將存儲的侵入半徑從in轉(zhuǎn)化為以m為單位
disp(Di);
h=h';

w=[h;Di];
fid=fopen('jieguo.xls','w');%將侵入半徑數(shù)據(jù)存儲到j(luò)ieguo.xls文件中
fprintf(fid,'%12.5f %12.5f\n',w);
fclose(fid);