- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一般区域二重、三重积分MATLAB计算方法
这里讨论的计算方法指的是利用现有的MATLAB函数来求解,而不是根据具体的数值计算方法来编写相应程序。目前最新版的2009a有关于一般区域二重积分的计算函数quad2d(详细介绍见 ),但没有一般区域三重积分的计算函数,而NIT工具箱似乎也没有一般区域三重积分的计算函数。 本文目的是介绍一种在7.X版本MATLAB(不一定是2009a)里求解一般区域二重三重积分的思路方法。需要说明的是,上述链接里已经讨论了一种求解一般区域二重三重积分的思路方法,就是将被积函数“延拓”到矩形或者长方体区域,但是这种方法不可避免引入很多乘0运算浪费时间。因此,新的思路将避免这些。由于是调用已有的MATLAB函数求解,在求一般区域二重积分时,效率和2009a的quad2d相比有一些差距,但是相对于延拓函数的做法,效率大大提高了。下面结合一些简单例子说明下计算方法。
diff(x) 如果x是向量,返回向量x的差分;如果x是矩阵,则按各列作差分
diff(x,k) k阶差分
q=polyder(p) 求得由向量p表示的多项式导函数的向量表示q
Fx=gradient(F,x) 返回向量F表示的一元函数沿x方向的导函数F(x),其中x是与F同维数的向量
z=trapz(x,y) x表示积分区间的离散化向量;y是与x同维数的向量,表示被积函数;z返回积分的近似值
z=guad(fun,a,b,tol) 自适应步长Simpson积分法求得Fun在区间[a,b]上的定积分,Fun为M文件函数句柄,tol为积分精度
z=dblquad(fun,a,b,c,d,tol,method) 求得二元函数Fun(x,y)的重积分
z=triplequad(fun,a,b,c,d,e,f,tol,method) 求得三元函数Fun(x,y,z)的重积分 譬如二元函数f(x,y) = x*y,y从sin(x)积分到cos(x),x从1积分到2,这个积分可以很容易用符号积分算出结果
syms x y
int(int(x*y,y,sin(x),cos(x)),1,2) 结果是 -1/2*cos(1)*sin(1)-1/4*cos(1)^2+cos(2)*sin(2)+1/4*cos(2)^2 = -0.635412702399943
如果你用的是2009a,你可以用
quad2d(@(x,y) x.*y,1,2,@(x)sin(x),@(x)cos(x),AbsTol,1e-12)
得到上述结果。 如果用的不是2009a,那么你可以利用NIT工具箱里的quad2dggen函数。 那么我们如果既没有NIT工具箱用的也不是2009a,怎么办呢?答案是我们可以利用两次quadl函数,注意到quadl函数要求积分表达式必须写成向量化形式,所以我们构造的函数必须能接受向量输入。见如下代码
function IntDemo
function f1 = myfun1(x)
? ???f1 = zeros(size(x));
? ???for k = 1:length(x)
? ?? ?? ?f1(k) = quadl(@(y) x(k)*y,sin(x(k)),cos(x(k)));
? ???end
? ???end
y = quadl(@myfun1,1,2)
end
myfun1函数就是构造的原始被积函数对y积分后的函数,这时候是关于x的函数,要能接受向量形式的x输入,所以构造这个函数的时候考虑到x是向量的情况。利用arrayfun函数(7.X后的版本都有这个函数,不了解这个函数的朋友可以查看帮助文档,或者搜索本版)可以将IntDemo函数精简成一句代码:
quadl(@(x) arrayfun(@(xx) quadl(@(y) xx*y,sin(xx),cos(xx)),x),1,2)
上面这行代码体现了用MATLAB7.X求一般区域二重积分的一般方法。可以这么理解这句代码:首先
@(x) arrayfun(@(xx) quadl(@(y) xx*y,sin(xx),cos(xx)),x)
定义了一个关于x的匿名函数,供quadl调用求最外重(x从1到2的)积分,这时候,x对于
arrayfun(@(xx) quadl(@(y) xx*y,sin(xx),cos(xx)),x)
就是已知的了。而
@(xx) quadl(@(y) xx*y,sin(xx),cos(xx))
定义的是对于给定的xx,求xx*y关于y的积分函数,这就相当于数学上积完第一重y的积分后得到一个关于xx的函数而
arrayfun(@(xx) quad
您可能关注的文档
最近下载
- baltur燃烧器TBG210P17690030中文使用维护手册.pdf VIP
- 审核员现场见证评价表参考实用文档.doc VIP
- 2023年军队文职人员招聘之军队文职公共科目押题练习试卷A卷附答案.docx VIP
- ai心理健康创业计划书.docx VIP
- 国开作业管理学基础-管理实训:第二章 查阅文献资料并写出评论参考(含答案)025.docx VIP
- 医学综合试题库及答案.docx VIP
- 通用工器具安全管理规范.docx VIP
- 2025届高考作文技巧之列提纲优化结构课件(共46张PPT).pptx VIP
- 光学装配与调校技术课件.pptx VIP
- 医学综合知识试题库+答案.pdf VIP
文档评论(0)