xihuangdongling
2011-07-31, 22:34
求助：请把下面的公式写成matlab语言
注意：其中K0（β）可以调用matlab里的K0=besselk(0,bata)
W可以根据下面的两个方法任选一个
方法一
方法二
我自己写的不正确如下：clear
mt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;
syms v
X=-50:50:250;Y=-50:50:50;
[x,y]=meshgrid(X,Y);
bata=sqrt(u^2*x.^2./(4*Dl^2)+u^2*y.^2./(4*Dl*Dt));
K0=besselk(0,bata)
miu=u^2*t/(4*Dl)
W=int((exp(-v-bata.^2./(4*v))/v,miu,inf)
C=mt*exp(x.*u./(2*Dl))*(2*K0-W)./(4*pi*M*n*sqrt(Dl*Dt));
figure(4);
surf(x,y,C);
注意：其中K0（β）可以调用matlab里的K0=besselk(0,bata)
W可以根据下面的两个方法任选一个
方法一
方法二
我自己写的不正确如下：clear
mt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;
syms v
X=-50:50:250;Y=-50:50:50;
[x,y]=meshgrid(X,Y);
bata=sqrt(u^2*x.^2./(4*Dl^2)+u^2*y.^2./(4*Dl*Dt));
K0=besselk(0,bata)
miu=u^2*t/(4*Dl)
W=int((exp(-v-bata.^2./(4*v))/v,miu,inf)
C=mt*exp(x.*u./(2*Dl))*(2*K0-W)./(4*pi*M*n*sqrt(Dl*Dt));
figure(4);
surf(x,y,C);