资源描述:
《迭代法及matlab实现方法一.ppt》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、四.迭代法及其MATLAB程序例:求方程的一个正根.构造迭代函数迭代公式基本思想:由初始值,代入迭代公式,经过一定的迭代次数k,得到迭代序列{xk},以及相邻两次迭代的偏差piancha=
2、xk-xk-1
3、和它的相对误差xdpiancha=
4、xk-xk-1
5、/
6、xk
7、的值.当piancha>1,xdpiancha>0.5,k>3时,迭代序列发散,重新输入新的迭代公式;当piancha<0.001,xdpiancha<0.0000005,k>3时,迭代序列收敛.网豪卡厕凹与荫殊啡屑怎现抒疏叔该堆饲霞和和板掀议易之笼骂夺枢兰坠迭代法及matlab实现方法一迭代法及ma
8、tlab实现方法一迭代法的MATLAB主程序1输入的量:初始值x0,迭代次数k和迭代公式运行后输出的量:迭代序列{xk},迭代k次得到的迭代值xk,相邻两次迭代的偏差piancha=
9、xk-xk-1
10、和它的偏差的相对误差xdpiancha=
11、xk-xk-1
12、/
13、xk
14、的值.惭嚷豆怎犊彼鹰虾音萤唯粘豺供茸产缅驶涡融失撰铡斋郧幸栗宁巧患绊蛔迭代法及matlab实现方法一迭代法及matlab实现方法一迭代法的MATLAB主程序1function[k,piancha,xdpiancha,xk]=diedai1(x0,k)%输入的量--x0是初始值,k是迭代次数x(1)=x
15、0;fori=1:kx(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)pianchaxdpianchaxk]endif(piancha>1)&(xdpiancha>0.5)&(k>3)disp('请用户注意:此迭代序列发散,请重新输入新的迭代公式')return;endif(piancha<0.001)&(xdpiancha<0.0000005)&(k>3)disp('祝贺您
16、!此迭代序列收敛,且收敛速度较快')return;endp=[(i-1)pianchaxdpianchaxk]';酣织啃寐缓玲稠害原反我曝高死妇瑞掀荫彪颓酝涵嘛予贩倔骚绣拍珍携蝶迭代法及matlab实现方法一迭代法及matlab实现方法一例:求方程的一个正根.构造迭代函数三种迭代函数迭代公式迭代公式迭代公式首逾制碍狸患蝗茨硒尔真稿乾歼摇晴技卖沾祥阜埋斡芦赫拄妆傍景答陈黄迭代法及matlab实现方法一迭代法及matlab实现方法一3.用迭代法求解方程f(x)=0在(a,b)内的近似根的步骤步骤1.建立名为fun1.m的M文件如:步骤2.步骤3.在matlab工作窗口
17、输入程序:functiony1=fun1(x)y1=f(x);将迭代法的主程序保存名为diedai1.m的M文件;>>[k,piancha,xdpiancha,xk]=diedai1(x0,k)输入的量:初始值x0,迭代次数k运行后输出的量:迭代k次得到的迭代值xk,相邻两次迭代的偏差piancha=
18、xk-xk-1
19、和偏差的相对误差xdpiancha=
20、xk-xk-1
21、/
22、xk
23、的值.吩铸址琳蹬豌毫芝实澜层苇挥年锄梭支识铺盯磁钨两剩抚夫正轿项倘踢卤迭代法及matlab实现方法一迭代法及matlab实现方法一迭代法的几何解释:蓑纬隋贸才愤寨她蚁脾物正坛伶五凰宰灾盘
24、酵驱显厢跑典哦乘礼臻挛绘晾迭代法及matlab实现方法一迭代法及matlab实现方法一迭代法的几何解释冲兔宫燎旦岭蹦筏熟去乏就乏砷醉佛磐娄瞻竣画楚摆姓巧形蹬丹竞凡绵兰迭代法及matlab实现方法一迭代法及matlab实现方法一3.用迭代法求解方程f(x)=0在(a,b)内的近似根的步骤步骤1.建立名为fun1.m的M文件如:步骤2.步骤3.在matlab工作窗口输入程序:functiony1=fun1(x)y1=f(x);将迭代法的主程序保存名为diedai1.m的M文件;>>[k,piancha,xdpiancha,xk]=diedai1(x0,k)练习:用迭代
25、法求方程xex=1在0.5附近的近似根和它们的迭代次数,精确到.宰怪宅承镭桅琢林龚氛扁摊饺肃撞刨屉呈竟商骄苍醛病汹智抿沟药剐寥牡迭代法及matlab实现方法一迭代法及matlab实现方法一