下面这个程序中Srthzt=AA*BB*EE,AA和BB都可以计算出来,但是EE这个积分计算得到的确是一个函数表达式,matlab为什么不能计算结果呢?
r=20;
thta=90;
z=10;
t=100;
ql=80;
alf=2;
ar=0.3;
az=0.3;
mu=0.004;
thtaw=45;
zw=5;
lw=2;
h=80;
syms zz n
AA=ql/(8*mu*sqrt((pi^3)*(ar^2)*az*(t^3))*sind(thtaw));
BB=exp(-alf*(r^2)/(4*az*t))/sqrt(t^3);
CC=exp(-(alf*(zz^2)*(cotd(thtaw))^2-2*alf*r*zz*cotd(thtaw)*cosd(thta)+(z-zw-zz)^2)/(4*az*t));
DD=exp(-(alf*(zz^2)*(cotd(thtaw))^2+2*alf*r*zz*cotd(thtaw)*cosd(thta)+(z+zw-zz)^2)/(4*az*t));
CD=CC+DD;
zz1=-lw*sind(thtaw)/2;
zz2=lw*sind(thtaw)/2;
EE=int(CD,zz,zz1,zz2);
Srthzt=AA*BB*EE
Srthzt =
-2902829188123145/590295810358705651712*pi^(1/2)*exp(-5/36)*10^(1/2)*erf(-1/40*2^(1/2)*10^(1/2)+1/12*10^(1/2))-2902829188123145/590295810358705651712*pi^(1/2)*exp(-5/4)*10^(1/2)*erf(-1/40*2^(1/2)*10^(1/2)+1/4*10^(1/2))+2902829188123145/590295810358705651712*pi^(1/2)*exp(-5/36)*10^(1/2)*erf(1/40*2^(1/2)*10^(1/2)+1/12*10^(1/2))+2902829188123145/590295810358705651712*pi^(1/2)*exp(-5/4)*10^(1/2)*erf(1/40*2^(1/2)*10^(1/2)+1/4*10^(1/2))
>> EE
EE =
-pi^(1/2)*exp(-5/36)*10^(1/2)*erf(-1/40*2^(1/2)*10^(1/2)+1/12*10^(1/2))-pi^(1/2)*exp(-5/4)*10^(1/2)*erf(-1/40*2^(1/2)*10^(1/2)+1/4*10^(1/2))+pi^(1/2)*exp(-5/36)*10^(1/2)*erf(1/40*2^(1/2)*10^(1/2)+1/12*10^(1/2))+pi^(1/2)*exp(-5/4)*10^(1/2)*erf(1/40*2^(1/2)*10^(1/2)+1/4*10^(1/2))