




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 共享洗车合同标准文本
- 个人租房经典合同标准文本
- 农村道路工程合同标准文本
- 企业培训课程开发与教学方法
- 些属于劳务合同标准文本
- 声波测井全波形反演及井旁介质参数估计
- 农村房契合同标准文本
- 中介临时用工合同标准文本
- 2024年广东公需课《百县千镇万村高质量发展工程与城乡区域协调发展》试题及答案
- 碱激发高钙体系早期产物形成与硬化调控研究
- 2025新人教版七下英语单词默写表
- 化工行业保安工作计划
- 人教版一年级数学上册《6-10的认识和加减法》同步练习题及答案
- 二十案例示轮回
- 老年营养示范化病房创建方案
- 设备安全操作培训
- 西方文化概论(第二版)课件全套 曹顺庆 第0-6章 绪论 西方文化的渊源与流变、西方文学 -西方社会生活与习俗
- 某地区现代有轨电车施工方案
- GB/T 6974.3-2024起重机术语第3部分:塔式起重机
- 物业品质巡查管理制度
- 高中物理-《互感与自感》课件-新人教版选修3
评论
0/150
提交评论