实验报告内容： 一：不动点迭代法解方程 二：牛顿插值法的MATLAB 实现 完成日期：2012 年6 月21 日星期四 数学实验报告一 日期：2012-6-21 班级 新能源11 姓名 王志红 学号 2110310021 实验 名称 用不动点迭代法解方程 【问题背景描述】 天文学中有一类著名的方程——开普勒方程 x=q*sinx+a （0> clear；clc; clf; >> f=inline('0.5*sin(x)+0.4'); >> g=inline('x') hold on >> fplot(g, [-2,2]) >> fplot(f, [-2,2]) >> hold off >> grid 输出结果如下所示： 所以，确定初值为x0=1 二：不断迭代 算法： 第一步：将 f(x0)赋值给 x1 第二步：确定x1-x0 的绝对值大小，若小于给定的误差值，则将 x1 当做方程的解，否则回到第一步 编写计算机程序： clear f=inline('0.5*sin(x)+0.4'); x0=1; x1=f(x0); k=1; while abs(x1-x0)>=1.0e-6 x0=x1; x1=f(x0); k=k+1; fprintf('k=%.0f,x0=%.9f,x1=%.9f\n',k,x0,x1) end 显示结果如下： k=2,x0=0.820735492,x1=0.765823700 k=3,x0=0.765823700,x1=0.746565483 k=4,x0=0.746565483,x1=0.739560873 k=5,x0=0.739560873,x1=0.736981783 k=6,x0=0.736981783,x1=0.736027993 k=7,x0=0.736027993,x1=0.735674699 k=8,x0=0.735674699,x1=0.735543758 k=9,x0=0.735543758,x1=0.735495216 k=10,x0=0.735495216,x1=0.735477220 k=11,x0=0.735477220,x1=0.735470548 k=12,x0=0.735470548,x1=0.735468074 k=13,x0=0.735468074,x1=0.735467157 >>。。。 ...
