




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2022年2月韶关市直机关遴选公务员面试真题带详解
- 2022年11月三明市直机关遴选公务员面试真题带详解
- 17护资试题及答案
- 2025翡翠首饰买卖合同模板
- 2025广告外包合同书范本
- 4-Bromotoluene-d7-生命科学试剂-MCE
- 2025合作合同模板合作伙伴协议范本
- 2025建筑材料供应调拨合同
- 2025房屋租赁合同范本正式版
- 2025购房租赁合同范本模板
- 复旦大学面试题及答案
- AI时代小学数学智慧课堂的构建与实践探索
- T-CECS 10400-2024 固废基胶凝材料
- 八年级语文上册第四单元整体公开课一等奖创新教学设计
- 智慧小区建设方案
- 2025年电潜螺杆泵项目可行性研究报告
- 肝门部胆管癌诊断和治疗指南(2025版)解读
- 新版统编版一年级道德与法治下册全册教案(完整版)教学设计含教学反思
- 2024年不动产登记代理人《地籍调查》考试题库大全(含真题、典型题)
- 财务服务协议书
- YC/Z 623-2024烟草商业企业卷烟物流应急作业指南
评论
0/150
提交评论