利用Matlab求线性方程组的通解

来源 :硅谷 | 被引量 : 0次 | 上传用户:yideng
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  [摘要]讨论线性方程组解的个数及求线性方程组的通解问题是线性代数中的常见问题,介绍利用Matlab软件解决这两种问题的方法。
  [关键词]Matlab 线性方程组 基础解系 通解
  中图分类号:O1-0 文献标识码:A 文章编号:1671-7597(2008)1220156-01
  
  一、引言
  
  线性代数是数学中的一个重要分支,很多理论问题和实际问题都需要借助于线性代数的理论工具来分析解决。学习线性代数有两大难点:一是概念、理论抽象,二是计算量大。由于这两个难点,初学者往往很难掌握好线性代数的知识理论。若能掌握好用于解数学问题的Matlab软件,则能轻松快捷的解决很多线性代数问题。
  (一)线性方程组有关定理
  对线性方程组 ,其中
  
  定理1:(1)若 ,则线性方程组无解;
  (2)若 ,则线性方程组存在唯一解;
  (3)若 ,则线性方程组有无穷多个解。
  定理2:对线性方程组的导出组 ,若 ,则导出组的基础解系含有个解向量,(1)若为导出组的基础解系,则 就是导出组的全部解,也称为通解;(2)若 是的一个特解,则就是
   的通解。
  
  二、线性方程组有关问题的Matllab解法
  
  (一)有关Matlab命令
  讨论线性方程组解的问题的有关Matlab命令见下表
  
  (二)利用左除命令给出线性方程组的特解
  对线性方程组,在Matlab软件里常用矩阵的左除命令“A\b”求出该线性方程组的一个解。若A为方阵,则A\b和inv(A)*b基本一致;若A不为方阵,A\b命令使Matlab软件自动选择适当的方法来求解。
  若 无解,则命令“A\b”将给出一个最小二乘意义上的近似解,即使 的长度达到最小。若 有无穷解,则命令“A\b”将给出一个具有最多零元素的特解。若 有唯一解,则命令“A\b”给出这个唯一解。
  (三)讨论线性方程组是否有解的Matlab方法
  对线性方程组 ,在Matlab软件里先给系数矩阵A和常数列矩阵b赋值,然后利用命令“rank(A),rank(A,b)”求出系数矩阵和增广矩阵的秩,再根据定理(1)所介绍的数学理论进行判断。
  (四)求线性方程组的通解的matlab方法
  为求线性方程组的通解,需先判断方程组是否有解。若有解则用“ A\b”求出一个特解,再用命令“null(A)”求出导出组的一个基础解系得出其通解,然后利用定理(2)的结论给出的通解。
  
  三、Matlab应用实例
  
  解:(1)先判断线性方程组解的情况,输入并运行以下代码:a=[1 1 1 1 1;3 2 1 1 -3;0 1 2 2 6;5 4 3 3 -1]; b=[7 -2 23 12]' ;r1=rank(a),r2=rank([a b])运行结果为r1=r2=2<5,说明该线性方程组有无穷多个解。(2)输入以下代码求出线性方程组的一个特解x0=a\b运行结果为:x0=(3.1667,0,0,0,3.8333) '(3)输入代码求出导出组aX=0的一个基础解系xx=null(a)运行结果是:
  Xx=
  0.75300.0176 -0.0000
  -0.4167 -0.7464 -0.0000
  -0.30430.4533 -0.7071
  -0.30430.45330.7071
  0.2723 -0.1778 -0.0000
  (4)线性方程组的通解为
  
  其中 为任意常数。
  
  参考文献:
  [1]胡良剑、孙晓君,《Matlab数学实验》[M],北京:高等教育出版社,2006.
  [2]赵静、但琦,《数学建模与数学实验》(第2版)[M],北京:高等教育出版社,2003.
  [3]同济大学应用数学系,《高等数学》(第5版)[M],北京:高等教育出版社,2004.
  
  作者简介:
  简绍勇,男,江西新余人,新余高等专科学校数学与信息科学系助教。
  
  注:“本文中所涉及到的图表、注解、公式等内容请以PDF格式阅读原文。”
其他文献
[摘要]《Flash动画制作》课程是职校计算机专业的一个必修课程。Flash教学要求教师结合学生实际情况,采用不同的教学方法,引入更为灵活的实践教学来培养学生的创作能力,让学生在实践中提高技能,学到知识。   [关键词]Flash教学特点 Flash教学方法  中图分类号:G42 文献标识码:A 文章编号:1671-7597(2008)1220138-01    Flash是一款矢量图形编辑和动画
期刊
[摘要]系统采用GIS技术,使用Microsoft Visual Basic 2005编程语言,利用TCP/IP通信技术,实现对防区的实时监控、报警处理、视频显示功能,并实现部分应急联动功能,初步构建了一个基于GIS技术的安防平台。相比过去同类平台,提出了引入GIS技术,从而使监控设备的地学信息更加直观。  [关键词]平安城市 组件式GIS Active X控件  中图分类号:TP3 文献标识码:
期刊
[摘要]建立齿轮传动约束优化数学模型;结合实例借助于数学软件MATLAB的优化工具箱进行寻优计算,给出优化设计程序,得到满足实际需要的最优化参数;通过对结果进行比较,证明这种方法功能强大,优化效果好,耗时短,无疑将成为机械优化设计领域中的重要工具。  [关键词]齿轮传动优化设计数学模型  中图分类号:TP2文献标识码:A 文章编号:1671-7597(2008)1220125-01    一、引言
期刊
[摘要]从十进制与二进制的特点入手,介绍一种把十进制数转换为二进制的新方法。  [关键词]十进制 二进制转换 定位减权  中图分类号:TP3 文献标识码:A 文章编号:1671-7597(2008)1220109-01    随着计算机的不断发展,让更多的人了解了计算机是按二进制的方式工作的,但是日常生活中,人们使用的是十进制,这就让人们在了解计算机的时候,不得不考虑十进制与二进制之间的相互转换。
期刊
[摘要]论述ID3算法,根据中医脏腑辨证数据的特点,简化经典的ID3算法,描述一个针对中医脏腑辨证数据的数据挖掘过程模型,并在这个模型的指导下,对临床病人的实时诊断症状和证候之间的关系进行分析、预测,在保证正确诊治率的前提下提高运行的效率。  [关键词]数据挖掘 决策树 ID3算法 中医辨证  中图分类号:R4 文献标识码:A 文章编号:1671-7597(2008)1220134-01    一
期刊
[摘要]通过对石化管道储运分公司ERP系统项目建过程分析,提出在油气长输管线上实施ERP项目的策略。即:企业最高领导层的重视;加强基础网络稳定性、安全性建设;还要重视培训,不断提升系统功能和加强系统运营维护等因素。阐述管道储运分公司SAP项目实施对策,并分析系统上线后取得的实际效果。  [关键词]SAP 管道 ERP系统 实施策略  中图分类号:TE9 文献标识码:A 文章编号:1671-7597
期刊
[摘要]3DS MAX软件作为制作动画的主流软件之一,自1996年由Kinetix推出3DS MAX1.0版本后不断更新,如今已经升级到3DS MAX2009版本。随着软件的不断升级和完善,3DS MAX软件已经被广泛运用到各大领域,在游戏、建筑、动画等领域尤为突出,软件的完善和市场的需求推动了动画产业的发展,增加了人才需求量。随着相关教育的普及,涌现了一大批毕业生充斥着整个人才市场,然而这些毕业
期刊
[摘要]弱势群体是指在社会性资源分配过程中,具有经济利益的贫困性,生活质量的低层次和随力上的脆弱性的特殊社会群体。在弱势群体中,青年占有很大很大部分的比例。弱势青年群体是一处于极其困难随力极为脆弱,极需社会各界支持与帮助的群体。社会支持网的规模、关系构成、紧密程度、趋同性、异质性等方面的特征,对个体的生存状况有着或多或少的影响,基于此,以青年弱势群体为对象,采用问卷调查研究的形式,对该群体的生存状
期刊
[摘要]介绍安全性评价中的LEC法,并将该方法引入电力生产进行安全隐患的风险分析,对贵州电网的电力生产运行现场工作安全隐患进行了LEC评价。  [关键词]安全性评价 LEC 电力生产安全隐患评价  中图分类号:TM7文献标识码:A文章编号:1671-7597 (2008) 0120008-01    一、供电企业安全性评价    安全性评价即风险评估,指综合运用安全系统工程的方法对系统的安全性进行
期刊
[摘要]微型计算机已进入到我们生活的方方面面,随着计算机技术的不断发展,其文件系统也由FAT16,到FAT32,进而到安全性更好的NTFS,而病毒,木马程序及计算机使用过程中的不当操作等的影响,可能使计算机系统崩溃,利用GHOST软件快速恢复系统则是我们经常采用的手段。为了系统安全性高,我们采用NTFS文件系统,而GHOST软件生成的是FAT32格式的文件,在系统被破坏时,要实现系统的快速恢复,就
期刊