高数上机作业答案.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
高数上机作业答案

1.(两个重要极限)计算下列函数的函数值并画出图形,观察两个重要极限值 (1)y=f(x)= 应用matlab编写程序做出图形 x=[0:0.01:8*pi]; ezplot(y-sin(x)/x);grid on 得到图形: 可明显观察到该函数趋近于0时极限值为1 当x趋向于正负无穷时极限值为0 (2)y=f(x)= 应用matlab编写程序做出图形 x=[0:0.001:8*pi]; ezplot(y-(1+x)^(1/x));grid on 得到图形: 可明显观察到该函数正向趋近于-1时极限值为正无穷, 负向趋近-1时极限为负无穷 当x趋向于正负无穷时极限值为1 2、人造地球卫星轨道可视为平面上的椭圆。我国的第一颗人造地球卫星近地点距离地球表面439km,远地点距地球表面2384km,地球半径为6371km,求该卫星的轨道长度。本问题可用椭圆参数方程: ???????????????????????? ???????????? ????????????????????????? ( 0 ≤ t ≤ 2 π ) ????????????? a , b 0 分别为椭圆的长短半轴, ??? 来描述人造地球卫星的轨道,该轨道的长度L就是如下的参数方程弧长积分: #includestdio.h #includemath.h #define p 3.1415926*2 #define nint main() { double a=7782.5, b=7721.5, c=0.0; int i; for (i=1;in;i++) c=c+(sqrt(a*a*sin((p*i)/n)*sin((p*i)/n)+b*b*cos((p*i)/n)*cos((p*i)/n)))*p/n; printf(%lf,c); return 0; } 答案:该卫星轨道的长度为48707.435320km。 3、(最大值问题)在一条宽20 m的道路两侧,分别安装了一只2 kW和一只3 kW的路灯,它们离地面的高度分别为5 m和6 m.在漆黑的夜晚,当两只路灯开启时,两只路灯连线的路面上最暗的点和最亮的点分别在哪里? K为比例系数,取1, P为光源功率; r为地面上一点M到光源的距离; a为光源到点M的光线与地面的夹角 两只路灯连线的路面上[xmin,fmin]=fminbnd(10/(25+x^2)^(3/2)+18/(436-40*x+x^2)^(3/2),0,20) xmin = 9.3383 fmin = 0.0182 [xmin,fmin]=fminbnd(-10/(25+x^2)^(3/2)-18/(436-40*x+x^2)^(3/2),0,20) xmin = 19.9767 fmin = -0.0845 加粗字体为输入程序,由于matlab中只能计算最小值,因而用-f(x)来找该函数的最大值,得出结论为,当两只路灯开启时,两只路灯连线的路面上最暗的点9.3383处,最亮的点19.9767处。 文件1jisi.m %Creat the function for jisi %Let x(1)=x, x(2)=y function dx=jisi(t,x,a,b,c) s=sqrt((c-x(1))^2+(a*t-x(2))^2); dx=[b*(c-x(1))/s;b*(a*t-x(2))/s]; 文件2seajisi.m clear,clf,shg %Set the definied time %ts=0:0.05:.5; % ts=0:0.1:1.6; n=length(ts); x0=[0 0]; %a=20;b=40;c=15; % a=35;b=40;c=15; % opt=odeset(RelTol,1e-6,AbsTol,1e-9); %[t,x]=ode45(@jisi,ts,x0,[],a,b,c); % [t,x]=ode45(@jisi,ts,x0,opt,a,b,c); %exact solution x1=c y1=a*t; %output t,x(t),y(t) and draw x(t),y(t) [t,x,y1] plot(t,x),grid,gtext(x(t),FontSize,16), gtext(y(t),FontSize,16),pause %draw y(x): the position of tatch jisi plot(x(:,1),x(:,2),r*),grid xlabel(x,FontSize,16),ylabel(y,FontSize,16)

文档评论(0)

yan698698 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档