毕业论文参考资料_非线性方程求根的迭代法.ppt

作业1： 求 ，要求精度为 . 要求（1）写出牛顿迭代格式，并分析其收敛速度（可仿照p97例4.3.1） （2）写出matlab实现程序（写清程序newtonqx(x0,tol,ftol,gxmax)中的各参数，另外自己写程序fnq,dfnq) 作业2 习题4.5 习题 4.8 弦截法 牛顿迭代法需要计算 ，有时候是件很麻烦的事 情，比如f(x)仅给出离散的形式，可用一阶差商 代替 ， 于是有弦截法迭代公式： 弦截法是超线性收敛的 例题 例 弦截法求方程 在区间（1，2）内的实根。 弦截法的matlab实现 function [k,piancha,xdpiancha,xk,yk]=gexian (x01,x02,tol,ftol,gxmax) x(1)=x01;x(2)=x02; for i=2: gxmax u(i)= fnq(x(i))*(x(i)-x(i-1)); v(i)= fnq(x(i))-fnq(x(i-1)); x(i+1)=x(i)- u(i)/( v(i)); piancha=abs(x(i+1)-x(i)); xdpiancha= piancha/( abs(x(i+1))+eps); i=i+1; xk= x(i); yk=fnq(x(i)); [(i-2) piancha xdpiancha xk yk] if (abs(yk)<ftol)&(( piancha <tol)|(xdpiancha< tol)) k=i-2; xk=x(i);yk=fnq(x(i)); [(i-2) piancha xdpiancha xk yk]； return; end end if i>gxmax disp('请注意：迭代次数超过给定的最大值gxmax.') k=i-2; xk=x(i);yk=fnq(x(i)); return; end * * * 例题 若取迭代函数 ， 不满足压缩映像原理，故不能肯定 收敛到方程的根。 简单迭代收敛情况的几何解释 是否取到合适的初值，是否构造合适的迭代格式，对于是否收敛是关键的。 对于初值，实际操作时，可以先画出函数图形，然后，观察根大概在什么地方。 对于迭代格式，可以对 求导，看看 是否小于1 迭代法收敛的阶 定义 设序列 收敛到 ，若有实数 和非零常数C，使得 其中， ，则称该序列是p 阶收敛的， 迭代法收敛的阶 当p=1时，称为线性收敛； 当p>1时，称为超线性收敛； 当p=2时，称为平方收敛或二次收敛。 误差估计 若 满足定理条件，则 下面定理给出判别迭代收敛阶的一个方法 定理： 记 是 的根， ，设 在 附近连续，若对 ，有 则基本迭代法 是P阶连续的。 基本迭代法的matlab实现 function [k,piancha,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)); i=i+1;xk=x(i);[(i-1) piancha xk] end Matlab中与或非，分别是：& | ~ 与或非 if (piancha >1)&(k>3) disp('请用户注意：此迭代序列发散，请重新输入新的迭代公式') return; end if (piancha < 0.001)&(k>3) disp('祝贺您！此迭代序列收敛，且收敛速度较快') return; end p=[(i-1) piancha xk]'; 关于程序里面的fun1,可以如下类似定义 function y1=fun1(x) y1=(10-x^2)/2; 作业： 1. 编程求方程