matlab二重积分求解
现在有这样一个函数Q=83;H=180;u=3;sigy=0.16.*x./(1+0.0001.*x).^0.5;sigz=0.12.*x;c=Q./(pi*sigy.*...
现在有这样一个函数
Q=83;
H=180;
u=3;
sigy=0.16.*x./(1+0.0001.*x).^0.5;
sigz=0.12.*x;
c=Q./(pi*sigy.*sigz*u+eps).*exp(-0.5*(y./(sigy+eps)).^2).*exp(-0.5*(H./(sigz+eps)).^2);
现在我想求解它在x[-9999，10000]和y[-9999,10000]平面上的积分。
函数三维模型如图
ps；我用a=int(int(c,x,-9999,10000),y,-9999,10000);求不出来啊，老是报错。
希望能给我最后的结果，以及方法。 展开
Q=83;
H=180;
u=3;
sigy=0.16.*x./(1+0.0001.*x).^0.5;
sigz=0.12.*x;
c=Q./(pi*sigy.*sigz*u+eps).*exp(-0.5*(y./(sigy+eps)).^2).*exp(-0.5*(H./(sigz+eps)).^2);
现在我想求解它在x[-9999，10000]和y[-9999,10000]平面上的积分。
函数三维模型如图
ps；我用a=int(int(c,x,-9999,10000),y,-9999,10000);求不出来啊，老是报错。
希望能给我最后的结果，以及方法。 展开
1个回答
展开全部
Q = 83;
H = 180;
u = 3;
sigy = @(x)0.16*x./sqrt(1+x/1e4);
sigz = @(x)0.12*x;
c = @(x, y)Q/(pi*sigy(x).*sigz(x)*u) ...
.*exp(-(y/(sigy(x))).^2/2) ...
.*exp(-(H/(sigz(x))).^2/2);
quadv(@(y)quadv(@(x)c(x,y), -1e4+1, 1e4), -1e4+1, 1e4)
ans =
717
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
