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