




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、、实验目的:1、学会判断非线性方程的根的存在性、根的分布范围;2、掌握迭代法求非线性方程的根。学会比较不同方法之间的优缺点;、实验题目:求非线性方程的根,精确到,给定方程为:(i): .、实验要求:A. 用你自己设计的一种线性收敛迭代法求方程(i)的根,然后再用斯蒂芬加速迭代计算。B. 用牛顿法求方程(i)的根,输出迭代初值,各次迭代值及迭代次数,并与(1)的结果比较。运用不动点迭代法:当时,将.化为x=(x) 因为x0,1,(x)0,1且(x)=在0,1内小于1,则(x)在0,1上存在唯一的不动点x*,使得x*=(x*);求f(x)的零点就等价于去(x)的不动点,将初始点=0.1代入(1)中
2、:得:=()反复迭代计算得:=(),(k=0,1,) 斯蒂芬森加速迭代法:埃特金方法不管原序列是怎样产生的,对进行加速计算,得到序列,如果把埃特金加速技巧与不动点迭代结合,则可得到如下的迭代法:,(*)称为斯蒂芬森(Steffensen)迭代法可将它写成另一种不动点迭代(*)其中(7.3.5)Newton迭代格式为: =-()/() k=0,1,2,实验过程源代码及运行结果Matlab程序:通过图像判断有无解:Clearformat longfor t=1:1000 x(t)=0.01*t; y1(t)=x(t)2-3*x(t)+2-exp(x(t);endy2=0;plot(x,y1,r*,
3、x,y2,bo)运行:不动点迭代法:迭代公式:x=(x) MATLAB程序:clearformat longk=1;x0(k)=0.5;x0(k+1)=(x0(k)2+2-exp(x0(k)/3;while abs(x0(k+1)-x0(k)0. k=k+1; x0(k+1)=(x0(k)2+2-exp(x0(k)/3;endx0k斯蒂芬森(Steffensen)迭代法:MATLAB程序:clearformat longx0=0.2;k=1;fi=inline(x.2+2-exp(x)/3,x);yk=fi(x0);zk=fi(yk);xk=x0-(yk-x0)2)/(zk-2*yk+x0);
4、ax(k)=xk;while abs(xk-x0)0. x0=ax(k); yk=fi(x0); zk=fi(yk); xk=x0-(yk-x0).2)/(zk-2*yk+x0); k=k+1; ax(k)=xk;endaxk运行后:牛顿迭代法:Matlab编程:clearformat longx0=0.5;k=1;x=x0-(x02-3*x0+2-exp(x0)/(2*x0-3-exp(x0);while abs(x-x0)0. x0=x; x=x0-(x02-3*x0+2-exp(x0)/(2*x0-3-exp(x0); k=k+1;endxk通过图像判断有无解:format longfo
5、r t=1:1000 x(t)=0.01*t; y1(t)=x(t)3+2*x(t)2+10*x(t)-20;endy2=0;plot(x,y1,r*,x,y2,bo)不动点迭代法:format longk=1;x0(k)=1.7;x0(k+1)=(20-x0(k)3-2*x0(k)2)/10;while abs(x0(k+1)-x0(k)0. x0=ax(k); yk=fi(x0); zk=fi(yk); xk=x0-(yk-x0).2)/(zk-2*yk+x0); k=k+1; ax(k)=xk;endaxk结果:ax = Columns 1 through 3 1.1475 1.7977 1.0470 Columns 4 through 5 1.0797 1.1373k = 5牛顿迭代法:format longx0=0.5;k=1;x=x0-(x03+2*x02+10*x0-20)/(3*x02+4*x0+10);while abs(x
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 决策权采购合同
- 车辆代步服务协议
- 资产管理合作协议内容与签订须知事项
- 光伏电站投资建设合作协议
- 农村劳务服务输出协议
- 配备司机汽车租赁合同
- 工厂车间承包生产合同书(5篇)
- 国际劳务合同格式(15篇)
- 海外院校录取与签证申请代理服务合同
- 医疗器械股权认购及合作开发合同范本
- 2025届浙江省杭州市建兰中学八年级英语第二学期期末综合测试试题含答案
- 小麦收割合同协议书模板
- 广东省惠州市惠阳区2023-2024学年一年级下学期语文期末随堂练习试卷(含答案)
- 起重吊装作业安全与操作规范培训
- 2025辽宁沈阳副食集团所属企业招聘25人笔试参考题库附带答案详解
- 2023电气装置安装工程 旋转电机施工及验收规范
- 施工项目资料管理试题及答案
- 2025年安全知识竞赛题库及答案(共150题)
- 第六单元 年、月、日 单元测试(含答案)2024-2025学年三年级下册数学人教版
- 国家开放大学2025年《机电控制工程基础》形考任务1-4答案
- 合资公司的组织结构与运营规范
评论
0/150
提交评论