Newton迭代法数值分析非线性方程组数值解法.docxVIP

信科0802 沈闻PAGE  PAGE 3 实验七、Newton迭代法 实验目的 了解非线性方程与方程组的数值解法；深入理解Newton迭代法，并用C++实现计算结果。 计算公式 将f(x)在点xn作Taylor展开: ——Taylor展开线性化 f(x)=0 近似于 f(xn)+ f′(xn)(x-xn)=0 (1) 从(1)解出x, 记为xn+1 ,则 它对应的迭代方程为 显然是f(x)=0的同解方程，故其迭代函数为 在 f(x)=0的根 x* 的某个邻域 内, 在 x* 的邻域R 内，对任意初值 ，应用公式（2）来解方程 三、程序代码 #includestdio.h #includemath.h #includeconio.h double f(double x); double ft(double x); main() { double X0,X1,tol; int n; int count=0; printf(请输入初始值：\n); scanf(%lf,X0); printf(请输入精度：\n); scanf(%lf,tol); printf(请输入最大循环次数\n); scanf(%d,n); X1=X0-f(X0)/ft(X0); while(fabs(X1-X0)=tol)/*在没有达到精度要求情况下就循环*/ { X0=X1; X1=X0-f(X0)/ft(X0); count++;/*循环次数计数*/ /*如果循环次数到用户输入的最大循环次数就退出，可以适当修改*/ if(count==n) { break; } } printf(方程的根为:%f\n,X1); getch(); } /************************* 函数功能：用于求f(x)的值 **************************/ double f(double x) { double y; y=x*x*x+2*x*x+10*x-20; return y; } /************************* 函数功能：用于求f(x)的导数值 **************************/ double ft(double x) { double y; y=3*x*x+4*x+10; return y; } 例题 用Newton迭代法求下面方程的一个正根,计算结果精确到7位小数.