www.ndwl.net > mAtlAB求二重积分

mAtlAB求二重积分

syms x y; f_1=0.2*exp(-0.5*((x-16)^2)+((y-10)^2)/0.64-((x-16)*(y-10))/2); xup=y/4+sqrt(6-1.5*y^2); xlow=y/4-sqrt(6-1.5*y^2); fs1=int(f_1,'x',xlow,xup); fs2=int(fs1,'y',-2,2); fs2 = int((exp(-(135*y)/4)*exp(675/4)*exp((27*y^2)/1...

用下列命令,可以求得其积分值。 >> syms x y >> s=int(2*y*exp(-y^2)*(int(exp(-x^2)/(x^2+y^2),x,-1,1))^2,0.2,1); >> vpa(s) 运行结果 ans =10.213463733952979925429035801331

程序: fun=sin(0.5*pi*x./y);%% a=int(int(fun,y,sqrt(x),x),x,1,2); b=simple(a) %化简 I=vpa(b,4) %得到4位近似解,也可以任意N位解 结果: I = 0.2719

没看明白你前面写的那一串是什么? 那个积分两条语句就可以算出来的: >> f=@(x,y)sin(x)./y;>> dblquad(f,pi/4,1,2,4)ans = 0.1156

假设y的范围为[-10,10],则获取F与y关系的代码及图像如下所示: >> fun = @(t,y,z)-20*sin(t).^2./((20*sin(t)).^2+(y-20*sin(t).^2).^2+(0.7-z).^2).^1.5; % 包含参数的被积函数>> y = -10:0.2:10; % 假设参数y的范围[-10,10]>> F = arrayfun(@...

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)an...

使用matlab的int函数可以方便的计算积分,以及多重积分。 设二重积分还是表达式为 z=z(x,y),积分域为下限 y1(x) 上限 y2(x),从 x1 到 x2,则二重积分代码为: int(int(z,y,y1,y2),x,x1,x2) 需要先定义符号变量 x,y,以及表达式 z,y1,y2 和...

从2012a开始,引入了几个计算数值积分的新函数integral、integral2、integral3,可以计算积分限为无穷大的积分。在2013a的文档中,建议使用这几个函数代替quad、quadl、quadv、dblquad、triplequad,后面几个函数在将来的版本中可能会被删除。 ...

这个可以用matlab的符号积分或者数值积分解决,下面提供4种方法: 1、直角坐标系符号积分 syms x yint_y = int(sin(pi*(x^2+y^2)),y,-sqrt(1-x^2),sqrt(1-x^2)); % 先对y积分I = vpa(int(int_y,x,-1,1)) I = 2.0 2、极坐标系符号积分 syms rI = ...

网站地图

All rights reserved Powered by www.ndwl.net

copyright ©right 2010-2021。
www.ndwl.net内容来自网络,如有侵犯请联系客服。zhit325@qq.com