第10章 代数方程组的MATLAB求解_第1页
第10章 代数方程组的MATLAB求解_第2页
第10章 代数方程组的MATLAB求解_第3页
第10章 代数方程组的MATLAB求解_第4页
第10章 代数方程组的MATLAB求解_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

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

文档简介

1、第10章 代数方程组的MATLAB求解编者 Outline10.1 线性方程组的求解10.2 多项式方程组的准解析解法10.3 超越方程组的求解10.1 超越方程组的求解1.克莱姆(Cramer)法则及其MATLAB实现 对于恰定线性方程组 若其系数行列式 , 则上述线性方程组有唯一解其中 是把D 中第j 列元素 对应地换成常数项 ,而其余各列保持不变所得到的行列式,这就是克莱姆(Cramer)法则。MATLAB中没有提供实现克莱姆法则的函数,这需要我们自行编写相应的函数文件 2.消去法及其MATLAB实现 上三角形方程组的求解 称形如 的方程组为上三角形方程组,写成矩阵形式为 其中 称 U为

2、上三角矩阵。若 即则上三角方程组有唯一解,且可从上式的最后一个方程解出 代入倒数第二个方程可得 一般地,设已求得 ,则由上述方程组的第i个方程可得 上述求解方程组的过程称为回代过程。 Gauss消去法 Gauss消去法是求解线性方程组较为有效的方法,它主要包括两个操作:消元和迭代。顺序Gauss消去法是指按行原先的位置进行消元的Gauss消去法。下面介绍顺序Gauss消去法的一般步骤。为叙述方便,将线性方程组 写成增广矩阵的形式其中3.矩阵分解法及其MATLAB实现 求解线性方程组除了消去法,还有矩阵分解方法。所谓矩阵分解方法,就是将矩阵 A 分解成两个或多个简单矩阵的乘积,由于分解后的矩阵具

3、有某种特殊性,因此便于方程组的求解。4.迭代法及其MATLAB实现 迭代法是求解线性方程组的一个重要的使用方法,特别适用于求解系数矩阵为稀疏矩阵的大型线性方程组。 Jacobi迭代法 Gauss-Seidel迭代法 逐次超松弛迭代法5.线性方程组的MATLAB函数求解齐次线性方程组的求解非齐次线性方程组的求解10.2 多项式方程组的准解析解法求下面的二元多项式方程组。对于本题,我们也可以使用roots函数求解首先将 代入 则得到关于 的一元多项式试求解下面的含自变量导数形式的多项式型方程组。10.3 超越方程组的求解1.牛顿法及其MATLAB实现牛顿法求解非线性方程组 的根首先将多元向量函数 在点 处展开其中, 是 的Jacobi矩阵。解方程得 上式写成迭代格式为:牛顿法的一般执行流程(右图) 图 Newton法执行流程2.超越方程组的MATLAB函数求解 MATLAB优化工具箱提供的fsolve函数是专门用来求解非线性方程组的实数根的函数。例:利用fsolve函数求函数 在 上的所有零点,其中参数 【分析】对于这样一个比较复杂的函数,我们不好判断它的零点的范围,只有先画出图形观察。执行如下语句,得到函数图形如图所示。 图 函数图

温馨提示

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

评论

0/150

提交评论