原创力文档- 1、本文档共35页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
问题背景和实验目的; 矩形法;矩形法;;解：; 理论值：;; 曲边小梯形的面积可以由直边小梯形的面积来近似; 如果我们 n 等分区间 [a,b]，即令：;解：; 2n 等分区间 [a,b] ，得; 设过以上三点的抛物线方程为：; 同理可得：; 整理后可得：;==>; 梯形法：trapz;前面的做法;quad(f,a,b,tol)
f = f(x) 为被积函数，[a,b] 为积分区间，tol 为计算精度;解：; 抛物线法计算二重积分： dblquad;>> f=inline('4*x*y+3*y^2');
>> I=dblquad(f,-1,1,0,2);例：计算二重积分;int(f,a,b)
计算 f 关于默认自变量 的定积分，积分区间为[a,b]。;>> syms x y;
>> f=y*sin(x);
>> int(f,x)
>> int(f,y)
>> int(f)
>> int('a+b'); 例：用 int 函数计算定积分：;double(a)
将 a 转化为双精度型，若 a 是字符，则取对应的 ASCII 码;>> x=1:0.001:2;
>> y=exp(x.^(-2));
>> trapz(x,y); 抛物线法：; 教材第 75 页：2、3、6;; 设过以上三点的抛物线方程为：; 整理后可得：; 抛物线法计算二重积分： dblquad;>> f=inline('4*x*y+3*y^2');
>> I=dblquad(f,-1,1,0,2);int(f,a,b)
计算 f 关于默认自变量 的定积分，积分区间为[a,b]。
文档评论（0）