牛顿迭代法求解非线性方程组的代码

上传人:夏** 文档编号:478854976 上传时间:2022-11-30 格式:DOC 页数:3 大小:42KB
返回 下载 相关 举报
牛顿迭代法求解非线性方程组的代码_第1页
第1页 / 共3页
牛顿迭代法求解非线性方程组的代码_第2页
第2页 / 共3页
牛顿迭代法求解非线性方程组的代码_第3页
第3页 / 共3页
亲,该文档总共3页,全部预览完了,如果喜欢就下载吧!
资源描述

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

1、牛顿迭代法求解非线性方程组非线性方程组如下:给定初值,要求求解精度达到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.编程牛顿迭代法解非线性方程

2、组,将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

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 建筑/环境 > 建筑资料

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号