原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
Q1.习题4第4题
用integral求取∫−||||的绝对精度为−的广义积分，并尝试用trapz及
−
符号计算此积分。（提示：注意integral指令相对误差控制对绝对精度的影响。本题假设
符号计算保留16位小数为真实值，trapz步长自定）
答：首先利用符号计算获取一个近似的真实值
symsxy(x)
y(x)=exp(-abs(x))*abs(sin(x));
si=vpa(int(y,-5*pi,10*pi),16)
si=
1.090331328569942
然后使用函数integral获得较精确的数值积分解
fx=@(x)exp(-abs(x)).*abs(sin(x));
formatlong;
s=integral(fx,-5*pi,10*pi,Abstol,1e-9)%如果要提高精度还需减小Restol
s=
1.090331540770801
formatshorte
err1=double(abs(s-si))
err1=
2.1220e-07
最后使用粗略的梯形数值积分函数trapz计算此积分，小步长速度慢，但误差更小
x=linspace(-5*pi,10*pi,1e5);
formatlong;
dx=x(2)-x(1);
st=trapz(exp(-abs(x)).*abs(sin(x)))*dx
st=
1.090331292961167
formatshorte
err2=double(abs(st-si))
err2=
3.5609e-08
实验结果表明，虽然integral在相同步长精度下往往误差更小，在相同误差情况下效率更
高，但是如果trapz函数设定了更小的间距或的分隔步骤，可以获得更高的精度。
Q2.的估计，如图，若设=,且圆心为(,)可生成×组符合(,)均匀分布的横坐
标与纵坐标。,)恰好落在正方形区域内。此时统计落在圆形区域(−)+
(−)≤内的点数,即可近似估计≈⋅。用实现这段代码。
答：本题只需将对应的算法实现即可，注意：无需使用for循环否则会降低运行的效率！
rngdefault;%为了结果可以再现，其实并不必要
x=rand(1,1e6);
y=rand(1,1e6);
r2=(x-1/2).^2+(y-1/2).^2;
Na=sum(r2=1/4);%ra=1/4为逻辑矩阵，符合条件的为1，求和即为计数
pi_app=4*Na/1e6
pi_app=
3.141648000000000
()()
Q3.使用diff与gradient函数，计算
文档评论（0）