第3讲+用MATLAB软件解线性方程组.ppt_第1页
第3讲+用MATLAB软件解线性方程组.ppt_第2页
第3讲+用MATLAB软件解线性方程组.ppt_第3页
第3讲+用MATLAB软件解线性方程组.ppt_第4页
第3讲+用MATLAB软件解线性方程组.ppt_第5页
已阅读5页,还剩4页未读 继续免费阅读

下载本文档

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

文档简介

1、第3讲 用MATLAB软件求逆矩阵和解线性方程组,制作: 江西科技师范学院 万重杰,3.1 矩阵函数,MATLAB提供了许多矩阵函数. 正是因为拥有了为数众多的、完善的矩阵函数,才使得MATLAB具有了强大的功能。,det 计算矩阵的行列式的值 inv 求矩阵的逆阵 rank 求矩阵的秩 V D=eig(A) 求矩阵A的特征值和特征向量 poly 求矩阵的特征多项式 rref 用初等变换将矩阵化成行阶梯形 null(A,r) 给出齐次线性方程组Ax=0 的基础解系 fliplr 矩阵左右翻转 flipud 矩阵上下翻转 trace 求矩阵的迹 diag 取得矩阵对角线元素,下面是几个常用的矩阵

2、函数:,在命令窗口运行帮助命令 :help elmat ,可以列举出大量的矩阵函数.,3.2 矩阵函数的应用,例31 设矩阵,解: A=3 -4 0; -1 5 2; 4 1 -6 det (A) %求矩阵的行列式的值 rank (A) %求矩阵的秩 inv (A) %求逆矩阵,求A的行列式、秩和逆矩阵。,3.3 解线性方程组,求线性方程组Ax=B的解,其中:,解法1 利用矩阵除法: X=AB,解法2 利用求逆矩阵函数 inv:X1=inv(A)*B,比较:解法1比解法2更简便, 解法1 的算法优于解法2 , 解法1可用于一般矩阵,而解法2只能用于非奇异的方阵 因此,只需运用解法1 .,例32

3、,3.3.1 求线性方程组的唯一解,3.3.2 求线性方程组的通解,求线性方程组Ax=B的通解。,例33 设,在命令窗口输入以下命令: (注意:这里给出的 A不 是方阵) A=1 1 -1 -1;2 -5 3 2;7 -7 3 1; B=5; -4; 7; format ratx1=AB %求得非齐次方程组Ax=B的一个特解x1Y=null(A,r) %求得齐次方程组Ax=0 的基础解系Y,则方程组Ax=B的通解为: x=x1+k1*Y(:,1)+k2*Y(:,2),x1 = 3 2 0 0,Y = 2/7 3/7 5/7 4/7 1 0 0 1,输出结果:,解法1: 利用除法 和 null 函数,在命令窗口输入以下命令:,format ratA=1 1 -1 -1;2 -5 3 2;7 -7 3 1; B=5; -4; 7;%用初等行变换将增广矩阵 A B 化成最简行阶梯形T T=rref(A B),于是可得方程组Ax=B的通解为:,解法2:利用 rref 函数,T = 1 0 -2/7 -3/7 3 0 1 -5/7 -4/7 2 0 0 0 0 0

温馨提示

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

评论

0/150

提交评论