优·线性代数实验 2.docVIP

一、 高斯消元法解方程实验1．Matlab简介 Matlab 是美国 The MathWorks公司出品的计算机科学计算软件，从1984年推出以来，受到广泛的推崇，在很多领域里，Matlab 已成为科技人员首选的计算机数学语言．它的语言简洁，功能强大，几乎涵盖了所有的数学计算内容，人机交互性能好，表达方式符合科技人员的思维习惯和书写习惯，短简的语句，便能完成许多复杂的计算．MATLAB是“矩阵实验室”（Matrix Laboratory）的缩写，它是一种以矩阵运算为基础的交互式程序语言，因此特别适于线性代数．线性代数是一门理论比较抽象、计算强度很大的数学学科，具有广泛的应用．在传统教材给出的线性代数的计算方法，如用手工计算，只能解决一些低阶、较少的问题，而在实际中出现的大量的线性问题，都是高阶的和变量很多的，用Matlab语言辅助线性代数的教学，近几年来已成为的模式．本将对Matlab语言作简单的介绍，章都安排用Matlab语言的实验，解决相应章节的计算问题．Matlab是科技工作者得力的科学计算工具，读者可查阅有关的书籍． 2．矩阵的表示Matlab程序后，会出现一个命令窗，Matlab语后键入 如要Matlab中输入一个矩阵A ＝在Matlab提示符〉〉后面键入 A=[1 2 3 4 5 6 7 8 9] 按回车A = 1 2 3 4 5 6 7 8 9 也可以键入： A=[1，2，3;4，5，6;7，8，9] 或A=[1 2 3;4 5 6;7 8 9]在程序中就代表所输入的矩阵。 3．线性方程组的高斯消元法 线性方程组的高斯消元法，等价于对相应的矩阵做行初等变换，化为型矩阵， 在Matlab语言中，调用一个函数rref()，可将矩阵化为型矩阵．化为阶梯形矩阵， 先键入矩阵A=[1 2 -1 4 7;1 2 3 4 3;-1 3 -2 4 9;6 -2 7 1 2] 屏幕显示： A = 1 2 -1 4 7 1 2 3 4 3 -1 3 -2 4 9 6 -2 7 1 2 再调用函数，将化为行最简型矩阵： rref(A) 结果显示为： ans = 1.0000 0 0 0 13.3333 0 1.0000 0 0 27.6667 0 0 1．0000 0 -1.0000 0 0 0 1.0000 -15.6667 例1．5　求解方程 键入方程矩阵 A = 1 2 3 9 2 -1 1 8 3 0 -1 3 （2）化为行阶梯矩阵 rref(A)ans = 1 0 0 2 0 1 0 -1 0 0 1 3 所以 x=2 ，y=-1， z=3例1：求解方程 A=[1 -2 1 3;2 -3 -1 7;5 -8 -1 20]; rref(A)ans = 1 0 -5 0 0 1 -3 0 0 0 0 1 对应的方程为显然方程无解． 例1．　求解方程 A=[3 4 -3 -6;-1 -1 2 4;1 2 1 2]; rref(A)ans = 1 0 -5 -10 0 1 3 6 0 0 0 0 对应的方程为： Z取任意值，得到的x，y，z都是方程的解，所以方程有无穷多个解． 2．用高斯消元法解下列方程组 （1） （2） （3） 二、 矩阵运算实验Matlab语言的基本计算对象是向量和矩阵，而把数看作一维向量． 向量的表示Matlab中，向量的赋值可用下列分式：