多曲线 提取标准曲线matlab,Theis直线图解法-特定条件下的标准曲线绘制

本帖最后由 zzbrte123123 于 2015-11-13 18:38 编辑

Theis直线图解法-特定条件下的标准曲线绘制，matlab新手，运行后报错  ??? Error using ==> loglog

Conversion to double from sym is not possible.

完全没有思绪，望各位大大帮忙看看！！！

>> u=[0.0000000001 0.0000000005 0.000000002 0.000000005 0.000000008 0.00000001 0.00000005 0.0000001 0.0000002 0.0000005 0.0000008 0.000001 0.000002 0.000005 0.000008 0.00001 0.00002 0.00005 0.00008 0.0001 0.0002 0.0005 0.0008 0.001 0.002 0.005 0.008 0.01 0.02 0.05 0.06 0.07 0.08 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 4 5 6 8 9 10 20 50];

n=length(u);

syms y;

C=exp(-y)/y;

j=1;

while(j

m1(j)=int(C,y,u(j),inf);

m2(j)=int(C,y,9*u(j),inf);

m3(j)=int(C,y,11*u(j),inf);

m4(j)=int(C,y,19*u(j),inf);

m5(j)=int(C,y,21*u(j),inf);

m6(j)=int(C,y,29*u(j),inf);

m7(j)=int(C,y,31*u(j),inf);

m8(j)=int(C,y,12.04159458*u(j),inf);

m9(j)=int(C,y,15*u(j),inf);

m10(j)=int(C,y,21.9317122*u(j),inf);

m11(j)=int(C,y,21.9317122*u(j),inf);

m12(j)=int(C,y,29.41088234*u(j),inf);

m13(j)=int(C,y,33.24154028*u(j),inf);

m14(j)=int(C,y,13.60147051*u(j),inf);

m15(j)=int(C,y,16.2788206*u(j),inf);

m16(j)=int(C,y,22.82542442*u(j),inf);

m17(j)=int(C,y,22.82542442*u(j),inf);

m18(j)=int(C,y,30.08321791*u(j),inf);

m19(j)=int(C,y,33.83784863*u(j),inf);

m20(j)=int(C,y,20.61552813*u(j),inf);

m21(j)=int(C,y,22.47220505*u(j),inf);

m22(j)=int(C,y,27.58622845*u(j),inf);

m23(j)=int(C,y,27.58622845*u(j),inf);

m24(j)=int(C,y,33.83784863*u(j),inf);

m25(j)=int(C,y,37.21558813*u(j),inf);

m26(j)=int(C,y,22.47220505*u(j),inf);

m27(j)=int(C,y,24.18677324*u(j),inf);

m28(j)=int(C,y,29*u(j),inf);

m29(j)=int(C,y,29*u(j),inf);

m30(j)=int(C,y,35*u(j),inf);

m31(j)=int(C,y,38.27531842*u(j),inf);

m32(j)=int(C,y,30.08321791*u(j),inf);

m33(j)=int(C,y,31.38470965*u(j),inf);

m34(j)=int(C,y,35.22782991*u(j),inf);

m35(j)=int(C,y,35.22782991*u(j),inf);

m36(j)=int(C,y,40.31128874*u(j),inf);

m37(j)=int(C,y,43.18564576*u(j),inf);

m38(j)=int(C,y,32.01562119*u(j),inf);

m39(j)=int(C,y,33.24154028*u(j),inf);

m40(j)=int(C,y,36.89173349*u(j),inf);

m41(j)=int(C,y,36.89173349*u(j),inf);

m42(j)=int(C,y,41.77319715*u(j),inf);

m43(j)=int(C,y,44.55333882*u(j),inf);

m44(j)=int(C,y,8.062257748*u(j),inf);

m45(j)=int(C,y,12.04159458*u(j),inf);

m46(j)=int(C,y,20.02498439*u(j),inf);

m47(j)=int(C,y,20.02498439*u(j),inf);

m48(j)=int(C,y,20.02498439*u(j),inf);

m49(j)=int(C,y,20.02498439*u(j),inf);

a1(j)=vpa((m1(j)),15);

a2(j)=vpa((m1(j)+m2(j)+m3(j)+m4(j)+m5(j)),15);

a3(j)=vpa((m1(j)+m2(j)+m3(j)+m4(j)+m5(j)+m6(j)+m7(j)+m8(j)+m9(j)+m10(j)+m11(j)+m12(j)+m13(j)),15);

a4(j)=vpa((m1(j)+m2(j)+m3(j)+m4(j)+m5(j)+m6(j)+m7(j)+m8(j)+m9(j)+m10(j)+m11(j)+m12(j)+m13(j)+m14(j)+m15(j)+m16(j)+m17(j)+m18(j)+m19(j)+m20(j)+m21(j)+m22(j)+m23(j)+m24(j)+m25(j)),15);

a5(j)=vpa((m1(j)+m2(j)+m3(j)+m4(j)+m5(j)+m6(j)+m7(j)+m8(j)+m9(j)+m10(j)+m11(j)+m12(j)+m13(j)+m14(j)+m15(j)+m16(j)+m17(j)+m18(j)+m19(j)+m20(j)+m21(j)+m22(j)+m23(j)+m24(j)+m25(j)+m26(j)+m27(j)+m28(j)+m29(j)+m30(j)+m31(j)+m32(j)+m33(j)+m34(j)+m35(j)+m36(j)+m37(j)),15);

a6(j)=vpa((m1(j)+m2(j)+m3(j)+m4(j)+m5(j)+m6(j)+m7(j)+m8(j)+m9(j)+m10(j)+m11(j)+m12(j)+m13(j)+m14(j)+m15(j)+m16(j)+m17(j)+m18(j)+m19(j)+m20(j)+m21(j)+m22(j)+m23(j)+m24(j)+m25(j)+m26(j)+m27(j)+m28(j)+m29(j)+m30(j)+m31(j)+m32(j)+m33(j)+m34(j)+m35(j)+m36(j)+m37(j)+m38(j)+m39(j)+m40(j)+m41(j)+m42(j)+m43(j)+m44(j)+m45(j)),15);

a7(j)=vpa((m1(j)+m2(j)+m3(j)+m4(j)+m5(j)+m6(j)+m7(j)+m8(j)+m9(j)+m10(j)+m11(j)+m12(j)+m13(j)+m14(j)+m15(j)+m16(j)+m17(j)+m18(j)+m19(j)+m20(j)+m21(j)+m22(j)+m23(j)+m24(j)+m25(j)+m26(j)+m27(j)+m28(j)+m29(j)+m30(j)+m31(j)+m32(j)+m33(j)+m34(j)+m35(j)+m36(j)+m37(j)+m38(j)+m39(j)+m40(j)+m41(j)+m42(j)+m43(j)+m44(j)+m45(j)+m46(j)+m47(j)+m48(j)+m49(j)),15);

j=j+1;

end

loglog(1./u,a1,1./u,a2,1./u,a3,1./u,a4,1./u,a5,1./u,a6,1./u,a7);

legend('a1','a2','a3','a4','a5','a6','a7',-1);

lh=findall(gca,'type','line');

xc=get(lh,'XData');

yc=get(lh,'YData');