- 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求解高等数学问题-1
姓名: 学号: 专业: 成绩:
一、实验名称:用MATLAB求解高等数学基本计算问题-1
二、实验目的:
1.学会利用MATLAB求解极限问题;
2.学会利用MATLAB求解导数、偏导数和高阶导数问题;
3. 学会利用MATLAB求解有关不定积分,定积分问题;
4. 学会利用MATLAB求解二重积分和三重积分问题;
5. 学会利用MATLAB求解极值和最值问题;
6. 学会利用MATLAB解决微分法在几何上的应用;
7. 学会利用MATLAB求解有关泰勒级数展开问题。
三、实验准备:
复习高等数学中相关主题的理论知识和方法。
四、实验内容
1.求极限;
2.设,求;
3.计算下列积分:
(1); (2);
4.计算三重积分,其中是由面上的直线及曲线所围成的平面图形绕轴旋转一周而成的空间闭区域。
注:需要画出被积区域的空间图形。
5.作出函数的图形,观察极限点、最值点的位置并求出,给出函数的单调区间;
6.求二元函数在点附近的极小值;
7.求圆柱螺旋线在对应点处的切线方程和法平面方程;
8.求球面在点处的切平面及法线方程;
9.用命令观测函数的麦克劳林展开式的第4,8,12,16项,然后在同一坐标系作出函数和它展开式的前几项构成的多项式函数的图形,观测这些多项式函数的图形向的图形逼近情况。
10.某个阿拉伯国家有一座著名的伊斯兰清真寺,它以中央大厅的金色巨大拱形圆顶名震遐迩。因年久失修,国王下令将清真寺顶部重新贴金箔装饰。据档案记载,大厅的顶部形状为半球面,其半径为30米。考虑到可能的损耗和其他技术因素,实际用量将会比清真寺顶部面积多1.5%。据此,国王的财政大臣拨出了可制造5800平方米有规定厚度金箔的黄金。建筑商人哈桑略通数学,他计算了一下,觉得黄金会有盈余。于是他以较低的承包价得到了这项装饰工程。但在施工前的测量中,工程师发现清真寺顶部实际上并非一个精确的半球面,而是一个半椭球面,其半立轴恰是30米,而半长轴和半短轴分别是30.6米和29.6米。这一来,哈桑犯了愁,他担心黄金是否还有盈余?甚至可能短缺,最后的结果究竟如何呢?
提示:若在直角坐标系下直接积分有困难,可考虑下面三种方法:
方法一:在极坐标系下进行积分,考虑广义极坐标法;
方法二:在直接坐标系下积分,考虑蒙特卡洛方法(随机投点法);
方法三:数值积分方法。
五、程序及实验结果
1. 程序:
syms x
limit(sqrt(x^2+x)-x,x,inf)
输出结果:
ans =1/2
2. 程序:
syms x y z zx zxx zxy
z=atan(x+y);
zx=diff(z,x);
zxx=diff(zx,x)
zxy=diff(zx,y)
输出结果:
zxx =-(2*x + 2*y)/((x + y)^2 + 1)^2
zxy =-(2*x + 2*y)/((x + y)^2 + 1)^2
3.(1)程序:
clear
syms x y
y=(sec(x))^3
f=int(y,x)
输出结果:f =log(tan(pi/4 + x/2))/2 + tan(x)/(2*cos(x))
(2)程序:
syms t r
s=int((r*cos(t)+r*sin(t)+1)*r,r,0,2*cos(t));
q=int(s,t,-pi/2,pi/2)
输出结果:q =2*pi
4. 程序:
clear;
z=linspace(0,2,40);
t=linspace(0,2*pi,40);
[z1,t1]=meshgrid(z,t);
x=sqrt(1+(z1-1).^2).*cos(t1);
y=sqrt(1+(z1-1).^2).*sin(t1);
surf(x,y,z1);
syms z
m=1/2*int(pi*(1+(z-1)^2)^2,z,0,2)
输出结果:m=(28*pi)/15
5. 程序:
syms x
ezplot(x^2*sin(x^2-x-2),[-2,2]);
[x,f]=fminbnd(x^2*sin(x^2-x-2),-2,2)
[x1,f1]=fminsearch(-x^2*sin(x^2-x-2),-1.554)
[x2,f2]=fminsearch(x^2*sin(x^2-x-2),-0.754)
[x3,f3]=fminsearch(-x^2*sin(x^2-x-2),0)
[x4,f4]=fminsearch(x^2*sin(x^2-x-2),1.62)
syms x
m=subs(x^2*sin(x^2-x-2),x,-2)
n=eval(m
原创力文档


文档评论(0)