matlab中Error: Missing variable or function.
clearmt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;X=-50:50:250;Y=-50:50:50;[x,y]=...
clear
mt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;
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(y.^-1*exp(-y.-bata^2/(4*y.)),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 =
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0042 Inf 0.0042 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
miu =
0.1005
??? W=int(y.^-1*exp(-y.-bata^2/(4*y.)),miu;inf)
|
Error: Missing variable or function. 展开
mt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;
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(y.^-1*exp(-y.-bata^2/(4*y.)),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 =
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0042 Inf 0.0042 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
miu =
0.1005
??? W=int(y.^-1*exp(-y.-bata^2/(4*y.)),miu;inf)
|
Error: Missing variable or function. 展开
1个回答
展开全部
Error: Missing variable or function 这个问题出在你在bata前的减号前还有一个点,只有点乘点除,没有点加点减得。还有(4*y.),y后面的点也有问题miu;inf中间的分号也有问题。
另外你使用的是积分变量y,但前面y已经被你使用了,是一个常数矩阵,所以这里积分没有变量,是错误的
求C的表达式也会出现逆上面提示的错误,因为W后面也多了个点。
我帮你把上面的错误更正了,但你的积分无法获得表达式解,可能是我没改对的要的公式吧
程序如下:
clear
mt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;
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)
syms t;
A = t.^-1.*exp(-t-bata.^2./(4*t));
W = int(A,'t',miu,100)
C=mt*exp(x.*u./(2*Dl)).*(2*K0-W)/(4*pi*M*n*sqrt(Dl*Dt));
figure(4);
surf(x,y,C);
另外你使用的是积分变量y,但前面y已经被你使用了,是一个常数矩阵,所以这里积分没有变量,是错误的
求C的表达式也会出现逆上面提示的错误,因为W后面也多了个点。
我帮你把上面的错误更正了,但你的积分无法获得表达式解,可能是我没改对的要的公式吧
程序如下:
clear
mt=3.734;u=0.0685;M=20;n=0.3;Dl=0.35;Dt=0.1;t=30;
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)
syms t;
A = t.^-1.*exp(-t-bata.^2./(4*t));
W = int(A,'t',miu,100)
C=mt*exp(x.*u./(2*Dl)).*(2*K0-W)/(4*pi*M*n*sqrt(Dl*Dt));
figure(4);
surf(x,y,C);
本回答被提问者和网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询
广告 您可能关注的内容 |