《计算方法课件:8非线性方程组》由会员分享,可在线阅读,更多相关《计算方法课件:8非线性方程组(8页珍藏版)》请在金锄头文库上搜索。
1、计计算算方方法法8 8非线性方程组的迭代法非线性方程组的迭代法计计算算方方法法直接推广Newton迭代为:单个方程的牛顿迭代形式:雅可比矩阵(Jacobi) 关键关键计计算算方方法法 (1) (1) 选定初值选定初值x0 ,计算计算f (x0) , f (x0) 计算步骤计算步骤:单个方程单个方程?(2) (2) 按公式按公式 迭代迭代 得新的近似值得新的近似值xk+1 (3) (3) 对于给定的允许精度对于给定的允许精度 ,如果如果 则终止迭代,取则终止迭代,取 ; ;否则否则k=k+1,再转再转 步骤步骤(2)(2)计算计算允许精度允许精度最大迭代最大迭代次数次数迭代信息迭代信息方程组呢方
2、程组呢?(1) (1) 选定初值选定初值X X0 ,计算计算F(X0) , J(X0) ( (2) 2) 按公式计算按公式计算X Xk+1k+1计计算算方方法法 (2) (2) 按公式按公式 迭代迭代 得新的近似值得新的近似值Xk+1 (3) (3) 对于给定的允许精度对于给定的允许精度 ,如果如果 则终止迭代,取则终止迭代,取 ; ;否则否则k=k+1,再转再转 步骤步骤(2)(2)计算计算(1) (1) 选定初值选定初值X X0 ,计算计算F(X0) , J(X0) 非线性方程组解法:非线性方程组解法:计计算算方方法法计计算算方方法法Matlab程序function P,iter,err=
3、newdim(FX,JF,P,delta,epsilon,max1)B=feval(FX,P);for k=1:max1 J=feval(JF,P); deltaX= uptrbk(J,-B); X=P+deltaX; err=norm(deltaX); relerr=err/(norm(X)+eps); P=X; B=feval(FX,X); iter=k; if (errdelta)|(relerrdelta)|(abs(B)epsilon) break endend计计算算方方法法调用方式:P=2.00 0.25;delta=0.000001;epsilon=0.000001;max1=100;P,iter,err=newdim(FX,JF,P,delta,epsilon,max1)B=feval(FX,P)err=norm(B)计计算算方方法法作业:H71、H72