迭代法及matlab实现方法一.ppt_第1页
迭代法及matlab实现方法一.ppt_第2页
迭代法及matlab实现方法一.ppt_第3页
迭代法及matlab实现方法一.ppt_第4页
迭代法及matlab实现方法一.ppt_第5页
已阅读5页，还剩3页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、四. 迭代法及其MATLAB程序,例:求方程 的一个正根.,构造迭代函数,迭代公式,基本思想: 由初始值, 代入迭代公式, 经过一定的迭代次数 k, 得到迭代序列xk, 以及相邻两次迭代的偏差piancha=|xk-xk-1|和它的相对误差xdpiancha=| xk-xk-1 |/|xk|的值.,当piancha1, xdpiancha0.5, k3时, 迭代序列发散, 重新输入新的迭代公式;,当piancha3时, 迭代序列收敛.,迭代法的MATLAB主程序1,输入的量: 初始值x0, 迭代次数 k 和迭代公式,运行后输出的量: 迭代序列 xk , 迭代k次得到的迭代值xk, 相邻两次迭代

2、的偏差 piancha=| xk-xk-1|和它的偏差的相对误差 xdpiancha=| xk-xk-1 |/|xk|的值.,迭代法的MATLAB主程序1,function k,piancha,xdpiancha,xk=diedai1(x0,k) % 输入的量-x0是初始值,k是迭代次数 x(1)=x0; for i=1:k x(i+1)=fun1(x(i);%程序中调用的fun1.m为函数y=(x) piancha= abs(x(i+1)-x(i); xdpiancha=piancha/( abs(x(i+1)+eps); i=i+1;xk=x(i);(i-1) piancha xdpian

3、cha xk end if (piancha 1),例:求方程 的一个正根.,构造迭代函数,三种迭代函数,迭代公式,迭代公式,迭代公式,3. 用迭代法求解方程f(x)=0在 (a,b)内的近似根的步骤,步骤1.,建立名为fun1.m的M文件如:,步骤2.,步骤3.,在matlab工作窗口输入程序:,function y1=fun1(x) y1=f(x);,将迭代法的主程序保存名为diedai1.m的M文件;, k,piancha,xdpiancha,xk=diedai1(x0,k),输入的量: 初始值x0, 迭代次数 k,运行后输出的量: 迭代k次得到的迭代值xk, 相邻两次迭代的偏差 piancha=| xk-xk-1|和偏差的相对误差 xdpiancha=| xk-xk-1 |/|xk|的值.,迭代法的几何解释:,迭代法的几何解释,3. 用迭代法求解方程f(x)=0在 (a,b)内的近似根的步骤,步骤1.,建立名为fun1.m的M文件如:,步骤2.,步骤3.,在matlab工作窗口输入程序:,function y1=fun1(x) y1=f(x);,将迭代法的主程序保存名为diedai1.m的M文件;, k,piancha,xdpiancha,x

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论