原创力文档- 1、本文档共35页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
数学软件与数学实验;问题背景和实验目的;矩形法;矩形法;步长;解：;理论值：;定积分几何意义;曲边小梯形的面积可以由直边小梯形的面积来近似;如果我们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]。;symsxy;
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);抛物线法：;用蒙特卡罗(MonteCarlo)法近似计算;蒙特卡罗法近似计算定积分-样本平均值法;特别地，假设a，b有限，可取g(x)为[a,b]上均匀分布．此时,设x1,…,xn是来自U(a,b)的随机数，那么?的一个估计为;例3计算定积分
事实上，其精确解为
样本平均值法求解：
注增加样本数目，可提高计算精度，但计算时间也会提高。
ybjzf(0,4,1000)result=7.1854
ybjzf(0,4,10000)result=7.2153
ybjzf(0,4,100000)result=7.2419;functionresult=ybjzf(a,b,mm)
%a是积分的下限
%b是积分的上限
%积分函数cos(x)+2
%mm是随机实验次数
sum=0;
xrandnum=unifrnd(a,b,1,mm);;
文档评论（0）