[理学]MATLAB的方程组解法.ppt

Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 非线性方程的解法 求解非线性方程的MATALB函数 fzero函数调用格式： [x,fval,exitflag]=fzero(…)—exitflag返回函数fzero退出的函数值 [x,fval,exitflag,output]=fzero(…)—output返回优化的信息 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 非线性方程的解法 求解非线性方程的MATALB函数 举例：求解f(x)=3x2-exp(x)+1=0的根 命令窗口输入： x0=4; —赋初始值 fzero(3*x^2-exp(x)+1,x0) x1=[4 10]; —区间未包含根 fzero(‘3*x^2-exp(x)+1’,x1) x1=[3 5]; —区间包含根 fzero(3*x^2-exp(x)+1,x1) Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 非线性方程组的解法 非线性方程组的求解 求解更为复杂、繁琐，一般要采用迭代法进行求解 不动点迭代法 牛顿迭代法 Broyden迭代法 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 非线性方程组的解法 不动点迭代法 思想：有n个未知数和n个方程的非线性方程组，F(x)=0，其迭代形式为：x=g(x)，不动点的迭代公式 x(k+1)=g(xk) 如果得到的序列{x(k)}满足limx(k)=x*，则x*为函数g(xk)的不动点 举例：求解非线性方程组x1-0.6sinx1-0.3cosx2=0和x2-0.6cosx1+0.3sinx2=0 1.编制文件BDD.m 2.编制函数文件gg.m 3.命令窗口输入：[y,n]=BDD([0.5 0.5]) Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 非线性方程组的解法 牛顿迭代法 思想：有n个未知数和n个方程的非线性方程组，F(x)=0，其迭代形式为：x=g(x)，不动点的迭代公式 x(k+1)=x(k)-(F(x(k))-1F(x(k)) 牛顿迭代法就是需要求解出线性方程组F(x(k))△x(k)=- F(x(k))的△x(k)，如果△x(k)足够小，则得到x(k+1)即为非线性方程组的解 举例：求解非线性方程组x12-10x1+x22+8=0和x1x22+x1-10x2+8=0 1.编制文件NewDD.m 2.编制函数文件fx1.m、dfx1.m 3.命令窗口输入：[y,n]=NewDD([0 0]) Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 非线性方程组的解法 Broyden迭代法 思想：牛顿迭代法具有较好的收敛性，但每步都要计算Fg(x(k))的值，用较简单的矩阵A(k)年来代替牛顿迭代法中的Fg(x(k)) ，就可以较为方便的进行迭代计算 举例：求解非线性方程组x1-0.7sinx1-0.2cosx2=0和x2-0.7cosx1+0.2sinx2=0 1.编制文件Broy.m 2.编制函数文件bf.m 3.命令窗口输入：x0=[0.5 0.5]; [y,n]=Broy(x0) Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY 西南交通大学摩擦学研究所 MATLAB的方程(组)解法 第六讲 王文健 wwj527@163.com 《MATLAB数据处理与应用》2011-2012学年选修课 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 主要内容 线性方程(组)的解法 非线性方程(组)解法 MATLAB统计分析 Tribology Research Institute SOUTHWEST JIAOTONG UNIVERSITY * 简介 线性方程(组) 直接法—在没有舍入的情况下，通过有限步四则运算求得方程组的准确解 迭代法—先给出一个解的初始值，然后按一定的法则逐步求出解的各个更准确的近似值的方法 非线性方程(组) 迭代法—不动点迭代法、Newton迭代法、Broyden Broyden迭代法 Tribology Research Instit