《迭代法求解开普列方程.doc》由会员分享，可在线阅读，更多相关《迭代法求解开普列方程.doc（11页珍藏版）》请在金锄头文库上搜索。
1、数值分析实验报告迭代法求解开普列方程作者：罗波 学号：2008203020081 课题题目开普列方程是天体力学中的一个重要方程，其中，x是某个行星的平均近点角，y是它的偏近点角，是它的轨道的偏心率。取=0.9，对区间0,中30个等分点的x计算出y的数值解。2 课题背景开普勒定律是开普勒发现的关于行星运动的定律。他于1609年在他出版的新天文学上发表了关于行星运动的两条定律，又于1618年，发现了第三条定律。开普勒很幸运地能够得到，著名的丹麦天文学家第谷布拉赫所观察与收集的，非常精确的天文资料。大约于1605年，根据布拉赫的行星位置资料，开普勒发现行星的移动遵守三条相当简单的定律。开普勒的定律给
2、予亚里士多德派与托勒密派在天文学与物理学上极大的挑战。他主张地球是不断地移动的；行星轨道不是周转圆的，而是椭圆形的；行星公转的速度不等恒。这些论点，大大地动摇了当时的天文学与物理学。经过了几乎一世纪披星戴月，废寝忘食的研究，物理学家终于能够用物理理论解释其中的道理。牛顿利用他的第二定律和万有引力定律，在数学上严格地证明开普勒定律，也让人们了解其中的物理意义。3 专业背景迭代算法是用计算机解决问题的一种基本方法。它利用计算机运算速度快、适合做重复性操作的特点，让计算机对一组指令（或一定步骤）进行重复执行，在每次执行这组指令（或这些步骤）时，都从变量的原值推出它的一个新值。牛顿迭代法（Newton
3、s method）又称为牛顿-拉夫逊方法（Newton-Raphson method），它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x)=0的根。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x)=0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。另外该方法广泛用于计算机编程中。对于不动眯迭代法，其迭代格式简单，通过对原方程变形构造通过重复计算能够得到满足一定精度的数值解，但由于其对于函数要求高，迭代
4、次数多，对于牛顿鬼迭代法和不动点迭代法都有改进，得到收敛性较好的迭代格式。4 问题分析对于超越方程无法求得其解析解，必须通过数值方法求其解，对于不同x、值，会造成方程有不同的收敛速度，经过Plot画图可以得出函数只有一种类型的拐点，对于牛顿迭代法其初值可以随便选择，通过画图选择其初始点为1.0，分别构造不动点格式与牛顿迭代格式，分别进行迭代，得到结构，比较其迭代次数，而得出牛顿迭代法与不动点迭代法其性能上差别。5 程序及注释不动点迭代法程序以及程序注释e=0.9; %给变量取值format long %设定数据格式n=29; %0 pi中取值个数x=0:pi/n:pi; %给出向量值k=0;
5、%记录迭代法数er=1e-8; %给出误差限ee=1; y0=1.5; %给出初值for i=1:1:30 %设定要计算的X值次数 while eeer y=x(i)+e*sin(y0); %不动点迭代格式 ee=abs(y-y0); %确定新误差值 y0=y; k=k+1; %失迭代次数记录 endm(i)=k; %记录X值迭代次数z(i)=y; %记录相应Y值ee=1; %复原值，进行下一个X值计算y0=1.5;k=0; end对于以上程序说明，由于题目给出值相对比较简单首先进行不动点迭代法计算，题目中要求给出30个0,值，语句x=0:pi/n:pi让x值为题目中需要计算的30个x值，用循
6、环变量i分别对于x值分别进行提取，使之嵌入到for循环中，对于迭代格式用while循环采用标准迭代格式进行，对于while循环的结果用向量m(1:30)、z(1:30)分别记录其迭代格式结果。牛顿迭代法程序及程序说明e=0.9; %给变量取值format longn=29; %0 pi中取值个数x=0:pi/n:pi; %给出向量值k=0; %记录迭代法数er=1e-8; %给出误差限ee=1; y0=1; %给出初值for i=1:1:30 %设定要计算的X值次数 while eeer y1=y0-(y0-e*sin(y0)-x(i)/(1-e*cos(y0) %牛顿迭代格式 ee=abs(
7、y0-y1); %确定新误差值 y0=y1; k=k+1; %失迭代次数记录 endm(i)=k; %记录X值迭代次数z(i)=y1; %记录相应Y值ee=1; %复原值，进行下一个X值计算y0=1;k=0;end程序与不动点迭代格式类似，仅牛顿迭代格式不一样，通过上面两程序运行，得到不同迭代格式算法下其收敛速度关系。6 数据及图形以下表格为不动点迭代格式数据记录X取值0.000000 0.108331 0.216662 0.324992 0.433323 0.541654 结果0.000000 0.660470 0.947341 1.144413 1.300692 1.433140 次数15
8、0482718139X取值0.649985 0.758315 0.866646 0.974977 1.083308 1.191639 结果1.549786 1.655118 1.751923 1.842065 1.926857 2.007264 次数6710131619X取值1.299969 1.408300 1.516631 1.624962 1.733292 1.841623 结果2.084019 2.157696 2.228751 2.297558 2.364428 2.429619 次数232731364249X取值1.949954 2.058285 2.166616 2.274946
9、 2.383277 2.491608 结果2.493357 2.555834 2.617219 2.677664 2.737305 2.796264 次数5665758699113X取值2.599939 2.708270 2.816600 2.924931 3.033262 3.141593 结果2.854654 2.912582 2.970147 3.027443 3.084562 3.141593 次数128142156168175178以下表格为牛顿失代法生成的数据X取值0.000000 0.108331 0.216662 0.324992 0.433323 0.541654 结果0.0
10、00000 0.660470 0.947341 1.144413 1.300692 1.433140 次数654555X取值0.649985 0.758315 0.866646 0.974977 1.083308 1.191639 结果1.549786 1.655118 1.751923 1.842065 1.926857 2.007264 次数666666X取值1.299969 1.408300 1.516631 1.624962 1.733292 1.841623 结果2.084019 2.157696 2.228751 2.297558 2.364428 2.429619 次数66566
11、6X取值1.949954 2.058285 2.166616 2.274946 2.383277 2.491608 结果2.493357 2.555834 2.617219 2.677664 2.737305 2.796264 次数778113716X取值2.599939 2.708270 2.816600 2.924931 3.033262 3.141593 结果2.854654 2.912582 2.970147 3.027443 3.084562 3.141593 次数814599937由上面结果可以看出来牛顿迭代法其收敛性远好于不动点迭代法，下面给出原函数图形，图中给x在0,上取值时函
12、数求解结果迭代次数的变化曲线图中红色曲线为牛顿迭代法迭代次数，蓝色线为不动点迭代法迭代次数随X值变化曲线，由图上可以看出来，在初值选定适当时，牛顿迭代法迭代次数少，算法稳定，不动点迭代法，迭代次数要多许多。7 实验总结通过实验验证了牛顿迭代法有二阶收敛性，只要初值选取适当，其收敛速度很快，对于函数变化，相比不动点迭代其也稳定得多，占用内存较小。在本实验中，对于编程开始忘记了对于每一次计算结果后，对于误差限、迭代次数以及初始点y0均要复原，开始没有察觉到此问题，程序没有调出。拉格郞日插值法逼近函数作者：罗波 学号：2008203020081 课题题目在区间5，5 上取 11 阶切比雪夫多项式的零点 k=0,1,10作插值结点，计算函数在结点处yk=f(xk)处值。构造插值多项式 其中取自变量点分别在点处L10(tk)，f(tk)并绘出两条曲线比较。 2 背景技术实践中常有这样的问题：由实验得到某一函数y = f (x)在一系列点x0, x1, xn处的值y0, yi, yn，其函数的解析表达式是未知的，需要构造一个简单函数P(x)作为y = f (x)的近似表达式；或者y = f (x)虽有解析式，但计算复杂，不便于使用，需要用一个比较简单且易于计算的函数P(x)去近似代替它；插值法是通过给定n+1个互异的插值节点，让你求一条n次代数
