《牛顿下山法实验报告》由会员分享，可在线阅读，更多相关《牛顿下山法实验报告（5页珍藏版）》请在金锄头文库上搜索。
1、实验一：牛顿下山法一：上机题目使牛顿下山法求解方程的解，已知方程为：x3 -x - 1 = 0 二：牛顿下山法运行步骤1 对原方程求导，得到导函数。2已知迭代公式，输入xO,求解x1.3对x1,x0，代入原函数判断其大小关系，根据下山条件继续求解4 判断所求是否满足精确度。5 最终输入 x0,x1。三：程序流程图四：源程序代码 #include stdio.h float f(float x) float f;f=x*x*x-x-1;return f;float g(float x) float g;g=3*x*x-1; return g;main()float a,b,c,d;c=1;pri
2、ntf( ”请输入初值:”)；scanf(%f,&a);for(int i=0；i=f(a)*f(a)c=0.5*c；a=b；elsea=b；printf(%f,b)；五：运行结果实验二：高斯赛德尔算法一：上机题目 使用高斯法解线性方程组：高斯赛德尔算法定义由雅可比迭代法可知，在计算Xk+1的过程中，采用的都是上一迭代步的结果Xk。 考察其计算过程，显然在计算新分量Xk+1时，已经计算得到了新的分量，iA+l A+2 A+1。有理由认为新计算出来的分量可能比上次迭代得到的分量有所改善。希望充分利用新计算出来的分量以提高迭代解法的效率，这就是高斯赛 德尔迭代法（简称 GS 迭代法） 对（64）式
3、进行改变可以得到 GS 迭代法的分量形式GS 迭代法的分量形式亦可表示为 也可写成矩阵形式。方程组的系数A在(58)式基础上还可进一步分解，若将 A继续分解为一个下三角阵Al和一个上三角阵Au。0 0 0三：程序代码#include stdio.hfloat max(float x,float y,float z)float max;max=x;if(ymax)max=y;if(zmax)max=z;return max;void main()float x3,a3,c3;int i;printf(请输入初值：)；for(i=0;i3;i+)scanf(%f,&xi)；dofor(int k=0;k3;k+)ak=xk;x0=-0.4*x1-0.2*x2-2.4;x1=0.25*x0-0.5*x2+5;x2=-0.2*x0+0.3*x1+0.3;for(int j=0;j0.0001);printf(%fn,x0);printf(%fn,x1);printf(%fn,x2);四：运行结果
