基matlab背景基础之三数值分析及其matlab实现.ppt

矩形域上计算二重积分的命令： dblquad(@fun,xmin,xmax,ymin,ymax,tol) 广义积分、二重和三重积分 长方体上计算三重积分的命令： triplequad(@fun,xmin,xmax,ymin,ymax, zmin,zmax,tol) 例 计算二重定积分 子函数： (1) 建立一个函数文件fxy.m： function f=fxy(x,y) f=exp(-x.^2/2).*sin(x.^2+y); (2) 调用dblquad函数求解 I=dblquad(fxy,-2,2,-1,1) I =1.57449318974494 匿名函数： f=@(x,y)(exp(-x.^2/2).*sin(x.^2+y)); I=dblquad(f,-2,2,-1,1) I =1.5745 例 计算三重定积分命令如下：fxyz=@(x,y,z)(4*x.*z.*exp(-z.*z.*y-x.*x));triplequad(fxyz,0,pi,0,pi,0,1,1e-7)ans= 1.7328 二重积分的计算 抛物线法计算二重积分： dblquad dblquad(f,a,b,c,d,tol) tol 为计算精度，若不指定，则缺省精度为 10-6 f(x,y) 可以由 inline 定义，或通过一个函数句柄传递 [a,b] 是第一积分变量的积分区间，[c,d] 是第二积分变量 的积分区间 按字母顺序，大写字母排在小写字母的前面 f=inline(4*x*y+3*y^2); I=dblquad(f, -1,1,0,2) f(x,y) 中关于第一自变量的运算是数组运算， 即把 x 看成是向量，y 看成是标量。 也可以全部采用数组运算 例：计算二重积分 dblquad(inline(4*x*y+3*x^2),-1,1,0,2) dblquad(inline(4*x*y+3*x.^2),-1,1,0,2) X 例：计算二重积分 dblquad 举例 例：计算二重积分 dblquad(@(x,y)4*x*y+3*x.^2 , -1,1, 0, 2) 指定 x、y 分别是第一和第二积分变量 dblquad(inline(4*x*y+3*x.^2), -1,1, 0, 2) 被积函数 f (x,y) 的另一种定义方法：匿名函数 dblquad(@(y,x)4*x*y+3*x.^2 , -1,1, 0, 2) 下面的命令运行结果和上面的一样吗？ dblquad 举例 其中 利用梯形法，先将[a,b]区间n 等分，hx=(b-a)/n, xi=a+ihx, i=0,1,…,n A a b d(x) c(x) 重积分的数值计算可通过 单积分组合计算 非矩形区域上的二重积分 非矩形区域重积分dblquad2.m function S=dblquad2(fun,a,b,clo,dhi,n) if nargin6,n=100;end x=linspace(a,b,n+1); S=0; for i=1:n S=S+dblquad(fun,x(i),x(i+1),feval(clo, … (x(i)+x(i+1))/2),feval(dhi,(x(i)+x(i+1))/2)); end M文件dblquad2.m给出二重积分计算法。 I=dblquad2(fun,a,b,clo,dhi, n) fun为被积函数f(x,y) clo和dhi是y的下限和上限函数c(x),d(x) a,b分别为x的下限和上限; n 为区间的等分数(默认100) 例 解答 调用函数dblquad2(@fun1,-1,1,@fun2,@fun3) function yy=fun1(x,y) yy=sqrt(1-x.^2); function y=fun2(x) y=-sqrt(1-x^2); function y=fun3(x) y=sqrt(1-x^2); I=dblquad2(fun, a, b, clo, dhi, n) 三重积分 求得三元函数Fun(x,y,z)的重积分, 格式类似dblquad。 z=triplequad(Fun,a,b,c,d,e,f) 例如：计算： triplequad(@(x,y,z)y*sin(x)+z*cos(x)… ,-1,1,0,1,0,pi) 练习 1，分别用quad,trapz,quadl积分法计算下面的定积分，并比较结果精度。 2、计算下列积分 作业 练习 3、计算重积分 其中D是由 所围成的区域 4