matlab二重积分.docVIP

一 使用两次一重积分 %%二重积分f= @(x,y)exp(sin(x))*ln(y)，y从5*x积分到x^2，x从10积分到20 1 (7.X后版本才有此函数quad2d) y1=quad2d(@(x,y) exp(sin(x)).*log(y),10,20,@(x)5*x,@(x)x.^2) 2 y2 =quadl(@(x) arrayfun(@(x) quadl(@(y)exp(sin(x)).*log(y),5*x,x.^2),x),10,20) 3 y3 = dblquad(@(x,y)exp(sin(x)).*log(y).*(y=5*x y=x.^2),10,20,50,400) 详细请看吴鹏老师的文章/viewthread.php 二 使用dblquad函数 q = dblquad(fun,xmin,xmax,ymin,ymax,tol,method) 该函数求f(x,y)在[a,b]×[c,d]区域上的二重定积分。参数tol，trace的用法与函数quad完全相同。 例8-5 计算二重定积分 (1) 建立一个函数文件fxy.m： function f=fxy(x,y) global ki; ki=ki+1; %ki用于统计被积函数的调用次数 f=exp(-x.^2/2).*sin(x.^2+y); (2) 调用dblquad函数求解。 global ki;ki=0; I=dblquad(fxy,-2,2,-1,1) ki I = 1.57449318974494 ki = 1038 来源 精通MATLAB科学计算 一书王正林,龚纯,何倩编写,电子工业出版社 三 复合辛普森公式(矩形积分区域) function q=DblSimpson(f,a,A,b,B,m,n) if(m==1 n==1) %辛普森公式 q=((B-b)*(A-a)/9)*(subs(sym(f),findsym(sym(f)),{a,b})+... subs(sym(f),findsym(sym(f)),{a,B})+... subs(sym(f),findsym(sym(f)),{A,b})+... subs(sym(f),findsym(sym(f)),{A,B})+... 4*subs(sym(f),findsym(sym(f)),{(A-a)/2,b})+... 4*subs(sym(f),findsym(sym(f)),{(A-a)/2,B})+... 4*subs(sym(f),findsym(sym(f)),{a,(B-b)/2})+... 4*subs(sym(f),findsym(sym(f)),{A,(B-b)/2})+... 16*subs(sym(f),findsym(sym(f)),{(A-a)/2,(B-b)/2})); else %复合辛普森公式 q=0; for i=0:n-1 for j=0:m-1 x=a+2*i*(A-a)/2/n; y=b+2*j*(B-b)/2/m; x1=a+(2*i+1)*(A-a)/2/n; y1=b+(2*j+1)*(B-b)/2/m; x2=a+2*(i+1)*(A-a)/2/n; y2=b+2*(j+1)*(B-b)/2/m; q=q+subs(sym(f),findsym(sym(f)),{x,y})+... subs(sym(f),findsym(sym(f)),{x,y2})+... subs(sym(f),findsym(sym(f)),{x2,y})+... subs(sym(f),findsym(sym(f)),{x2,y2})+... 4*subs(sym(f),findsym(sym(f)),{x,y1})+... 4*subs(sym(f),findsym(sym(f)),{x2,y1})+... 4*subs(sym(f),findsym(sym(f)),{x1,y})+... 4*subs(sym(f),findsym(sym(f)),{x1,y2})+... 16*subs(sym(f),findsym(sym(f)),{x1,y1}); end end end q=((B-b)*(A-a)/36/m/n)*q; 四 复合梯形公式(矩形积分区域) function q=DblTraprl(f,a,A,b,B,m,n) if(m==1 n==1) %梯形公式 q=((B-b)*(A-a)/4)*(subs(sym(f),findsym(sym(f)),{a,b})+... subs(sym(f),findsym(sym(