全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
matlab实现牛顿迭代法求解非线性方程组已知非线性方程组如下3*x1-cos(x2*x3)-1/2=0x12-81*(x2+0.1)2+sin(x3)+1.06=0exp(-x1*x2)+20*x3+(10*pi-3)/3=0求解要求精度达到0.00001首先建立函数fun储存方程组编程如下将fun.m保存到工作路径中:functionf=fun(x);%定义非线性方程组如下%变量x1x2x3%函数f1f2f3symsx1x2x3f1=3*x1-cos(x2*x3)-1/2;f2=x12-81*(x2+0.1)2+sin(x3)+1.06;f3=exp(-x1*x2)+20*x3+(10*pi-3)/3;f=f1f2f3;建立函数dfun用来求方程组的雅克比矩阵将dfun.m保存到工作路径中:functiondf=dfun(x);%用来求解方程组的雅克比矩阵储存在dfun中f=fun(x);df=diff(f,x1);diff(f,x2);diff(f,x3);df=conj(df);编程牛顿法求解非线性方程组将newton.m保存到工作路径中:functionx=newton(x0,eps,N);con=0;%其中x0为迭代初值eps为精度要求N为最大迭代步数con用来记录结果是否收敛fori=1:N;f=subs(fun(x0),x1x2x3,x0(1)x0(2)x0(3);df=subs(dfun(x0),x1x2x3,x0(1)x0(2)x0(3);x=x0-f/df;forj=1: length(x0);il(i,j)=x(j);endifnorm(x-x0)epsr=subs(F,findsym(F),x0);%迭代公式tol=norm(r-x0);%注意矩阵的误差求法,norm为矩阵的欧几里德范数n=n+1;x0=r;if(n100000)%迭代步数控制disp(迭代步数太多,可能不收敛!);return;endendx0=0 0 0;r,n,data=budong(x0); disp(不动点计算结果为) x1=1 1 1; x2=2 2 2; x,n,data=new_ton(x0); disp(初始值为0,牛顿法计算结果为:) x,n,data=new_ton(x1); disp(初始值为1,牛顿法计算结果为:) x,n,data=new_ton(x2); disp (初始值为2,牛顿法计算结果为:) budong.m functionr,n,data=budong(x0, tol) if nargin=-1 tol=1e-3: end x1=budong fun(x0); n=1; while(norm(x1-x0)tol)&(n500) x0=x1; x1=budong_fun(x0); n=n+1: data(:,n)=x1; end r=x1: new_ton.m function x,n,data=new_ton(x0, tol) if nargin=-1 tol=1e-8; end x1=x0-budong_fun(x0)/df1(x0); n=1; while (norm(x1-x0)tol) x0=x1; x1=x0-budong_fun(x0)/df1(x0); n=n+1; data(:,n)=x1; end x=x1; budong_fun.m function f=budong_fun(x) f(1)=3* x(1)-cos(x(2)*x(3)-1/2; f(2)=x(1)2-81*(x(2)+0.1)2+sin(x(3)+1.06; f(3)=exp(-x(1)*x(2)+20* x(3)+10* pi/3-1; f=f(1)*f(2)*f(3); df1.m function f=df1(x) f=3sin(x(2)*x(3)*x(3) sin(x(2)*x(3)*x(2) 2* x(1)-162*(x(2)+0.1)cos(x(3) exp(-x(1)*x(2)*(-x(2)exp(-x(1)*x(2)*(-x(1)20; 结果: 不动点计算结果为 r= 1.0e+012* NaN -Inf 5.6541 初始值为0,牛顿法计算结果为: x= 0.5000 -0.
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 携带手机用腕带项目成效分析报告
- 2023年高三物理备课组工作总结报告三篇
- 两轮机动车用轮胎市场分析及投资价值研究报告
- 闪光信号灯项目成效分析报告
- 流程意识与流程方法课件
- 市场营销方案策划书
- 2023年学校常规管理情况汇报材料
- 办公用淀粉浆糊产品相关项目实施方案
- 湖南省常德市石门县新铺乡凉井教学点2023-2024学年二年级下学期期中数学试题
- 第三单元解决问题的策略(单元测试)-2023-2024学年六年级下册数学苏教版
- 煤矿瓦斯检查工安全培训课件
- 宫外孕失血性休克应急预案演练记录
- 配电房数据采集解决方案
- 蜀绣文化国外现状研究报告
- 《药品销售渠道》课件
- 三维五轴联动激光加工设备全球前14强生产商排名及市场份额分析报告
- 医院充电宝服务为患者提供便利
- 预防火灾网盘
- 北师大六年级数学下册比例单元整理和复习
- 粤教版(新)一年级下第2单元 我来做小船7 我的小船
- 通勤车租赁投标方案(技术标)
评论
0/150
提交评论