整理matlab二重积分.pptxVIP

一 使用两次一重积分 %%二重积分 f= @(x,y)exp(sin(x))*ln(y)，y 从 5*x 积分到x^2，x 从 10 积分到 20 (7.X 后版本才有此函数quad2d) y1=quad2d(@(x,y) exp(sin(x)).*log(y),10,20,@(x)5*x,@(x)x.^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 计算二重定积分 建立一个函数文件 fxy.m： function f=fxy(x,y) global ki; ki=ki+1; %ki 用于统计被积函数的调用次数 f=exp(-x.^2/2).*sin(x.^2+y); 调用dblquad 函数求解。 global ki;ki=0; I=dblquad(fxy,-2,2,-1,1) ki I = 1.57449318974494 ki = 1038 来源 精通MATLAB 科学计算 一书王正林,龚纯,何倩编写,电子工业出版社;三 复合辛普森公式(矩形积分区域);3;4