MATLAB在线性代数课程中应用.docVIP

MATLAB在线性代数课程中应用 　　【摘要】 结合线性代数课程特点，引入计算软件MATLAB辅助教学，探讨MATLAB在线性代数中的矩阵运算、行列式计算、向量组的线性相关性、线性方程组求解以及特征值和特征向量等若干问题方面的应用，以期提高教学质量，改进教学效果. 　　【关键词】 矩阵；线性方程组；特征值；MATLAB 　　线性代数课程是高等学校的重要数学基础课，是用数学知识解决实际问题的一个强有力的工具，广泛地应用于信号处理、系统控制、电子通信等学科领域. 线性代数课程除了培养学生的基本运算能力、抽象思维能力、逻辑推理能力外，还应注重培养学生的数学建模能力与数值计算能力，使学生学会用数学方法解决实际问题，并能利用计算机进行一定的科学计算. 　　目前传统的线性代数教材仍然是以理论为主导，偏重理论体系的系统性、完整性，过多强调证明和推导，却忽视了概念、原理和模型的实际背景，缺乏解决实际问题的训练. 再加上课程本身所固有的抽象性和计算的繁琐性，学习兴趣不高，针对线性代数课程的这些特点，必须对现有的教学模式进行改革. 利用计算机软件辅助教学，既能提高教学效率，减轻学生计算负担，促进线性代数教学理论走向应用，又能锻炼和培养学生的数值计算能力，使得学生逐渐喜欢线性代数课程. 本文主要通过实例介绍MATLAB在线性代数教学中的应用. 　　1. 矩阵的基本运算 　　（1）矩阵的初等变换 　　在线性代数中，矩阵的初等变换是最基本也是最重要的一种运算，求矩阵的逆、矩阵的秩、判断向量组的线性相关性、解线性方程组等问题都离不开初等变换，而矩阵的初等变换相对比较繁琐，故教学过程中老师会省略很多步骤，这样不利于学生的学习，引入MATLAB可以方便解决此类问题. 　　初等变换包括行变换和列变换，下面重点介绍MATLAB在初等行变换中的命令，列变换命令类似. （I）交换矩阵A第i行和第j行，A（[j，i]，：）=A（[i，j]，：）；（II）矩阵A第i行乘以k倍，A（i，：）=k*A（i，：）；（III）将矩阵A第i行的k倍加到第j行上去，A（j，：）=k*A（i，：）+A（j，：） . 　　例1 设A = 0 -2 13 0 -2-2 3 0，用初等变换求A-1. 　　解 将（A，E）通过初等行变换化成（E，A-1），求得逆矩阵A-1. MATLAB程序和运行结果： 　　 A=[0 -2 1；3 0 -2；-2 3 0]； 　　 B=[A，eye（3）]； 　　B（2，：）=1*B（3，：）+B（2，：）； B（3，：）=2*B（2，：）+B（3，：）； B（3，：）=4*B（1，：）+B（3，：）； 　　B（1，：）=2*B（3，：）+B（1，：）； B（2，：）=（-3）*B（3，：）+B（2，：）； B（2，：）=2*B（1，：）+B（2，：）； 　　 B（[1，2]，：）=B（[2，1]，：）； B（[3，2]，：）=B（[2，3]，：） 　　执行后，输出结果为： 　　B = 1 0 0 6 3 4 　　0 1 0 4 2 3 　　0 0 1 9 4 6 　　即求得A-1 = 6 3 44 2 39 4 6 　　（2）矩阵的其他典型计算 　　矩阵A的典型计算还包括：矩阵的转置A′；方阵行列式det（A）；矩阵逆inv（A）；矩阵的秩rank（A）；矩阵的行最简型rref（A）. 例如例题1中的矩阵A，MATLAB命令det（A）、inv（A）分别得到运行结果是1和3. 　　2. 行列式与方程组求解 　　根据克拉默法则，若线性方程组系数行列式D ≠ 0，则方程组有唯一解. 当线性方程组为AX = B，则用矩阵左除X=A/B，即X = A-1B；当线性方程组为XA = B，则用矩阵右除X=B/A，即X = BA-1. 　　例2 当a取何值时，线性方程组 　　（1 - a）x1 - 2x2 + 4x3 = 12x1 + （3 - a）x2 + x3 = 4x1 + x2 + （1 - a）x3 = 2有唯一解？ 　　解 当系数行列式不等于零时，方程组有唯一解. 因此，先求解a使得系数行列式等于零. MATLAB计算程序和运行结果： 　　syms a 　　format rat 　　A=[1-a，-2，4；2，3-a，1；1，1，1-a]； b=[1，4，2]； D=det（A）； 　　a0=solve（D） %系数行列式等于零的根 　　X=A＼b %利用左除求出方程组的符号解 　　X=subs（X，a，-1） %特别的当时，方程组的根 　　执行后，输出结果为： 　　a0 = 　　[ 0， 2， 3] 　　X = 　　（- a^2 + 4*a + 2）/（a*（a^2 - 5*a