原创力文档- 1、本文档共13页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
Matlab实验报告
西安交通大学
车辆01班
教师：阮小娥
郝杰
2011.06.22
一、实验问题：
开普勒方程近似解与方程求根
天文学中，有一类著名的方程—开普勒方程,x=qsinx+a(0q1,a
为常数)，开普勒方程是用来确定行星在其运动轨道上的位置的。如
何求解方程并使其解达到一定的精度要求呢？
开普勒方程近似解与方程求根。
二、问题的分析、
（1）理论知识：
开普勒方程是一个超越方程，很难得出严格的分析解，但是，
已经证明这个方程存在唯一解。如果已知某一作椭圆运动的天体
的轨道要素，利用二体问题的关系式可以得到任意给定时刻t时
的平近点角M，而后采用图解法、数值法或近似迭代法求解开普
勒方程得出偏近点角E，再利用二体问题的其他积分而得到t时
刻天体在轨道上的坐标和速度。对于抛物线轨道和双曲线轨道也
有相应的开普勒方程。首先为了求解开普勒方程，通常将求函数
f（x）的零点或方程f（x）=0的根的问题。其次，在平面坐标系
中绘制函数f（x）的曲线，大致了解f（x）的零点情况和位置，
如果有多个零点，则要分区逐个求解。二体问题运动方程的一个
积分。它反映天体在其轨道上的位置与时间t的函数关系。对于
椭圆轨道，开普勒方程可以表示为E－esinE＝M，式中E为偏近
点角，M为平近点角，都是从椭圆轨道的近日点开始起算，沿逆
时针方向为正，E和M都是确定天体在椭圆轨道上的运动和位置
的基本量。如果定义天体在轨道上运动的平均角速度为n，天体
过近日点的时刻为τ，则对任一给定时刻t，天体从近日点出发
所走过的角度就是平近点角M＝n（t－τ）。这样，开普勒方程给
出了天体在轨道上运动的位置与时间t的关系。
三、实验目的：
通过对开普勒方程求根问题的解决，熟悉非线性方程的几种数
值方法，并进行试验。
数学建模及程序设计及结果展示：
Number1.用二分法求方程x=0.5sinx+1的近似根（误差小于
0.00001）
程序设计：
f=inline(x-0.5*sin(x)-1);
a=1;
b=2;
dlt=1.0e-5;
k=1;
whileabs(b-a)dlt
c=(a+b)/2;
iff(c)==0
break;
elseiff(c)*f(b)0
a=c;
else
b=c;
end
fprintf(k=%d,x=%.5f
,k,c);
k=k+1;
end
实验结果：
Number2：用“切线法”求方程x=0.5sinx+1的近似根（误差小于
0.00001）
.
程序设计：
f=inline();
x-0.5*sin(x)-1
df=inline(1-0.5*cos(x));
d2f=inline(0.5*sin(x));
a=1;
b=2;
dlt=1.0e-5;
iff(a)*d2f(a)0
x0=a;
else
x0=b;
end
m=min(
文档评论（0）