迭代法求解开普列方程.docVIP

数值剖析实验报告 迭代法求解开普列方程 作者：罗波 学号：200820302008 课题题目 开普列方程是天体力学中的一个重要方程xysiny，此中，x是某个行 星的均匀近点角，y是它的偏近点角，ε是它的轨道的偏爱率。取ε=0.9，对区间[0,]中30个平分点的x计算出y的数值解。 课题背景 开普勒定律是开普勒发现的对于行星运动的定律。他于1609年在他第一版的 《新天文学》上发布了对于行星运动的两条定律，又于1618年，发现了第三条 定律。 开普勒很好运地能够获得，有名的丹麦天文学家第谷·布拉赫所察看与采集的，特别精准的天文资料。大概于1605年，依据布拉赫的行星地点资料，开普勒发现行星的挪动恪守三条相当简单的定律。 开普勒的定律赐予亚里士多德派与托勒密派在天文学与物理学上极大的挑战。他主张地球是不停地挪动的；行星轨道不是周转圆的，而是椭圆形的；行星公转的速度不等恒。这些论点，大大地摇动了当时的天文学与物理学。经过了几乎一世纪披星带月，忘餐废寝的研究，物理学家终于能够用物理理论解说此中的道理。牛顿利用他的第二定律和万有引力定律，在数学上严格地证明开普勒定律，也让人们认识此中的物理意义。 专业背景 迭代算法是用计算机解决问题的一种基本方法。它利用计算机运算速度快、 合适做重复性操作的特色，让计算机对一组指令（或必定步骤）进行重复履行， 在每次履行这组指令（或这些步骤）时，都从变量的原值推出它的一个新值。牛 顿迭代法（Newtonsmethod）又称为牛顿-拉夫逊方法（Newton-Raphsonmethod）， 它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多半 1 作者:罗波学号:200820302008 方程不存在求根公式，所以求精准根特别困难，甚至不行能，进而找寻方程的近 似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来找寻方程f(x)=0的根。牛顿迭代法是求方程根的重要方法之一，其最大长处是在方程f(x)=0的单根邻近拥有平方收敛，并且该法还能够用来求方程的重根、复根。此外该方法宽泛用于计算机编程中。对于不动眯迭代法，其迭代格式简单，经过对原方程变形结构x(xk),经过重复计算能够获得知足必定精度的数值解，但因为其对 于函数要求高，迭代次数多，对于牛顿鬼迭代法和不动点迭代法都有改良，获得收敛性较好的迭代格式。 问题剖析 对于超越方程xysiny没法求得其分析解，一定经过数值方法求其解， 对于不同x、ε值，会造成方程有不同的收敛速度，经过Plot绘图能够得出函数只有一种种类的拐点，对于牛顿迭代法其初值能够随意选择，经过绘图选择其初始点为1.0，分别结构不动点格式与牛顿迭代格式，分别进行迭代，获得结构，比较其迭代次数，而得出牛顿迭代法与不动点迭代法其性能上差异。 程序及说明 不动点迭代法程序以及程序说明 e=0.9;%给变量取值 formatlong%设定数据格式 n=29;%[0pi]中取值个数 x=0:pi/n:pi;%给出向量值 k=0;%记录迭代法数 er=1e-8;%给出偏差限 ee=1; y0=1.5;%给出初值 fori=1:1:30%设定要计算的X值次数 whileeeer y=x(i)+e*sin(y0);%不动点迭代格式 ee=abs(y-y0);%确立新偏差值 y0=y; k=k+1;%失迭代次数记录 end 2 数值剖析实验报告 m(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值，用循环变量i分别对于x值分别进行提取，使之嵌入到for循环中，对于迭代格式用while循环采纳标准迭代格式进行，对于while循环的结果用向 量m(1:30)、z(1:30)分别记录其迭代格式结果。 牛顿迭代法程序及程序说明 e=0.9;%给变量取值 formatlong n=29;%[0pi]中取值个数 x=0:pi/n:pi;%给出向量值 k=0;%记录迭代法数 er=1e-8;%给出偏差限 ee=1; y0=1;%给出初值 fori=1:1:30%设定要计算的X值次数 whileeeer y1=y0-(y0-e*sin(y0)-x(i))/(1-e*cos(y0))%牛顿迭代格式 ee=abs(y0-y1);%确立新偏差值 y0=y1; k=k+1;%失迭代次数记录 end m(i)=k;%记录X值迭代次数 z(i)=y1;%记录相应Y值 ee=1;%还原值，进行下一个X值计算 y0=1; k=0; end 程序与不动点迭代格式近似，仅