资源预览内容
第1页 / 共3页
第2页 / 共3页
第3页 / 共3页
亲,该文档总共3页全部预览完了,如果喜欢就下载吧!
资源描述
牛顿迭代法求解非线性方程组非线性方程组如下:给定初值,要求求解精度达到0.000011.首先建立函数,方程编程如下,将F.m保存到工作路径中:function f=F(x) f(1)=x(1)2-10*x(1)+x(2)2+8; f(2)=x(1)*x(2)2+x(1)-10*x(2)+8; f=f(1),f(2) ;2.建立函数,用于求方程的jacobi矩阵,将DF.m保存到工作路径中:function df=DF(x) df=2*x(1)-10,2*x(2);x(2)2+1,2*x(1)*x(2)-10; %jacobi矩阵是一阶偏导数以一定方式排列成的矩阵。3.编程牛顿迭代法解非线性方程组,将newton.m保存在工作路径中:clear,clc; x=0,0; f=F(x); df=DF(x); fprintf(%d %.7f %.7fn,0,x(1),x(2);N=4; for i=1:N y=dff; x=x-y; f=F(x); df=DF(x); fprintf(%d %.7f %.7fn,i,x(1),x(2); if norm(y)0.0000001 break; else end endezplot(x2-10*x+y2+8,-6,6,-6,6);hold onezplot(x*y2+x-10*y+8,-6,6,-6,6);运行结果如下:0 0.0000000 0.00000001 0.8000000 0.88000002 0.9917872 0.99171173 0.9999752 0.99996854 1.0000000 1.0000000
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号