nm_w_i_k.m

基于概率积分模型的下沉、倾斜、曲率计算

clear;
clc;
 syms x s z qy z2 xx yy tl zx tt;
 syms q d t qx_jf;
%走向各点间距
nm=load('nmdist_bf.txt');
%以开切眼为原点。
nm=nm';nm0=0;
[md,nd]=size(nm);
   for i=1:1:nd
    nm(i)=nm0+nm(i);
    nm0=nm(i);
   end
   nm=nm-39.5;
kc=[115,217,290];
d3=[45 71 115]
%吴侃曲线  
 syms xx yy;
%以开切眼为原点。
w0=4000;r=100;
v=[2.5 3 2.5];
%v=2.5;
yw=110;h=200;lw=8;gw=v*1.8/r;
%rq=163
y_jf=w0/r*int(exp(-pi*(yw-yy+8)^2/r^2),yy,0,300);y_jf=double(y_jf);
% for i=1:1:nd
%     for j=1:1:3 
%      w(i,j)=y_jf/r*int(exp(-pi*(nm(i)-xx)^2/r^2)*(1-exp(-gw(j)*(d3(j)-xx/v(j)))),xx,0,kc(j));     
%      vw(i,j)=gw(j)*y_jf/r*int(exp(-pi*(xw(i)-xx)^2/r^2)*exp(-gw(j)*(d_j(j)-xx/v(j))),xx,0,fy);
%    end
% end
% w=double(w);
% onekc=ones(21,14);
% for i=1:1:3
%     plot(nm,-w(:,i),'b+-')
%     text(200,-30,num2str(i))
%     figure
% end
% for i=1:1:nd
%     plot(d3,vw(i,:),'b+-')
%     text(0,0,num2str(nm(i)))
%     figure
% end

%速度,增加时间次数，求1033点的速度曲线
%kcv=[115 217 277 278.7 285.2 291.2 298 299.5 306 312];
d3=[0 5 10 15 20 25 33 40 45 50 53 55 59 60 63 65 68 70 72 75 78 80 82 85 88 90 92 95 103 104 105 106 107 108 109 110];
[md3 nd3]=size(d3);
kc=2.8*d3;
%vv=[3.5 3.7 2.6 2.8 2.8 2.8 2.8 2.8 2.8 2.8];
v=2.8;gw=v*1.8/r;
%gwv=vv*1.8/r;
for i=1:1:nd
    for j=1:1:nd3 
     w(i,j)=y_jf/r*int(exp(-pi*(nm(i)-xx)^2/r^2)*(1-exp(-gw*(d3(j)-xx/v))),xx,0,kc(j));
     %vw(i,j)=gwv(j)*y_jf/r*int(exp(-pi*(nm(i)-xx)^2/r^2)*exp(-gwv(j)*(d3v(j)-xx/vv(j))),xx,0,kcv(j));
     vw(i,j)=gw*y_jf/r*int(exp(-pi*(nm(i)-xx)^2/r^2)*exp(-gw*(d3(j)-xx/v)),xx,0,kc(j));
   end
end
vw=double(vw);w=double(w);
for i=10:1:36
    plot(nm,-w(:,i)/100,'y+-',nm,vw(:,i),'b+-')
    text(0,0,num2str(kc(i)))
    figure
end
for i=1:1:nd
    plot(d3,vw(i,:),'b+-')
    text(0,0,num2str(nm(i)))
    figure
end
% %下沉以距离为横轴
% %下沉实测、概率、吴侃比较，以时间为横轴
% for i=1:1:21
%     plot(d3,-wg(i,:),'ro-',d3,-wl(i,:),'b+-',d3,-ww(i,:),'y*-')
%     text(100,-30,num2str(i));figure
% end