牛顿迭代法解线性方程matlab程序,牛顿迭代法matlab程序(解线性方程组)

《牛顿迭代法matlab程序(解线性方程组)》由会员分享，可在线阅读，更多相关《牛顿迭代法matlab程序(解线性方程组)(4页珍藏版)》请在金锄头文库上搜索。

1、牛顿迭代法matlab程序(解线性方程组)作者：佚名来源：转载发布时间： 2009-3-7 16:55:53 减小字体增大字体1功能本程序采用牛顿法，求实系数高次代数方程f(x)=a0xn+a1xn-1+an-1x+an=(an )(1)的在初始值 x0附近的一个根。2.使用说明(1)函数语句Y=NEWTON_1(A,N,X0,NN, EPS1) 调用 M 文件 newton_1.m。(2)参数说明A n+1 元素的一维实数组，输入参数，按升幂存放方程系数。N整变量，输入参数，方程阶数。X0 实变量，输入参数，初始迭代值。NN整变量，输入参数，允许的最大迭代次数。EPS1实变量，输入参数，控制。

2、根的精度。3方法简介解非线性方程 f(x)=0 的牛顿法是把非线性方程线性化的一种近似方法。把f(x) 在 x0点附近展开成泰勒级数f(x)=f(x0)+(x x0)f(x0)+(x x0)2+取其线性部分，作为非线性方程f(x)=0 的近似方程，则有f(x0)+f(x0)(xx0)=0设 f(x0)0 则其解为x1=x0f(x0)/f(x0) 再把 f(x) 在 x1附近展开成泰勒级数，也取其线性部分作f(x)=0 的近似方程。若 f(x1)0 ，则得x2=x1f(x1)/f(x1) 这样，得到牛顿法的一个迭代序列xn+1=xnf(xn)/f(xn) 4 newton_1.m程 序funct。

3、ion y=newton_1(a,n,x0,nn,eps1 ) x(1)=x0；b=1；i=1；while(abs(b)eps1*x(i) i=i+1 ；x(i)=x(i-1 )-n_f(a,n,x(i-1)/n_df(a,n,x(i-1) ；b=x(i)-x(i-1) ；if(inn)error( nn is full )；return；end end y=x(i) ；i 程序中调用的 n_f.m 和 n_df.m 文件如下：function y=n_df(a,n,x)%方程一阶导数的函数y=0.0;for i=1:n y=y+a(i)*(n +1-i)*x(n-i);end functio。

4、n y=n_df(a,n,x) y=0.0；for i=1:n y=y+a(i)*(n+-i)*x?(n-i)；end 5程序附注(1)程序中调用 n_f.m 和 n_df.m 文件。 n_f.m 是待求根的实数代数方程的函数,n_d f.m 是方程一阶导数的函数。由使用者自己编写。(2)牛顿迭代法的收敛速度：如果f(x)在零点附近存在连续的二阶微商， 是 f(x)的一个重零点，且初始值x0充分接近于 ，那么牛顿迭代是收敛的，其收敛速度是二阶的，即平方收敛速度。6例题用牛顿法求下面方程的根f(x)=x3+2x2+10x-20 y =y+a(i)*(n+1-i)*x(n-i);7运行结果a=1,。

5、2,10,-20 ；n=3 ；x0=1；nn=1000；eps1=1e-8；y=newton_1(a,n,x0,nn,eps1) y= 1.368808107821373e+000 i= 6 function fp = newton_interpolation(x,y,p) % Script for Newtons Interpolation. % Muhammad Rafiullah Arain % Mathematics a(1) = y(1); for k = 1 : n - 1 d(k, 1) = (y(k+1) - y(k)/(x(k+1) - x(k); end for j = 2 : n - 1 for k = 1 : n - j d(k, j) = (d(k+1, j - 1) - d(k, j - 1)/(x(k+j) - x(k); end end d for j = 2 : n a(j) = d(1, j-1); end Df(1) = 1; c(1) = a(1); for j = 2 : n Df(j)=(p - x(j-1) .* Df(j-1); c(j) = a(j) .* Df(j); end fp=sum(c); 。